Subscript indices must either be real positive integers or logicals.

2 views (last 30 days)
function T=baryeval2(tau,s,c,belta,x)
u=zeros(size(x))
d=zeros(size(x))
T=zeros(size(x))
for k=1:5
u=u+s(k).*belta(k).*(x-tau(k)).^(-1);
d=d+c(k).*belta(k).*(x-tau(k)).^(-1);
end
T=u./d;
end

Answers (4)

Azzi Abdelmalek
Azzi Abdelmalek on 27 Feb 2016
Edited: Azzi Abdelmalek on 27 Feb 2016
You cant use an index equal to zero, Matlab allows only positive integer indices, or logical indices
  4 Comments
Yintong Zheng
Yintong Zheng on 27 Feb 2016
I changed my code and now I get answers. But apparently something is wrong! Please help me to fix it.
Yintong Zheng
Yintong Zheng on 27 Feb 2016
I insert all the values of tau, s, c, belta at command window. And all those '-Inf' are what I got.

Sign in to comment.


the cyclist
the cyclist on 27 Feb 2016
MATLAB has 1-based indexing. When you enter
tau(k)
for k = 0, there is no such element.

Walter Roberson
Walter Roberson on 27 Feb 2016
You have, in part, (x-tau(k)).^(-1) . If there is any x equal to tau(k) then the subtraction would give 0 and raising that to power -1 would give infinity. Infinity times any non-zero value gives infinity (infinity times 0 gives nan), so you would be adding infinity to the value, resulting in infinity. And then since you are effectively totaling values, you are going to get an infinite result.
  2 Comments
Walter Roberson
Walter Roberson on 27 Feb 2016
If you do have an x equal to some tau, then inf or nan is the correct answer for the formulae you have written. The situation is exactly like asking to total 1/x when x is allowed to be 0.

Sign in to comment.


Image Analyst
Image Analyst on 27 Feb 2016

Categories

Find more on Startup and Shutdown in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!