数字信号处理3.docx
- 文档编号:24707925
- 上传时间:2023-05-31
- 格式:DOCX
- 页数:17
- 大小:284.33KB
数字信号处理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);
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理