matlab 數組橫向縱向拼接

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];