24小时热门版块排行榜    

查看: 1177  |  回复: 3

邪恶好小子

银虫 (小有名气)

[求助] 传递函数的求解方程怎么弄?

syms s
l=0.1;
b=0.020;
h1=0.008;
h2=0.005;
h3=0.003;
ro1=7800;
ro2=2600;
ro3=7800;
E1=2.06e11;
E3=2.06e11;
G2=9.8e9*(1+0.1i);
c=1+(h1+h3)/(2*h2);
I1=b*h1^3/12;
I3=b*h3^3/12;
D=E1*I1+E3*I3;
A1=b*h1;
A2=b*h2;
A3=b*h3;
m=ro1*A1+ro2*A2+ro2*A2;
K=E1*A1*E3*A3/(E1*A1+E3*A3);
Y=K*c^2/D;
g=(G2*A2)/(K*h2^2);
F=[0,1,0,0,0,0;0,0,1,0,0,0;0,0,0,1,0,0;0,0,0,0,1,0;0,0,0,0,0,g;m*s^2/D,0,-m*s^2/(g*D),0,(1+Y),0];
M=[1,0,0,0,0,0;0,1,0,0,0,0;0,m*s^2/(g*D),0,-Y/g,0,1/g;0,0,0,0,0,0;0,0,0,0,0,0;0,0,0,0,0,0];
N=[0,0,0,0,0,0;0,0,0,0,0,0;0,0,0,0,0,0;0,m*s^2/g,0,-D*(1+Y),0,D;m*s^2/g,0,-D*(1+Y),0,D/g,0;m*s^2/g,0,-D*Y,0,D/g,0];
efs=expm(F);
si=det(M+N*efs);
solve(si,s)
这样可以求解S么?我试过,计算机直接死机,请问有什么更好的解非线性方程的?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

杨海

新虫 (初入文坛)

题主这个问题解决了吗?我现在也遇到这个问题了,不知怎么弄。
2楼2015-04-25 06:24:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)

是不是那个e^At形式的传递函数??
如果是的话,拉普拉斯变换
3楼2015-04-25 16:29:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

杨海

新虫 (初入文坛)

题主这个问题解决了吗?怎么解决的?
4楼2015-04-25 21:03:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 邪恶好小子 的主题更新
信息提示
请填处理意见