how can i plot a graph using the data
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
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 commentaires
dpb
le 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?
Réponses (0)
Voir également
Catégories
En savoir plus sur Time Series dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!