直升机旋翼微动效应仿真

基于微多普勒效应给出的关于直升机旋翼产生微动效应的仿真结果，对雷达目标识别和信号处理具有实际的参考意义。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 
% Radar returns from rectangular rotor blades
%
% Radar: C-band wavelength = 0.06 m
%        Range resolution 0.5 m
%        Observation time: T = 1 sec
%        Number of pulses 10240
%        Location: X = 500 m; Y = 0 m; Z = 500 m
% Rotor: Center location X = 0; Y = 0; Z = 0 m
%        Blade lengh L = 6 m （L1:0.5 - L2:6.5） wide W = 1 m
%        Rotation rate: 4 r/s
%
% By V.C. Chen
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all

c = 2.99792458e8;
j = sqrt（-1）;

% radar paramerters
radarloc = [5000500];  % radar location
xr = radarloc（1）;
yr = radarloc（2）;
zr = radarloc（3）;
rangeres = 0.5;  % （m）: designed range resolution
% total number of range bins
nr = floor（2*norm（radarloc）/rangeres）;

% time smapling
T = 1; % time duration
nt = 10240; % number of time samples
ts = linspace（0Tnt）; % time span

f0 = 5e9;
lambda = c/f0;

% rotor blades
Nb = 2; % number of blades （2 or 3）
L = 6; % blade length （m）
W = 1; % blade wide （m）
L1 = 0.5; % one end （root） of blade （m）
L2 = 6.5; % other end （tip） of blade （m）
a = （L2-L1）/2; % rectangular blade parameter
b = W/2; % rectangular blade parameter

rotorloc = [000]; % rotor center location
% direction of radar with respect to rotor center
radar_direction = atan2（radarloc（2）-rotorloc（2）radarloc（1）-rotorloc（1））;
Omega = 4*2*pi; %blade otation rate

% rotor blades‘ rotation center
x0 = 0;
y0 = 0;
z0 = 0;
rotor_center = [x0;y0;z0];

for k=1:nt
   psi（k） = 0;
   theta（k） = 0;
   phi（k） = Omega*ts（k）;
end

data = zeros（nrnt）;
for k = 1:nt
    % blade 1
    Rzxz = XConvention（phi（k）theta（k）psi（k））;  
    CM1 = Rzxz*[6.5;0;0];
    xcm1 = CM1（1）;
    ycm1 = CM1（2）;
    zcm1 = CM1（3）;
    % distance
    distance1 = sqrt（（xr-xcm1）^2+（yr-ycm1）^2+（zr-zcm1）^2）;
    % theta and phi angles
    Theta1 = atan2（sqrt（（xr-xcm1）^2+（yr-ycm1）^2）zr-zcm1）;
    Phi1 = -atan2（yr-ycm1xr-xcm1）;
    % rectangular blade shape
    rcs1 = rcs_rect（abTheta1 Phi1f0）;       
    amp1 = sqrt（rcs1）;
    PHs1 = amp1*（exp（j*4*pi*distance1/lambda））;
    data（floor（distance1/rangeres）k） = ...
        data（floor（distance1/rangeres）k）+ PHs1;           
         
    % blade 2
    Rzxz = XConvention（phi（k）+2*pi/Nbtheta（k）psi（k））;  
    CM2 = Rzxz*[6.5;0;0];
    xcm2 = CM2（1）;
    ycm2 = CM2（2）;
    zcm2 = CM2（3）;    
     % distance
    distance2 = sqrt（（xr-xcm2）^2+（yr-ycm2）^2+（zr-zcm2）^2）;
    % theta and phi angles
    Theta2 = atan2（sqrt（（xr-xcm2）^2+（yr-ycm2）^2）zr-zcm2）;
    Phi2 = -atan2（yr-ycm2xr-xcm2）;
    % rectangular balde shape
    rcs2 = rcs_rect（abTheta2 Phi2f0）;       
    amp2 = sqrt（rcs2）;
    PHs2 = amp2*（exp（j*4*pi*distance2/lambda））;
    data（floor（distance2/rangeres）k） = ...
        data（floor（distance2/rangeres）k）+ PHs2;           
         
    if Nb == 3
        % blade 3
        Rzxz = X

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6446  2018-01-09 17:38  程序\RadarRectBladeReturns.m

     文件        783  2010-09-15 04:39  程序\rcs_rect.m

     文件       3064  2018-01-22 21:36  程序\RotorBladesTheory.m

     文件       1127  2010-09-14 00:47  程序\stft.m

     文件       1138  2010-09-14 03:30  程序\XConvention.m

     目录          0  2018-11-24 16:27  程序

----------- ---------  ---------- -----  ----

                12558                    6