實驗二數字濾波器的設計與實現報告班級姓名學號時間一實驗目的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