| 3 | 1/1 | 返回列表 |
| 查看: 2890 | 回復(fù): 2 | |||
njutsjh木蟲 (小有名氣)
|
[求助]
關(guān)于用Broyden法解非線性方程組的問題 已有1人參與
|
|
我用Broyden法解非線性方程組的時(shí)候,Matlab需要計(jì)算很長時(shí)間才能出結(jié)果,是不是代碼有不合理的地方,該怎么修改,謝謝大家! 代碼如下: function =Broyden(func,u,TOL,m) % 應(yīng)用Broyden方法求非線性方程組法f(x)=0。 % func為待求解的非線性方程組。 % u為初始近似向量。 % TOL為誤差容限。 % m為最大迭代次數(shù)。 N=size(func); n=N(1); var=findsym(func); A=subs(jacobian(func),var,u); v=subs(func,var,u); H=inv(A); s=-H*v; u=u+s; for k=1:m w=v; v=subs(func,var,u); y=v-w; z=-H*y; p=-s'*z; if(p==0) disp('Method failed!'); return; end C=p*eye(n)+(s+z)*s'; H=C*H/p; s=-H*v; u=u+s; disp(vpa(u',6)); if(max(abs(s))<TOL) break; end end u=vpa(u,6); |
榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
金蟲 (正式寫手)
| 試試這個(gè),你源程序里面太多符號(hào)運(yùn)算,還有就是你的Jacobi矩陣是調(diào)用http://www.pcs.cnu.edu/~bbradie/matlab/systems/broyden.m |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +3 | Amber00 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 309求調(diào)劑 +17 | 誰不是少年 2026-03-29 | 17/850 |
|
|
[考研] 0856材料化工調(diào)劑 總分330 +17 | zhubinhao 2026-03-27 | 17/850 |
|
|
[考研] 085600 一志愿9 總分351 求調(diào)劑學(xué)校 +7 | czhcz 2026-03-31 | 7/350 |
|
|
[考研] 304求調(diào)劑 +8 | 素年祭語 2026-03-31 | 11/550 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +15 | 101次希望 2026-03-28 | 15/750 |
|
|
[考研] 08開頭看過來。! +3 | wwwwffffff 2026-03-31 | 5/250 |
|
|
[考研] 330分求調(diào)劑 +6 | qzenlc 2026-03-29 | 6/300 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-29 | 6/300 |
|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +7 | 里子木yy 2026-03-29 | 7/350 |
|
|
[考研] 313求調(diào)劑 +6 | 賣個(gè)關(guān)子吧 2026-03-31 | 6/300 |
|
|
[考研]
|
Gymno 2026-03-30 | 6/300 |
|
|
[考研]
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
[考研] 材料專碩 085600求調(diào)劑 +7 | BBQ233 2026-03-30 | 7/350 |
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|