随机信号上机作业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，N t/T 97.8474,所以N取小于97的整数 snr -15;%信噪比 % 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; % 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 ; figure;plot m_hb ; % for r 1:N for h 1:m s_hb2 h,r hbb2 r-1 *m+h ; end end figure;mesh 1:N,1:511,s_hb2 ; % for h 1:m r_fft2 h,: abs fft s_hb2 h,: ; end figure; mesh 1:N,1:511,r_fft2 ; % % end 二：仿真线性调频连续波雷达的信号处理。设线性调频带宽为5MHz，时宽为102.2us，雷达载频为10GHz，输入输出为高斯白噪声，视频输入信噪比为-15dB，相干积累总时宽不大于10ms，给出回波视频表达式，脉压和FFT后的表达式；仿真给出脉压和FFT后的输出图形；通过仿真说明各级处理的增益，与各级时宽和带宽的关系；仿真说明脉压时多卜勒敏感现象。 线性调频信号：LFM信号的数