soundread.m

数字信号处理技术 Digital Signal Processing

[y,fs,bits]=WAVREAD('dajiahao.wav');
sound(y,fs,bits);
Y=fft(y,4096);figure;showim(fftshift(Y))

fb=2000;fc=2200;as=100;ap=1;fs=22050;  % lowpass
wc=2*fc/fs;wb=2*fb/fs;
[n,wn]=ellipord(wb,wc,ap,as);
[b,a]=ellip(n,ap,as,wn);


% fb1=2400;fb2=4800;fc1=2200;fc2=5000;as=100;ap=1;fs=22050;  % bandpass
% wc1=2*fc1/fs;wc2=2*fc2/fs;wb1=2*fb1/fs;wb2=2*fb2/fs;
% [n,wn]=ellipord([wb1,wb2],[wc1,wc2],ap,as);
% [b,a]=ellip(n,ap,as,wn);

 
% fb=5500;fc=5000;as=100;ap=1;fs=22050;      % highpass
% wc=2*fc/fs;wb=2*fb/fs;
% [n,wn]=ellipord(wb,wc,ap,as);
% [b,a]=ellip(n,ap,as,wn,'high');


figure;
freqz(b,a,512,fs);
x=filter(b,a,y);
X=fft(x,4096);
figure;showim(fftshift(X));

% sum(abs(x).^2)/sum(abs(y).^2)