| 查看: 328 | 回复: 1 | |||
[交流]
很简单的循环求解，为什么解不出啊？
|
|
a1=0;a2=0;b1=1; syms('r') i=1; for j=0:pi/8:pi/2 theta(i)=j; sigmaxx(i)=0.5*a1*r^(-0.5)*(1.5*cos(0.5*theta(i))+0.5*cos(2.5*theta(i)))+0.5*b1*r^(-0.5)*(3.5*sin(0.5*theta(i))+0.5*sin(2.5*theta(i)))+a2*4; sigmayy(i)=0.5*a1*r^(-0.5)*(2.5*cos(0.5*theta(i))-0.5*cos(2.5*theta(i)))+0.5*b1*r^(-0.5)*(0.5*sin(0.5*theta(i))-0.5*sin(2.5*theta(i))); sigmaxy(i)=0.5*a1*r^(-0.5)*(-0.5*sin(0.5*theta(i))+0.5*sin(2.5*theta(i)))+0.5*b1*r^(-0.5)*(1.5*cos(0.5*theta(i))+0.5*cos(2.5*theta(i))); sigma1(i)=0.5*(sigmaxx(i)+sigmayy(i))+sqrt((0.5*(sigmaxx(i)-sigmayy(i)))^2+sigmaxy(i)*sigmaxy(i)); sigma2(i)=0.5*(sigmaxx(i)+sigmayy(i))-sqrt((0.5*(sigmaxx(i)-sigmayy(i)))^2+sigmaxy(i)*sigmaxy(i)); sigma3(i)=0.3*(sigma1(i)+sigma2(i)); z(i)=(sigma1(i)-sigma2(i))^2+(sigma2(i)-sigma3(i))^2+(sigma3(i)-sigma1(i))^2-2*26*26; p(i)=solve('z(i)','r') i=i+1; end 显示： Warning: Explicit solution could not be found. > In solve at 140 ??? Subscripted assignment dimension mismatch. >> |
回复此楼» 猜你喜欢
2026年国自然面上资助率
已经有20人回复
-
一个有机合成实验室都需要哪些设备？
已经有5人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有22人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有8人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
你们的时间戳变了吗
已经有4人回复
-
产物和副产物价值比较
已经有5人回复
-
时间戳他又来了
已经有17人回复
2楼2016-02-28 17:21:13