一般ニュートンメソッド関数ファイルMatlab

Question

私はMatlabに関する知識が非常に限られています。一般的なニュートンラフソン関数を作成しようとしていますが、入力引数が足りないというエラーが表示されるたびに発生します。私のコードは、f（関数）、c0（初期推測）、n（ステップ数）の3つの入力が必要です。私は、コマンド・ウィンドウにこれを入れて、私はこの関数は、スクリプトファイルに

g=@(x)(sin((pi.*x)/2)+(1/x)-(10.*x)); 

を書かれている

function [c] = Q3(f,c0,n) 


for i=1:n 

c(0)=c0; 
f=f(x); 
fp=diff(f,x); 

c(i)=c(i-1)-subs(f,x,c(i-1))/subs(fp,c(i-1)); 

end 


disp(c) 

[C] = Q3（G（x）は、1：これは、これまでの私のコードです、n）それはうまくいくと思っていますが、明らかに私のコードは動作する必要があります。 おかげ

Answer 1

これはトリックを行う必要があり、あなたが述べたように、関数gが定義されている：

g=@(x)(sin((pi.*x)/2)+(1/x)-(10.*x));

をも、あなたが例えば、nとc0を定義する必要があります。

clearvars 

g=@(x)(sin((pi.*x)/2)+(1/x)-(10.*x)); 

n=10 
c0=1 

c=Q3(g,c0,n) 

をと別のファイルでは、NRの関数を記述します：

function [c] = Q3(f,c0,n) 

h=1e-4; 
c(1)=c0; 

for i=1:n 
    g=f(c(i)); 
    gp=(f(c(i)+h)-f(c(i)-h))/(2*h) 
    c(i+1)=c(i)-g/gp 
end 

disp(c) 

この場合、数値微分近似のためにh=1e-4を選択しますが、変更することができます。私はh<1e-2を提案する。