基于MATLAB语音信号采集与分析加上高频噪声并去噪声源程序仿真.pdfVIP

语音信号的综合仿真分析

%仿真说明：有些人仿真不出来并不是源代码不对，而是自己没有操作好，下面是一些操作

注意事项，希望能帮助到各位！

1、源代码中的（3.wav）音频文件是用windows系统自带录音软件录制的，在附件--娱乐

--录音机，保存文件时一定要保存二进制格式.wav。

2、.wav文件必须和源程序代码M文件保存在一起，即同一路径下，最好是根目录下，比

如都保存在D:\，仿真时将源代码中wavread(3.wav)改成wavread(D:\3.wav)即

可。本仿真是保存在桌面实现的！

3、本仿真是自己的录音，若要实现仿真必须自己录制.wav音频文件。

一、综合仿真源程序代码：

[y,fs,bits]=wavread(3.wav);

%sound(y,fs,bits);

Y=fft(y,40000);%采样点数40000

figure(1);

subplot(211);plot(y);title(原始信号波形);

subplot(212);plot(abs(Y));title(原始信号频谱);

fc=20000;%载波频率

y1=modulate(y,fc,fs,fm);%对原语音信号调制

Y1=fft(y1,40000);%采样点数40000

figure(2);

subplot(211);plot(y1);title(调制后信号波形);

subplot(212);plot(abs(Y1));title(调制后信号频谱);

%sound(y1);

%调制信号+原始信号波形

y2=y+y1;

Y2=fft(y2,40000);%采样点数40000

figure(3);

subplot(211);plot(y2);title(调制信号+原始信号波形);

subplot(212);plot(abs(Y2));title(调制信号+原始信号频谱);

%sound(y2);

figure(4);

%用双线性变换法设计的低通滤波器

fp=1000;fc=1200;As=100;Ap=1;fs1=22050;

wc=2*fc/fs1;wp=2*fp/fs1;

[n,wn]=ellipord(wp,wc,Ap,As);

[b,a]=ellip(n,Ap,As,wn);

freqz(b,a,512,fs1);

x=filter(b,a,y2);%对要进行低通滤波

X=fft(x,40000);

figure(5);

subplot(211);plot(x);title(低通滤波后信号波形);

subplot(212);plot(abs(X));title(低通滤波后信号频谱);

sound(x,40000);

二、仿真结果说明

①录制一段自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波

形和频谱图；

原始信号波形

1

0.5

0

-0.5

00.511.522.5

x105

原始信号频谱

500

400

300

200

100

0

00.511.522.533.54