| 3 | 1/1 | 返回列表 |
| 查看: 2888 | 回復(fù): 2 | |||
njutsjh木蟲(chóng) (小有名氣)
|
[求助]
關(guān)于用Broyden法解非線(xiàn)性方程組的問(wèn)題 已有1人參與
|
|
我用Broyden法解非線(xiàn)性方程組的時(shí)候,Matlab需要計(jì)算很長(zhǎng)時(shí)間才能出結(jié)果,是不是代碼有不合理的地方,該怎么修改,謝謝大家! 代碼如下: function =Broyden(func,u,TOL,m) % 應(yīng)用Broyden方法求非線(xiàn)性方程組法f(x)=0。 % func為待求解的非線(xiàn)性方程組。 % 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ù)版主 (文壇精英)
![]() |
專(zhuān)家經(jīng)驗(yàn): +518 |
金蟲(chóng) (正式寫(xiě)手)
| 試試這個(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 | 家佳佳佳佳佳 2026-03-29 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085601材料工程找調(diào)劑 +14 | oatmealR 2026-03-29 | 15/750 |
|
|
[考研] 329求調(diào)劑,一志愿西北工業(yè)大學(xué),材料工程(085601) +4 | 小小機(jī)靈蟲(chóng) 2026-03-29 | 10/500 |
|
|
[考研] 求調(diào)劑 +4 | 青春裁為三截 2026-03-29 | 4/200 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 343求調(diào)劑 +5 | 愛(ài)羈絆 2026-03-28 | 5/250 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
[考研] 320分,材料與化工專(zhuān)業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 環(huán)境專(zhuān)碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 化工專(zhuān)碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開(kāi)開(kāi)心心沒(méi)煩惱 2026-03-23 | 4/200 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問(wèn)津 2026-03-23 | 5/250 |
|