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

实验二数字滤波器的设计与实现报告班级姓名学号时间一实验目的1.掌握采样定理的应用技术2.掌握数字滤波器的设计方法3.熟练应用频谱分析技术4.掌握matlab下低通滤波器应用程序的设计方法二实验结果=1\*GB3①学号和20134081=19，f1=19*100=1900Hz,f2=19*300=5700Hz=2\*GB3②滤波三附录程序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);p