通信原理编码程序.docVIP

一、实验目的 1、熟悉Matlab环境并编写函数实现功能； 2、掌握数字基带传输系统基本知识。 二、实验内容 1、产生一个长度为1000的二进制随机序列，“0”的概率为0.7,”1”的概率为0.3； 下面是产生的随机序列前20个符号对应的波形： 2、对上述数据进行归零AMI编码，脉冲宽度为符号宽度的50%，波形采样率为符号率的8倍，话出前20个符号对应的波形（同时给出前20位信源序列）； 3、改用HDB3码，画出前20个符号对应的波形； 4、改用密勒码，画出前20个符号对应的波形； 5、分别对上述1000个符号的波形进行功率谱估计，画出功率谱； 6、改变信源“0”的概率，观察AMI码的功率谱变化情况。 原序列全为1时，功率谱为 可以看到原序列中”0”的个数越少，在1/2码速率处出现的尖峰会越高，也就说能量会更集中。 三、实验代码 clear; fs=8; %采样率倍数 X=ceil(rand(1,1000)-0.5) %长度为1000的序列,0的概率为0.7,1的概率为0.3 L=length(X); %AMI j=-1; %-1和1交替变换 for i=1:L if X(i)==1 y1(i)=-j; j=y1(i); else y1(i)=X(i); end end %HDB3 count=0 %计算连0的个数 y2=y1; sign=0; %标志V码极性交替 V=zeros(1,L); %存储破坏脉冲V的位置 B=zeros(1,L); %存储调节脉冲B的位置 for i=1:L if y1(i)==0 count=count+1; if count==4 if i4 y2(i)=y2(i-4); else y2(i)=1; end V(i)=y2(i); if y2(i)==sign y2(i)=-y2(i); y2(i-3)=y2(i); B(i-3)=y2(i-3); V(i)=y2(i); for k=i+1:L y2(k)=-y2(k); end end sign=y2(i); end else count=0; end end %Miler y3=zeros(1,2*L); sign1=0; %10和01交替表示 sign0=1; %计算连0的个数,1个0 sign0=1,多个0 sign0=0 for i=1:L if X(i)==1 y3(2*i-1)=sign1; y3(2*i)=abs(1-sign1); sign1=abs(1-sign1); sign0=1; else if sign0==1 y3(2*i-1)=sign1; y3(2*i)=sign1; sign0=0; else sign1=abs(1-sign1); y3(2*i-1)=sign1; y3(2*i)=sign1; sign0=0; end end end %波形采样 y=zeros(1,L*fs); X_out=zeros(1,L*fs); y1_out=zeros(1,L*fs); y2_out=zeros(1,L*fs); y3_out=zeros(1,L*fs); for i=1:L*fs if mod(i-1,8)=5 mod(i-1,8)=2 % 50%脉冲宽度 y(i)=1; else y(i)=0; end X_out(i)=X(ceil(i/8)); y1_out(i)=y(i)*y1(ceil(i/8)); %归零码 y2_out