Coding a function to find zeros using a variant of the bisection/secant method

Question

Lavorizia Vaughn 2021-9-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1462689-coding-a-function-to-find-zeros-using-a-variant-of-the-bisection-secant-method

编辑： Walter Roberson 2023-5-23

programming_assignment_1.pdf

Hello, I have a programming assignment where I have to implement a matlab function that is a variant of the bisection and secant method. Please see attachment for exact details. I am having problems with the code. Will someone help me with the necessary fixes? My code is as follows: —————————————————————————————————————-

if true
  % code
function p = bisection(f, a, b, tol)
w=1;
for i=1:100
  p= a +(w*f(a)*(a-b)/f(b)-w*f(a));
  fprintf(a,b,p,f(p));
  if f(p)*f(b)>0
      w=0.5;
  else
      w=1;
      a=b;
  end
  b=p;
  if abs(b-a)<tol
      abs(f(p))<tol
      break;
  end
end.

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

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2021-9-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1462689-coding-a-function-to-find-zeros-using-a-variant-of-the-bisection-secant-method#answer_797194

在 MATLAB Online 中打开

No, no error about needing to use () . However,

output = bisection(@(x) sin(x)+3*cos(x).^2, -pi/2, pi/2, 1e-5)
Error using fprintf
Invalid file identifier.  Use fopen to generate a valid file identifier.

Error in solution>bisection (line 7)
    fprintf(a,b,p,f(p));
function p = bisection(f, a, b, tol)
w=1;
for i=1:100
    p= a +(w*f(a)*(a-b)/f(b)-w*f(a));
    fprintf(a,b,p,f(p));
    if f(p)*f(b)>0
        w=0.5;
    else
        w=1;
        a=b;
    end
    b=p;
    if abs(b-a)<tol
        abs(f(p))<tol
        break;
    end
end
end

This is because the first parameter to fprintf must be one of:

A file identifier
A format character vector or format string scalar
A device object for some kinds of devices such as serialport objects

The first parameter you are passing is a, which is numeric but is not a valid file identifier.

6 个评论
显示 4更早的评论隐藏 4更早的评论

Lavorizia Vaughn 2021-9-29

Yeah dude, I’m still kind of lost. You mind telling me exactly what to change?

Walter Roberson 2021-9-29

在 MATLAB Online 中打开

    fprintf('Lovorizia, the value for a was %g .\nAnd b was %g .\n, Then there was p, which was %g .\n And lastly, if you still care, f(p) was %g .\n', a, b, p, f(p));
    

请先登录，再进行评论。