幸儿李
2010-04-16, 15:44
各位大虾,利用下面的程序求出矩阵c之后,若想通过“c(2,1)=0”来求解W(也就是系统的各阶固有频率),该采用什么样的方法?新手上路,请大虾们指教一二
% function f=cdjz(n)
w=sym('w');
c=ones(2);
% f1=sym('[x;m]');
j=[46474 879.69 462.3 445.15 528.63 193.45 5.376 12.124 2.912 11.611 4.935];
k=[0 36.186 0.0297 4.1489 22.21 0.2329 19.529 0.1338 0.0117 0.0146 0.1717];
for i=2:11
c=c*[1 1/k(i);-j(i)*w^2 1-j(i)*w^2/k(i)];
end
% f=c*f1;
% function f=cdjz(n)
w=sym('w');
c=ones(2);
% f1=sym('[x;m]');
j=[46474 879.69 462.3 445.15 528.63 193.45 5.376 12.124 2.912 11.611 4.935];
k=[0 36.186 0.0297 4.1489 22.21 0.2329 19.529 0.1338 0.0117 0.0146 0.1717];
for i=2:11
c=c*[1 1/k(i);-j(i)*w^2 1-j(i)*w^2/k(i)];
end
% f=c*f1;