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

yuxinglong · 发表于 2011-4-10 23:59

回复 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 · 发表于 2011-4-11 00:01

回复 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 · 发表于 2011-4-11 00:11

本帖最后由 meiyongyuandeze 于 2011-4-11 08:12 编辑

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

yuxinglong · 发表于 2011-4-11 00:16

回复 18 # meiyongyuandeze 的帖子

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

ChaChing · 发表于 2011-4-11 00:16

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

yuxinglong · 发表于 2011-4-11 00:20

回复 20 # ChaChing 的帖子

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

wangrong631 · 发表于 2011-4-11 15:17

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

tangtao7874 · 发表于 2011-4-13 21:13

回复 15 # yuxinglong 的帖子

路过，顶一下

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 求救，求救 求固有频率和主振型

[编程技巧] 求救，求救求固有频率和主振型