求助：在matlab中怎样用fft求非线性方程数值解后的频谱特性？

yannanr

向大家请教：
一个非线性系统，振动方程是达芬方程D2x+uDx+w^2*x+a*x^3=f(wt)
用龙格－库塔法求解它的数值解后，在画出它的频谱特性，fft可以用吗？
数值求解之后得到的x是一系列的数值，fft可以对这个数值向量进行频谱变换吗？