How can i modify this code so x can = 0

1 次查看（过去 30 天）

Muhammad Choudhury 2022-2-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1648595-how-can-i-modify-this-code-so-x-can-0

评论： Image Analyst 2022-2-13

function qw=fNR(x)
qw=zeros(100,1)
qw(2)=x
i=2
while abs(qw(i)-qw(i-1))>1e-6
    qw(i+1)=qw(i)-ff1(qw(i))/fdf1(qw(i))
    i=i+1
end
qw=qw(2:i-1)
[qw(end) i] 
plot(qw)
end 

how can i model this functio so it can take fNR(0)

fdf1

function y=fdf1(x)
y=3*x.^2-12*x+11;
end 

ff1

function y=ff1(x)
y=x.^3-6*x.^2+11*x-5;
end 

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Image Analyst 2022-2-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1648595-how-can-i-modify-this-code-so-x-can-0#answer_894780

在 MATLAB Online 中打开

Try this:

function qw=fNR(x)
qw=zeros(100,1)
qw(2)=x
i=2
maxIterations = 1000;
loopCounter = 0;
while (abs(qw(i)-qw(i-1))>1e-6) && loopCounter < maxIterations
    qw(i+1)=qw(i)-ff1(qw(i))/fdf1(qw(i))
    i=i+1
    loopCounter = loopCounter + 1;
end
qw=qw(2:i-1)
if isempty(qw)
    message = sprintf('qw is empty')
    uiwait(warndlg(message));
    return;
end
[qw(end) i]
plot(qw)
end

5 个评论
显示 3更早的评论隐藏 3更早的评论

DGM 2022-2-12

Oops. I didn't realize I was commenting on the answer instead of on the question.

Image Analyst 2022-2-13

@Muhammad Choudhury I know that was your error with the original code. Obviously you never ran my code because my code throws up an alert that the array is empty rather than your error.

请先登录，再进行评论。

请先登录，再回答此问题。