likefeng 发表于 2008-12-8 19:20

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

在此类微分方程组中，参数均含有正余旋函数。怎样才能得到附件的幅频图？请大家提示。微分方程组如下：（图中横坐标中的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 发表于 2008-12-8 19:23

微分方程组如下：

咕噜噜 发表于 2008-12-8 19:38

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

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

likefeng 发表于 2008-12-8 19:44

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

likefeng 发表于 2008-12-9 09:20

大家给点意见吧，如果右端只是简单的函数值，我可令x1=asinwt,x2=bsinwt,带入方程求出幅值与频率关系，但当右端为复杂函数关系时，如何做出附件的幅频图，使横坐标为w1/w0?
页: [1]
查看完整版本: 请问：如何做出此类方程组的幅频曲线？