实验一 Matlab软件的使用..

由天下分享时间：2025/2/28 19:23:14 加入收藏我要投稿点赞

=1;

h([find((n>=0)&(n<=8))]) =0.5;

subplot(3,1,1); stem(x); title('x(n)'); axis([0,30,0,2]); subplot(3,1,2); stem(h); title('h(n)'); axis([0,30,0,2]); X=fft(x); H=fft(h); Y=X.*H; y=ifft(Y); subplot(3,1,3);

实验结果： 1、

2、

实验四 IIR数字滤波器的设计和实现

1、实验目的和要求

㈠掌握冲击响应不变法和双线形变换法设计IIR数字滤波器的原理和方法；㈡通过观察对实际信号的滤波，获得对数字滤波的感性认识。 2、实验内容

㈠巴特沃思滤波器设计实现模拟低通滤波器；

㈡利用冲击响应不变法和双线形变换法对模拟滤波器进行A/D转换； 3、实验原理

用数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲激响应等于

ha(t),让h(n)正好

ha(t)的采样值，即h(n)?ha(nT)。

4、主要仪器设备计算机，Matlab软件。

实验步骤：

设计一个巴特沃思数字低通滤波器，设计指标如下：通带内大于1dB；阻带现数字滤波。

function shiyan3() fp=200*pi; fs=300*pi; wp=0.2*pi; ws=0.3*pi; Rp=1; As=15; T=1; Fs=1000;

[N,fc]=buttord(fp,fs,Rp,As,'s'); [b,a]=butter(N,fc,'s'); w=[0:1000*2*pi];

[hf,w]=freqs(b,a,1000); [d,c]=impinvar(b,a,Fs); wd=[0:512]*pi/512; hw1=freqz(d,c,wd); [f,e]=bilinear(b,a,Fs); hw2=freqz(f,e,wd); OmegaP=(2/T)*tan(wp/2); OmegaS=(2/T)*tan(ws/2); ep=sqrt(10^(Rp/10)-1); Ripple=sqrt(1/(1+ep.^2));

???p?0.2?幅度衰减不

???s?0.35?幅度衰减不小于15dB。编制实现该数字滤波器程序并且实

Attn=1/10^(As/20); subplot(3,1,1);

plot(wd/pi,abs(hw1)/abs(hw1(1))); title('幅度响应(冲击响应不变法)'); xlabel('w(pi)'); ylabel('H'); ylabel('H'); axis([0,1,0,1.1]);

set(gca,'XTickmode','manual','XTick',[0,0.2,0.35,1.1]); set(gca,'YTickmode','manual','YTick',[0,Attn,Ripple,1]); grid;

subplot(3,1,2);

plot(wd/pi,20*log10((abs(hw2)+eps)/abs(hw2(1)))); title('幅度响应(双线性变换法(dB))'); xlabel('w(pi)'); ylabel('H'); axis([0,1,-40,5]);

set(gca,'XTickmode','manual','XTick',[0,0.2,0.35,1.1]); set(gca,'YTickmode','manual','YTick',[-50,-15,-1,0]); grid;

subplot(3,1,3);

plot(wd/pi,-angle(hw2)); title('相位响应'); xlabel('w(pi)'); ylabel('pi unit'); %axis([0,1,0,1.1]);

set(gca,'XTickmode','manual','XTick',[0,0.2,0.35,1.1]); set(gca,'YTickmode','manual','YTick',[-1,0,1]); grid; 实验结果：

实验五 FIR数字滤波器的设计和实现