数字信号处理课程设计离散时间信号处理.docx

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