24小时热门版块排行榜

返回列表

tshgth

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 911.1
散金: 802
红花: 7
帖子: 211
在线: 259.1小时
虫号: 77696
注册: 2005-07-01
专业: 岩土与基础工程

[求助] 解符号方程的错误提示问题

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);
不知什么原因，请大家帮忙给指点一下。

回复此楼

» 猜你喜欢

石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有233人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有1人回复
投稿chemical physical letters不送审？已经有2人回复
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生已经有4人回复

1楼 2012-05-05 10:15:06

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 tshgth 的主题更新

返回列表