Index exceeds the number of array elements (10), How do I solve this error??

1 view (last 30 days)
Devesh Singh
Devesh Singh on 15 Jun 2019
Answered: Star Strider on 15 Jun 2019
Error in Main_Program (line 101)
Tg(xx)=(((a*Tw(xx)^2-b*Tw(xx)+c)*Tw(xx)+(Ar*Ta(xx)*hw)+Ar*Ts(xx)*(0.048*Ta(xx)-9))/((a*Tw(xx)^2-b*Tw(xx)+c)...
global N;
global ee;
dt=1;
tf=24;
nm=12;
t(1)=0;
np=(tf-t(1))/dt;
dt1=(t(1)+1)/dt;
Tw0(1)=0;
l=0.05;
a=0.02612;
b=15.76;
c=2392;
d=0.048;
e=3.8213;
w=0.0;
s=5.67*10^-8;
Ta1=[15 16.5 18 22.5 24 26 27 27 26.5 25.5 24];
Ta2=[19 23 25.5 28.5 29.5 30 32 32 31.5 31 30];
Ta3=[25.5 29.5 32 35 37 37 37.5 38 37 36 34];
Ta4=[29.5 33 35 37 38.5 38.5 38.5 39 38.5 38.5 38];
Ta5=[30 32.5 34.5 37 38.5 39.5 41 41.5 41 40 39.5];
Ta6=[33 35.5 37 39.5 41 41 42 42.5 42 41.5 40.5];
Ta7=[29.5 32.5 34 35.5 36.5 37 38 38 33 33 32.5];
Ta8=[23 25.5 28.5 30 31 32 33 33 33.5 33 32.5];
Ta9=[29 29 30 31 34 35.5 36.5 37 38 38 37];
Ta10=[33 34 35 38 39.5 40 40 39 38 37];
Ta11=[23.5 25 28 29.5 31 32 34 35 34 32.5 31];
Ta12=[20 23.5 25 26 29 31 31.5 30 30 30.5 30];
RR=[15;46;74;105;135;166;181;212;243;273;304;349];
sw=[20 23 24 25 28 28 25 25 25 25 23 20];
Z=1;
Z1=10;
V1=4;
Z0=0.03;
V=(log2(Z/Z0)/log2(Z1/Z0))*V1;
hw=2.8+3*V;
Aw=1;
Ag=1.46;
Ar=Ag/Aw;
for L=1:nm;
N=RR(L);
GTOTAL=g_t();
if L==1
disp(['Variation of temperature productivity and solar intensity'])
disp(['of January'])
Ta=Ta1';
elseif L==2
disp(['Variation of temperature productivity and solar intensity'])
disp(['of February'])
Ta=Ta2';
elseif L==3
disp(['Variation of temperature productivity and solar intensity'])
disp(['of March'])
Ta=Ta3';
elseif L==4
disp(['Variation of temperature productivity and solar intensity'])
disp(['of April'])
Ta=Ta4';
elseif L==5
disp(['Variation of temperature productivity and solar intensity'])
disp(['of May'])
Ta=Ta5';
elseif L==6
disp(['Variation of temperature productivity and solar intensity'])
disp(['of June'])
Ta=Ta6';
elseif L==7
disp(['Variation of temperature productivity and solar intensity'])
disp(['of July'])
Ta=Ta7';
elseif L==8
disp(['Variation of temperature productivity and solar intensity'])
disp(['of August'])
Ta=Ta8';
elseif L==9
disp(['Variation of temperature productivity and solar intensity'])
disp(['of September'])
Ta=Ta9';
elseif L==10
disp(['Variation of temperature productivity and solar intensity'])
disp(['of October'])
Ta=Ta10';
elseif L==11
disp(['Variation of temperature productivity and solar intensity'])
disp(['of November'])
Ta=Ta11';
elseif L==12
disp(['Variation of temperature productivity and solar intensity'])
disp(['of December'])
Ta=Ta12';
end
ee=sw(L);
Tw(1)=ee;
Ts=0.0552*Ta.^(1.5);
for xx=1:11
kk(xx)=xx+7;
q(xx)=xx;
%Ts=[10 10 9.9 9.5 8.6 8.8 9.5 10.4 11.6 12.4 13.5 15.6 15.7 15.5 14.2
%14 13.1 12.2 11.5 11.2 11 10.3 10 10];
Tg(xx)=(((a*Tw(xx)^2-b*Tw(xx)+c)*Tw(xx)+(Ar*Ta(xx)*hw)+Ar*Ts(xx)*(0.048*Ta(xx)-9))/((a*Tw(xx)^2-b*Tw(xx)+c)...
+(Ar*hw)+Ar*(0.048*Ta(xx)-9))-e);
Ti(xx)=Tw(xx)/2.0+Tg(xx)/2.0;
K(xx)=0.0244+0.7673*10^-4*Ti(xx);
Cw(xx)=999.2+0.1343*Ti(xx)+0.01*10^-4*Ti(xx)^2-6.758*10^-8*Ti(xx)^3;
Ub(xx)=K(xx)/1;
Eg=0.98;
x=647.27-(Tw(xx)+273.15);
a1=3.2437814;
b1=5.86826*10^-3;
c1=1.1702379*10^-8;
d1=2.1878462*10^-3;
xg(xx)=647.27-(Tg(xx)+273.15);
Pw(xx)=165960.72*10^-(x*(a1+b1*x+c1*x^3)/((Tw(xx)+273.15)*(1+d1*x)));
Pg(xx)=165960.72*10^(-(xg(xx)*(a1+b1*x+c1*xg(xx)^3))/((Tg(xx)+273.15)*(1+d1*x)));
Pw1(xx)=(101300/760)*Pw(xx);
Pg1(xx)=(101300/760)*Pg(xx);
y=Tw();
hfg(xx)=3044205.5-1679.1109*(Tw(xx)+273)-1.1425*(Tw(xx)+273)^2;
hrwg(xx)=0.9*s*(Tw(xx)^2+Tg(xx)^2)*(Tw(xx)+Tg(xx));
hcwg(xx)=0.884*((Tw(xx)-Tg(xx))+((Pw1(xx)-Pg1(xx))/(268900-Pw1(xx)))*Tw(xx))^(1/3);
hewg(xx)=(9.15*10^-7*hcwg(xx)*(Pw1(xx)-Pg1(xx))*hfg(xx))/(Tw(xx)-Tg(xx));
hrgs(xx)=Eg*s*(Tg(xx)^2+Ts(xx)^2)*(Tg(xx)+Ts(xx));
hcgs(xx)=hw*(Tg(xx)-Ta(xx))/(Tg(xx)-Ts(xx));
Ui(xx)=hrwg(xx)+hcwg(xx)+hewg(xx);
Uo(xx)=hw*(Tg(xx)-Ta(xx))/(Tg(xx)-Ts(xx))+hrgs(xx);
Ut1(xx)=(1/(Ui(xx)+1/(Ar*Uo(xx))));
Ut(xx)=1/Ut1(xx);
Eue(xx)=(hewg(xx)/Ut(xx))*Ui(xx)*(Tw(xx)-Tg(xx));
Mw(xx)=(Eue(xx))*10^8/hfg(xx);
p(xx)=Eue(xx);
Eff(xx)=Mw(xx)*hfg(xx)*10^-6/GTOTAL(xx);
Tw(xx+1)=(10+(Ta(xx)+Tw(xx))/2)+(Tw(xx)+(Eff(xx)*GTOTAL(xx))+(Ut(xx)*Ts(xx)))+(Ub(xx)*Ta(xx))/((Cw(xx)/Aw)+Ut(xx)+Ub(xx))*w;
z11(xx)=Tw(xx);
t(xx+1)=t(xx)+dt
dx1b(xx)=Aw/Cw(xx)*(Eff(xx)*GTOTAL(xx)-Ut(xx)*(Tw(xx)-Ts(xx))-Ub(xx)*(Tw(xx)-Ta(xx)));
Tw(xx+1)=Tw(xx)+dx1b(xx)*dt1;
Tw1(xx)=(Tw(xx)*Tw0(1)+y(xx))*dt1;
dx1e(xx)=Aw/Cw(xx)*(Eff(xx)*GTOTAL(xx)-Ut(xx)*(Tw(xx+1)-Ts(xx))-Ub(xx)*(Tw(xx+1)+Ta(xx)));
dTw(xx)=(dx1b(xx)+dx1e(xx))/2;
Tw(xx+1)=Tw(xx)+dTw(xx)*dt;
Tw(xx+1)=Tw1(xx);
end
format short;
Tg1=Tg'
Tw2=Tw1';
dd=kk';
Mw1=Mw'*10^-5;
disp([' Time(hr) WatTemp GlasTemp ambienttemp SkyTemp Productivity x10^2 ']);
disp([dd Tw2 Tg1 Ta Ts Mw1]);
if L==1
figure(11),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(12),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(13),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(14),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(15),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==2
figure(21),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(22),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(23),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(24),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(25),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==3
figure(31),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(32),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(33),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(34),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(35),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==4
figure(41),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(42),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(43),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(44),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(45),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==5
figure(51),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(52),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(53),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(54),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(55),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==6
figure(61),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(62),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(63),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(64),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(65),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==7
figure(71),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(72),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(73),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(74),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(75),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==8
figure(81),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(82),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(83),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(84),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(85),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==9
figure(91),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(92),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(93),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(94),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(95),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==10
figure(101),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(102),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(103),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(104),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(105),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==11
figure(111),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(112),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(113),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(114),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(115),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
elseif L==12
figure(121),plot(dd,Tw2,'-kd',dd,Tg1,'-bs',dd,Ta,'-ro')
figure(122),plot(Mw1,Ta,'-kd',Mw1,Tw2,'-bs')
figure(123),plot(dd,GTOTAL,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(124),plot(dd,Eue,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
figure(125),plot(dd,Mw1,'-kd','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',5)
end;
end;

Sign in to answer this question.

Accepted Answer

Star Strider
Star Strider on 15 Jun 2019
Five lines above it, you initialise the for loop as:
for xx=1:11
Perhaps:
for xx=1:10
will eliminate the error

More Answers (0)

Categories

Find more on Solar Power in Help Center and File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!