实验二数字滤波器的设计与实现报告.docxVIP

实验二 数字滤波器的设计与实现报告 班级 姓名 学号 时间 一 实验目的1.掌握采样定理的应用技术2.掌握数字滤波器的设计方法3.熟练应用频谱分析技术4.掌握matlab下低通滤波器应用程序的设计方法二 实验结果①学号19，f1=19*100=1900Hz,f2=19*300=5700Hz②滤波三 附录程序SID=[2,0,1,3,4,0,8,1];%我的学号累加和SS=sum(SID);f1=SS*100;f2=SS*300;fs=15000;N=256;n=0:N-1;t=n/fs;x1=sin(2*pi*f1*t);%figure(1);%plot(t,x1);%作正弦信号的时域波形;xlabel(t);ylabel(y);%title(正弦信号1900Hz时域波形);%grid;%进行FFT变换并做频谱图y=fft(x1,N);%进行fft变换mag=abs(y);%求幅值f=(0:length(y)-1)*fs/length(y);%进行对应的频率转换figure(1);subplot(3,1,1);plot(f,mag);%做频谱图set(gca,XTick,0:2500:15000);ylabel(低频幅值);title(正弦信号1900HZ幅频谱图);x2=sin(2*pi*f2*t);%figure(3);%plot(t,x2);%作正弦信号的时域波形;set(gca,XTick,0:2500:15000);%title(正弦信号5700Hz时域波形);%grid;%进行FFT变换并做频谱图y=fft(x2,N);%进行fft变换mag=abs(y);%求幅值f=(0:length(y)-1)*fs/length(y);%进行对应的频率转换figure(1);subplot(3,1,2);plot(f,mag);%做频谱图set(gca,XTick,0:2500:15000);xlabel(频率(Hz));ylabel(高频幅值);title(正弦信号5700Hz幅频谱图);x=x1+0.5*x2; figure(2);subplot(2,2,1);plot(t,x);%作正弦信号的时域波形;title(混合信号X时域波形);grid;%进行FFT变换并做频谱图y=fft(x,N);%进行fft变换mag=abs(y);%求幅值f=(0:length(y)-1)*fs/length(y);%进行对应的频率转换figure(1);subplot(3,1,3);plot(f,mag);%做频谱图set(gca,XTick,0:2500:15000);xlabel(频率(Hz));ylabel(幅值);title(混合信号X幅频谱图);% 设计滤波器进行滤波%进行FFT变换并做频谱图y=fft(x,N);%进行fft变换mag=abs(y);%求幅值f=(0:length(y)-1)*fs/length(y);%进行对应的频率转换figure(2);subplot(2,2,2);plot(f,mag);%做频谱图set(gca,XTick,0:2500:15000);xlabel(频率(Hz));ylabel(幅值);title(混合信号X幅频谱图);Wn = (2*2500)/15000;[B,A] = butter(5,Wn,low);y = filter(B,A,x);figure(2);subplot(2,2,3);plot(y);title(滤波后信号);C=fft(y,N);%进行fft变换mag=abs(C);%求幅值f=(0:length(C)-1)*fs/length(C);%进行对应的频率转换figure(2);subplot(2,2,4);plot(f,mag);%做频谱图set(gca,XTick,0:2500:15000);xlabel(频率(Hz));ylabel(幅值);title(滤波后幅频谱图);