时变微分方程组的幅频曲线如何画？

IceHill

我浏览了许多关于时频曲线画法的帖子，但还是不得要领。特别是我还不会使用maple。
哪位高手能给我一个matlab算法的具体思路（最好有个例子什么的）。
问题如下：AX‘’+B(t)X‘+C(t)=0,其中X为4*1矩阵，A为4*4矩阵，B(t)和C(t)都是含有sin(t)的矩阵，
我现在想画出X中两个变量的时频曲线，该怎么办？
初次拜访，请多指教！

[ 本帖最后由 花如月 于 2007-9-7 12:56 编辑 ]

无水1324

X中的两个变量是什么意思？
是什么变量？

IceHill

不好意思，没说清楚。
就是说X是4*1列向量，2个代表位移，两个代表受力
我想知道头两个变量的幅频曲线。
我查了一些资料，用谐波平衡法求解，可是我所知道的谐波法
都是一些1元或者2元的，都是先手算，然后在代入matlab
现在是四元的，而且是变系数的，手工明显不可以了，所以
我想知道matlab里面有没有一些库函数可以直接算。
谢谢大家。

fineshang

我也在这个问题上打转呢，我是12自由度的二阶微分方程组，刚度矩阵是时变的，方程组矩阵形式为：
[M]X''+CX'+K(t)X=F(t)，现在只是刚刚解出时域数值解，我还想得出频域解，做出各变量随频率的变化情况曲线图。请做过这方面研究的朋友，指点一下，非常着急，非常感谢！！！

花如月

请楼上的和楼主把问题描述清楚，新来的先阅读一下置顶帖子了解一下版规。这样可以节省你们的时间，也可以节省我们的时间。不要让我们带着疑问去探讨你们的问题。

最后说明一下，本版区关于函数绘图的例子和资料非常丰富。应该可以解决你们的问题，从你们的描述来看解都有了。就是一个数据绘图的问题，能用简单的例子描述清楚就好了

fineshang

方程形如：[m]{X''}+[C]{X'}+[K(t)]{X}={P}
其中[m]=diag{m1,m2,m3,m4,m5,m6,m7,m8,m9,m10,m11,m12}
    {X}=（x1，x2，x3，x4，x5，x6，x7，x8，x9，x0，x11，x12）
    [C]为12×12阶的常数矩阵
    [K]为12×12阶的时变矩阵，其中每个元素均是时间t的函数，Kij=Kij(t)
    {P}={p1(t),p2(t),p3(t),p4(t),p5(t),p6(t),p7(t),p8(t),p9(t),p10(t),p11(t),p12(t)}’
用数值解法（四阶龙格库塔法），借助matlab进行求解，解出了每个自变量不同时刻的解，但是怎么得到频域解，输出频率特性曲线呢？非常感谢，请明示！

[ 本帖最后由 eight 于 2007-9-7 13:51 编辑 ]