求救，求救 求固有频率和主振型

yuxinglong

回复 12 # meiyongyuandeze 的帖子

clc
clear
syms W I k;
I1=I;
I2=2*I;
I3=4*I;
k2=k;
k3=k;
%W=w^2
I0=[I1 0 0;0 I2 0;0 0 I3];
K=[0 0 0;0 k2 0;0 0 k3];
C4=1
for i=2:3
    C1=[1 0;-W*I0(i,i) 1];
    C2=[1 1/K(i,i);0 1];
    C3=C1*C2;
    C4=C3*C4;
end
C=C4*[1 0;-W*I0(1,1) 1];
%C3=[1 1/k3;-W*I3 1-W*I3/k3];
%C2=[1 1/k2;-W*I2 1-W*I2/k2];
%C1=[1 0;-W*I1 1];
%C=C3*C2*C1;
f=C(2,1);
f
W=solve(f,W);
w=sqrt(W);
w
这个程序哪有问题 怎么改善能像我说的那样呢

yuxinglong

回复 12 # meiyongyuandeze 的帖子

clc
clear
syms W I k;
I1=I;
I2=2*I;
I3=4*I;
k2=k;
k3=k;
%W=w^2
I0=[I1 0 0;0 I2 0;0 0 I3];
K=[0 0 0;0 k2 0;0 0 k3];
C4=1
for i=2:3
    C1=[1 0;-W*I0(i,i) 1];
    C2=[1 1/K(i,i);0 1];
    C3=C1*C2;
    C4=C3*C4;
end
C=C4*[1 0;-W*I0(1,1) 1];
%C3=[1 1/k3;-W*I3 1-W*I3/k3];
%C2=[1 1/k2;-W*I2 1-W*I2/k2];
%C1=[1 0;-W*I1 1];
%C=C3*C2*C1;
f=C(2,1);
f
W=solve(f,W);
w=sqrt(W);
w
这个程序哪有问题   怎么改善才能像我说的那样呢

meiyongyuandeze

回复 15 # yuxinglong 的帖
你自己先运行下程序以后，把错误代码发给我们啊！

yuxinglong

回复 18 # meiyongyuandeze 的帖子

不明白意思啊   大神  你加一下我Q   我有问题 714208227

ChaChing

"程序哪有问题?"
matlab的好处不就是互动式!? 直接试不就知道哪有问题
静下心搜下, 消化下老帖, 相信LZ行的

yuxinglong

回复 20 # ChaChing 的帖子

你对传递矩阵法熟悉不？ 就这道题  怎么用传递矩阵法编啊！

wangrong631

感觉ansys做这个挺容易的~！

tangtao7874

回复 15 # yuxinglong 的帖子

路过，顶一下