Index exceeds the number of array elements?

2020 5 27
1 个回答
回答已采纳
更新时间：2020 5 27
5 次查看（30 天）

0 个投票

Hello,
I am trying to extract data from Xsteam.I am trying to retrive data with respect to temperature. the code is as follows.
p = 10;
a = 1000;
b = 0.1;
rho_s = zeros(1,((a-1)/b));
L = length(rho_s);
for i = 1:L
T(i)= 1+b*i;
rho_s(i) =XSteam('rho_pT', p, T(i));
drho_s_dT = (rho_s(i+1)-rho_s(i))/(T(i+1)-T(i));
end
when i run this code I am getting an error as "Index exceeds the number of array elements". Can some one check and help me how can I do this.

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

John Marriott
John Marriott 2020-5-27
Hi
on this line - drho_s_dT = (rho_s(i+1)-rho_s(i))/(T(i+1)-T(i)); you are trying to take the ith+1 element of rho_s without having created it first. E.g on the first iteration you are trying to subtract the first value in rho_s from the value generated in the second iteration of the loop but this hasnt been created yet
perhaps instead you could use
for i = 1:length(rho_s)
T(i)= 1+b*i;
rho_s(i) =XSteam('rho_pT', p, T(i));
if i >=2
drho_s_dT(i) = (rho_s(i)-rho_s(i-1))/(T(i)-T(i-1));
end
end
also i added a loop index to the variable drho_s_dT becasue it didnt make sense doing pressumably you want to save each time its performing this calculation or else what is the point?
Rik
Rik 2020-5-27
Do you mean the XSteam function published here?
vishnuvardhan naidu tanga
编辑：Rik 2020-5-27
@Rik yes the Xsteam data mention there
@John Marriott Thank you so much for the clarification I will have a look into it.
@John Marriott
I have run the code. It works but after 3729 the column turns to NaN. can you please help me in that.
John Marriott
John Marriott 2020-5-27
Ok - the columns of which variable outputs NaN?- it doesnt return NaN when I run it - are you using the XSteam function posted by Rik?
yes and I have changed the code a little. I have consider only temperature instead of pressure.
a = 1000;
b = 0.1;
rho_s = zeros(1,((a-1)/b));
L = length(rho_s);
for i = 1:L
T(i)= 1+b*i;
rho_s(i) =XSteam('rhoV_T', T(i));
if i >=2
drho_s_dT(i) = (rho_s(i)-rho_s(i-1))/(T(i)-T(i-1));
end
end

请先登录，再进行评论。

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

 采纳的回答

Ameer Hamza
Ameer Hamza 2020-5-27

0 个投票

It seems that you are trying to find a derivative of rho_s w.r.t. T, you can use gradient() function
p = 10;
a = 1000;
b = 0.1;
rho_s = zeros(1,((a-1)/b));
L = length(rho_s);
T = zeros(1, L);
for i = 1:L
T(i)= 1+b*i;
rho_s(i) =XSteam('rho_pT', p, T(i));
end
drho_s_dT = gradient(rho_s, T);

2 个评论
显示 无 隐藏 无

Hello Ameer Hamza
Thank you for the suggestion it works. But when i consider interms of temperature only again from 3729 the column turns to NaN. is there any suggestion for that. the code is
a = 1000;
b = 0.1;
rho_s = zeros(1,((a-1)/b));
L = length(rho_s);
for i = 1:L
T(i)= 1+b*i;
rho_s(i) = XSteam('rhoV_T', T(i));
end
drho_s_dT = gradient(rho_s, T);
Ameer Hamza
Ameer Hamza 2020-5-27
This is probably caused by the thermodynamical equation implemented in XSteam() and not a MATLAB specific issue. You will need to understand the equations and the implementation of this function to see what might be going wrong. Also, you aren't passing variable 'p' in your code, which might be causing the issue.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Matrix Indexing 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by