clc;clear all; %/ ***** ***** ***** ***** ***** 讀入音頻文件 ***** ***** ***** ***** ***** ***** */ long=input(’想處理的音樂的長度（推薦值200 000以下，太長會很慢）:’); disp(’ ***** ***** 請欣賞發(fā)送的的音樂 ***** *** *’) [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);%用巴特沃斯濾波器進(jìn)行濾波 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濾波后信號的頻譜’); %/ ***** ***** ***** ***** ***** 對音頻信號進(jìn)行pcm編碼 ***** ***** ***** ***** ***** / L=length(y); pcmy=pcmcoder(y,L,long); %/ ***** ***** ***** ***** ***** 對編碼后的信號進(jìn)行信道編碼-漢明編碼 ***** ***** **/ s=pcmy; L=length(s); h=zeros(1,L+L/4*3);%產(chǎn)生編碼序列 N=L+L/4*3; h=hanmincoder(s,L); %/ ***** ***** ***** ***** ***** 對編碼后的信號進(jìn)行qpsk調(diào)制 ***** ***** ***** ***** / sig=qpskcoder(N,h); %/ ***** ***** ***** ***** ***** 發(fā)送信號進(jìn)入AWGN信道 ***** ***** ***** ***** ***** */ %加入噪聲 n=noise(N,L); r=sig+n;%檢測器的輸入模塊 %/ ***** ***** ***** ***** ***** 對接收到的信號進(jìn)行qpsk解調(diào) ***** ***** ***** ***** / rr=qpskdecoder(N,r); %/ ***** ***** ***** ***** ***** 信道譯碼-漢明糾錯和譯碼 ***** ***** ***** ***** *** / y=hanmingdecoder(rr,L) ; %/ ***** ***** ***** ***** ***** 對音頻信號進(jìn)行pcm譯碼 ***** ***** ***** ***** *** */ pcmyout=pcmdecoder(y,long); disp(’ ***** ***** 請欣賞接收到的音樂 ***** *** *’) sound(pcmyout,fs); disp(’ ***** ***** 程序結(jié)束，謝謝 ***** *** *’)