随机信号上机作业1分析.doc

随机信号分析上机作业 姓名：王霞 学号：111040709 学院：电光学院 指导老师：顾红 一：仿真伪相位编码连续波雷达的信号处理。设码频为 5MHz，伪码周期内码长为 511，雷达载频为10GHz，输入噪声为高斯白噪声，视频输入信噪比为-15dB，相干积累总时宽不大于 10ms，给出回波视频表达式，脉压和 FFT后的表达式；仿真给出脉压和 FFT 后的输出图形；通过仿真说明各级处理的增益，与各级时宽和带宽的关系；仿真说明脉压时多卜勒敏感现象和多卜勒容限及其性能损失（脉压主旁比与多卜勒的曲线） 回波视频表达式： 脉压后的表达式： FFT后的表达式： 仿真： 脉压后的输出图形 脉压增益为：G1=511 与时宽带宽的关系：脉压增益=时宽带宽积 FFT后的输出图形 FFT的增益为G2=N/2=97/2=48.5 脉压多普勒敏感现象 ，脉压输出 d）小多普勒时，脉压输出 e）脉压长度为多普勒整周期时，脉压输出 由图可以看出：小多普勒的时候，脉压比小于511，脉压长度为多普勒整周期时，脉压输出基本为0。 多普勒容限 脉压长度刚好为多普勒信号一个周期一半时 , 脉压长度刚好为多普勒信号一个周期一半时，脉压输出损失2.4db。 改变多普勒频率，记录主旁瓣比，如表一所示： 表一 Fd/KHz 0 1 2 3 4 6 7 8 9 10 主旁瓣比/db 52 39 34 29 4 24 20 17 10 5 g）主旁瓣比与多普勒频率关系 附MATLAB程序： clc; clear; fm=5e6;%码频 fs=30e6; fc=10e9;%载频 t=10e-3;%相干积累总时宽 T=1/fm*511;%周期 N=97; %积累时间不大于10ms，Nt/T=97.8474,所以N取小于97的整数 snr=-15;%信噪比 %----------------511位m序列产生---------------------% fbconnection=[0 0 0 1 0 0 0 0 1]; n = length(fbconnection); m = 2^n-1; register = [1 0 1 0 1 0 1 0 1];%定义移位寄存器的初始状态 mseqmatrix(1)= register(n); for i = 2:m newregister(1)= mod(sum(fbconnection.*register),2); for j = 2:n newregister(j)= register(j-1); end register = newregister; mseqmatrix(i) = register(n); end m_511=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 r=1:N for h=1:m s_hb1(h,r)=hbb1((r-1)*m+h); end end figure;mesh(1:N,1:511,s_hb1); % %------------------------------------------------ for h=1:m r_fft(h,:)=abs(fft(s_hb1(h,:))); end figure; mesh(1:N,1:511,r_fft); %------------------------------------------------- % for a=1:10 % fd=a*1e3; fd=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); fi