| 5 | 1/1 | 返回列表 |
| 查看: 1605 | 回復: 4 | ||
世紀乘風銀蟲 (小有名氣)
|
[求助]
使用matlab 中的fsolve 老提示出錯
|
|
我最近想用Matlab 中fsolve求解一組非線性方程。編寫程序如下: function f= myfun(t) syms x a b t0 t4; a=5000;b=10000;t0=0;t4=1; t1=t(1);t2=t(2);t3=t(3); f(1)=-int(a+b*x,t1,t2)+(t1-t0)*(a+b*t1); f(2)=-int(a+b*x,t2,t3)+(t2-t1)*(a+b*t2); f(3)=-int(a+b*x,t3,t4)+(t3-t2)*(a+b*t3); end 但在調(diào)用fsolve(@myfun,[0.2 0.5 0.7]) 計算時,老報錯,如下: Undefined function or method 'isfinite' for input arguments of type 'sym'. Error in ==> trustnleqn at 27 if any(~isfinite(Fvec)) Error in ==> fsolve at 366 [x,FVAL,JACOB,EXITFLAG,OUTPUT,msgData]=... 如果不用M函數(shù)文件,以上是有結果的,但不知編成M函數(shù)文件時始終提示不對。 望高手幫忙,萬分感謝。 |
木蟲 (職業(yè)作家)
中國特色
|
改成下面的形式: function f=myfun(t) syms x a b t0 t4 t1=t(1);t2=t(2);t3=t(3); f1=-int(a+b*x,t1,t2)+(t1-t0)*(a+b*t1); f2=-int(a+b*x,t2,t3)+(t2-t1)*(a+b*t2); f3=-int(a+b*x,t3,t4)+(t3-t2)*(a+b*t3); y1=subs(f1,{a,b,t0,t4},{5000 10000 0 1}); y2=subs(f2,{a,b,t0,t4},{5000 10000 0 1}); y3=subs(f3,{a,b,t0,t4},{5000 10000 0 1}); f=[y1,y2,y3]; end |

銀蟲 (小有名氣)
銀蟲 (小有名氣)

| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑(0805) +6 | 煙汐憶海 2026-03-16 | 16/800 |
|
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學,080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 321求調(diào)劑 +3 | 何潤采123 2026-03-18 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學666 2026-03-16 | 6/300 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 一志愿,福州大學材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
[考研] 學碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
[考研] 321求調(diào)劑(食品/專碩) +3 | xc321 2026-03-12 | 6/300 |
|