英雄論壇

不是都已經po在板上了 ..= ="

我....也要做類似的東西ㄋ......
阿來....把程式PO 上來讓我改一改吧....
cccccc

阿~已經寫完ㄌ
為了慶祝我寫完了

還打了好久的坦克呢

看~懶趴熊~

阿來不要偷懶。

..........

function [xv,it]=broyden(x,f,n,tol)
fr=zeros(n,1); it=0; xv=x;
Br=eye(n);
fr=feval(f, xv);
while norm(fr)>tol
it=it+1;
pr=-Br*fr;
tau=1;
xv1=xv+tau*pr; xv=xv1;
oldfr=fr; fr=feval(f,xv);
y=fr-oldfr; oldBr=Br;
oyp=oldBr*y-pr;pB=pr'*oldBr;
for i=1:n
for j=1:n
M(i,j)=oyp(i)*pB(j);
end;
end;
Br=oldBr-M./(pr'*oldBr*y);
end;

<font size=-1>[ 這篇文章被編輯過: LAI 在 2001-11-18 11:46 ]</font>

在 2001-11-17 23:37, onitsuka[Fantasy Ranger] 寫了:
外星人

o_o?外星人?!在哪阿

感謝幫忙挑錯
現在已經有結果摟
fr =

-11.5000
-2.8660


rootvals =

0.9397
0.3420


iter =

12

(2)Broyden's method
還在膠著中= =

喔喔喔....
ㄚ來真是數值分析女王ㄚ....

...............
誇模~~
這是亂碼ㄇ?

外星人

可以算但是算不出結果o_o!

in command window
>> [rootvals,iter]=newtonmv([1,2]','f308','f309',2,0.005)

function jf=f309(v)
x=v(1) ;y=v(2);
jf=zeros(2,2);
jf(1,:)=[3*x^2-3*y -6*x*y];
jf(2,:)=[6*x*y 3*x^2-3*y^2];

<font size=-1>[ 這篇文章被編輯過: LAI 在 2001-11-18 11:10 ]</font>

function f=f308(v)
x=v(1);y=v(2);
f=zeros(2,1)
f(1)=x^3-3*x*y^2-0.5;
f(2)=3*(x^2)*y-y^3-(3^.5)/2;