how can i plot a graph using the data
2 views (last 30 days)
Show older comments
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 Comments
dpb
on 22 May 2014
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?
Answers (0)
See Also
Categories
Find more on Time Series in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!