matlab 数组横向纵向拼接

最新推荐文章于 2024-11-23 23:29:46 发布

蜡木浦

最新推荐文章于 2024-11-23 23:29:46 发布

阅读量2.2k

点赞数

CC 4.0 BY-SA版权

分类专栏： matlab

本文链接：https://blog.youkuaiyun.com/luoluo19550418/article/details/99934962

matlab 专栏收录该内容

1 篇文章

订阅专栏

本文介绍了一个使用MATLAB进行月相计算和观测数据处理的程序。该程序读取特定目录下所有.dat文件的名称，从中提取观测时间、方位角和仰角等信息，并计算相对于给定时刻的月相角。此外，程序还整合了从.asc文件中读取的数据，最终形成了包含观测时间、月相角、方位角、仰角和USB数据的完整数据集。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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