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代码问题可私信交流。
部分理论引用网络文献，若有侵权联系博主删除。
文章出处登录后可见！
已经登录？立即刷新
