应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助关于matlab解决一个复杂方程问题
81/1返回列表
查看: 1218  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

傻帽1004

铁虫 (初入文坛)

[求助] 求助关于matlab解决一个复杂方程问题 已有2人参与

求助Sk=69,Ac=35.19,Nc=1.3,qb=2.15e5,Ag=15350,qf=2.13e5,l3=30,As=825,Ns=1.3,a=1.38e-3,b=2.0e-5,T=540,t=4000,qbt=1.75e5,
Dm=349,p=2.1,n2=1.78,l0=450,abt=1.405e-5,OT=510,A11=-5.345e-9,Ab=16054,n1=1.608,m2=1.02,A2=5.83e-9,qft=1.71e5,Di=300,l1=60,
l2=40,aft=1.455e-5,tf=68,Dk=1.679,qpt=1.0e5
syms,Sg
solve('(Sk/Ac)^(1/Nc)+qb*Ag*Sk+2*qf*Ag*Sk*l3=(Sg/As)^(1/Ns)+Dk*(a+b*T)*log(t)+qbt*(Ag*Sg+pi*Dm^2*p/4)+(A11*Ab^(-n1)*(Ag*Sg+pi*Dm^2*p/4)^n1*t^m2)+A2*Ab^(-n2)*(Ag*Sg+pi*Dm^2*p/4)^n2*t)*l0+abt*OT*l0+2*qft*(pi*Di^2*p*l1/4+pi*(Dm^2-Di^2)*p*l2/4+Ag*Sg*l3)+2*qpt*p-2*aft*tf*OT','Sg')
??? Maximum recursion limit of 500 reached. Use set(0,'RecursionLimit',N)
to change the limit.  Be aware that exceeding your available stack space can
crash MATLAB and/or your computer.
请问哪位大神指导我这个求解Sg的方程错在哪里了，应该怎么改正！很急 ！！！感激不尽啊！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-07-22 10:20:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
用数值解试试：fsolve
赞一下 回复此楼
2楼2015-07-22 11:13:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

傻帽1004

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by dingd at 2015-07-22 11:13:35
用数值解试试：fsolve

能具体点不，不胜感激。用数值解怎么写程序
赞一下 回复此楼
3楼2015-07-22 11:17:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

左右括号都不匹配吧，看Matlab的fsolve帮助
赞一下 回复此楼
4楼2015-07-22 12:56:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

傻帽1004

铁虫 (初入文坛)
引用回帖:
4楼: Originally posted by dingd at 2015-07-22 12:56:03
左右括号都不匹配吧，看Matlab的fsolve帮助

哥们儿  能具体帮我调出正确的解来吗
赞一下 回复此楼
5楼2015-07-22 13:39:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
傻帽1004: 金币+20, ★★★很有帮助 2015-07-23 15:39:46
CODE:
clear all;clc
format long

x0=10;

Sk=69,Ac=35.19,Nc=1.3,qb=2.15e5,Ag=15350,qf=2.13e5,l3=30,As=825,Ns=1.3,a=1.38e-3,b=2.0e-5,T=540,t=4000,qbt=1.75e5,
Dm=349,p=2.1,n2=1.78,l0=450,abt=1.405e-5,OT=510,A11=-5.345e-9,Ab=16054,n1=1.608,m2=1.02,A2=5.83e-9,qft=1.71e5,Di=300,l1=60,
l2=40,aft=1.455e-5,tf=68,Dk=1.679,qpt=1.0e5;

[x fval]=fsolve(@(Sg)     (Sk/Ac)^(1/Nc)+qb*Ag*Sk+2*qf*Ag*Sk*l3-((Sg/As)^(1/Ns)+Dk*(a+b*T)*log(t)+qbt*(Ag*Sg+pi*Dm^2*p/4)+(A11*Ab^(-n1)*(Ag*Sg+pi*Dm^2*p/4)^n1*t^m2+A2*Ab^(-n2)*(Ag*Sg+pi*Dm^2*p/4)^n2*t)*l0+abt*OT*l0+2*qft*(pi*Di^2*p*l1/4+pi*(Dm^2-Di^2)*p*l2/4+Ag*Sg*l3)+2*qpt*p-2*aft*tf*OT),x0);
x=x
fval=fval

计算结果：
x =

  62.212209181227031


方程残差fval =

     0

你给的方程太长，请检查我上述代码中方程形式和你的是否完全一致。注意，代码中方程形式是按照f(x)=0的形式，即等式右边等于0来改写的。
赞一下(2人) 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
6楼2015-07-23 12:45:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

傻帽1004

铁虫 (初入文坛)
引用回帖:
6楼: Originally posted by 月只蓝 at 2015-07-23 12:45:56
clear all;clc
format long

x0=10;

Sk=69,Ac=35.19,Nc=1.3,qb=2.15e5,Ag=15350,qf=2.13e5,l3=30,As=825,Ns=1.3,a=1.38e-3,b=2.0e-5,T=540,t=4000,qbt=1.75e5,
Dm=349,p=2.1,n2=1.78,l0=450,abt=1.405e- ...

仁兄你太狠了，能留个联系方式不
赞一下 回复此楼
7楼2015-07-23 15:40:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星空里轨迹

木虫 (正式写手)
引用回帖:
6楼: Originally posted by 月只蓝 at 2015-07-23 12:45:56
clear all;clc
format long

x0=10;

Sk=69,Ac=35.19,Nc=1.3,qb=2.15e5,Ag=15350,qf=2.13e5,l3=30,As=825,Ns=1.3,a=1.38e-3,b=2.0e-5,T=540,t=4000,qbt=1.75e5,
Dm=349,p=2.1,n2=1.78,l0=450,abt=1.405e- ...

联系方式方便给不，本人大四狗，做论文画图需要matlab希望得到你的帮助，我会给金币的，谢谢
赞一下 回复此楼
8楼2015-08-16 18:37:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 傻帽1004 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭