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

《数字信号处理课程设计--离散时间信号处理》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、离散时间信号处理课程设计设计要求：某雷达系统接收机框架如图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接收信号参数1.请设计ADC的采样率fs，画出x[n]的频谱X(ejw)2.请推导正交解调器的输出S_I[n]和S_Q[n]的表达式，画出频谱示意图；若要将信号

2、谱搬移到零中频，请确定NCO的频率W03.请设计一个FIR线性相位数字LPF对正交解调的输出进行处理，要求杂散抑制比超过50dBc，确定滤波器设计指标、给出设计过程和结果。4.请确定采样率变换模块的参数(抽取或内插系数)，画出输出信号谱。设计过程1.设计ADC的采样率fs，画出x[n]的频谱X(ejw)Adc采样率fs>=2*f0，fs取100MHZ。φ(t)为相位调制函数(调制带宽B=2MHz)，取φ(t)=cos(2π**t)作为单音调制信号，则,,,程序如下：clear;clc;clf;n0=20000;%点数sprate=;%采样频率100Mt0=n0/spra

3、te;%总时长fca=;%载波频率fsi=;%信号频率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');gridon;legend('modulatedsignal');subplot(3,1,2);w=linspace(-1,1,n0);plot(w,fftshift(abs(fft(Dsb))));title('幅频');xlabel('w'

4、);ylabel('

5、X(ejw)

6、');gridon;legend('modulatedsignal');subplot(3,1,3);plot(w,fftshift(angle(fft(Dsb))));title('相频');xlabel('w');ylabel('∠X(ejw)');gridon;legend('modulatedsignal');上图中频域图0.8代表0.8*/2=40M2.推导正交解调器的输出S_I[n]和S_Q[n]的表达式，画出频谱示意图；若要将信号谱搬移到零中频，请确定NCO的频率W0W0==4π/5S_I[n]=x[n]*cos[w0n

7、]=cos[4πn/5+cos[πn/50]]*cos[4πn/5]Sl=Dsb.*cos(2*pi*fca*t);subplot(3,1,1);plot(t(1:100),Sl(1:100));title('时域');xlabel('t');ylabel('voltage');gridon;legend('modulatedsignal');subplot(3,1,2);plot(w,fftshift(abs(fft(Sl))));title('幅频');xlabel('w');ylabel('

8、X(ejw)

9、');gridon;legend('modulatedsi

10、gnal');subplot(3,1,3);plot(w,fftshift(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(3,1,1);plot(t(1:100),Sq(1:100));title('时域');xlabel('t');ylabel('voltage')

11、;gridon;legend('modulatedsignal');subplot(3,1,2);plot(w,fftshift(abs(fft(Sq))));title('幅频');xlabel('w');ylabel('

12、X(ejw)

13、');gridon;legend('modulatedsignal');subplot(3,1,3);plot(w,fftshift(angle(fft(Sq))));title('相频');xlabel('w');ylabel('∠X(ejw)');gridon;legend('modulatedsi