%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?
%?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