通信原理编码程序

一、实验目的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.3L=length(X);%AMIj=-1;%-1和1交替变换fori=1:LifX(i)==1y1(i)=-j;j=y1(i);elsey1(i)=X(i);endend%HDB3count=0%计算连"0"的个数y2=y1;sign=0;%标志V码极性交替V=zeros(1,L);%存储破坏脉冲V的位置B=zeros(1,L);%存储调节脉冲B的位置fori=1:Lify1(i)==0count=count+1;ifcount==4ifi>4y2(i)=y2(i-4);elsey2(i)=1;endV(i)=y2(i);ify2(i)==signy2(i)=-y2(i);y2(i-3)=y2(i);B(i-3)=y2(i-3);V(i)=y2(i);fork=i+1:Ly2(k)=-y2(k);endendsign=y2(i);endelsecount=0;endend%Milery3=zeros(1,2*L);sign1=0;%"10"和"01"交替表示sign0=1;%计算连"0"的个数,1个"0"sign0=1,多个"0"sign0=0fori=1:LifX(i)==1y3(2*i-1)=sign1;y3(2*i)=abs(1-sign1);sign1=abs(1-sign1);sign0=1;elseifsign0==1y3(2*i-1)=sign1;y3(2*i)=sign1;sign0=0;elsesign1=abs(1-sign1);y3(2*i-1)=sign1;y3(2*i)=sign1;sign0=0;endendend%波形采样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);fori=1:L*fsifmod(i-1,8)<=5&mod(i-1,8)>=2%50%脉冲宽度y(i)=1;elsey(i)=0;endX_out(i)=X(ceil(i/8));y1_out(i)=y(i)*y1(ceil(i/8));%归零码y2_out(i)=y(i)*y2(ceil(i/8));y3_out(i)=y3(ceil(i/4));end%画出波形figure(1)plot(X_out(1:160));axis([1,160,-0.2,1.2]);title('原始序列');figure(2)subplot(311);plot(y1_out(1:160));axis([1,160,-1.2,1.2]);title('归零AMI码');subplot(312);plot(y2_out(1:160));axis([1,160,-1.2,1.2]);title('HDB3码')subplot(313);plot(y3_out(1:160));axis([1,160,-0.2,1.2]);title('Miler码');%功率谱密度C_y1=xcorr(y1_out);%求自相关函数C_y2=xcorr(y2_out);C_y3=xcorr(y3_out);S_y1=abs(fft(C_y1));%自相关函数作傅里叶变换S_y2=abs(fft(C_y2));S_y3=abs(fft(C_y3));figure(3)subplot(311);plot(10*log10(S_y1));ti