multi-beam.m

一个多波束合成算法的程序

year=2004;%年份
day=24;%日期
hour=7;%钟点
minute=0;%分钟
s=hour+minute/60;%时间
n=31+29+31+30+31+30+day;%七月积日
if 2<=s<4.4%时刻修正
    n=n-0.2;
elseif 4.4<=s<6.8
    n=n-0.1;
elseif 9.2<=s<11.6
    n=n+0.1;
elseif 11.6<=s<14
    n=n+0.2;
elseif 14<=s<16.4
    n=n+0.3;
elseif 16.4<=s<18.8
    n=n+0.4;
elseif 18.8<=s<21.2
    n=n+0.5;
else
    n=n;
end;
n=n+0.2734-0.3;%年度和经度修正
no=79.6764+0.2422*(year-1985)-fix((year-1985)/4);%标准日
dayd=n-no;%日差
sa=2*pi*dayd/365.2422;%日角
ed=0.3723+23.2567*sin(sa)+0.1149*sin(2*sa)-0.1712*sin(3*sa)-0.758*cos(sa)+0.3656*cos(2*sa)+0.0201*cos(3*sa);%赤纬角
et=0.0028-1.9857*sin(sa)+9.9059*sin(2*sa)-7.0924*cos(sa)-0.6882*cos(2*sa);%时差
long=117+14.150/60;%合肥经度
lat=31+50.718/60;%合肥纬度
sd=s+(long-120)/15;%合肥地方时
so=sd+et/60;%真太阳时
tau=(so-12)*15;%太阳时角
%tau=15*s+long-300;%太阳时角
h=asin(sin(pi*ed/180)*sin(pi*lat/180)+cos(pi*ed/180)*cos(pi*lat/180)*cos(pi*tau/180))*180/pi;%太阳高度角
h1=(h-fix(h))*60;%分
if so>=12
    a=acos((sin(pi*h/180)*sin(pi*lat/180)-sin(ed*pi/180))/(cos(h*pi/180)*cos(pi*lat/180)))*180/pi;%真太阳时12点以后太阳方位角
else
    a=360-acos((sin(pi*h/180)*sin(pi*lat/180)-sin(ed*pi/180))/(cos(h*pi/180)*cos(pi*lat/180)))*180/pi;%真太阳时12点以前太阳方位角
end 
%a=asin(cos(pi*ed/180)*sin(tau*pi/180)/cos(h*pi/180))*180/pi;%太阳方位角
a1=(a-fix(a))*60;%分
h
h1
a
a1