| 查看: 382 | 回复: 0 | ||
[求助]
解符号方程的错误提示问题
|
|
for theta=(0:10:90)*3.1415926/180; syms Rp b1=1/(2-yita)*(1+rou)*alpha./Rp.*(1-kezei*Rp./(Rp-R0))+alpha*kezei*(1+rou)./((3-yita).*(Rp-R0)); b2=kezei*rou./((1-yita).*(Rp-R0))+1/yita.*(1-kezei*Rp./(Rp-R0))*rou./Rp; b3=1/(2-yita)*(1+rou)*alpha./Rp.^2.*(1-kezei.*Rp./(Rp-R0))*R0+alpha*kezei*(1+rou)./((3-yita).*(Rp-R0))*R0^2./Rp.^2; b4=kezei*rou./((1-yita).*(Rp-R0))+1/yita.*(1-kezei.*Rp/(Rp-R0))*rou/R0; myfun=thegma_c0.*(b1-b2)+(P0-thegma_c0.*(b3-b4)).*(R0./Rp).^(1-yita)-thegma_pe; ezplot(myfun); hold on end Rp=fzero(@(Rp)(thegma_c0.*(b1-b2)+(P0-thegma_c0.*(b3-b4)).*(R0./Rp).^(1-yita)-thegma_pe),4); 本人程序如上，其中Rp为待求量，其余均为常数，程序可以运行到ezplot(myfun)，运行到最后一行解方程时就提示 ??? Function 'isfinite' is not defined for values of class 'sym'. Error in ==> isfinite at 19 [varargout{1:nargout}] = builtin('isfinite', varargin{:}); Error in ==> fzero at 309 elseif ~isfinite(fx) || ~isreal(fx) Error in ==> plasarea_onlyruanhua at 48 Rp=fzero(@(Rp)(thegma_c0.*(b1-b2)+(P0-thegma_c0.*(b3-b4)).*(R0./Rp).^(1-yita)-thegma_pe),4); 不知什么原因，请大家帮忙给指点一下。 |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有273人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
找到一些相关的精华帖子，希望有用哦~
求带积分项（积分项除了积分变量还有一个符号常量）的非线性方程
已经有7人回复
- 求matlab 解方程的问题
已经有5人回复
- 求助Matlab解方程问题
已经有6人回复
- 求助：非线性方程系数回归问题？
已经有10人回复
- spss多元逐步线性回归方程的问题，求解答！
已经有4人回复
- 【求助】mathematic解方程问题
已经有4人回复
- 【求助】请教一个关于解方程的问题
已经有8人回复
- 【求助】解方程组的错误
已经有7人回复
- 【求助】解含时薛定谔方程的问题
已经有5人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
