x_fianl.m

大学期间

%随机数字基带信号
x=round(rand(1,30));%随机产生的1行100列数字基带信号
subplot(2,3,1);
impulscx(x,2);                       %调用impulscx(x,ts) ts=2函数
axis([0,30,0,1.1]);
title('随机数字基带信号')         %为所画图命名

%数字基带信号功率谱密度
x1=xcorr(x);                                         %求x的自相关函数x1
x1=fft(x1);                                         %求x的一阶傅立叶变换
subplot(2,3,2);     
plot(abs(x1));
title('数字基带信号功率谱密度')              %数字基带信号画出功率谱密度图

%双相编码
y=[x,x];
for i=1:length(x)
    if x(i)==1  
            y(2*i-1)=1;        %遇到“1”则转换成10
            y(2*i)=0;
    else 
            y(2*i-1)=0;        %遇到“1”则转换成01
            y(2*i)=1;
         
    end
end
subplot(2,3,3);
impulscx(y,1);                       %调用impulscx(x,ts) ts=2函数
axis([0,30,0,1.1]);
title('双相编码');

%双相编码后功率谱密度
yp=xcorr(y);                                         %求y的自相关函数yp
yp=fft(yp);                                         %求y的一阶傅立叶变换
subplot(2,3,4);  
plot(abs(yp));
title('双相编码后功率谱密度');              %数字基带信号画出功率谱密度图

%双相解码
xd=[x];
for i=1:length(y)/2
    if (y(2*i-1)==1 && y(2*i)==0) xd(i)=1
    elseif (y(2*i-1)==0 && y(2*i)==1) xd(i)=0
    end
end
subplot(2,3,5);
impulscx(xd,2);
axis([0,30,0,1.1]);
title('双相解码') ;                        %画出解码后的图