1 简介

1. 样本音频信号的频谱分析

2. 对音频信号的滤波

3. 对音频信号进行PCM编码

4. 对编码后的信号进行汉明编码

5. 对编码后的信号进行qpsk调制

6. 发送信号进入AWGN信道

7. 对接收到的信号进行qpsk解调

8. 信道译码-汉明纠错和译码

9. 对音频信号进行pcm译码并播放

2 部分代码

clc;clear all;
%/*************************读入音频文件*******************************/
long=input('想处理的音乐的长度(推荐值200 000以下,太长会很慢):');
disp('**********请欣赏发送的的音乐*********')
[x,fs]=audioread('高山流水.wav',[1 long]);
sound(x,fs);
X=fft(x,long);
magX=abs(X);
angX=angle(X);
figure;%画图
subplot(321);plot(x);title('原始信号波形');
subplot(322);plot(abs(X)); title('原始信号频谱');
%/*************************对音频信号滤波*******************************/
N=5;wc=4000/6000;
[b,a]=butter(N,wc);%用巴特沃斯滤波器进行滤波
X=fft(x);
subplot(323);plot(x);title('滤波前信号的波形');
subplot(324);plot(abs(X));title('滤波前信号的频谱');
y=filter(b,a,x);
Y=fft(y);
subplot(325);plot(y);title('IIR滤波后信号的波形');
subplot(326);plot(abs(Y));title('IIR滤波后信号的频谱');
%/*************************对音频信号进行pcm编码*************************/
L=length(y);
pcmy=pcmcoder(y,L,long);
%/*************************对编码后的信号进行信道编码-汉明编码************/
s=pcmy;
L=length(s);
h=zeros(1,L+L/4*3);%产生编码序列
N=L+L/4*3;
h=hanmincoder(s,L);
%/*************************对编码后的信号进行qpsk调制********************/
sig=qpskcoder(N,h);
%/*************************发送信号进入AWGN信道**************************/
%加入噪声
n=noise(N,L);
r=sig+n;%检测器的输入模块
%/*************************对接收到的信号进行qpsk解调********************/
rr=qpskdecoder(N,r);
%/*************************信道译码-汉明纠错和译码***********************/
y=hanmingdecoder(rr,L) ;
%/*************************对音频信号进行pcm译码************************/
pcmyout=pcmdecoder(y,long);
disp('**********请欣赏接收到的音乐*********')
sound(pcmyout,fs);
disp('**********程序结束,谢谢*********')

3 仿真结果

【滤波器】基于matlab实现语音通信仿真附报告_路径规划

【滤波器】基于matlab实现语音通信仿真附报告_matlab代码_02编辑

4 参考文献

[1]刘洲洲. 基于原型滤波器的语音信号滤波分析与仿真实现[J]. 电子设计工程, 2013, 21(5):4.

博主简介:擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真,相关matlab代码问题可私信交流。

部分理论引用网络文献,若有侵权联系博主删除。

【滤波器】基于matlab实现语音通信仿真附报告_路径规划_03

【滤波器】基于matlab实现语音通信仿真附报告_路径规划_04编辑