数字信号处理3.docx
《数字信号处理3.docx》由会员分享，可在线阅读，更多相关《数字信号处理3.docx（17页珍藏版）》请在冰点文库上搜索。
数字信号处理3
低通clc;
%音频信号生成的.wav文件导入MATLAB
x=wavread（'2.wav'）;%读一个给定的音频信号
N=length（x）;
n=1:
N;
%画出音频信号的时域波形
plot（n,x）;
title（'导入的音频信号的时域波形'）;
xlabel（'数据个数'）;
ylabel（'量化幅度'）;
t=3.93;
fs=N/t;%采样频率
sound（x,fs）;
pause
figure;
%用FFT计算此音频信号的频谱
X=fft（x）;
f=0:
fs/N:
fs/2-1/N;
subplot（3,1,1）;
plot（f,abs（X（1:
（N+1）/2）））;
title（'原音频信号的频谱'）;
xlabel（'实际频率（HZ）'）;
ylabel（'幅度响应'）;
%设计一个低通的IIR滤波器
Wp1=500;%通带的截频，单位HZ
Ws1=1000;%阻带的截频，单位HZ
Rp1=0.5;%通带衰减，单位DB
Rs1=40;%通带衰减，单位DB
[n1,Wn1]=cheb1ord（Wp1/（fs/2）,Ws1/（fs/2）,Rp1,Rs1）;%获得切比雪夫滤波器滤波阶数n1，和截止频率Wn1
[b1,a1]=cheby1（n1,Rp1,Wn1）;%确定切比雪夫滤波器的系统函数H_lp（s）的分子多项式b1，和分母多项式a1
[H1,W1]=freqz（b1,a1）;%求系统的频率响应H（wjw）;
%画出该低通的IIR滤波器幅频响应
subplot（3,1,2）;
plot（W1*fs/（2*pi）,abs（H1））;
xlabel（'实际频率（Hz）'）;
ylabel（'滤波器的幅度响应'）;
y1=filter（b1,a1,x）;%用此低通滤波器对音频信号进行滤波
Y1=fft（y1）;
subplot（3,1,3）;
plot（f,abs（Y1（1:
（N+1）/2）））;
xlabel（'实际频率（HZ）'）;
ylabel（'滤波后的音频信号幅度响应'）;
pause
figure;
subplot（2,1,1）;
plot（n,x）;
title（'导入的音频信号的时域波形'）;
xlabel（'数据个数N'）;
ylabel（'量化幅度'）;
sound（x,fs）;
pause
subplot（2,1,2）;
plot（n,y1）;
title（'经过低通IIR型滤波器滤波后的音频信号的时域波形'）;
xlabel（'数据个数N'）;
ylabel（'量化幅度'）;
sound（y1,fs）;
高通
clc;
%音频信号生成的.wav文件导入MATLAB
x=wavread（'2.wav'）;%读一个给定的音频信号
N=length（x）;
n=1:
N;
%画出音频信号的时域波形
plot（n,x）;
title（'导入的音频信号的时域波形'）;
xlabel（'数据个数'）;
ylabel（'量化幅度'）;
t=3.93;
fs=N/t;%采样频率
sound（x,fs）;
pause
figure;
%用FFT计算此音频信号的频谱
X=fft（x）;
f=0:
fs/N:
fs/2-1/N;
subplot（3,1,1）;
plot（f,abs（X（1:
（N+1）/2）））;
title（'原音频信号的频谱'）;
xlabel（'实际频率（HZ）'）;
ylabel（'幅度响应'）;
%设计一个低通的IIR滤波器
Wp1=500;%通带的截频，单位HZ
Ws1=1000;%阻带的截频，单位HZ
Rp1=0.5;%通带衰减，单位DB
Rs1=40;%通带衰减，单位DB
[n1,Wn1]=cheb1ord（Wp1/（fs/2）,Ws1/（fs/2）,Rp1,Rs1）;%获得切比雪夫滤波器滤波阶数n1，和截止频率Wn1
[b1,a1]=cheby1（n1,Rp1,Wn1,'high'）;%确定切比雪夫滤波器的系统函数H_lp（s）的分子多项式b1，和分母多项式a1
[H1,W1]=freqz（b1,a1）;%求系统的频率响应H（wjw）;
%画出该低通的IIR滤波器幅频响应
subplot（3,1,2）;
plot（W1*fs/（2*pi）,abs（H1））;
xlabel（'实际频率（Hz）'）;
ylabel（'滤波器的幅度响应'）;
y1=filter（b1,a1,x）;%用此低通滤波器对音频信号进行滤波
Y1=fft（y1）;
subplot（3,1,3）;
plot（f,abs（Y1（1:
（N+1）/2）））;
xlabel（'实际频率（HZ）'）;
ylabel（'滤波后的音频信号幅度响应'）;
pause
figure;
subplot（2,1,1）;
plot（n,x）;
title（'导入的音频信号的时域波形'）;
xlabel（'数据个数N'）;
ylabel（'量化幅度'）;
sound（x,fs）;
pause
subplot（2,1,2）;
plot（n,y1）;
title（'经过低通IIR型滤波器滤波后的音频信号的时域波形'）;
xlabel（'数据个数N'）;
ylabel（'量化幅度'）;
sound（y1,fs）;
带通
clc;
%音频信号生成的.wav文件导入MATLAB
x=wavread（'2.wav'）;%读一个给定的音频信号
N=length（x）;
n=1:
N;
%画出音频信号的时域波形
plot（n,x）;
title（'导入的音频信号的时域波形'）;
xlabel（'数据个数'）;
ylabel（'量化幅度'）;
t=3.93;
fs=N/t;%采样频率
sound（x,fs）;
pause
figure;
%用FFT计算此音频信号的频谱
X=fft（x）;
f=0:
fs/N:
fs/2-1/N;
subplot（3,1,1）;
plot（f,abs（X（1:
（N+1）/2）））;
title（'原音频信号的频谱'）;
xlabel（'实际频率（HZ）'）;
ylabel（'幅度响应'）;
%设计一个低通的IIR滤波器
Wp1=[940,1560];%通带的截频，单位HZ
Ws1=[575,2000];%阻带的截频，单位HZ
Rp1=0.1;%通带衰减，单位DB
Rs1=40;%通带衰减，单位DB
[n1,Wn1]=cheb1ord（Wp1/（fs/2）,Ws1/（fs/2）,Rp1,Rs1）;%获得切比雪夫滤波器滤波阶数n1，和截止频率Wn1
[b1,a1]=cheby1（n1,Rp1,Wn1）;%确定切比雪夫滤波器的系统函数H_lp（s）的分子多项式b1，和分母多项式a1
[H1,W1]=freqz（b1,a1）;%求系统的频率响应H（wjw）;
%画出该低通的IIR滤波器幅频响应
subplot（3,1,2）;
plot（W1*fs/（2*pi）,abs（H1））;
xlabel（'实际频率（Hz）'）;
ylabel（'滤波器的幅度响应'）;
y1=filter（b1,a1,x）;%用此低通滤波器对音频信号进行滤波
Y1=fft（y1）;
subplot（3,1,3）;
plot（f,abs（Y1（1:
（N+1）/2）））;
xlabel（'实际频率（HZ）'）;
ylabel（'滤波后的音频信号幅度响应'）;
pause
figure;
subplot（2,1,1）;
plot（n,x）;
title（'导入的音频信号的时域波形'）;
xlabel（'数据个数N'）;
ylabel（'量化幅度'）;
sound（x,fs）;
pause
subplot（2,1,2）;
plot（n,y1）;
title（'经过低通IIR型滤波器滤波后的音频信号的时域波形'）;
xlabel（'数据个数N'）;
ylabel（'量化幅度'）;
sound（y1,fs）;
带阻
clc;
%音频信号生成的.wav文件导入MATLAB
x=wavread（'2.wav'）;%读一个给定的音频信号
N=length（x）;
n=1:
N;
%画出音频信号的时域波形
plot（n,x）;
title（'导入的音频信号的时域波形'）;
xlabel（'数据个数'）;
ylabel（'量化幅度'）;
t=3.93;
fs=N/t;%采样频率
sound（x,fs）;
pause
figure;
%用FFT计算此音频信号的频谱
X=fft（x）;
f=0:
fs/N:
fs/2-1/N;
subplot（3,1,1）;
plot（f,abs（X（1:
（N+1）/2）））;
title（'原音频信号的频谱'）;
xlabel（'实际频率（HZ）'）;
ylabel（'幅度响应'）;
%设计一个低通的IIR滤波器
Wp1=[940,1560];%通带的截频，单位HZ
Ws1=[575,2000];%阻带的截频，单位HZ
Rp1=0.1;%通带衰减，单位DB
Rs1=40;%通带衰减，单位DB
[n1,Wn1]=cheb1ord（Wp1/（fs/2）,Ws1/（fs/2）,Rp1,Rs1）;%获得切比雪夫滤波器滤波阶数n1，和截止频率Wn1
[b1,a1]=cheby1（n1,Rp1,Wn1,'stop'）;%确定切比雪夫滤波器的系统函数H_lp（s）的分子多项式b1，和分母多项式a1
[H1,W1]=freqz（b1,a1）;%求系统的频率响应H（wjw）;
%画出该低通的IIR滤波器幅频响应
subplot（3,1,2）;
plot（W1*fs/（2*pi）,abs（H1））;
xlabel（'实际频率（Hz）'）;
ylabel（'滤波器的幅度响应'）;
y1=filter（b1,a1,x）;%用此低通滤波器对音频信号进行滤波
Y1=fft（y1）;
subplot（3,1,3）;
plot（f,abs（Y1（1:
（N+1）/2）））;
xlabel（'实际频率（HZ）'）;
ylabel（'滤波后的音频信号幅度响应'）;
pause
figure;
subplot（2,1,1）;
plot（n,x）;
title（'导入的音频信号的时域波形'）;
xlabel（'数据个数N'）;
ylabel（'量化幅度'）;
sound（x,fs）;
pause
subplot（2,1,2）;
plot（n,y1）;
title（'经过低通IIR型滤波器滤波后的音频信号的时域波形'）;
xlabel（'数据个数N'）;
ylabel（'量化幅度'）;
sound（y1,fs）;
升级会员 