taiyingri_l.m

来自「平时练习时的一些仿真程序」· M 代码 · 共 39 行

H=AnyChooseY([93 11 5 00 00],[93 11 30 23 00]);
L=length(H);
%28*24*60;
sq=zeros(1440,1);
for n=1:1440:L-1439
    sq=sq+H(n:n+1439);
end
ndate=L/1440;
sq=sq/ndate;

figure(2)
subplot(311)
plot(H) 
% grid on 
% set(gca,'XMinorGrid','on','YMinorGrid','on')
% hold on

subplot(312)
plot(sq)     %月均值曲线

H2=H-repmat(sq,ndate,1);
ML=zeros(1490,1);
NL=0;
for n=1:1490:L-1489
    ML=ML+H2(n:n+1489);
    NL=NL+1;
end
ML=ML/NL;
subplot(313)
plot(ML)     %月平均太阴日变化曲线L