Matlab基本知识和信号处理工具箱.ppt

LOGO 函数firls设计出的线性相位FIR滤波器，可使指定频段内的理想分段线性函数与滤波器幅度响应之间的误差平方和最小。 b=firls(n,f,m)可设计出n阶FIR滤波器，其幅频特性匹配于由f和m给出的特性。 LOGO 例3 设计一24阶的fir滤波器，使其具有分段线性通带。画出实际和理想的频率响应。其程序如下： f=[0 0.3 0.4 0.6 0.7 0.9]; m=[0 1 0 0 0.5 0.5]; b=firls(24,f,m); for i=1:2:6 plot([f(i) f(i+1)],[m(i),m(i+1)],’—’), hold on end [h,f]=freqz(b,1,512,2); plot(f,abs(h)); grid on LOGO 实验一 上一次试验中，我们设计了如下的IIR滤波器： 示例三 设计一低通切比雪夫I型数字滤波器，通带范围0-100Hz，通带波纹3dB，阻带衰减 -30dB，系统采样率为1000Hz。 请用fir1函数完成相同的设计要求。（提示：请参考本次实验的例1） 思考：如果fir滤波器的设计采用与iir滤波器相同的阶次，请对比两者的幅频响应，思考为什么。 LOGO 实验二 两个频率为0.2pi和0.4pi的正弦信号叠加在一起，请用fir2函数设计一FIR滤波器将频率为0.2pi的正弦信号提取出来。 LOGO n=0:1:99; x=cos(0.2*pi*n)+cos(0.3*pi*n); subplot(211);stem(n,x);title(长度为100的信号x(n)); xlabel(n); axis([0,99,-2.5,2.5]); X=fft(x);magX=abs(X(1:51)); k=0:50;w=2*pi/100*k; subplot(212);plot(w,magX);title(幅频响应); xlabel(frequency in pi); axis([0,pi,0,60]); LOGO f=[0 0.25 0.3 1] m=[1 1 0 0] b=fir2(30,f,m) [h,w]=freqz(b,1,128); figure(2) plot(f,m,w/pi,abs(h)); y=conv(x,b) figure(3) stem(y) LOGO LOGO [s,fs]=wavread(I:\anjian1); figure(1) plot(s) axis([0, 5e5, -1.5,1.5]) s1=s(38e4:39e4); fs1=fft(s1); k=(0:floor(length(fs1)/2)); kf=k*fs/length(fs1); figure(2) plot(kf,abs(fs1(1:length(kf)))) shg LOGO LOGO 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * LOGO 其中 DFT与IDFT变换： 一、用Matlab实现DFT 与IDFT变换 LOGO 用Matlab实现上面的变换和反变换，程序如下： function Xk=dft(xn,N) n=[0:1:N-1]; k=[0:1:N-1]; WN=exp(-j*2*pi/N); nk=n*k; WNnk=WN.^nk; Xk=xn*WNnk; LOGO function xn=idft(Xk,N) n=[0:1:N-1]; k=[0:1:N-1]; WN=exp(-j*2*pi/N); nk=n*k; WNnk=WN.^(-nk); xn=(Xk*WNnk)/N; LOGO 二、计算分辨率和物理分辨力 的分析 调用上述dft函数给出下列序列的DFT变换X(k)： 再调用idft函数，由X(k)变换回x(n)，看是否一致。 LOGO 利用Matlab来分别讨论频率的计算分辨率和物理分辨力 为了说明两者之间的差异，请思考下面的例子： Matlab实现程序如下所示： LOGO n=0:1:99; x=cos(0.48*pi*n)+cos(0.52*pi*n); figure(1); n1=0:1:9;y1=x(1:10); subplot(211);stem(n1,y1); title(长度为10的