| 9 | 1/1 | 返回列表 |
| 查看: 2574 | 回復: 8 | ||
z65668396金蟲 (正式寫手)
|
[求助]
Matlab用牛頓法求解非線性方程組問題
|
|
我的問題是 比如這個方程組 x1.^2-10x1+x2.^2+A=0 x1*x2.^2+x1-10x2+8=0 其中A是在一定范圍變化的 例如從-5~5; 用牛頓法解這個方程組,繪制方程組的解隨A變化的曲線 如果A是常量 牛頓法很簡單 比如A=8 程序如下: 牛頓迭代M文件 function s=NewtonIterate(x,eps) %Newton迭代法求非線性方程組 %x為迭代初值,eps為允許誤差值 if nargin==1 eps=1.0e-6; elseif nargin<1 error return end x1=fx1(x);%非線性方程組 x2=-dfx1(x);%非線性方程組導數 x3=inv(x2); x0=x3*x1'; while norm(x0)>=eps %循環(huán)迭代 x=x0'+x; x1=fx1(x); x2=-dfx1(x); x3=inv(x2); x0=x3*x1'; end s=x0'+x; return 非線性方程組M文件 function y=fx1(x) y(1)=x(1)*x(1)-10*x(1)+x(2)*x(2)+8; y(2)=x(1)*x(2)*x(2)+x(1)-10*x(2)+8; y=[y(1) y(2)]; 偏導M文件 function y=dfx1(x) y(1)=2*x(1)-10; y(2)=2*x(2); y(3)=x(2)*x(2)+1; y(4)=2*x(1)*x(2)-10; y=[y(1) y(2);y(3) y(4)]; 如果A是變量,就要改動非線性方程組的M函數, function y=fx1(x) A=-5:5; y(1)=x(1)*x(1)-10*x(1)+x(2)*x(2)+A; y(2)=x(1)*x(2)*x(2)+x(1)-10*x(2)+8; y=[y(1) y(2)]; 運行報錯 不知道哪位大神能夠解答一下 非常感謝 |

木蟲 (職業(yè)作家)
中國特色

木蟲 (職業(yè)作家)
中國特色

木蟲 (職業(yè)作家)
中國特色

金蟲 (正式寫手)

木蟲 (職業(yè)作家)
中國特色

金蟲 (正式寫手)

鐵桿木蟲 (職業(yè)作家)
金蟲 (正式寫手)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0854電子信息求調劑 +3 | α____ 2026-03-22 | 3/150 |
|
|
[考研] 北科281學碩材料求調劑 +6 | tcxiaoxx 2026-03-20 | 6/300 |
|
|
[考研] 275求調劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 17/850 |
|
|
[考研] 291 求調劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
[考研] 354求調劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考研] 286分人工智能專業(yè)請求調劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 311求調劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 求調劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 一志愿中國石油大學(華東) 本科齊魯工業(yè)大學 +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 一志愿武理材料305分求調劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境專碩,總分308求調劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 295求調劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 288求調劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 086500 325 求調劑 +3 | 領帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 材料學碩318求調劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 085601專碩,總分342求調劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|