can anyone help me for loop

I have question for example w=15*(ts-12), w start from some value not from 1 to some value and the number of values 360 , how can make a loop start from the first value of w to last value
for w=15*(ts-12)
H(w)=sind(w+1)-sind(w)
end

Respuestas (1)

Azzi Abdelmalek
Azzi Abdelmalek el 15 de Jul. de 2016
Editada: Azzi Abdelmalek el 15 de Jul. de 2016

0 votos

w=linspace(some_value1,some_value2,360)
H=sind(w+1)-sind(w)

2 comentarios

Jose Louis
Jose Louis el 15 de Jul. de 2016
Editada: Azzi Abdelmalek el 15 de Jul. de 2016
thank you, Mr. Azzi, but the value of w not constant, anyway this my code the problem with the loop i don't know how to write the loop
n=1:1:365 ; %the day of the year [a number 1 through 365]
fi=19.9450; %is the latitude
GMT=2 ;% Time Zone for example +2
LT=5; % local time (LT).
lambda=50.067;% longitude
Zc=15*GMT;%Local Standard Time Meridian
B=(360/365)*(n-81)*pi/180;% in degrees and d is the number of days since the start of the year
E=9.87*sin(2*B)-7.53*cos(B)-1.5*sin(B); % (in minutes) is an empirical equation that corrects for the eccentricity of the Earth's orbit and the Earth's axial tilt.
tc=4*(lambda-Zc)+E;%Time Correction Factor
ts=LT+(tc/60);%Local Solar Time (LST) can be found by using the previous two corrections to adjust the local time (LT).
w=15*(ts-12);%Hour Angle converts the local solar time (w) into the number of degrees which the sun moves across the sky
segma=23.45*sin((284+n)/36.25); %The solar declination
Gsc=1.367;%solar constant [1.367 kW/m2],
Gon=Gsc*(1+0.033*cos(360*n/365));
for J=1:length(w)
H(w)=sind(w+1)-sind(w);
%B=diff(sind(1:365))
end
Azzi Abdelmalek
Azzi Abdelmalek el 15 de Jul. de 2016
Editada: Azzi Abdelmalek el 15 de Jul. de 2016
fi=19.9450; %is the latitude
GMT=2 ;% Time Zone for example +2
LT=5; % local time (LT).
lambda=50.067;% longitude
Zc=15*GMT;%Local Standard Time Meridian
for n=1:365 %the day of the year [a number 1 through 365]
B=(360/365)*(n-81)*pi/180;% in degrees and d is the number of days since the start of the year
E=9.87*sin(2*B)-7.53*cos(B)-1.5*sin(B); % (in minutes) is an empirical equation that corrects for the eccentricity of the Earth's orbit and the Earth's axial tilt.
tc=4*(lambda-Zc)+E;%Time Correction Factor
ts=LT+(tc/60);%Local Solar Time (LST) can be found by using the previous two corrections to adjust the local time (LT).
w(n)=15*(ts-12);%Hour Angle converts the local solar time (w) into the number of degrees which the sun moves across the sky
segma=23.45*sin((284+n)/36.25); %The solar declination
Gsc=1.367;%solar constant [1.367 kW/m2],
Gon=Gsc*(1+0.033*cos(360*n/365));
end
H=sind(w+1)-sind(w)
%B=diff(sind(1:365))

Iniciar sesión para comentar.

Categorías

Más información sobre Earth and Planetary Science en Centro de ayuda y File Exchange.

Etiquetas

Preguntada:

el 15 de Jul. de 2016

Editada:

el 15 de Jul. de 2016

Community Treasure Hunt

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

Start Hunting!

Translated by