数字课程设计.docVIP

数字信号处理 课程设计 学院：信息工程学院 班级：信息08—1班 姓名：### 学号：### 课程设计目的： 通过此次课程设计使同学们加深对数字信号课程所学的知识的理解，巩固课堂上所学习的知识并熟悉其应用，并通过实际操作熟练掌握Matlab软件的应用，掌握用Matlab实现信号仿真及计算的方法及各种函数。 课程设计内容： 一.用 Matlab产生下列信号并与人工分析结果进行比较 x(t)=cos(2πt)+ cos(20πt) 0t5 二.用 DFT计算下列信号的频谱： x(t)=5cos(πt/8+π/4) 三.已知某LTI系统的微分方程为 2 y``(t)+y`(t)+8 y(t)=x(t) , 用Matlab求该系统的单位冲激响应和阶跃响应，并与理论值进行比较。 四.设计一个带通数字滤波器，使之满足： 通带：0.25π= w =0.35π ; wl=0.2π ; wu=0.4π; Ap=0.8dB ; Ar=60dB ; 五.某随机信号由两余弦信号与噪声组成： x[k]=cos(20πk)+ cos(40πk)+s[k] , s[k]为均值为0 ，方差为1的高斯白噪声 绘制此随机信号的时域波形 ； 分别用周期图法、平均周期图法和Welch分析该序列的功率谱估计。 课程设计步骤及结果： 一、程序代码： t=-2:0.01:2; x=cos(2*π*t)+cos(20*π*t); plot(t,x) 图形： 二.程序代码： T0=16;N=20;T=T0/N; t=0:T:T0; 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(Frequency(Hz));ylabel(magnitude); subplot(2,1,2);stem(f,angle(fftshift(X))); xlabel(Frequency(Hz));ylabel(phase); 图形： 三.程序代码： b=[1]; a=[2 1 8]; t=0:0.1:10; y=impulse(b,a,t); y1=step(b,a,t); subplot(2,1,1);plot(t,y);title(冲激响应); subplot(2,1,2);plot(t,y1);title(阶跃响应); 图形： 四.程序代码： %（1）给定数字频率指标 设计带通滤波器 omega0=0.5*pi; %通带中心频率 omegac1=0.25*pi; omegac2=0.35*pi; Apmax=0.8; %通带截频 omegar1=0.2*pi; omegar2=0.4*pi; Armin=60; %阻带截频 T=5*pi*10^(-6); %取样间隔 %（2）用双线性变换法设计滤波器 % 1）按给定数字频率求得相应的模拟角频率：频率预畸变 wc1=(2/T)*tan(omegac1/2); wc2=(2/T)*tan(omegac2/2); wr1=(2/T)*tan(omegar1/2); wr2=(2/T)*tan(omegar2/2); w0=(2/T)*tan(omega0/2); B=wc2-wc1; % 2) 将指标转换为归一化模拟低通频率指标 norm_wr1=(((wr1^2)-(w0^2))/(B*wr1)); norm_wr2=(((wr2^2)-(w0^2))/(B*wr2)); norm_wc1=(((wc1^2)-(w0^2))/(B*wc1)); norm_wc2=(((wc2^2)-(w0^2))/(B*wc2)); % 3) 设计归一化模拟低通滤波器确定归一化截频： if (abs(norm_wr1)-abs(norm_wr2))0 norm_wr=abs(norm_wr2); else norm_wr=abs(norm_wr1); end norm_wc=1; % 确定归一化模拟低通滤波器阶数： N=buttord(norm_wc,norm_wr,Apmax,Armin,s); % 设计归一化模拟低通滤波器： [b_LP,a_LP]=butter(N,norm_wc,s); % 4) 把归一化模拟低通滤波器转换为模拟带通滤波器 [b_BP,a_BP]=lp2bp(b_LP,a_LP,w0,B); % 5) 转换为数字滤波器 Fs=1/T; [b,a]=bilinear(b