We present a new computer code for the random phase approximation calculation with Skyrme density functional, in which the Fourier-series expansion method is used. With this code, we investigate rotational effects on low-frequency octupole vibrations emerging from superdeformed yrast states in 40 Ca and 44 Ti. We predict that the octupole instabilities toward a reflection-asymmetric non-axial shape occur at high spin.