clear all;
clc;
ri0=1;shi0=11;fen0=12;
shuo=ri0+shi0/24+fen0/(24*60);
filepath='E:\夏季\观测\效率\';
dat=[filepath '*.dat']; dirfile=dir(dat); %包括 name date bytes isdir datenum
ri=[],shi=[],fen=[],az=[],el=[]; %观测时间
for i=1:length(dirfile),
ri=[ri;dirfile(i).name(12:13)];shi=[shi;dirfile(i).name(14:15)];fen=[fen;dirfile(i).name(16:17)];az=[az;dirfile(i).name(21:23)];el=[el;dirfile(i).name(25:26)];
end
ri=str2num(ri);shi=str2num(shi);fen=str2num(fen);az=str2num(az);el=str2num(el);
phase=[ri,shi,fen]; %横向拼接
period=29.530589;day=[];deltaday=[];chi=[];
for i=1:length(phase(:,1)),
daytmp=phase(i,1)+phase(i,2)/24+phase(i,3)/(24*60);
deltadaytmp=daytmp-shuo;
chitmp=(deltadaytmp/period)*360-180; %月相
day=[day;daytmp];deltaday=[deltaday;deltadaytmp],chi=[chi;chitmp]; %纵向拼接
end
phase=[phase,day,deltaday,chi];
beam1=[az,el];
ascfile=[filepath 'filepar.asc'];
data=load(ascfile);
usb1=data(:,1);
beam1=[beam1,usb1(7:12),usb1(1:6),chi];