matlab哈明窗带阻,MATLAB数字滤波器程序 Hamming窗带通滤波器

%hamming带通

%产生语音信号

clc;

Fs=48000;

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

x=x(:,1);

sound(x,FS);

%频谱分析

n=length(x);

n=2^nextpow2(n); %选取变换的点数

t=(0:(n-1))/FS;%计算音频信号的长度

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

figure(1);

subplot(2,1,1); plot(t,x); axis tight; title('语音信号时域图'); xlabel('t/s');ylabel('幅度')

y=fft(x,n);          %对n点进行傅里叶变换到频域

f=FS*(0:n/2-1)/n;             % 对应点的频率

subplot(2,1,2);

plot(f,abs(y(1:n/2)));axis tight;xlabel('f/Hz');ylabel('幅度');title('语音信号频域图');

pause(2.5);

%###########################################################################################################

%产生噪声信号

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

Nois