Victor C.Chen 《雷达中的微多普勒效应》英文原版代码

国内电子工业那版书里面的附录代码删除了部分，很多缺函数无法运行，英文原版代码很全

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% 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个文件信息