%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%?micro-Doppler?signature?of?rotating?ALCM
%?
%?ALCM:?length?of?6.4?m?and?wingspan?of?3.4?m
%???????rotating?（yaw）:?0-359?degree;
%?Radar?Signal:?3?cm?wavelength?
%???????????????1?us?short?pulse?waveform?
%???????????????pulse?duration?0.0533?us
%?Pulse?repetition?frequency?（PRF）:?15000?Hz
%?ALCM?translational?velocity:?0
%?Number?of?range?samples:?100
%?Range?dimension:?20?m
%?Number?of?realization?（pulses）:?8192
%?Total?rotating?time?period:?0.55?s
%?Data?type:?complex?I&Q
%
%?By?V.C.Chen
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear?all

%?read?the?ALCM?file:

[filename]?=?‘RotatingALCM.dat‘;
%?read?data?header?and?data
[npt1t2t3nrt4]=textread（filename‘%d?%s?%s?%s?%d?%s‘1）;
[c1c2]=textread（filename‘%n%n%*[^\n]‘‘headerlines‘1‘delimiter‘‘‘）;
for?k=1:nr
????Cdata（k1:np）=c1（k+[0:np-1]*nr1）+j*c2（k+[0:np-1]*nr1）;
end

rngpro?=?fftshift（fft（Cdata）1）;

Cdata?=?Cdata/max（max（abs（Cdata）））;

%?figure（1）
%?colormap（jet）
%?imagesc（[1?np][1?nr]20*log10（abs（Cdata（:1:10:end）+eps）））;
%?grid?on;?box?on;?zoom?on;
%?axis?xy;
%?xlabel（‘Pulses‘）;
%?ylabel（‘Range?cells‘）;
%?title（‘Range?Profiles?（dBm^2）‘）
%?clim?=?get（gca‘CLim‘）;
%?set（gca‘CLim‘clim（2）?+?[-60?-10]）;
%?colorbar;
%?drawnow

%?Micro-Doppler?signature
PRF?=?15000;
PRI?=?1/PRF;
T?=?PRI*np;
F?=?1/PRI;

%?analyzed?signal
x?=?sum（Cdata）;

%?divide?long?data?into?a?number?of?shorter?（512）?data?segments
wd?=?512;
wdd2?=?wd/2;
wdd8?=?wd/8;
ns?=?np/wd;?%?total?number?of?segments

%?calculate?short-time?Fourier?transform?for?each?short?data?segment
disp（‘Calculating?STFT?for?segments?...‘）
for?k?=?1:ns
????disp（strcat（‘Progress:?segment?no.‘num2str（k）‘/‘...
????????num2str（round（ns））））
????sig（1:wd1）?=?x（1（k-1）*wd+1:（k-1）*wd+wd）;
????TMP?=?stft（sig24）;
????TF2（:（k-1）*wdd8+1:（k-1）*wdd8+wdd8）?=?TMP（:1:8:wd）;
end
TF?=?TF2;
disp（‘Calculating?STFT?for?time-shifted?segments?...‘）
TF1?=?zeros（size（TF））;
for?k?=?1:ns-1
????disp（strcat（‘Progress:??shifted?no.‘num2str（k）‘/‘...
????????num2str（round（ns-1））））
????sig（1:wd1）?=?x（1（k-1）*wd+1+wdd2:（k-1）*wd+wd+wdd2）;
????TMP?=?stft（sig24）;
????TF1（:（k-1）*wdd8+1:（k-1）*wdd8+wdd8）?=?TMP（:1:8:wd）;
end
disp（‘Removing?the?edge?effect?...‘）
for?k?=?1:ns-1
????TF（:k*wdd8-8:k*wdd8+8）?=?...
????????TF1（:（k-1）*wdd8+wdd8/2-8:（k-1）*wdd8+wdd8/2+8）;
end

%?display?final?time-frequency?micro-Doppler?signature
figure（2）
colormap（jet）
imagesc（[0T][-F/2F/2]20*log10（fftshift（abs（TF）1）+eps））
xlabel（‘Time?（s）‘）
ylabel（‘Doppler?（Hz）‘）
title（‘Micro-Doppler?Signature‘）
axis?xy
clim?=?get（gca‘CLim‘）;
set（gca‘CLim‘clim（2）?+?[-80?-5]）;
colorbar
drawnow

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-10-18?09:14??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\
?????文件????????1898??2010-09-28?01:08??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\readme.txt
?????目錄???????????0??2018-10-18?09:13??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第一章\
?????文件????????2878??2010-09-14?00:53??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第一章\ALCMSignature.m
?????文件????51274356??2005-08-25?05:37??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第一章\RotatingALCM.dat
?????文件????????1127??2010-09-14?00:47??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第一章\stft.m
?????目錄???????????0??2018-10-18?09:13??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\
?????目錄???????????0??2018-10-18?09:13??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\HeavyTop?Precession\
?????文件????????1508??2010-09-14?09:25??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\HeavyTop?Precession\circularplate.m
?????文件????????1455??2010-09-16?10:12??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\HeavyTop?Precession\cylinder3.m
?????文件?????????459??2010-09-14?06:38??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\HeavyTop?Precession\EulerTransfMatrix.m
?????文件????????7515??2014-05-07?14:36??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\HeavyTop?Precession\RadarTopReturns.m
?????文件????????1238??2010-09-15?04:42??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\HeavyTop?Precession\rcscircplate.m
?????文件????????2661??2010-09-15?04:43??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\HeavyTop?Precession\rcsfrustum.m
?????文件????????1127??2010-09-14?00:47??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\HeavyTop?Precession\stft.m
?????文件????????1829??2010-09-14?06:39??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\HeavyTop?Precession\TopDynamic.m
?????文件????????1138??2010-09-14?03:30??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\HeavyTop?Precession\XConvention.m
?????目錄???????????0??2018-10-18?09:13??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Helicopter?Rotor?Blades\
?????文件????????5319??2010-09-14?03:47??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Helicopter?Rotor?Blades\RadarRectBladeReturns.m
?????文件?????????783??2010-09-15?04:39??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Helicopter?Rotor?Blades\rcs_rect.m
?????文件????????3064??2010-09-14?03:50??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Helicopter?Rotor?Blades\RotorBladesTheory.m
?????文件????????1127??2010-09-14?00:47??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Helicopter?Rotor?Blades\stft.m
?????文件????????1138??2010-09-14?03:30??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Helicopter?Rotor?Blades\XConvention.m
?????目錄???????????0??2018-10-18?09:13??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Pendulum\
?????文件????????1356??2010-09-14?02:00??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Pendulum\Pendulum.m
?????文件????????1342??2010-09-14?02:12??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Pendulum\PendulumDamping.m
?????文件????????1361??2010-09-14?02:19??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Pendulum\PendulumDampingDriving.m
?????文件????????1144??2010-09-14?02:14??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Pendulum\PendulumDampingDrivingDynamic.m
?????文件????????1053??2010-09-14?02:14??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Pendulum\PendulumDampingDynamic.m
?????文件????????1266??2010-09-14?02:15??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Pendulum\PendulumDriving.m
?????文件????????1062??2010-09-14?02:15??Victior?c.chen《雷達中的微多普勒效應》英文原版代碼\第三章\Pendulum\PendulumDrivingDynamic.m
............此處省略27個文件信息