參與評論

評論

登錄后參與評論

全部評論(1條)

• 

  qxxheart 2017-12-15 08:20:38

  以前我做的一個完整的信號濾波前后的程序，可能數(shù)據(jù)有些不一樣，但用MATLAB做濾波器的大體思路都差不多，希望對你有用。（完全自創(chuàng)，要給分哦） f=8*10^3; %信號頻率 f1=30*10^3; %噪聲頻率 l=100; %信號的長度 fs=100*10^3; %采樣頻率 %%%%%%%%%%%%%%%%生成信號 t=(0:l-1)/fs; s=0.5*sin(2*pi*f*t)+cos(2*pi*f1*t); figure(1); subplot(1，2，1); plot(fs*t(1:60)，s(1:60)); title('時域波形'); xlabel('時間'); ylabel('幅度'); grid on; nf=2^nextpow2(l); Y=fftshift(fft(s，nf)); %k=-fs/2+(0:nf-1)*fs/nf; k=fs*linspace(0，1，nf)-fs/2; subplot(1，2，2); plot(k，abs(Y)); %信號的頻譜圖 title('信號頻譜'); xlabel('頻率'); ylabel('幅度'); grid on; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%濾波器設定 ws=10*10^3; %通帶截止頻率 ws1=20*10^3; %阻帶起始頻率 ds=50; %阻帶衰減 wt=2*pi*ws/fs; %經采樣后的通帶截止角頻率 wz=2*pi*ws1/fs; %阻帶的 wc=(wt+wz)/2; %歸一化后的濾波器截止頻率 n=ceil(6.6*pi/(wz-wt)); %t=(n-1)/2; q=fir1(n-1，wc/pi，hanning(n));%濾波器時域函數(shù) w=linspace(0，pi，512); h1=freqz(q，1，512); %進行512個點的傅里葉變換 figure(2); plot(w/pi，20*log10(abs(h1))); title('濾波器頻譜圖'); xlabel('頻率'); ylabel('幅度'); grid ; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%經濾波器后的信號 signal=conv(s，q); nx=(0:n+l-2)/fs; figure(3); subplot(1，2，1); plot(nx(n:60)*fs，signal(n:60)); title('經濾波后的信號時域圖'); xlabel('時間'); ylabel('幅度'); grid on; lx=2^nextpow2(nx); Y=fftshift(fft(signal，lx)); k=-fs/2+(0:lx-1)*fs/lx; % subplot(1，2，2); plot(k，abs(Y)); %信號的頻譜圖 title('濾波后信號頻譜'); xlabel('頻率'); ylabel('幅度'); grid on;

  贊(4)

  回復(0)

  評論

  評論

  登錄后參與評論