线性调频信号回波处理Word文档下载推荐.docx

1、 %带宽pri=1e-3; %脉冲重复周期fs=70e6; %采样率R=100e3; %距离v=0; %速度Pt=10*log10（1e6）; %功率lamda=10*log10（0.1）;%波长G=30; %天线增益dBL=6; %损耗dBFn=3; %噪声系数dBRCS=10*log10（2）; %RCS%*信号参数计算*%Pr=Pt+2*G+2*lamda+RCS-3*10*log10（4*pi）-4*10*log10（R）-L; %回波功率计算A=sqrt（10（Pr/10）; %回波信号幅度tou=2*R/3e8; %收发时间fd=2*v*fr/3e8; %回波多普勒频率k=bw/p

2、w; %调频带宽频率变化率%*发射信号产生*%t=（0:1/fs:pri）.; %时间sampleTbegin=ceil（0*fs+1）; %采样时间起始点sampleTEnd=ceil（pw*fs+1）; %采样时间完毕点sTrans=zeros（length（t）,1）; %产生发射信号sTrans（sampleTbegin:sampleTEnd,1）=. sin（2*pi*（f0*t（sampleTbegin:sampleTEnd,1）+k*t（sampleTbegin:sampleTEnd,1）.2/2）;%图像展示figure;plot（t*1e6,sTrans）; %发射信号时域波

3、形 xlabel（时间/us）;ylabel（幅度freTrans=fftshift（fft（sTrans）; %发射信号频域波形f=linspace（-fs/2,fs/2,length（t）;plot（f/1e6,abs（freTrans）;频率/MHz%*回波信号产生*%beginsample=ceil（tou*fs）; %采样起始点endsample=ceil（tou+pw）*fs）; %采样完毕点sEcho=zeros（length（t）,1）; %产生回波信号sEcho（beginsample:endsample,1）=. A*sin（2*pi*（f0+fd）*（t（beginsam

4、ple:endsample,1）-tou）+k*（t（beginsample:endsample,1）-tou）.2/2）;plot（t*1e6,sEcho）; %回波信号时域波形 freEcho=fftshift（fft（sEcho）;plot（f/1e6,abs（freEcho）; %回波信号频域波形%*叠加接收机噪声*%kn=1.38e-23; T=290;Bn=bw;NoiseDev=kn*T*fs*10（Fn/10）; %噪声功率，因后续要通过滤波器，%此处带宽用采样带宽，保证噪声功率密度noise=sqrt（NoiseDev）*randn（length（sEcho）,1）;sEch

5、o=sEcho+noise; %信号混叠噪声Wn=fir1（50,（f0-bw）/（fs/2） （f0+2*bw）/（fs/2）; %带通滤波器设计sEcho=filter（Wn,1,sEcho）; %带通滤波处理%*正交鉴相*%LocalI=sin（2*pi*f0*t）;LocalQ=cos（2*pi*f0*t）;sI=sEcho.*LocalI;sQ=sEcho.*LocalQ;filtWn=fir1（50,0.2）; %低通滤波器设计filsI=filter（filtWn,1,sI）; %低通滤波处理filsQ=filter（filtWn,1,sQ）;fils=filsI+sqrt（-1

6、）.*filsQ; %组成复信号plot（t*1e6,fils）; %零中频信号frefilsI=fftshift（fft（fils）;plot（f,abs（frefilsI）;%*匹配滤波*%tlocal=（0:pw）.slocal=exp（j*2*pi*（k*tlocal.2/2）;%窗是对称的，因此窗长为奇数if rem（length（slocal）,2）=0 slocal=slocal;0;endwrect=window（rectwin,length（slocal）;srect=slocal.*wrect;srectproc=ifft（conj（fft（srect,length（sEc

7、ho）.*fft（fils）;plot（t*3e8/2/1e3,20*log10（abs（srectproc）-max（20*log10（abs（srectproc）,b距离/Km幅度/dB%*加窗处理*%whamming=window（hamming,length（slocal）;shamming=slocal.*whamming;shammingproc=ifft（conj（fft（shamming,length（sEcho）.*fft（fils）;wblackman=window（blackman,length（slocal）;sblackman=slocal.*wblackman;sb

8、lackmanproc=ifft（conj（fft（sblackman,length（sEcho）.*fft（fils）;plot（wrect,hold on;plot（whamming,rplot（wblackman,k采样点plot（t*3e8/2/1e3,20*log10（abs（shammingproc）-max（20*log10（abs（shammingproc）,plot（t*3e8/2/1e3,20*log10（abs（sblackmanproc）-max（20*log10（abs（sblackmanproc）,Figure1为发射信号时域波形Figure2为发射信号幅频特性曲线

9、Figure3为回波信号时域波形图Figure4为回波信号幅频特性曲线Figure5为叠加高斯噪声后的回波信号时域波形Figure6为叠加噪声后的回波信号幅频特性曲线Figure7为接收信号下变频至零中频信号时的波形图Figure8为零中频信号的幅频特性曲线Figure9为回波信号在接收机处经匹配滤波并作归一化处理后的显示结果。蓝、红、黑线分别为不加窗、汉明窗、布莱克曼窗的窗形。蓝、红、黑线分别为匹配滤波后的信号不加窗、汉明窗、布莱克曼窗的结果。二、思考题因为在MATLAB软件仿真的过程当中，时间、频率、采样点是离散的，且函数ceil的特性是+1取整，因此影响到fd，进而影响到最后的距离，从而出现了仿真之后的距离比设定距离稍多的情况。