孙利栋3100403032.doc

信号与系统大作业 题 目 滤波器 姓 名 孙利栋 　 学 号 3100403032 　　 专业班级 通信101 　 　 指导教师 朱勇 　　 学 院 信息学院 　　 完成日期 2012年5 月 31 日 【实验目的】 加深对matlab的学习 学习滤波器的工作原理 【实验原理】 输入信号由许多信号组合而成，需要把他们分开。首先需要把原信号x=cos(2*pi*200*t)+cos(2*pi*400*t)+cos(2*pi*800*t)通过FFT傅里叶快速变换成频域上的波形。然后通过低通，高通，带通，把需要的波形滤出来。 由上式可知，我们在频域上滤出的波形是200Hz,400Hz,800Hz的量。低通滤出200Hz的波，带通滤出400Hz的，高通滤出800Hz的。我选择了巴特沃斯滤波器。 巴特沃斯滤波器的特点：通带处幅值特性平坦切比雪夫一般而言，阶数越高，则逼近越精确，但计算代价也随之上升，所以性能与代价总需要寻求一个平衡点滤波器在低通频带里 ，滤波器对低频信号的增益是平稳的，就是放大倍数一样，所有的低频信号得到一样的放大倍数， 但是对不同频率的信号相位有些偏移，这样的结果是可能使通过滤波器的信号产生失真。但是当误差允许的范伟内，可以认识相位失真可以不考虑。 [B,A]=butter(N,wc); 原始式子[B,A] = butter(N, wc,S) N为设计滤波器的阶数?N越大，通带和阻带的近似性越好，过渡带也越陡。 S ?为字符串，表明设计的滤波器类型，low低通/high高通/stop带阻 为要求的标准化截至频率clear;clear clf; f1=200;f2=400;f3=800;fs=2400;t=(0:1000)/fs%设置采样频率 x=cos(2*pi*f1*t)+cos(2*pi*f2*t)+cos(2*pi*f3*t);%原信号 L=length(x);N=2^(nextpow2(L));%nextpow2：靠的最近的2的指数 。比如数1000，那么靠的最近的（且比他大的）数就是1024=2^10，这个函数一般用在需要数组长度为2的指数的情况下，比如fft hw=fft(x,N)%快速傅里叶变换 figure(1); subplot(2,1,1);plot(t,x);grid on;title(滤波前信号x);xlabel(时间/s);% 原始信号 subplot(2,1,2);plot((0:N-1)*fs/L,abs(hw)/max(abs(hw)));grid on;% 查看信号频谱 title(滤波前信号频谱图);xlabel(频率/Hz);ylabel(振幅|H(e^jw)|); %巴特沃斯低通 Ap=1;As=50;% 定义通带及阻带衰减 wp=2*220/fs ; %通带频率归一化后Matlab使用归一化的频率参数（临界频率除以采样频率的1/2）。这样临界频率参数的取值范围在0和1之间1代表Fs/2（用角频率表示的时候对应π）ws=2*370/fs; %归一化 [N,wc]=buttord(wp,ws,Ap,As);%估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc [B,A]=butter(N,wc);%设计Butterworth低通滤波器 [H,W]=freqz(B,A);%求数字低通滤波器的频率响应 y1=filter(B,A,x);%余弦波叠加Y = FILTER(B,A,X)，输入X为滤波前序列，Y为滤波结果序列，B/A 提供滤波器系数，B为分子，A为分母整个滤波过程是通过下面差分方程实现的：a(1)*y(n) = b(1)*x(n) + b(2)*x(n-1) + ... + b(nb+1)*x(n-nb) L=length(y1);N=2^(nextpow2(L));hw_1=fft(y1,N) figure(2); subplot(2,1,1);plot(t,y1);title(经低通滤波后信号x);xlabel(时间t/s); subplot(2,1,2);plot((0:N-1)*fs/L,abs(hw_1));title(经