zhs4276689

木蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 3293
散金: 30
紅花: 2
帖子: 345
在線: 171.2小時
蟲號: 334368
注冊: 2007-03-30
專業(yè): 金屬基復(fù)合材料

[求助] 用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

結(jié)果：

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

請問是什么問題？如何解決？謝謝！

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有199人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有27人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有47人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

matlab求解非線性方程組已經(jīng)有16人回復(fù)
matlab用多元線性最小二乘法求解多元線性方程的問題已經(jīng)有11人回復(fù)
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復(fù)
求matlab求解一個方程組的問題已經(jīng)有6人回復(fù)
求各位前輩幫忙！matlab解復(fù)數(shù)方程，出現(xiàn)圖中的警告，怎來修改呢？謝謝了啊已經(jīng)有5人回復(fù)
matlab 求解隱式常微分方程時用ode15i 函數(shù)出現(xiàn)的問題求助已經(jīng)有3人回復(fù)
求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【求助】求助用matlab解偏微分方程已經(jīng)有15人回復(fù)
【求助】請教蟲友關(guān)于matlab差分方程的求解和畫圖程序已經(jīng)有4人回復(fù)

1樓 2012-03-05 15:07:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
zhs4276689(金幣+5): 謝謝，繼續(xù)幫我解答哦，呵呵 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
我覺得，沒什么大的問題，計(jì)算的值 -2.3283e-010距離0已經(jīng)足夠接近了，你如果還不滿足，那需要調(diào)節(jié)一些參數(shù)。

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2樓2012-03-05 15:38:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhs4276689

木蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 3293
散金: 30
紅花: 2
帖子: 345
在線: 171.2小時
蟲號: 334368
注冊: 2007-03-30
專業(yè): 金屬基復(fù)合材料

引用回帖:

樓: Originally posted by dbb627 at 2012-03-05 15:38:56:
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.3 ...

就是中間出現(xiàn)提示，覺得不舒服。
另外我在用另一種方法進(jìn)行求解時，也老出現(xiàn)這樣的問題，或者說是solve 提前結(jié)束，但計(jì)算出的值距離0也很近，這樣解出的結(jié)果可信嗎？會不會根據(jù)賦值的不同，還會出現(xiàn)一個更適合實(shí)際的結(jié)果呢？

贊一下

回復(fù)此樓

3樓2012-03-05 16:31:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhs4276689(金幣+15): 謝謝！ 2012-03-05 17:46:47

中間出現(xiàn)提示可以關(guān)閉，在optimset中設(shè)置一下，這個看你初值給的怎么樣，一般來說你要求的精度不是非常高的話，是可以接受的
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. 可以認(rèn)為已經(jīng)收斂了

  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.

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

4樓2012-03-05 17:11:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	6/300	2026-03-18 13:10 by zhq0425
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無際的草原
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 本人考085602 化學(xué)工程專碩 +16	不知道叫什么！ 2026-03-15	18/900	2026-03-17 17:05 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實(shí)習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 311求調(diào)劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky