数字信号课程设计.docxVIP

PAGE PAGE # PAGE PAGE # 数字信号处理课程设计 学院：信息工程学院 班级：信息 1 班 姓名： XXX 学号： XXXXXXXXXX 课程设计目的： 通过此次课程设计使同学们加深对数字信号课程所学的知识的理解， 巩固课 堂上所学习的知识并熟悉其应用， 并通过实际操作熟练掌握 Matlab 软件的应用， 掌握用 Matlab 实现信号仿真及计算的方法及各种函数。 课程设计内容： 用 Matlab 产生下列信号并与人工分析结果进行比较 x(t)=co s(2 n t) + cos(20 n t) 0t5 用DFT计算下列信号的频谱： x(t)=5cos( n t/8+ n /4) 已知某LTI系统的微分方程为2 y(t)+y(t)+8 y(t)=x(t), 用 Matlab 求该系统的单位冲激响应和阶跃响应，并与理论值进行比较。 . 设计一个带通数字滤波器，使之满足： 通带：0.25 n = w =0.35 n ; wl=0.2 n ; wu=0.4 n ; Ap=0.8dB ; Ar=60dB ; . 某随机信号由两余弦信号与噪声组成： x[k]=cos(20 n k)+ cos(40 n k)+s[k], s[k] 为均值为 0 ，方差为 1 的高斯白噪声 (1) 绘制此随机信号的时域波形 ； (2)分别用周期图法、平均周期图法和 Welch分析该序列的功率谱估计。 课程设计步骤及结果 ： 、程序代码： t=-2:0.01:2; x=cos(2* n *t)+cos(20* n *t); plot(t,x) 图形: 二.程序代码： T0=16;N=20;T=T0/N; t=O:T:TO; x=5*cos(pi/8*t); X=1/N*fft(x,N); f=1/T/N*(-N/2:(N/2-1));%N,f=1/T/N*(-(N-1)/2:(N-1)/2); subplot(2,1,1);stem(f,abs(fftshift(X))); xlabel(Freque ncy(Hz));ylabel(mag nitude); subplot(2,1,2);stem(f,a ngle(fftshift(X))); xlabel(Freque ncy(Hz));ylabel(phase); 图形: ■4 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 Freque ncy(Hz) e 2 -unm 1 -0 --0.8esanp-0.6 -0.4 -0.20 0.2 e 2 - u n m 1 - 0 - -0.8 esanp -0.6 -0.4 -0.2 0 0.2 0.4 0.6 Freque ncy(Hz) 2 0 -2 三.程序代码： b=[1]； a=[2 1 8]; t=0:0.1:10; y=impulse(b,a,t); y仁 step(b,a,t); subplot(2,1,1);plot(t,y);title(冲激响应); subplot(2,1,2);plot(t,y1);title(阶跃响应); 图形： 0.6 冲激响应 0.4 0.2 0 -0.2 0 1 2 3 4 5 6 7 8 9 10 阶跃响应 0.4 0.3 0.2 0.1 0 0 1 2 3 4 5 6 7 8 9 10 程序代码： % ( 1)给定数字频率指标设计带通滤波器 omega0=0.5*pi; %通带中心频率 omegac 仁0.25*pi; omegac2=0.35*pi; Apmax=0.8; %通带截频 omegar1=0.2*pi; omegar2=0.4*pi; Armin=60; % 阻带截频 T=5*pi*10A(-6); % 取样间隔 %( 2)用双线性变换法设计滤波器 % 1)按给定数字频率求得相应的模拟角频率：频率预畸变 wc1=(2/T)*ta n(omegac1/2); wc2=(2/T)*ta n(omegac2/2); wr1=(2/T)*tan(omegar1/2); wr2=(2/T)*tan(omegar2/2); w0=(2/T)*ta n(omega0/2); B=wc2-wc1; % 2)将指标转换为归一化模拟低通频率指标 n orm_wr1=(((wr1A2)-(wOA2))/(B*wr1)); n orm_wr2=(((wr2A2)-(wOA2))/(B*w ⑵)； n orm_wc1=(((wc1A2)-(wOA2))/(B*wc1)); norm_wc2=(((wc2A2)-(w0A2))/(B*wc2)); % 3) 设计归一化模拟低通滤波器确定归一化截频： if (abs(norm_wr1)-abs(norm_wr2