how can i plot a graph using the data
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
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 Kommentare
dpb
am 22 Mai 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?
Antworten (0)
Siehe auch
Kategorien
Mehr zu Time Series finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!