標(biāo)題: MATLAB數(shù)字濾波器程序 Hamming窗帶通濾波器 [打印本頁]
作者: 小許家的豬    時間: 2020-1-2 00:48
標(biāo)題: MATLAB數(shù)字濾波器程序 Hamming窗帶通濾波器

Hamming窗帶通濾波器

%Hamming窗帶通濾波器的設(shè)計

fp1=200;fc1=100;

fpu=8000;fcu=10000;

wlp=2*pi*fp1/FS;wls=2*pi*fc1/FS;

wup=2*pi*fpu/FS;wus=2*pi*fcu/FS;

Bt=wlp-wls;            

N0=ceil(6.6*pi/Bt);

N=N0+mod(N0+1,2);

wc=[(wls+wlp)/2/pi,(wus+wup)/2/pi];   

hn=fir1(N-1,wc,hamming(N));

[h1,w1] = freqz(hn,1,512,FS);

figure(4);

subplot(2,1,1);plot(hn);axis ([700 900 -0.5 1]);xlabel('n');ylabel('幅度');title('hamming窗帶通時域圖');

subplot(2,1,2);plot(w1,abs(h1));title('hamming窗帶通頻譜圖');axis tight;xlabel('f/Hz');ylabel('耗損(dB)');

pause(1)；

%對濾波后的信號進(jìn)行分析變換

X=conv(hn,x);

n2=length(X);      

n2=2^nextpow2(n2);

f=FS*(0:n2/2-1)/n2;  

figure(5);subplot(2,1,1);plot(X);title('濾波后的信號時域圖');axis tight;xlabel('f/Hz');ylabel('幅度');

X1=fft(X,n2);subplot(2,1,2);plot(f,abs(X1(1:n2/2))); axis tight;xlabel('f/Hz');ylabel('幅度');title('濾波后的信號頻譜');

sound(X,FS);

圖1  語音信號

圖2  加噪信號

圖3  濾波信號

Hamming窗帶通濾波器程序


  1. %hamming帶通

  3. %產(chǎn)生語音信號

  5. clc;

  7. Fs=48000;

  9. [x,FS]=audioread('C:\Users\SDHH\Documents\錄音\luyin.m4a');

  11. x=x(:,1);

  13. sound(x,FS);

  15. %頻譜分析

  17. n=length(x);      

  19. n=2^nextpow2(n); %選取變換的點(diǎn)數(shù)

  21. t=(0:(n-1))/FS;%計算音頻信號的長度

  23. x=[x',zeros(1,n-length(x))]';

  25. figure(1);

  27. subplot(2,1,1); plot(t,x); axis tight; title('語音信號時域圖'); xlabel('t/s');ylabel('幅度')

  29. y=fft(x,n);          %對n點(diǎn)進(jìn)行傅里葉變換到頻域

  31. f=FS*(0:n/2-1)/n;             % 對應(yīng)點(diǎn)的頻率

  33. subplot(2,1,2);

  35. plot(f,abs(y(1:n/2)));axis tight;xlabel('f/Hz');ylabel('幅度');title('語音信號頻域圖');

  37. pause(2.5);

  39. %###########################################################################################################

  41. %產(chǎn)生噪聲信號

  43. noise=1*sin(2*pi*20000*t)+1*sin(2*pi*200*t);

  45. Noise=fft(noise,n);%對n點(diǎn)進(jìn)行傅里葉變換到頻域

  47. figure(2);

  49. subplot(2,1,1);plot(t,noise);axis tight;xlabel('t/s');ylabel('幅度');title('加噪聲信號時域圖');

  51. subplot(2,1,2)

  53. plot(f,abs(Noise(1:n/2)));     %加噪語音信號的頻譜圖

  55. axis axis([1 40000 1 50000] );xlabel('f/Hz');ylabel('幅度');title('加噪語音信號頻譜圖');

  57. pause(1);

  59. %#############################################################################################################

  61. x1=x+noise';  %將兩個信號疊加成一個新的信號——加噪聲處理

  63. sound(x1,FS);

  65. %加噪后頻譜分析

  67. noise_1=fft(x1,n);%對n點(diǎn)進(jìn)行傅里葉變換到頻域

  69. figure(3);

  71. subplot(2,1,1)

  73. plot(t,x1);axis tight;xlabel('t/s');ylabel('幅度');title('加噪聲信號時域圖');

  75. subplot(2,1,2)

  77. plot(f,abs(noise_1(1:n/2)));     %加噪語音信號的頻譜圖

  79. axis ([1 30000 1 4000] ;xlabel('f/Hz');ylabel('幅度');title('加噪語音信號頻譜圖');

  81. pause(2.5);

  83. %#################################################################################################################

  85. %hamming窗帶通濾波器設(shè)計

  87. fp1=200;fc1=100;

  89. fpu=8000;fcu=10000;

  91. wlp=2*pi*fp1/FS;wls=2*pi*fc1/FS;

  93. wup=2*pi*fpu/FS;wus=2*pi*fcu/FS;

  95. Bt=wlp-wls;            

  97. N0=ceil(6.6*pi/Bt);

  99. N=N0+mod(N0+1,2);

  101. wc=[(wls+wlp)/2/pi,(wus+wup)/2/pi];   

  103. hn=fir1(N-1,wc,hamming(N));

  105. [h1,w1] = freqz(hn,1,512,FS);

  107. figure(4);

  109. subplot(2,1,1);plot(hn);axis tight;xlabel('n');ylabel('幅度');title('hamming窗帶通時域圖');

  111. subplot(2,1,2);plot(w1,abs(h1));title('hamming窗帶通頻譜圖');axis ([700 900 -0.5 1]);xlabel('f/Hz');ylabel('耗損(dB)');

  113. pause(1);

  115. %#################################################################################################################

  117. %對濾波后的信號進(jìn)行分析變換

  119. X=conv(hn,x);

  121. n2=length(X);      

  123. n2=2^nextpow2(n2);

  125. f=FS*(0:n2/2-1)/n2;  

  127. figure(5);subplot(2,1,1);plot(X);title('濾波后的信號時域圖');axis tight;xlabel('f/Hz');ylabel('幅度');

  129. X1=fft(X,n2);subplot(2,1,2);plot(f,abs(X1(1:n2/2))); axis tight;xlabel('f/Hz');ylabel('幅度');title('濾波后的信號頻譜');

  131. sound(X,FS);
復(fù)制代碼

完整的Word格式文檔51黑下載地址：
Hamming窗帶通濾波器.docx (82.7 KB, 下載次數(shù): 20)





歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1