[求助]IDL 下 雷达地面信号回波的带通滤波

我现在在做一个小课题（用IDL Matlab太慢了） 遇到一点麻烦 我是第一次用Idl 不清楚怎么做滤波就是在下面的程序中 模拟的地面回波信号 老师要求再仿真雷达的扫描过程 以及将同时落入同一波束的两个点目标的前一个滤出来，点目标的坐标都是自己放的，只要能够有滤波功能就行了 程序附在下面俯仰角15。； 锐化比60；波束扫描范围72O 108。；主波束方位向宽3．5。；波束扫描速度l8。／s； 发射信号带宽80 MHz；脉冲重复频率2000 Hz；脉冲宽度40 s ;参数设置c=double(3e8)lamd=double(0.03)B=double(80.e6)fs=double(100.e6)Tp=double(10.e-6)va=double(200.)H=double(6000.)gama_a=double(3.*!pi/180.)OC=double(40000.);APC到成像区中心距离short_R=double(sqrt((OC^2-H^2)));航路捷径;基本关系rou_r=double(c/2./B);15mpixel_r=double(c/2./fs);7.5k=B/Tp;调频斜率PRF=double(2000.)M=65.;相干处理点数Tr_points=floor(fs*Tp) 1.Ta_points=Mt_r=(findgen(Tr_points)-ceil(Tr_points/2))/fst_a=(findgen(M)-floor(M/2))/PRFazimuth_samples=Ta_pointsrange_samples=double(2048.)range_center=ceil(range_samples/2)range_start=range_center-ceil(tr_points/2);存放中心点回波数据的起点;几何坐标设置APC_start_x=double(0.)APC_start_y=double(0.)APC_start_z=double(H)APC_x=double(0.)APC_y=va*t_aAPC_z=double(H)dot_x=[short_R,short_R 50,short_R-50,short_R,short_R]dot_y=[0.,0.,0.,500.,-500.]dot_z=[0.,0.,0.,0.,0.];回波数据储存 ；【应该就是在这下面加 每64个处理点做一次滤波 然后储存 但是我是第一次用Idl 不清楚怎么做滤波】echo=complexarr(azimuth_samples,range_samples)echo_tmp=complexarr(Tr_points)Rt=sqrt((APC_x-dot_x)^2 (APC_y-dot_y)^2 (APC_z-dot_z)^2)for ii=0,azimuth_samples-1 do beginfor n=0,4 do begin Rt=sqrt((APC_x-dot_x(n))^2 (APC_y(ii)-dot_y(n))^2 (APC_z-dot_z(n))^2) ref_t=2*OC/c delay=2*Rt/c delay_relative=ceil((delay-ref_t)*fs); echo_tmp=exp(complex(0,k*!pi*t_r^2))*exp(complex(0,-4.*!pi/lamd*Rt)) echo(ii,range_start delay_relative:range_start delay_relative tr_points-1)=echo(ii,range_start delay_relative:range_start delay_relative tr_points-1) echo_tmp echo_tmp=complexarr(Tr_points)endforendforwindow,1contour,echo;距离向脉冲压缩;match_t=exp(complex(0,-!pi*k*t_r^2));match_f=fft(match_t,range_samples)fs_in_range=-fs/2 (findgen(range_samples))*fs/range_samplesmatch_f=exp(complex(0,fs_in_range^2.*!pi/k))edge_in_range=ceil((fs-B)/2./fs*range_samples)match_f(0:edge_in_range-1.)=0.match_f(range_samples-edge_in_range:range_samples-1.)=0.;match_tmp=complexarr(range_samples);match_t=exp(complex(0,k*!pi*t_r^2));match_tmp(0:tr_points-1)=match_t .. Oy(fh%k#  
&m{SWV+  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册