数字信号处理实验报告-离散时间信号通过线性时不变系统精选.docxVIP

数字信号处理实验报告实验名称：离散时间信号通过线性时不变系统姓名：专业：年级：学号：指导教师：一、综合设计目的掌握信号通过IIR系统的实现方法。了解直接型和级联型网络结构实现IIR滤波时的量化效应。实验原理(1)根据滤波器技术指标计算数字滤波器阶数。(2)根据滤波器技术和通/阻带截止频率计算数字滤波器的系统函数。(3)系统量化对不同结构的敏感程度不同。三、实验过程①对数据一进行处理代码：closeall;clc;clearall;qq=16;[y,fs,bits]=wavread(C:\Users\ET\Desktop\×?o?éè??\êy?Y\′??éè?μ?ò?à?D?o?IIR\music1-1.wav);Wp=2*900;Ws=2*1500;Rp=0.1;Rs=60;Fs=32000;%éè????2¨?÷2?êywp=Wp/Fs;ws=Ws/Fs;%[N,Wpo]=cheb2ord(wp,ws,Rp,Rs);%?ó?D±è??·ò2Dí1éò??ˉ?￡?aμíí¨??2¨?÷μ??×êy%[b,a]=cheby2(N,Rs,Wpo);%éè????2¨?÷%[N,Wpo]=cheb1ord(wp,ws,Rp,Rs);%?ó?D±è??·ò1Dí1éò??ˉ?￡?aμíí¨??2¨?÷μ??×êy%[b,a]=cheby1(N,Rp,Wpo);%éè????2¨?÷% [N,Wpo]=buttord(wp,ws,Rp,Rs);%?ó°íì????1Dí1éò??ˉ?￡?aμíí¨??2¨?÷μ??×êy% [b,a]=butter(N,Wpo);%éè????2¨?÷[N,Wpo]=ellipord(wp,ws,Rp,Rs);%?óí??2Dí1éò??ˉ?￡?aμíí¨??2¨?÷μ??×êy[b,a]=ellip(N,Rp,Rs,Wpo);%éè??í??2??2¨?÷[h,w]=freqz(b,a,512);%?ó??2¨?÷μ?·ù?μ?ìó|h=20*log10(abs(h)/max(abs(h)));bqr=intbR(b,qq);aqr=intbR(a,qq); [hr,wr]=freqz(bqr,aqr,512);hr=20*log10(abs(hr)/max(abs(hr))); [sos,g]=tf2sos(b,a); b1=sos(1,1:3);a1=sos(1,4:6); bqr1=intbR(b1,qq); aqr1=intbR(a1,qq); [h1,w1]=freqz(bqr1,aqr1,512); sos1(1,1:6)=[bqr1 aqr1]; b2=sos(2,1:3);a2=sos(2,4:6); bqr2=intbR(b2,qq); aqr2=intbR(a2,qq); [h2,w1]=freqz(bqr2,aqr2,512); sos1(2,1:6)=[bqr2 aqr2]; b3=sos(3,1:3);a3=sos(3,4:6); bqr3=intbR(b3,qq); aqr3=intbR(a3,qq); [h3,w1]=freqz(bqr3,aqr3,512); sos1(3,1:6)=[bqr3 aqr3];% b4=sos(4,1:3);a4=sos(4,4:6); bqr4=intbR(b4,qq); aqr4=intbR(a4,qq); [h4,w1]=freqz(bqr4,aqr4,512); sos1(4,1:6)=[bqr4 aqr4];% b5=sos(5,1:3);a5=sos(5,4:6); bqr5=intbR(b5,qq); aqr5=intbR(a5,qq); [h5,w1]=freqz(bqr5,aqr5,512); sos1(5,1:6)=[bqr5 aqr5];% b6=sos(6,1:3);a6=sos(6,4:6); bqr6=intbR(b6,qq); aqr6=intbR(a6,qq); [h6,w1]=freqz(bqr6,aqr6,512); sos1(6,1:6)=[bqr6 aqr6];%b7=sos(7,1:3);a7=sos(7,4:6); bqr7=intbR(b7,qq); aqr7=intbR(a7,qq); [h7,w1]=freqz(bqr7,aqr7,512); sos1(7,1:6)=[bqr7 aqr7];%b8=sos(8,1:3);a8=sos(8,4:6); bqr8=intbR(b8,qq); aqr8=intbR(a8,qq); [h8,w1]=freqz(bqr8,aqr8,512); sos1(8,1:6)=[bqr8