《数字信号处理》课程设计.docx
- 文档编号:28038431
- 上传时间:2023-07-07
- 格式:DOCX
- 页数:15
- 大小:682.69KB
《数字信号处理》课程设计.docx
《《数字信号处理》课程设计.docx》由会员分享,可在线阅读,更多相关《《数字信号处理》课程设计.docx(15页珍藏版)》请在冰豆网上搜索。
《数字信号处理》课程设计
输入声音、
Fn=3000;Fs=44100;W3DB=0.06;
y=wavread('C:
\DocumentsandSettings\Administrator\桌面\my.wav');
[y,fs,nbits]=wavread('my.wav');
sound(y,Fs);pause;
Fn=3000;Fs=44100;W3DB=0.06;
y=wavread('C:
\DocumentsandSettings\Administrator\桌面\my.wav');
>>Fn=3000;Fs=44100;W3DB=0.06
y=wavread('C:
\DocumentsandSettings\Administrator\桌面\my.wav');
W3DB=
0.0600
>>sound(y,Fs);pause;
>>
Fn=3000;Fs=44100;W3DB=0.06
y=wavread('C:
\Users\Administrator\Desktop\my.wav');
y=y(:
1);
X=fft(y,3000);
magY=abs(y);
angY=angle(y);
subplot(221);plot(y);title('原始信号波形');
subplot(222);plot(y);title('原始信号频谱');
subplot(223);plot(magY);title('原始信号幅值');
subplot(224);plot(angY);title('原始信号相位');
加噪声
N=length(y);%求出语音信号的长度
t=0:
5*pi
noise=5*cos(2*pi*3000*t);%噪声信号的函数
z=fft(noise);%快速傅里叶变换
subplot(2,1,1),plot(noise);title('噪声信号波形');
subplot(2,1,2),plot(abs(z));title('噪声信号频谱');
axis([0,10,0,5]);%xmin,xmax,ymin,ymax
叠加
《1》>>fs=44100;
x1=wavread('C:
\Users\Administrator\Desktop\my.wav');
>>t=(0:
length(y)-1)/8000;
f=fs*(0:
1023)/2048;
noise=[5*cos(2*pi*3800*t)];
>>y=y(:
1);
>>z=y'+noise;
>>subplot(2,1,1),plot(noise);title('噪声信号波形');
subplot(2,1,2),plot(abs(z));title('噪声信号频谱');
axis([0,10,0,5]);%xmin,xmax,ymin,ymax
>>sound(noise,44100);
>>sound(noise,44100);
>>sound(y,44100);
>>sound(y'+noise);
《2>>
y=wavread('C:
\Users\Administrator\Desktop\my.wav');
N=length(y);
t=(0:
length(y)-1)/8000;
f=fs*(0:
1023)/2048;
noise=[5*cos(2*pi*3800*t)];
y=y(:
1);
k=y'+noise;
Subplot(2,1,1);
plot(s);title('叠加时域波形');
S=fft(s);
Subplot(2,1,2);
plot(abs(S));title('叠加频域波形');
axis([0,10,0,5]);
sound(y'+noise,44100);
IIR低通滤波器
Ft=8000;
Fb=1000;
Fc=1200;
wb=2*pi*Fb/Ft;
wc=2*pi*Fc/Ft;
fb=2*Ft*tan(wb/2);
fc=2*Fc*tan(wc/2);
[N,wd]=buttord(wb,wc,1,100,'s');
[b,a]=butter(N,wd,'s');
[bz,az]=bilinear(b,a,0.5);
[h,w]=freqz(bz,az);
plot(w*8000*0.5/pi,abs(h));
Legend('用butter设计');
title('IIR低通滤波器');
grid;
2.4FIR数字滤波器设计
Ft=8000;
Fb=2800;
Fc=3000;
wp=2*Fb/Ft;
ws=2*Fc/Ft;
Ap=1;
As=100;
p=1-10.^(-Ap/20);
s=10.^(-As/20);
fpts=[wpws];
mag=[10];
dev=[ps];
[N,wc,beta,ftype]=kaiserord(fpts,mag,dev);
b=fir1(N,wc,kaiser(N+1,beta));
[h,w]=freqz(b,1);
plot(w/pi,abs(h));
title('FIR低通滤波器');
grid;
2.5滤波器对信号滤波
IIR低通滤波:
Fn=3000;Fs=44100;W3DB=0.06;
y=wavread('C:
\Users\Administrator\Desktop\my.wav');
N=length(y);
t=(0:
length(y)-1)/8000;
f=fs*(0:
1023)/2048;
noise=[5*cos(2*pi*3800*t)];
y=y(:
1);
k=y'+noise;
Subplot(2,2,1);plot(k);title('滤波前的时域波形');
K=fft(k);
Subplot(2,2,2);plot(abs(K));title('滤波前的频域波形');
axis([0,100,0,50]);
Ft=8000;
Fp=1000;
Fs=1200;
wp=2*pi*Fp/Ft;
ws=2*pi*Fs/Ft;
fp=2*Ft*tan(wp/2);
fs=2*Ft*tan(ws/2);
[N,wc]=buttord(wp,ws,1,50,'s');
[b,a]=butter(N,wc,'s');
[bz,az]=bilinear(b,a,0.5);
[h,w]=freqz(bz,az);
z11=filter(bz,az,k);
m11=fft(z11);
subplot(2,2,3);plot(z11,'g');title('滤波后的时域波形');
subplot(2,2,4);plot(abs(m11),'r');title('滤波后的时域波形');
axis([0,100,0,50]);
sound(z11,44100);
FIR低通滤波:
y=wavread('my.wav');
N=length(y);
t=(0:
length(y)-1)/8000;
f=fs*(0:
1023)/2048;
noise=[5*cos(2*pi*3800*t)];
y=y(:
1);
k=y'+noise;
Subplot(2,2,1);plot(k);title('滤波前的时域波形');
K=fft(k);
Subplot(2,2,2);plot(abs(K));title('滤波前的频域波形');
Ft=8000;
Fb=2800;
Fc=3000;
wb=2*Fb/Ft;
wc=2*Fc/Ft;
Ap=1;
As=100;
p=1-10.^(-Ap/20);
s=10.^(-As/20);
fpts=[wbwc];
mag=[10];
dev=[ps];
[N,wc,beta,ftype]=kaiserord(fpts,mag,dev);
b=fir1(N,wc,kaiser(N+1,beta));
[h,w]=freqz(b,1);
Ft=8000;
Fb=2800;
Fc=3000;
wb=2*Fb/Ft;
wc=2*Fc/Ft;
Ap=1;
As=100;
p=1-10.^(-Ap/20);
s=10.^(-As/20);
fpts=[wpws];
mag=[10];
dev=[ps];
[N,wc,beta,ftype]=kaiserord(fpts,mag,dev);
b=fir1(N,wc,kaiser(N+1,beta));
[h,w]=freqz(b,1);
z21=fftfilt(b,k);
subplot(2,2,3);plot(z11,'g');title('滤波后的时域波形');
subplot(2,2,4);plot(abs(m11),'r');title('滤波后的时域波形');
axis([0,100,0,50]);
sound(z21,44100);
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号处理 数字信号 处理 课程设计
![提示](https://static.bdocx.com/images/bang_tan.gif)