Question

我需要使用牛顿方法找到函数的根。我从键盘输入间隔和准确度。这是我的代码

disp('Newton method')
fx=@(g) 5*sin(g.^3-2*g.^2-1);
fx1=@(g) 5*g*(3*g-4)*cos(-g.^3+2*g.^2+1);
fx2=inline('-5*((4-6*g)*cos(-g.^3+2*g.^2+1)-(4*g-3*g.^2).^2*sin(-g.^3+2*g.^2+1))'); 
e=input ('Enter accuracy:');
a=input ('enter a:');
b=input ('enter b:');
x0=a:e:b;
y= 5*sin(x0.^3-2*x0.^2-1);
y2= -5*((4-6*x0)*cos(-x0.^3+2*x0.^2+1)-(4*x0-3*x0.^2).^2*sin(-x0.^3+2*x0.^2+1));
plot (x0,y),grid
xlabel('x'),ylabel('y')
fa=fx(a);
n=0;
if (fa*y2>0)
    x1=a;
else
    x1=b;

end;
 while(abs(fx(x1))>e)
       n=n+1;
    x1=x1-(fx(x1))/(fx1(x1));

 end;
 disp(sprintf('Answer:%g',x1))
    disp(sprintf('Number of iterations:%g',n))

当我编译时，它说：

Error using  * 
Inner matrix dimensions must agree.

Error in Untitled3 (line 10)
y2= -5*((4-6*x0)*cos(-x0.^3+2*x0.^2+1)-(4*x0-3*x0.^2).^2*sin(-x0.^3+2*x0.^2+1));

Answer 1

您将两个1xn向量相乘，这是不可能的。这种乘法会导致错误：

y2= -5*((4-6*x0)*cos(-x0.^3+2*x0.^2+1)-(4*x0-3*x0.^2).^2*sin(-x0.^3+2*x0.^2+1));
                ^                                       ^
                |                                       |

使用逐元素乘法.*可能是解决方案，但我不知道您尝试实施的是什么。

矩阵维必须一致，牛顿法

1 个答案: