線性調(diào)頻信號(hào)回波管理

.-%程序：線性調(diào)頻信號(hào)回波處理%LFMEchoProc.mclc;clear;closeall;

%清楚命令窗口信息%清除工作區(qū)的變量%清除繪圖%****************參數(shù)設(shè)置*************%感謝閱讀fr=3e9; %射頻f0=10e6; %中頻pw=256e-6; %脈寬bw=2e6; %帶寬pri=1e-3; %脈沖重復(fù)周期fs=70e6; %采樣率R=100e3; %距離v=0; %速度Pt=10*log10(1e6); %功率lamda=10*log10(0.1);%波長G=30; %天線增益dBL=6; %損耗dBFn=3; %噪聲系數(shù)dBRCS=10*log10(2); %RCS%***************信號(hào)參數(shù)計(jì)算*************%感謝閱讀Pr=Pt+2*G+2*lamda+RCS-3*10*log10(4*pi)-4*10*log10(R)-L;%回波功率計(jì)算謝謝閱讀A=sqrt(10^(Pr/10));%回波信號(hào)幅度tou=2*R/3e8;%收發(fā)時(shí)間fd=2*v*fr/3e8;%回波多普勒頻率k=bw/pw;%調(diào)頻帶寬頻率變感謝閱讀化率%************************發(fā)射信號(hào)產(chǎn)生*******************%感謝閱讀t=(0:1/fs:pri).';sampleTbegin=ceil(0*fs+1);謝謝閱讀sampleTEnd=ceil(pw*fs+1);sTrans=zeros(length(t),1);謝謝閱讀sTrans(sampleTbegin:sampleTEnd,1)=...感謝閱讀

%時(shí)間%采樣時(shí)間起始點(diǎn)%采樣時(shí)間結(jié)束點(diǎn)%產(chǎn)生發(fā)射信號(hào).-sin(2*pi*(f0*t(sampleTbegin:sampleTEnd,1)+k*t(sampleTbegin:sampleTEnd精品文檔放心下載,1).^2/2));%圖像展示figure;plot(t*1e6,sTrans); %發(fā)射信號(hào)時(shí)域波精品文檔放心下載形xlabel('時(shí)間/us');ylabel('幅度');freTrans=fftshift(fft(sTrans)); %發(fā)射信號(hào)頻域波謝謝閱讀形f=linspace(-fs/2,fs/2,length(t));精品文檔放心下載figure;plot(f/1e6,abs(freTrans));精品文檔放心下載xlabel('頻率/MHz');ylabel('幅度');%***********************回波信號(hào)產(chǎn)生*******************%精品文檔放心下載beginsample=ceil(tou*fs); %采樣起始點(diǎn)感謝閱讀endsample=ceil((tou+pw)*fs); %采樣結(jié)束點(diǎn)感謝閱讀sEcho=zeros(length(t),1); %產(chǎn)生回波信號(hào)感謝閱讀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); %回波信號(hào)時(shí)域波精品文檔放心下載形xlabel('時(shí)間/us');ylabel('幅度');freEcho=fftshift(fft(sEcho));精品文檔放心下載f=linspace(-fs/2,fs/2,length(t));精品文檔放心下載figure;plot(f/1e6,abs(freEcho)); %回波信號(hào)頻域波謝謝閱讀形xlabel('頻率/MHz');ylabel('幅度');%**********************疊加接收機(jī)噪聲*******************%感謝閱讀kn=1.38e-23;.-T=290;Bn=bw;NoiseDev=kn*T*fs*10^(Fn/10); %噪聲功率，因后續(xù)要通過濾波器，謝謝閱讀%此處帶寬用采樣帶寬，保證噪聲功率密度noise=sqrt(NoiseDev)*randn(length(sEcho),1);精品文檔放心下載sEcho=sEcho+noise;%信號(hào)混疊噪聲Wn=fir1(50,[(f0-bw)/(fs/2)(f0+2*bw)/(fs/2)]);%帶通濾波器設(shè)計(jì)sEcho=filter(Wn,1,sEcho);%帶通濾波處理謝謝閱讀figure;plot(t*1e6,sEcho);xlabel('時(shí)間/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); %低通濾波器設(shè)計(jì)精品文檔放心下載filsI=filter(filtWn,1,sI); %低通濾波處理精品文檔放心下載filsQ=filter(filtWn,1,sQ); %低通濾波處理謝謝閱讀fils=filsI+sqrt(-1).*filsQ; %組成復(fù)信號(hào)感謝閱讀figure;plot(t*1e6,fils); %零中頻信號(hào)xlabel('時(shí)間/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));謝謝閱讀%窗是對稱的，因此窗長為奇數(shù)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('采樣點(diǎn)');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為發(fā)射信號(hào)時(shí)域波形.-Figure2為發(fā)射信號(hào)幅頻特性曲線.-Figure3為回波信號(hào)時(shí)域波形圖.-Figure4為回波信號(hào)幅頻特性曲線.-Figure5為疊加高斯噪聲后的回波信號(hào)時(shí)域波形.-Figure6為疊加噪聲后的回波信號(hào)幅頻特性曲線.-Figure7為接收信號(hào)下變頻至零中頻信號(hào)時(shí)的波形圖感謝閱讀.-Figure8為零中頻信號(hào)的幅頻特性曲線.-Figure9為回波信號(hào)