发布时间 : 星期日 文章南京邮电大学DSP实验报告更新完毕开始阅读f8e455a19e314332396893e9
反三角波频谱图变密，两种图频谱图没有相似处。
(5) 用FFT分别实现xa(n)（p＝8，q＝2）和xb(n)（a＝0.1，f＝0.0625）的16点循环卷积和线性卷积。 解： clear n=0:1:15; p=8; q=2;
xa=exp(-(n-p).^2/q);
a=0.1; f=0.0625;
xb = exp(-a*n).*sin(2*pi*f*n);
N = 16;
Xa=fft(xa,N); Xb=fft(xb,N); Y=Xa.*Xb; y=ifft(Y,N); m=0:1:N-1 stem(m,y) xlabel('m') ylabel('y')
title('循环卷积')
循环卷积21.510.5y0-0.5-105m1015
clear n=0:1:15; p=8; q=2;
xa=exp(-(n-p).^2/q);
a=0.1; f=0.0625;
xb = exp(-a*n).*sin(2*pi*f*n);
N =31;
Xa=fft(xa,N); Xb=fft(xb,N); Y=Xa.*Xb; y=ifft(Y,N); m=0:1:N-1 stem(m,y) xlabel('m') ylabel('y')
title('线性卷积')
线性卷积21.510.5y0-0.5-1051015m202530
实验三 IIR数字滤波器的设计
一、实验目的
(1) 掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。 (2) 观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。
(3) 熟悉巴特沃思滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。
二、实验内容 （1）P162 例4.4
设采样周期T=250?s(采样频率fs=4kHz)，分别用脉冲响应不变法和双线性变换法设计一个三阶巴特沃思低通滤波器，其3dB边界频率为fc=1kHz。
解
%脉冲响应不变法
fc=1000;
fs=4000;
[B,A]=butter(3, 2*pi*fc,'s'); [num1,den1]=impinvar(B,A,fs); [h1,w]=freqz(num1,den1); f = w/pi*fs/2;
%双线性变换法 fc=1000; fs=4000;
[B,A]=butter(3, 2*fs*tan(pi*fc/fs),'s'); [num2,den2]=bilinear(B,A,fs); [h2,w]=freqz(num2,den2); f = w/pi*fs/2;
plot(f,abs(h1),'-.',f,abs(h2),'-'); grid on
legend('脉冲响应不变法','双线性变换法') xlabel('频率/Hz') ylabel('幅度')
(选做) (2)fc=0.2kHz,?=1dB,fr=0.3kHz,At=25dB,T=1ms；分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。
每一小题均给出实验过程与结果（含实验程序、运行的数据结果和图形）； 解： clear; fc = 200; Rp = 1; fr = 300; Rs = 25; T = 0.001; fs = 1/T;