线性调频信号回波管理

.-%程序：线性调频信号回波处理%LFMEchoProc.mclc;clear;closeall;

%清楚命令窗口信息%清除工作区的变量%清除绘图%****************参数设置*************%感谢阅读fr=3e9; %射频f0=10e6; %中频pw=256e-6; %脉宽bw=2e6; %带宽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/pw;%调频带宽频率变感谢阅读化率%************************发射信号产生*******************%感谢阅读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); %发射信号时域波精品文档放心下载形xlabel('时间/us');ylabel('幅度');freTrans=fftshift(fft(sTrans)); %发射信号频域波谢谢阅读形f=linspace(-fs/2,fs/2,length(t));精品文档放心下载figure;plot(f/1e6,abs(freTrans));精品文档放心下载xlabel('频率/MHz');ylabel('幅度');%***********************回波信号产生*******************%精品文档放心下载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(beginsample:endsample,1)-tou)+k*(t(beginsample:endsample,1)-tou).^2/2));谢谢阅读figure;plot(t*1e6,sEcho); %回波信号时域波精品文档放心下载形xlabel('时间/us');ylabel('幅度');freEcho=fftshift(fft(sEcho));精品文档放心下载f=linspace(-fs/2,fs/2,length(t));精品文档放心下载figure;plot(f/1e6,abs(freEcho)); %回波信号频域波谢谢阅读形xlabel('频率/MHz');ylabel('幅度');%**********************叠加接收机噪声*******************%感谢阅读kn=1.38e-23;.-T=290;Bn=bw;NoiseDev=kn*T*fs*10^(Fn/10); %噪声功率，因后续要通过滤波器，谢谢阅读%此处带宽用采样带宽，保证噪声功率密度noise=sqrt(NoiseDev)*randn(length(sEcho),1);精品文档放心下载sEcho=sEcho+noise;%信号混叠噪声Wn=fir1(50,[(f0-bw)/(fs/2)(f0+2*bw)/(fs/2)]);%带通滤波器设计sEcho=filter(Wn,1,sEcho);%带通滤波处理谢谢阅读figure;plot(t*1e6,sEcho);xlabel('时间/us');ylabel('幅度');freEcho=fftshift(fft(sEcho));谢谢阅读f=linspace(-fs/2,fs/2,length(t));谢谢阅读figure;plot(f/1e6,abs(freEcho));xlabel('频率/MHz');ylabel('幅度');%*******************正交鉴相****************%感谢阅读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).*filsQ; %组成复信号感谢阅读figure;plot(t*1e6,fils); %零中频信号xlabel('时间/us');ylabel('幅度');frefilsI=fftshift(fft(fils));感谢阅读f=linspace(-fs/2,fs/2,length(t));感谢阅读figure;plot(f,abs(frefilsI));xlabel('频率/MHz');ylabel('幅度');%************匹配滤波**************%感谢阅读.-tlocal=(0:1/fs:pw).';slocal=exp(j*2*pi*(k*tlocal.^2/2));谢谢阅读%窗是对称的，因此窗长为奇数ifrem(length(slocal),2)==0感谢阅读slocal=[slocal;0];endwrect=window(@rectwin,length(slocal));感谢阅读srect=slocal.*wrect;srectproc=ifft(conj(fft(srect,length(sEcho))).*fft(fils));谢谢阅读figure;plot(t*3e8/2/1e3,20*log10(abs(srectproc))-max(20*log10(abs(srectproc)精品文档放心下载)),'b');xlabel('距离/Km');ylabel('幅度/dB');%****************加窗处理***************%感谢阅读wrect=window(@rectwin,length(slocal));精品文档放心下载srect=slocal.*wrect;srectproc=ifft(conj(fft(srect,length(sEcho))).*fft(fils));精品文档放心下载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;感谢阅读sblackmanproc=ifft(conj(fft(sblackman,length(sEcho))).*fft(fils));谢谢阅读figure;plot(wrect,'b');holdon;plot(whamming,'r');holdon;plot(wblackman,'k');xlabel('采样点');ylabel('幅度');figure;plot(t*3e8/2/1e3,20*log10(abs(srectproc))-max(20*log10(abs(srectproc)精品文档放心下载)),'b');holdon;plot(t*3e8/2/1e3,20*log10(abs(shammingproc))-max(20*log10(abs(shammin谢谢阅读gproc))),'r');holdon;.-plot(t*3e8/2/1e3,20*log10(abs(sblackmanproc))max(20*log10(abs(sblack谢谢阅读manproc))),'k')xlabel('距离/Km')ylabel('幅度/dB')Figure1为发射信号时域波形.-Figure2为发射信号幅频特性曲线.-Figure3为回波信号时域波形图.-Figure4为回波信号幅频特性曲线.-Figure5为叠加高斯噪声后的回波信号时域波形.-Figure6为叠加噪声后的回波信号幅频特性曲线.-Figure7为接收信号下变频至零中频信号时的波形图感谢阅读.-Figure8为零中频信号的幅频特性曲线.-Figure9为回波信号