请问：如何做出此类方程组的幅频曲线？

likefeng

在此类微分方程组中，参数均含有正余旋函数。怎样才能得到附件的幅频图？请大家提示。微分方程组如下：（图中横坐标中的w0就是公式中的w ,w2=2w1）

F1=2+4*cos(1.6*t) -0.3*w^2*cos(w*t).*cos(1.6*t;

F2=9*cos(1.6*t)+0.09*w^(2)*sin(w*t).*sin(w*t).*sin(1.6*t).*cos(1.6*t)

C12=1078.62*(-0.3)*sin(1.6*t);

C21=1078.62*0.3 *sin(1.6*t);

K11=0.09*w^(2)*sin(w*t).*sin(w*t) -0.3*w*sin(w*t).*cos(s*t)*1.6);

K12=0.3*w^(2)*cos(w*t).*sin(s*t)-0.3*w*sin(w*t).*cos(1.6*t)*1.6)； K21=0.09*w^(2)*sin(w*t).*sin(w*t)； K22=4；

likefeng

微分方程组如下：

咕噜噜

用数值法求解吧，饭后做FFT

[ 本帖最后由 咕噜噜 于 2008-12-8 19:39 编辑 ]

likefeng

谢谢，我用数值法做了，但前提是要先设定w1,w0的值，可得出时域内的解，但我不知道FFT变换后，横坐标表示什么？怎么得到横坐标是w1/w0?

likefeng

大家给点意见吧，如果右端只是简单的函数值，我可令x1=asinwt,x2=bsinwt,带入方程求出幅值与频率关系，但当右端为复杂函数关系时，如何做出附件的幅频图，使横坐标为w1/w0?