通信系统课群综合训练与设计-PCM、CMI、汉明码、ASK、衰落信道

断更了很久，一方面是因为之后开始准备期末考试和课程设计，另一方面也是因为自己因为放假堕落了（好吧，这个是主要原因）。
课程设计做了挺久，整合了很多，参考了不少内容才做出来。服务学弟学妹了，就把我的这个发出来好了。
其中，信源为自己构造的时间函数，数字化方式为PCM编码，基带码为CMI码，信道码为汉明码，调制方式为ASK调制与解调，信道类型为衰落信道，接收端与发送端对应。

1.main函数
%empty the processing
clc;
close all;
clear all;
 
%original signal
clear;
fs=8000;%抽样频率
fc=100;
dt=1/fs;
Ts=dt;
t=0:dt:0.02-dt;
y=0.3*cos(2*pi*fc*t)+0.5*sin(4*pi*fc*t)+0.2*sin(8*pi*fc*t);%输入模拟信号
figure(1);
plot(t,y);
ylim([-2,2]);
title('源信号');
 
%PCM modulation
n=length(t);
y_PCM=PCMENC2(y,n);
figure(2);
stairs(y_PCM);
%zoom xon;
xlim([0 100]);
ylim([-0.2,1.2]);
title('PCM编码信号');
 
%CMI modulation
y_CMI = CMI_coding(y_PCM) ;
t2 = 0:Ts:41*Ts;
figure(3)
stairs(y_CMI);
xlim([0 200]);
ylim([-0.2 1.2]);
%zoom xon;axis([0 length(y_CMI) -0.2 1.2]);
title('CMI编码信号');
 
%Hamming modulation
y_Hamming=encode(y_CMI,7,4,'hamming') ;
figure(4)              
stairs(y_Hamming);
xlim([0 200]);
ylim([-0.2 1.2]);
%zoom xon;axis([0 length(y_Hamming) -0.2 1.2]);
title('Hamming编码');
 
%ASK modulaiton and rayleigh channel
y_DASK = ASK(y_Hamming);
 
%Hamming demodulation
y_DHamming=decode(y_DASK,7,4,'hamming'); 
figure(8)
stairs(y_DHamming);
xlim([0 200]);
ylim([-0.2 1.2]);
%zoom xon;axis([0 length(y_DHamming) -0.2 1.2]);
title('Hammming译码输出信号（CMI）');
 
%CMI demodulation
y_DCMI = CMI_decoding(y_DHamming);
figure(9)
stairs(y_DCMI);
xlim