通用牛顿法函数文件 Matlab答案

【问题标题】：General Newton Method Function File Matlab通用牛顿法函数文件 Matlab
【发布时间】：2017-03-10 03:46:09
【问题描述】：

我对 Matlab 的了解非常有限，我正在尝试制作一个通用的 Newton Raphson 函数，但每次都会出现错误，提示输入参数不足。我的代码需要三个输入，f（函数）、c0（初始猜测）和 n（步数）。到目前为止，这是我的代码：

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)

我在脚本文件中写了这个函数

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

然后我将其放入命令窗口 [c]=Q3(g(x),1,n) 希望它可以工作，但显然我的代码需要工作。谢谢

【问题讨论】：

【解决方案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。

【讨论】：