【滤波器】基于matlab实现语音通信仿真附报告

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 仿真结果

4 参考文献

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

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

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