欢迎来到天天文库
浏览记录
ID:33030408
大小:344.50 KB
页数:14页
时间:2019-02-19
《随机信号上机作业》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、随机信号分析上机作业姓名:王霞学号:111040709学院:电光学院指导老师:顾红一:仿真伪相位编码连续波雷达的信号处理。设码频为5MHz,伪码周期内码长为511,雷达载频为10GHz,输入噪声为高斯白噪声,视频输入信噪比为-15dB,相干积累总时宽不大于10ms,给出回波视频表达式,脉压和FFT后的表达式;仿真给出脉压和FFT后的输出图形;通过仿真说明各级处理的增益,与各级时宽和带宽的关系;仿真说明脉压时多卜勒敏感现象和多卜勒容限及其性能损失(脉压主旁比与多卜勒的曲线)回波视频表达式:脉压后的表达式:FFT后的表达式:仿真:a)
2、脉压后的输出图形(有噪声)脉压增益为:G1=511与时宽带宽的关系:脉压增益=时宽带宽积b)FFT后的输出图形(有噪声)FFT的增益为:G2=N/2=97/2=48.5脉压多普勒敏感现象c)无多普勒时,脉压输出图形d)小多普勒时,脉压输出图形e)脉压长度为多普勒整周期时,脉压输出图形由图可以看出:无多普勒时,脉压比为511,小多普勒的时候,脉压比小于511,脉压长度为多普勒整周期时,脉压输出基本为0。多普勒容限f)脉压长度刚好为多普勒信号一个周期一半时,脉压输出图形,脉压长度刚好为多普勒信号一个周期一半时,脉压输出损失2.94db
3、。性能损失改变多普勒频率,记录主旁瓣比,如表一所示:表一Fd/KHz012345678910主旁瓣比/db80523934294242017105g)主旁瓣比与多普勒频率关系附MATLAB程序:clc;clear;fm=5e6;%码频fs=30e6;fc=10e9;%载频t=10e-3;%相干积累总时宽T=1/fm*511;%周期N=97;%积累时间不大于10ms,N4、%fbconnection=[000100001];n=length(fbconnection);m=2^n-1;register=[101010101];%定义移位寄存器的初始状态mseqmatrix(1)=register(n);fori=2:mnewregister(1)=mod(sum(fbconnection.*register),2);forj=2:nnewregister(j)=register(j-1);endregister=newregister;mseqmatrix(i)=register(n);endm_55、11=2*mseqmatrix-1;%------------------0延时,0多普率的回波--------------------------%hb=repmat(m_511,1,N);hb=awgn(hb,snr,'measured');hb1=hb;%------------------脉压————————————%pipei=fliplr(m_511);hbb1=conv(pipei,hb1);figure;plot(hbb1);%-----------------距离门重排------------------%for6、r=1:Nforh=1:ms_hb1(h,r)=hbb1((r-1)*m+h);endendfigure;mesh(1:N,1:511,s_hb1);%%------------------------------------------------forh=1:mr_fft(h,:)=abs(fft(s_hb1(h,:)));endfigure;mesh(1:N,1:511,r_fft);%-------------------------------------------------%fora=1:10%fd=a*1e3;f7、d=5e6/511/4;i=1:511*N;duopule=cos(2*pi*fd/fm*i);hb2=hb.*duopule;%------------------脉压————————————%pipei=fliplr(m_511);hbb2=conv(pipei,hb2);m_hb1=abs(hbb2);m_hb=abs(hbb2)/max(m_hb1);m_hb=10*log(m_hb);figure;plot(m_hb);%-----------------距离门重排------------------%forr=1:Nf8、orh=1:ms_hb2(h,r)=hbb2((r-1)*m+h);endendfigure;mesh(1:N,1:511,s_hb2);%-------------------------------------------forh=1:mr_f
4、%fbconnection=[000100001];n=length(fbconnection);m=2^n-1;register=[101010101];%定义移位寄存器的初始状态mseqmatrix(1)=register(n);fori=2:mnewregister(1)=mod(sum(fbconnection.*register),2);forj=2:nnewregister(j)=register(j-1);endregister=newregister;mseqmatrix(i)=register(n);endm_5
5、11=2*mseqmatrix-1;%------------------0延时,0多普率的回波--------------------------%hb=repmat(m_511,1,N);hb=awgn(hb,snr,'measured');hb1=hb;%------------------脉压————————————%pipei=fliplr(m_511);hbb1=conv(pipei,hb1);figure;plot(hbb1);%-----------------距离门重排------------------%for
6、r=1:Nforh=1:ms_hb1(h,r)=hbb1((r-1)*m+h);endendfigure;mesh(1:N,1:511,s_hb1);%%------------------------------------------------forh=1:mr_fft(h,:)=abs(fft(s_hb1(h,:)));endfigure;mesh(1:N,1:511,r_fft);%-------------------------------------------------%fora=1:10%fd=a*1e3;f
7、d=5e6/511/4;i=1:511*N;duopule=cos(2*pi*fd/fm*i);hb2=hb.*duopule;%------------------脉压————————————%pipei=fliplr(m_511);hbb2=conv(pipei,hb2);m_hb1=abs(hbb2);m_hb=abs(hbb2)/max(m_hb1);m_hb=10*log(m_hb);figure;plot(m_hb);%-----------------距离门重排------------------%forr=1:Nf
8、orh=1:ms_hb2(h,r)=hbb2((r-1)*m+h);endendfigure;mesh(1:N,1:511,s_hb2);%-------------------------------------------forh=1:mr_f
此文档下载收益归作者所有