how can i plot a graph using the data
2 次查看(过去 30 天)
显示 更早的评论
this is a raw data file
clear all;
close all;
clc;
tic;
fp1=fopen('09JL2008.r2','r');
Startheader= struct('npar',fread(fp1,1,'int16'),...
'nwds',fread(fp1,1,'int16'),...
'nrgb',fread(fp1,1,'int16'),...
'nfft',fread(fp1,1,'int16'),...
'ncoh',fread(fp1,1,'int16'),...
'nicoh',fread(fp1,1,'int16'),...
'ipp',fread(fp1,1,'int16'),...
'pwd',fread(fp1,1,'int16'),...
'cflag',fread(fp1,1,'int16'),...
'nwin',fread(fp1,1,'int16'),...
'w1start',fread(fp1,1,'int16'),...
'w1len',fread(fp1,1,'int16'),...
'w2start',fread(fp1,1,'int16'),...
'w2len',fread(fp1,1,'int16'),...
'year',fread(fp1,1,'int16'),...
'month',fread(fp1,1,'int16'),...
'day',fread(fp1,1,'int16'),...
'hour',fread(fp1,1,'int16'),...
'min',fread(fp1,1,'int16'),...
'sec',fread(fp1,1,'int16'),...
'nbeam',fread(fp1,1,'int16'),...
'beam',fread(fp1,1,'int16'),...
'scan',fread(fp1,1,'int16'));
timestart=sprintf('%d:%d:%d.',Startheader.hour,Startheader.min,Startheader.sec);
secstart=(Startheader.hour*3600+Startheader.min*60+Startheader.sec);
fseek(fp1,-(8*Startheader.nrgb*Startheader.nfft+128),1);
%for be=1:5
header= struct('npar',fread(fp1,1,'int16'),...
'nwds',fread(fp1,1,'int16'),...
'nrgb',fread(fp1,1,'int16'),...
'nfft',fread(fp1,1,'int16'),...
'ncoh',fread(fp1,1,'int16'),...
'nicoh',fread(fp1,1,'int16'),...
'ipp',fread(fp1,1,'int16'),...
'pwd',fread(fp1,1,'int16'),...
'cflag',fread(fp1,1,'int16'),...
'nwin',fread(fp1,1,'int16'),...
'w1start',fread(fp1,1,'int16'),...
'w1len',fread(fp1,1,'int16'),...
'w2start',fread(fp1,1,'int16'),...
'w2len',fread(fp1,1,'int16'),...
'year',fread(fp1,1,'int16'),...
'month',fread(fp1,1,'int16'),...
'day',fread(fp1,1,'int16'),...
'hour',fread(fp1,1,'int16'),...
'min',fread(fp1,1,'int16'),...
'sec',fread(fp1,1,'int16'),...
'nbeam',fread(fp1,1,'int16'),...
'beam',fread(fp1,1,'int16'),...
'scan',fread(fp1,1,'int16'));
nbeam=header.nbeam;
nrgb=header.nrgb;
nfft=header.nfft;
ncoh=header.ncoh;
nicoh=header.nicoh;
ipp=header.ipp;
pwd=header.pwd;
scan=header.scan;
header.day;
header.month;
header.year;
header.hour;
header.min;
header.sec;
beamdate = sprintf('%2d-%d-%4d.', header.day,header.month,header.year);
timeend = sprintf('%d:%d:%d.',header.hour,header.min,header.sec);
secend=(header.hour*3600+header.min*60+header.sec);
%%%data read
nrgb=150;nfft=512;
for rb=1:nrgb
for nft=1:nfft
datI(rb,nft)=fread(fp1,1,'float32');
datQ(rb,nft)=fread(fp1,1,'float32');
end
end
how to load time with the above information and how i can plot a graph in time series using datI and datQ separately i want a datI in timeseries separately and datQ separately with respect to timeseries
2 个评论
dpb
2014-5-22
Seems peculiar request -- aren't the dat arrays a fft at a time, not timeseries data? Or are you meaning something like a waterfall plot of the spectra or somesuch?
What is the time step in the series?
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Time Series 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!