数字信号处理Test讲解.ppt

s1=[0 0 2 3 5 2 ]; s2=[2 3 5 2 1 -1]; XY=fft(s2,32).*conj(fft(s1,32)); R12=real(ifft(XY)); L=round(length(R12)/2); r12=[R12(L+2:2*L) R12(1:L)]; % figure(3) ,stem(-L+1:L-1,r12),grid figure(4), stem(R12),grid 32点圆周相关 其中哪些点等于线性相关？？ 线性相关 线性相关的取值范围？ Test1 设计切比雪夫高通数字滤波器 fp=300; fr=200; T=.001; ap=0.8; as=20; wp=2*fp*T; ws=2*fr*T; % 以fs/2(或pi)归一化的通带、阻带截止频率，可为矢量。 [n,wn]=cheb1ord(wp,ws,ap,as); % 用双线性变换完成AF DF ； % ap：通带内波纹， as：阻带最大衰减（分贝）； % n：阶数； wn：3dB截止频率； Ch3 [b,a]=cheby1(n,ap,wn,high); % n: 阶数; % ap ：通带内波纹（分贝）； % wn为3dB截止频率，0 wn 1； 如果wn=[w1 w2]，默认为带通，‘stop’带阻。 % high：高通，默认：低通， % b,a：系统函数的分子、分母多项式系数（矢量）。 [h,w]=freqz(b,a,100,whole); % h: 在单位圆上采样100点的频率响应 % w：频率矢量 c=20*log10(abs(h)); (Hz) figure(1), plot(w(1:50), c(1:50)); title(切比雪夫高通数字滤波器); axis([0,500,-80,0]); grid xlabel ( 频率/Hz ); ylabel ( 幅度/dB ); Test2 fs=1000; wp=2*pi*200;ws=2*pi*300 [n,wn]=buttord(wp,ws,1,25,s); %设计Butterworth低通模拟滤波器 通带衰减小于1分贝、阻带（0.6pi）衰减大于25分贝 n: 阶数 wn:低通模拟滤波器的3分贝通带截止频率 [B,A] = BUTTER(n,wn,s); % 设计Butterworth低通模拟滤波器 B,A：模拟滤波器传递函数分子、分母多项式系数 [b,a]=impinvar(B,A,fs); % 用脉冲响应不变法将模拟低通滤波器转换为数字滤波器b,a：分子、分母多项式系数 [b,a]=bilinear(B,A,fs); Test 3 [n,wn]=ellipord(wp,ws,ap,as); [b,a] = ellip(n,ap,as,wn); Help ellip Test4 w1=4*pi/30; w2=6*pi/30; %计算 w3=3*pi/30; w4=12*pi/30; cw0=sin(w1+w2)/(sin(w1)+sin(w2)); %计算 wp=abs((cw0-cos(w1))/sin(w1)); %计算 ws1=abs((cw0-cos(w3))/sin(w3)); % 计算 ws2=abs((cw0-cos(w4))/sin(w4)); % 计算 [n1,wn1]=buttord(wp,ws1,3,20,s); % n1:对应阻带1及其衰减的滤波器阶数 [n2,wn2]=buttord(wp,ws2,3,5,s); % n2:对应阻带2及其衰减的滤波器阶数 n=max(n1,n2); % 取最大阶数以满足所有指标 [b,a]=butter(n,[2/15 3/15]); % [2/15 3/15] 通带截止频率fl=2kc,fh=3kc用fs/2作归一化 ch5 Test2 n=15; b=fir1(n-1,wn,hanning(n)); % 阶数：n-1, 窗长：n。 % 0wn1; fs/2归一化后的频率 % 带通、带阻：wn=[w1 w2]; wn=[0.3 0.5]; % 多带：wn=[w1 w2 …….wN]; 缺省：Hamming窗 % fir1(N,Wn,‘h