Creating a Function That Finds Zeros

5 次查看（过去 30 天）

显示更早的评论

Lavorizia Vaughn 2021-9-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1463669-creating-a-function-that-finds-zeros

评论： Lavorizia Vaughn 2021-9-30

在 MATLAB Online 中打开

Hello,

I am trying to make a function of the form

function p = findmanyzeros(f, a, b, n, tol)

Which finds zeros in the interval [a, b] using the following strategy:

1. Compute n+1 equidistant points xk for k=0,...,n, between a and b

2. For k = 1,...,n, if f(xk) and f(xk−1) have different signs, compute a zero using findzero

3. The output vector p should contain all the computed zeros

I believe the code i have (which is below) is terribly wrong, but its my best shot. I would appreciate some help. thanks.

function p = findmanyzeros(f, a, b, n)
x = linspace(a,b,n+1);
for k = 1:n
    if f(x(k))*f(x(k+1))<0 
        p(k)=findzero(f);
        break;
    end

11 个评论
显示 9更早的评论隐藏 9更早的评论

Walter Roberson 2021-9-30

在 MATLAB Online 中打开

for k = 0:n

okay, that is a range of n+1 values.

Note that n will be recorded at the time the for loop is started, and changes to n after that will not affect how many iterations the loop does.

n=n+1;

Each of the n+1 times through the loop, you make n larger by 1 (which does not affect how many loops are done.) At the end, n will have become one more than twice as large as it was originally.

x = x(k);

x has not been defined yet, so that is an error. x has certainly not been defined as a function of one parameter, or as an array to be indexed by k. If it were an array, then you would have a problem because k starts at 0, and indexing from 0 is not permitted. If it were an array, then because you are assigning to the complete variable x, after the first x(k) was extracted, the entire array x would be overwritten with the single value extracted. For example if x started as [3 1 4 2] and k were 1, then x(1) would be 3, and that 3 would be written over top of the entire variable x, making x be just 3 . And then the next iteration, x would no longer have multiple items, so x(k) on the right hand side would be an index-out-of-range.

The equidistant points are supposed to be in the range a to b, but nothing in that loop uses a or b.

Walter Roberson 2021-9-30

在 MATLAB Online 中打开

You posted very similar code in another Question. In that other place, you noted that it warned about p = [p,z] and it gave you an error about findzero not being found. Is that the same issues as here?

MATLAB does not supply any findzero() function, so you would have to supply your own.

Hint: compare

N = 20;
x = 1:N;
z = [];
for K = 1 : N
   if isprime(K); z = [z, K]; end
end
z
z = 1×8
     2     3     5     7    11    13    17    19

N = 20;
x = 1:N;
P = false(1,N);
for K = 1 : N
   if isprime(K); P(K) = true; end
end
x(P)
ans = 1×8
     2     3     5     7    11    13    17    19

Notice that the second of these does not grow any arrays dynamically.

Lavorizia Vaughn 2021-9-30

yes i simply needed to rename a file to findzero thank you

请先登录，再进行评论。

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

回答（0 个）

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

类别

MATLAB Mathematics

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

产品

MATLAB

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by

Creating a Function That Finds Zeros

11 个评论
显示 9更早的评论隐藏 9更早的评论

回答（0 个）

另请参阅

类别

标签

产品

Community Treasure Hunt

Creating a Function That Finds Zeros

11 个评论 显示 9更早的评论隐藏 9更早的评论

回答（0 个）

另请参阅

类别

标签

产品

Community Treasure Hunt

WeChat

11 个评论
显示 9更早的评论隐藏 9更早的评论