求频响函数和响应均方值

guozheng029

拜托各位帮忙！
    我要求解多自由度系统随机激励下振动的频响函数和响应均方值，已求得系统得固有频率和固有振型， 用MATLAB编程如下：（假设三个自由度）
clear;
syms x l ;
fai=sym('[sin(pi*x/l);sin(2*pi*x/l);sin(3*pi*x/l)]');
fait=fai.';
fft=fai*fait;
m1=int(fft,x,0,l);
faix1=subs(fai,{x l},{1/4,1});
faix2=subs(fai,{x l},{3/4,1});
m2=faix1*faix1.';
m3=faix2*faix2.';
M=m1+m2+m3;
C=int(0.2*fft,x,0,l);
k1=diff(fai,x,1);
K=int(1000*k1*k1.',x,0,l);
M=subs(M,{l},{10});
K=subs(K,{l},{10});
[V,D]=eig(K,M);
D=diag(D);
DD=sqrt(D);
上述程序中fai为假设模态，M，C，K分别为质量，阻尼和刚度矩阵，DD为系统固有频率
对应第 i 阶（i=1,2,3）固有频率的固有振型为 xi=假设模态的转置乘以Vj（j=1,2,3）。
接下来我该怎么求解呢？matlab程序怎么实现？

guozheng029

上面程序中倒数第四行少了一个赋值语句
C=subs(C,{l},{10});
然后进行主坐标变换将振动方程化成解耦方程，但是程序运行结果C不是一个对角矩阵，我不知道问题出哪了？？？

空山长风

自由度很少,且线性状况,你直接求解析解就可以了