数字信号处理课程设计.docx

离散时间信号处理课程设计学号：2012029070019姓名：谢超任课老师：唐斌设计要求：某雷达系统接收机框架如图1所示。接收机输入信号其中中心频率=40MHz为相位调制函数(调制带宽B=2MHz)信号谱如图2所示。输出基带信号I[n]、Q[n]数据率(采样率)为2.5MHz。LPFADCNCOLPF采样率变换采样率变换数字信号处理单元Xc(t)cos[w0n]sin[w0n]fsX[n]S_I[n]S_Q[n]图1接收机原理框架BfXc(f)f0-f0B图2接收信号参数请设计ADC的采样率fs画出x[n]的频谱X(ejw)请推导正交解调器的输出S_I[n]和S_Q[n]的表达式画出频谱示意图；若要将信号谱搬移到零中频请确定NCO的频率W0请设计一个FIR线性相位数字LPF对正交解调的输出进行处理要求杂散抑制比超过50dBc确定滤波器设计指标、给出设计过程和结果。请确定采样率变换模块的参数(抽取或内插系数)画出输出信号谱。设计过程1.设计ADC的采样率fs画出x[n]的频谱X(ejw)Adc采样率fs>=2*f0fs取100MHZ。φ(t)为相位调制函数(调制带宽B=2MHz)取φ(t)=cos(2π*1000000*t)作为单音调制信号则程序如下：clear;clc;clf;n0=20000;%点数sprate=100000000;%采样频率100Mt0=n0/sprate;%总时长fca=40000000;%载波频率fsi=1000000;%信号频率t=linspace(0t0n0);Dsb=cos(2*pi*fca*t+cos(fsi*2*pi*t));subplot(311);plot(t(1:100)Dsb(1:100));title('时域');xlabel('t');ylabel('voltage');gridon;legend('modulatedsignal');subplot(312);w=linspace(-11n0);plot(wfftshift(abs(fft(Dsb))));title('幅频');xlabel('w');ylabel('|X(ejw)|');gridon;legend('modulatedsignal');subplot(313);plot(wfftshift(angle(fft(Dsb))));title('相频');xlabel('w');ylabel('∠X(ejw)');gridon;legend('modulatedsignal');上图中频域图0.8代表0.8*100000000/2=40M2.推导正交解调器的输出S_I[n]和S_Q[n]的表达式画出频谱示意图；若要将信号谱搬移到零中频请确定NCO的频率W0W0==4π/5S_I[n]=x[n]*cos[w0n]=cos[4πn/5+cos[πn/50]]*cos[4πn/5]Sl=Dsb.*cos(2*pi*fca*t);subplot(311);plot(t(1:100)Sl(1:100));title('时域');xlabel('t');ylabel('voltage');gridon;legend('modulatedsignal');subplot(312);plot(wfftshift(abs(fft(Sl))));title('幅频');xlabel('w');ylabel('|X(ejw)|');gridon;legend('modulatedsignal');subplot(313);plot(wfftshift(angle(fft(Sl))));title('相频');xlabel('w');ylabel('∠X(ejw)');gridon;legend('modulatedsignal');S_Q[n]=x[n]*sin[w0n]=cos[4πn/5+cos[πn/50]]*sin[4πn/5]Sq=Dsb.*sin(2*pi*fca*t);subplot(311);plot(t(1:100)Sq(1:100));title('时域');xlabel('t');ylabel('voltage');gridon;legend('modulatedsignal');subplot(312);plot(wfftshift(abs(fft(Sq))));title('幅频');xlabel('w');ylabel('|X(ejw)|');gridon;legend('modulatedsignal');subplot(313);plot(wfftshift(angle(