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

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(’ ***** ***** 程序结束，谢谢 ***** *** *’)