| 查看: 1438 | 回复: 3 | ||||
zhs4276689木虫 (正式写手)
|
[求助]
用matlab求解方程出问题，请帮忙看看
|
|
我求解一个方程，中间虽然解出来了，但还是有些问题，请帮忙看下 程序： clc; Ts=800; Cs=1134;ks=218; Ps=2591; Tm=933; Cl=921; z=0.65 as=ks/Cs/Ps t=ks*(Tm-Ts)*exp(-z^2)/sqrt(as)/(1+erf(z)) f31=@(x)(t-(60.3+0.033*x)*(x-Tm)*exp(-z^2*as/(2.39e-5+1.9e-8*x))/sqrt(2.39e-5+1.9e-8*x)*erf(z*sqrt(as/(2.39e-5+1.9e-8*x)))) x=fsolve(f31,1100) Ts f(1)=t-(60.3+0.033*x)*(x-Tm)*exp(-z^2*as/(2.39e-5+1.9e-8*x))/sqrt(2.39e-5+1.9e-8*x)*erf(z*sqrt(as/(2.39e-5+1.9e-8*x))); f 结果： Equation solved, fsolve stalled. fsolve stopped because the relative size of the current step is less than the default value of the step size tolerance squared and the vector of function values is near zero as measured by the default value of the function tolerance. x = 1.1734e+003 Ts = 800 f = -2.3283e-010 请问是什么问题？如何解决？谢谢！ |
回复此楼» 猜你喜欢
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有195人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab求解非线性方程组
已经有16人回复
- matlab用多元线性最小二乘法求解多元线性方程的问题
已经有11人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- matlab 求解隐式常微分方程时用ode15i 函数出现的问题求助
已经有3人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】求助用matlab解偏微分方程
已经有15人回复
- 【求助】请教虫友关于matlab差分方程的求解和画图程序
已经有4人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhs4276689(金币+5): 谢谢，继续帮我解答哦，呵呵 2012-03-05 16:28:18
感谢参与，应助指数 +1
zhs4276689(金币+5): 谢谢，继续帮我解答哦，呵呵 2012-03-05 16:28:18
|
clc; Ts=800; Cs=1134;ks=218; Ps=2591; Tm=933; Cl=921; z=0.65 as=ks/Cs/Ps t=ks*(Tm-Ts)*exp(-z^2)/sqrt(as)/(1+erf(z)) f31=@(x)(t-(60.3+0.033*x)*(x-Tm)*exp(-z^2*as/(2.39e-5+1.9e-8*x))/sqrt(2.39e-5+1.9e-8*x)*erf(z*sqrt(as/(2.39e-5+1.9e-8*x)))) [x,feval,exflag]=fsolve(f31,10,optimset('TolFun',1e-8)) [x,feval,exflag]=fsolve(f31,10,optimset('TolFun',1e-8)) Equation solved, fsolve stalled. fsolve stopped because the relative size of the current step is less than the default value of the step size tolerance squared and the vector of function values is near zero as measured by the selected value of the function tolerance. x = 1.1734e+003 feval = -2.3283e-010 exflag = 2 我觉得，没什么大的问题，计算的值 -2.3283e-010距离0已经足够接近了，你如果还不满足，那需要调节一些参数。 |
2楼2012-03-05 15:38:56
zhs4276689
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3293
- 散金: 30
- 红花: 2
- 帖子: 345
- 在线: 171.2小时
- 虫号: 334368
- 注册: 2007-03-30
- 专业: 金属基复合材料
3楼2012-03-05 16:31:37
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhs4276689(金币+15): 谢谢！ 2012-03-05 17:46:47
zhs4276689(金币+15): 谢谢！ 2012-03-05 17:46:47
|
中间出现提示可以关闭，在optimset中设置一下，这个看你初值给的怎么样，一般来说你要求的精度不是非常高的话，是可以接受的 clc; Ts=800; Cs=1134;ks=218; Ps=2591; Tm=933; Cl=921; z=0.65 as=ks/Cs/Ps t=ks*(Tm-Ts)*exp(-z^2)/sqrt(as)/(1+erf(z)) f31=@(x)(t-(60.3+0.033*x)*(x-Tm)*exp(-z^2*as/(2.39e-5+1.9e-8*x))/sqrt(2.39e-5+1.9e-8*x)*erf(z*sqrt(as/(2.39e-5+1.9e-8*x)))) [x,feval,exflag]=fsolve(f31,10,optimset( 'Display','off',TolFun',1e-8)) x = 1.1734e+003 feval = -2.3283e-010 exflag = 2 fsolve的exflag = 2 表示 2 Change in X too small. 可以认为已经收敛了 1 FSOLVE converged to a root. 2 Change in X too small. 3 Change in residual norm too small. 4 Computed search direction too small. 0 Too many function evaluations or iterations. -1 Stopped by output/plot function. -2 Converged to a point that is not a root. -3 Trust region radius too small (Trust-region-dogleg) or Regularization parameter too large (Levenberg-Marquardt). -4 Line search failed. |
4楼2012-03-05 17:11:12
10