matlab求复杂代数方程的解的问题-经管之家官网!

人大经济论坛-经管之家 收藏本站
您当前的位置> 软件培训>>

Matlab软件培训

>>

matlab求复杂代数方程的解的问题

matlab求复杂代数方程的解的问题

发布:cj56748938 | 分类:Matlab软件培训

关于本站

人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!

经管之家新媒体交易平台

提供"微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯"等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

提供微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

clearall;E=30;v=0.4;kappa=3-4*v;mju=E/(2*(1+v));E11=80;E22=12;E33=20;v12=0.28;v13=0.24;v23=0.08;G12=6;v21=(E22/E11)*v12;v31=(E33/E11)*v13;v32=(E33/E22)*v23;w=(1-v13*v31)*(1-v23*v32)-(v21+v23*v31)*(v12 ...
免费学术公开课,扫码加入


clear all;
E=30;
v=0.4;
kappa=3-4*v;
mju=E/(2*(1+v));
E11=80;
E22=12;
E33=20;
v12=0.28;
v13=0.24;
v23=0.08;
G12=6;
v21=(E22/E11)*v12;
v31=(E33/E11)*v13;
v32=(E33/E22)*v23;
w=(1-v13*v31)*(1-v23*v32)-(v21+v23*v31)*(v12+v32*v13);
c11=E11*(1-v23*v32)/w;
c12=E11*(v21+v23*v31)/w;
c22=E22*(1-v13*v31)/w;
c66=G12;
a=c22*c66;
b=c11*c22-c12^2-2*c12*c66;
c=c11*c66;
s1=sqrt((b+sqrt(b^2-4*a*c))/(2*a));
s2=sqrt((b-sqrt(b^2-4*a*c))/(2*a));
alpha1=(c11-c66*s1^2)/(c12*s1+c66*s1);
alpha2=(c11-c66*s2^2)/(c12*s2+c66*s2);
lmd=0.0:0.01:1.0;
N=length(lmd);
for m=1:N
lambda=lmd(N);
a1=2*mju*(lambda+1)*(s1^lambda);
a2=2*mju*(lambda+1)*(s2^lambda);
b1=-(lambda+1)*c66*(1+alpha1/s1)*(s1^lambda);
b2=-(lambda+1)*c66*(1+alpha2/s2)*(s2^lambda);
t=tan(lambda*pi/2);
g1=(((lambda+1)*(1-kappa)*alpha1-(lambda-1)*(1+kappa))*a1+((lambda-kappa)*(1-kappa)*s1+(lambda+kappa)*(1+kappa)*b1))/((1-kappa^2)*(t+1/t));
g2=(((lambda+1)*(1-kappa)*alpha2-(lambda-1)*(1+kappa))*a2+((lambda-kappa)*(1-kappa)*s2+(lambda+kappa)*(1+kappa)*b2))/((1-kappa^2)*(t+1/t));
l1=(a1*(alpha1+1)+b1*(s1-1))/((1+kappa)*(t+1/t));
l2=(a2*(alpha2+1)+b2*(s2-1))/((1+kappa)*(t+1/t));
x1=(-((lambda-1)*(1+kappa)+(lambda+1)*(1-kappa)*t^2*alpha1)*a1+((lambda+kappa)*(1+kappa)-(lambda-kappa)*(1-kappa)*t^2*s1)*b1)/((1-kappa^2)*(1+t^2));
x2=(-((lambda-1)*(1+kappa)-(lambda+1)*(1-kappa)*t^2*alpha2)*a2+((lambda+kappa)*(1+kappa)-(lambda-kappa)*(1-kappa)*t^2*s2)*b2)/((1-kappa^2)*(1+t^2));
y1=-(a1*(t^2*alpha1-1)+b1*(t^2*s1+1))/((1+kappa)*(1+t^2));
y2=-(a2*(t^2*alpha2-1)+b2*(t^2*s2+1))/((1+kappa)*(1+t^2));
k1=(g1+(lambda+1)*l1)*sin(lambda*pi)+(x1+(lambda+1)*y1)*cos(lambda*pi);
k2=(g2+(lambda+1)*l2)*sin(lambda*pi)+(x2+(lambda+1)*y2)*cos(lambda*pi);
k3=(x1+(lambda-1)*y1)*sin(lambda*pi)+(g1+(lambda-1)*l1)*cos(lambda*pi);
k4=(x2+(lambda-1)*y2)*sin(lambda*pi)+(g2+(lambda-1)*l2)*cos(lambda*pi);
eq(m)=k1*k4-k2*k3;
solve('eq(m)','lambda')
end
请问有没有哪位大侠能够看出程序有没有什么错误的地方?或者求出lambda的值~谢啦
「经管之家」APP:经管人学习、答疑、交友,就上经管之家!
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
本文关键词:

本文论坛网址:https://bbs.pinggu.org/thread-1652033-1-1.html

人气文章

1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责;
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
经管之家 人大经济论坛 大学 专业 手机版