通信原理matalab实验汇总——二三元码的调制与解调

已于 2024-01-07 11:39:35 修改
阅读量449 收藏 10
点赞数 7
分类专栏： # 通信原理 文章标签： 通信原理 二元码 三元码 调制与解调
于 2024-01-05 09:47:26 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_63831368/article/details/135400499
版权

通信原理matalab仿真实验汇总——二三元码的调制与解调

本文主要涉及通信原理中重要的二三原码的调制与解调的matalab代码的实现。
关于通信原理还有其他文章可参考：

1、信息量、码元、比特、码元速率、信息速率详细解析——实例分析
2、模拟系统的AM信号的调制与频域分析
3、AM解调、DSB、VSB的调制与解调
4、模拟非线性调制——相位调制PM、频率调制FM、NBFM窄带调频
5、低通抽样与带通抽样
6、简单二元码的编码与功率谱
7、数字双向码、密勒码、传号反转（CMI）码、AMI、HDB3
8、ASK（幅度移键控）、FSK（频率移键控）、PSK（相位移键控）的matalab仿真

文章目录

一、 二元码的编码

1.1 所有二元码的matalab代码实现

%-----二元码的编码----
% 定义二进制信码
code = [1 1 1 0 1 0 0 0 1 1 0 0 1 0];
% 创建时间轴数据
t = 0:0.5:(length(code)*100-1)*0.5;
 
length_t = length(t);
%查看t的长度
% disp(length_t);
 
% 单极性非归零码：1的时候为1，0的时候为0
NRZ = [];%存储单极性非归零码
for i = 1:length(code)
    if code(i) == 1
        %使用方括号拼接两个向量[]
        NRZ = [NRZ ones(1, 100)];%为了和t长度一样，一次的跨度也要为100
    else
        NRZ = [NRZ zeros(1, 100)];
    end
end
 
% 双极性非归零码:1的时候为1，0的时候为-1
polar_NRZ = [];
for i = 1:length(code)
    if code(i) == 1
        polar_NRZ = [polar_NRZ ones(1, 100)];
    else
        polar_NRZ = [polar_NRZ -1*ones(1, 100)];
    end
end
 
% 单极性归零码：1的时候分一半为0（归零），0的时候为0
RZ = [];
for i = 1:length(code)
    if code(i) == 1
        RZ = [RZ ones(1, 50) zeros(1, 50)];%从高电平到低电平过渡
    else
        RZ = [RZ zeros(1, 100)];
    end
end
 
% 双极性归零：在双极性非归零码的基础上，给1和-1均加入过渡值,即一半归零
polar_RZ = [];
for i = 1:length(code)
    if code(i) == 1
        polar_RZ = [polar_RZ ones(1, 50) zeros(1, 50)];
    else
        polar_RZ = [polar_RZ -1*ones(1, 50) zeros(1, 50)];
    end
end
 
%传号差分码：波形跳变取值为1
NRZ_M = [];
previous_bit = 0; % 初始值为0，用于表示前一个比特的数值
for i = 1:length(code)
    if code(i) == 1
        bit_value = xor(previous_bit, 1); % 当前位与前一位异或运算,相同为0，不同为1
    else
        bit_value = previous_bit;
    end
    NRZ_M = [NRZ_M ones(1, 100)*bit_value];
    
    previous_bit = bit_value; % 更新前一位比特的数值
end
 
 
% 空号差分码：波形跳变取值为0
NRZ_S = [];
previous_bit = 0; % 初始值为0，用于表示前一个比特的数值
for i = 1:length(code)
    if code(i) == 1
        bit_value = previous_bit;
    else   
        bit_value = xor(previous_bit, 1); % 当前位与前一位异或运算
    end
    NRZ_S = [NRZ_S ones(1, 100)*bit_value];
    
    previous_bit = bit_value; % 更新前一位比特的数值
end
 
%数字双向码
%0：用10表示；1：用01表示
manchester_code  = [];
for i = 1:length(code)
    if code(i) == 1
        manchester_code = [manchester_code -1*ones(1,50) ones(1,50)];
    else
        manchester_code = [manchester_code ones(1,50) -1*ones(1,50)];
    end
end
 
 
%密勒码
miler_code = [];
flag = 1;
count_0 = 0;%存储连0的个数
for i=1:length(code)
    if code(i) == 1
        count_0 = 0;%清空连0个数
        miler_code = [miler_code -1*flag*ones(1,50) flag*ones(1,50)]
        if flag == 1
            flag = -1;
        else
            flag = 1;
        end
    else  %等于0
        count_0 = count_0+1;
        if(count_0 < 2)%不连零的情况
            miler_code = [miler_code -flag*ones(1,100)]
        else%连零
            miler_code = [miler_code flag*ones(1,100)]
            flag = -flag;
        end
    end
end
            
 
%传号反转码
%1交替的用00和11表示，0固定用01表示
CMI=[];
flag = 1; % 初始极性为正
for i = 1:length(code)
    if code(i) == 1
        CMI = [CMI ones(1, 100)*flag];
       if flag==1
           flag=-1;
       else
           flag=1;
       end
    else
        CMI = [CMI ones(1, 50)*-1 ones(1, 50)]; % 0：交替极性
    end
end
 
figure(1)
subplot(2,1,1)
plot(t, NRZ, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -0.5 1.5])
title('单极性非归零码')
xlabel('时间')
ylabel('幅度')
grid on;
 
subplot(2,1,2)
plot(t, polar_NRZ, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -1.5 1.5])
title('双极性非归零码')
xlabel('时间')
ylabel('幅度')
grid on;

figure(2)
subplot(2,1,1)
plot(t, RZ, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -0.5 1.5])
title('单极性归零码')
xlabel('时间')
ylabel('幅度')
grid on;

subplot(2,1,2)
plot(t, polar_RZ, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -1.5 1.5])
title('双极性归零码')
xlabel('时间')
ylabel('幅度')
grid on;
 
figure(3)
subplot(2,1,1)
plot(t, NRZ_M, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -0.5 1.5])
title('传号差分码')
xlabel('时间')
ylabel('幅度')
grid on;
 
subplot(2,1,2)
plot(t, NRZ_S, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -1 1.5])
title('空号差分码')
xlabel('时间')
ylabel('幅度')
grid on;
 
figure(3)
subplot(2,1,1)
plot(t, manchester_code, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -1 1.5])
title('数字双相码')
xlabel('时间')
ylabel('幅度')
grid on;
 
subplot(2,1,2)
plot(t, miler_code, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -1 1.5])
title('密勒码')
xlabel('时间')
ylabel('幅度')
grid on;
 
figure(4)
subplot(1,1,1)
plot(t, CMI, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -1 1.5])
title('传号反转码')
xlabel('时间')
ylabel('幅度')
grid on;

1.2实现效果图：

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

1.3 二元码编码

  1. 初始化：

    • 定义了一个二进制信码code。
    • 创建了时间轴数据t。

  2. 单极性非归零码（NRZ）：

    • 当码值为1时，信号为高电平；当码值为0时，信号为低电平。

  3. 双极性非归零码（polar_NRZ）：

    • 与单极性非归零码类似，但当码值为0时，信号为负电平。

  4. 单极性归零码（RZ）：

    • 当码值为1时，信号从高电平到低电平过渡；当码值为0时，信号为低电平。

  5. 双极性归零码（polar_RZ）：

    • 在双极性非归零码的基础上，对1和-1的信号添加了过渡值，即一半为低电平。

  6. 传号差分码（NRZ_M）：

    • 与前一个比特位的值相反的比特位会导致信号波形的跳变。

  7. 空号差分码（NRZ_S）：

    • 当码值为0时，信号波形保持不变；当码值为1时，信号波形跳变。

  8. 数字双向码（manchester_code）：

    • 当码值为1时，信号波形从负电平跳变到正电平；当码值为0时，信号波形从正电平跳变到负电平。

  9. 密勒码（miler_code）：

    • 当连续两个或多个码值为0时，信号波形会有特定的跳变。

  10. 传号反转码（CMI）：

  • 对于连续的1，信号保持连续的电平；对于0，信号波形交替变化。
  1. 绘图：
  • 使用MATLAB的plot函数绘制了每种码型的波形。

二、三元码的编码

2.1 三元码编码matalab代码

% 三元码的编码AMI码，HDB3码，B6ZS码
clc
clear
% 定义二进制信码
code = [1 1 0 1 0 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0 0 1];
% 创建时间轴数据
t = 0:0.5:(length(code)*100-1)*0.5;
 
 
% 单极性非归零码：1的时候为1，0的时候为0
NRZ3 = [];%存储单极性非归零码
for i = 1:length(code)
    if code(i) == 1
        %使用方括号拼接两个向量[]
        NRZ3 = [NRZ3 ones(1, 100)];%为了和t长度一样，一次的跨度也要为100
    else
        NRZ3 = [NRZ3 zeros(1, 100)];
    end
end
 
%AMI码
%二进制码0用0电平表示
%二进制码1交替地用＋1和－1的半占空归零码表示
AMI=[];
flag = 1; % 初始极性为正
for i = 1:length(code)
    if code(i) == 1
        AMI = [AMI ones(1, 50)*flag zeros(1,50)]
        if flag==1
           flag=-1;
        else
           flag=1;
        end
    else
        AMI = [AMI zeros(1,100)]
    end
end
 
%HDB3码
HDB3=[];
flag = 1; % 初始极性为正
count_0 = 0;%记录连0个数
count2V_1 = 0;%计算两个V之间1的个数，用于确定是B00V型还是000V型
count_1 = 0; %计算总的1的个数，用于确定B的极性
 
for i = 1:length(code)
    if code(i) == 1   %当不连0时，与AMI一样
         HDB3 = [HDB3 ones(1, 50)*flag zeros(1,50)]
         count_0 = 0; %连0个数清0,开启新一轮计数
         count_1 = count_1+1;%总的1的个数加1
         count2V_1 = count2V_1 + 1;%两个V之间1的个数加1
        if flag==1
           flag=-1;
        else
           flag=1;
        end
    else  %出现了0
        count_0 = count_0+1;
        if(count_0 < 4)  %仍然和AMI码一样
            HDB3 = [HDB3 zeros(1,100)]
        else   %出现了四个连零
            if mod(count2V_1 ,2 ) == 1 %奇数个，为000V类型,V的极性跟着上一个B走(10)
                HDB3 = [HDB3 ones(1,50) zeros(1,50)]
                count_0 = 0;  %连零个数清零
                count2V_1 = 0 ;%两个V之间1的个数清零，总的1的个数不清零
            else %偶数个，为B00V型,V的极性跟着上一个B走(-10)
                     if mod(count_1 ,2) == 0 %偶数个1，B的极性轮到1，V跟着B走
                         HDB3( (i-3)*100-99:(i-3)*100-50 ) = 1*ones(1,50)  %B
                         HDB3 = [HDB3 ones(1,50) zeros(1,50)]              %V
                         flag = -flag;
                         count_0 = 0;  %连零个数清零，开启新一轮连0计数
                         count2V_1 = 0 %两个V之间1的个数清零，总的1的个数不清零
                         count_1 = count_1 + 1;
                     else %奇数个1，B的极性轮到-1，V跟着B走
                         HDB3( (i-3)*100-99:(i-3)*100-50 ) = -1*ones(1,50)  %B
                         HDB3 = [HDB3 -1*ones(1,50) zeros(1,50) ]           %V
                         flag = -flag;
                         count_0 = 0;
                         count2V_1 = 0;
                         count_1 = count_1 +1;
                 end
            end
        end
    end
end
 
%BNZS码
BNZS = []
count_0=0;
count_1=0;
flag =1;
for i = 1:length(code)
    if code(i) == 1  %当不连0时，与AMI一样
         BNZS = [BNZS ones(1, 50)*flag zeros(1,50)]
         count_0 = 0;        %连0个数清0
         count_1 = count_1+1;%总的1的个数加1
        if flag==1
           flag=-1;
        else
           flag=1;
        end
    else  %出现了0
        count_0 = count_0+1;
        if(count_0 < 3)
            BNZS = [BNZS zeros(1,100)]
        else %出现了3个连零，用0vb，接下来判断VB的极性
            if mod(count_1,2)==1 %奇数个1，V的极性跟着前一个B走，为正，00->10
                BNZS( (i-1)*100-99:(i-1)*100-50 ) = 1*ones(1,50)%V
                BNZS = [BNZS -1*ones(1,50) zeros(1,50)]         %B
                count_1 = count_1 +1;
                count_0=0;
            else%偶数个1
                BNZS( (i-1)*100-99:(i-1)*100-50 ) = -1*ones(1,50)
                BNZS = [BNZS ones(1,50) zeros(1,50)]
                count_1 = count_1 +1;
                count_0=0;
            end
        end
    end
end
 
figure(1)
subplot(2,1,1)
plot(t, NRZ3, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -1 1.5])
title('单极性非归零码')
xlabel('时间')
ylabel('幅度')
grid on;
 
subplot(2,1,2)
plot(t, AMI, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -1 1.5])
title('AMI码')
xlabel('时间')
ylabel('幅度')
grid on;
 
figure(2)
subplot(2,1,1)
plot(t, HDB3, 'LineWidth', 2)
%stairs(0:length(HDB3)-1,HDB3);
axis([0 (length(code)*100-1)*0.5 -1 1.5])
title('HDB3码')
xlabel('时间')
ylabel('幅度')
grid on;
 
 
subplot(2,1,2)
plot(t, BNZS, 'LineWidth', 2)
axis([0 (length(code)*100-1)*0.5 -1 1.5])
title('BNZS码')
xlabel('时间')
ylabel('幅度')
grid on;

2.2 实现效果

在这里插入图片描述
在这里插入图片描述

2.3 三元码编步骤

用于生成和绘制三种三元码的波形：AMI码、HDB3码和BNZS码。

  1. 初始化：

    • 定义了一个三元信码code。
    • 创建了时间轴数据t。

  2. 单极性非归零码（NRZ3）：

    • 当码值为1时，信号为高电平；当码值为0时，信号为低电平。

  3. AMI码：

    • 二进制码0用0电平表示。
    • 二进制码1交替地用+1和-1的半占空归零码表示。

  4. HDB3码：

    • HDB3是一种维护信号的直流偏置为零的编码方案，它在传输时保持相对恒定的电平。
    • 当连续出现4个零时，HDB3码会替换前面出现的3个零。

  5. BNZS码：

    • BNZS码是一种与HDB3类似的码，但其处理连续的零的方式略有不同。

  6. 绘图：

    • 使用MATLAB的plot函数绘制了每种码型的波形。

三、二元码的解码

3.1 二元码解码matalab代码


% 二进制解码
clc
clear
% 定义二进制信码
code = [1 1 1 0 1 0 0 0 1 1 0 0 1 0];
% 创建时间轴数据
t = 0:0.5:(length(code)*100-1)*0.5;
 
%% 单极性非归零码：1的时候为1，0的时候为0
NRZ = [];%存储单极性非归零码
for i = 1:length(code)
    if code(i) == 1
        %使用方括号拼接两个向量[]
        NRZ = [NRZ ones(1, 100)];%为了和t长度一样，一次的跨度也要为100
    else
        NRZ = [NRZ zeros(1, 100)];
    end   
end
%解码：1解码为1，0解码为0
decode_NRZ = [];
for i = 1:100:length(NRZ)
    if NRZ(i) == 0
        decode_NRZ = [decode_NRZ 0];
    else
        decode_NRZ = [decode_NRZ 1];
    end
end
 
%% 双极性非归零码:1的时候为1，0的时候为-1
polar_NRZ = [];
for i = 1:length(code)
    if code(i) == 1
        polar_NRZ = [polar_NRZ ones(1, 100)];
    else
        polar_NRZ = [polar_NRZ -1*ones(1, 100)];
    end
end
%解码:1解码为1，-1解码为0
decode_polar_NRZ = [];
for i = 1:100:length(polar_NRZ)
    if polar_NRZ(i) == -1
        decode_polar_NRZ = [decode_polar_NRZ 0];
    else
        decode_polar_NRZ = [decode_polar_NRZ 1];
    end
end
 
%% 单极性归零码：1的时候分一半为0（归零），0的时候为0
RZ = [];
for i = 1:length(code)
    if code(i) == 1
        RZ = [RZ ones(1, 50) zeros(1, 50)];%从高电平到低电平过渡
    else
        RZ = [RZ zeros(1, 100)];
    end
end
%解码,不管是否归零，以100为跨度，只看第一个
decode_RZ = [];
for i=1:100:length(RZ)  %设置步长，只看第一个
    if RZ(i) == 1
        decode_RZ = [decode_RZ 1]
    else
        decode_RZ = [decode_RZ 0]
    end
end
 
%% 双极性归零：在双极性非归零码的基础上，给1和-1均加入过渡值,即一半归零
polar_RZ = [];
for i = 1:length(code)
    if code(i) == 1
        polar_RZ = [polar_RZ ones(1, 50) zeros(1, 50)];
    else
        polar_RZ = [polar_RZ -1*ones(1, 50) zeros(1, 50)];
    end
end
%解码
decode_polar_RZ = []
for i= 1:100:length(polar_RZ)
    if polar_RZ(i) == 1
        decode_polar_RZ = [decode_polar_RZ 1]
    else
        decode_polar_RZ = [decode_polar_RZ 0]
    end
end
 
%% 传号差分码：波形跳变取值为1
NRZ_M = [];
previous_bit = 0; % 初始值为0，用于表示前一个比特的数值
for i = 1:length(code)
    if code(i) == previous_bit
        bit_value = 0; % 当前位与前一位异或运算,相同为0，不同为1
    else
        bit_value = 1;
    end
    NRZ_M = [NRZ_M ones(1, 100)*bit_value];
    
    previous_bit = bit_value; % 更新前一位比特的数值
end
%解码
decoded_NRZ_M = [];
previous_bit = 0; % 初始值为0，用于表示前一个比特的数值
for i=1:100:length(NRZ_M)
    if  NRZ_M(i) == previous_bit
        decoded_NRZ_M = [decoded_NRZ_M 0]
    else
        decoded_NRZ_M = [decoded_NRZ_M 1]
    end
    previous_bit = NRZ_M( i )
end
 
 
%% 空号差分码：波形跳变取值为0
NRZ_S = [];
previous_bit = 0; % 初始值为0，用于表示前一个比特的数值
for i = 1:length(code)
    if code(i) == previous_bit
        bit_value = 1;
    else   
        bit_value = 0; % 当前位与前一位异或运算
    end
    NRZ_S = [NRZ_S ones(1, 100)*bit_value];
    previous_bit = bit_value; % 更新前一位比特的数值
end
%解码
decoded_NRZ_S = [];
previous_bit = 0; % 初始值为0，用于表示前一个比特的数值
for i=1:100:length(NRZ_S)
    if  NRZ_S(i) == previous_bit
        decoded_NRZ_S = [decoded_NRZ_S 1]
    else
        decoded_NRZ_S = [decoded_NRZ_S 0]
    end
    previous_bit = NRZ_S( i )
end
 
%% 数字双向码
%0：用10表示；1：用01表示
manchester_code  = [];
for i = 1:length(code)
    if code(i) == 1
        manchester_code = [manchester_code -1*ones(1,50) ones(1,50)];
    else
        manchester_code = [manchester_code ones(1,50) -1*ones(1,50)];
    end
end
%解码 1--0，0--1
decode_manchester_code  = [];
for i = 1:100:length(manchester_code)
    if manchester_code(i) == 1
        decode_manchester_code = [decode_manchester_code 0]
    else
        decode_manchester_code =[decode_manchester_code 1]
    end
end
 
 
%% 密勒码
miler_code = [];
flag = 1;
count_0 = 0;%存储连0的个数
for i=1:length(code)
    if code(i) == 1
        count_0 = 0;%清空连0个数
        miler_code = [miler_code -1*flag*ones(1,50) flag*ones(1,50)]
        if flag == 1
            flag = -1;
        else
            flag = 1;
        end
    else  %等于0
        count_0 = count_0+1;
        if(count_0 < 2)%不连零的情况
            miler_code = [miler_code -flag*ones(1,100)]
        else%连零
            miler_code = [miler_code flag*ones(1,100)]
            flag = -flag;
        end
    end
end
            
%密勒码解码49==99-----0（连0）
decode_miler_code = []
for i=1:100:length(miler_code)
    if miler_code(i+49) == miler_code(i+99)
        decode_miler_code = [decode_miler_code 0]
    else
        decode_miler_code = [decode_miler_code 1]
    end
end
 
%% 传号反转码
%1交替的用00和11表示，0固定用01表示
CMI=[];
flag = 1; % 初始极性为正
for i = 1:length(code)
    if code(i) == 1
        CMI = [CMI ones(1, 100)*flag];
       if flag==1
           flag=-1;
       else
           flag=1;
       end
    else
        CMI = [CMI ones(1, 50)*-1 ones(1, 50)]; % 0：交替极性
    end
end
 
%传号反转码解码
decode_CMI = []
for i = 1:100:length(CMI)
    if CMI(i+49) * CMI(i+99) == -1
        decode_CMI = [decode_CMI 0]
    else
        decode_CMI = [decode_CMI 1]
    end
end
 
disp("源码：")
disp(code)
disp("单极性非归零码解码")
disp(decode_NRZ)
disp("双极性非归零码解码")
disp(decode_polar_NRZ)
disp("单极性归零码解码")
disp(decode_RZ)
disp("双极性归零码解码")
disp(decode_polar_RZ)
disp("传号差分码解码")
disp(decoded_NRZ_M);
disp("空号差分码解码")
disp(decoded_NRZ_S)
disp("数字双相码解码")
disp(decode_manchester_code)
disp("密勒码解码")
disp(decode_miler_code)
disp("传号反转码解码")
disp(decode_CMI)

四、三元码的解码

4.1 三元码的解码代码实现

%三元码的解码
clc
clear
% 定义二进制信码
code = [1 1 0 1 0 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0 0 1];
%创建时间轴数据
t = 0:0.5:(length(code)*100-1)*0.5;
 
%% AMI码编码
AMI=[];
flag = 1; % 初始极性为正
for i = 1:length(code)
    if code(i) == 1
        AMI = [AMI ones(1, 50)*flag zeros(1,50)]
        if flag==1
           flag=-1;
        else
           flag=1;
        end
    else
        AMI = [AMI zeros(1,100)]
    end
end
%解码
decode_AMI = []
for i=1:100:length(AMI)
    if AMI(i) == 0
        decode_AMI = [decode_AMI 0]
    else
        decode_AMI = [decode_AMI 1]
    end
end
 

%% HDB3码编码
HDB3=[];
flag = 1; % 初始极性为正
count_0 = 0;%记录连0个数
count2V_1 = 0;%计算两个V之间1的个数，用于确定是B00V型还是000V型
count_1 = 0; %计算总的1的个数，用于确定B的极性
 
for i = 1:length(code)
    if code(i) == 1   %当不连0时，与AMI一样
         HDB3 = [HDB3 ones(1, 50)*flag zeros(1,50)]
         count_0 = 0; %连0个数清0,开启新一轮计数
         count_1 = count_1+1;%总的1的个数加1
         count2V_1 = count2V_1 + 1;%两个V之间1的个数加1
        if flag==1
           flag=-1;
        else
           flag=1;
        end
    else  %出现了0
        count_0 = count_0+1;
        if(count_0 < 4)  %仍然和AMI码一样
            HDB3 = [HDB3 zeros(1,100)]
        else   %出现了四个连零
            if mod(count2V_1 ,2 ) == 1 %奇数个，为000V类型,V的极性跟着上一个B走(10)
                HDB3 = [HDB3 ones(1,50) zeros(1,50)]
                count_0 = 0;  %连零个数清零
                count2V_1 = 0 %两个V之间1的个数清零，总的1的个数不清零
            else %偶数个，为B00V型,V的极性跟着上一个B走(-10)
                     if mod(count_1 ,2) == 0 %偶数个1，B的极性轮到1，V跟着B走
                         HDB3( (i-3)*100-99:(i-3)*100-50 ) = 1*ones(1,50)  %B
                         HDB3 = [HDB3 ones(1,50) zeros(1,50)]              %V
                         flag = -flag;
                         count_0 = 0;  %连零个数清零，开启新一轮连0计数
                         count2V_1 = 0 %两个V之间1的个数清零，总的1的个数不清零
                         count_1 = count_1 + 1;
                     else %奇数个1，B的极性轮到-1，V跟着B走
                         HDB3( (i-3)*100-99:(i-3)*100-50 ) = -1*ones(1,50)  %B
                         HDB3 = [HDB3 -1*ones(1,50) zeros(1,50) ]           %V
                         flag = -flag;
                         count_0 = 0;
                         count2V_1 = 0;
                         count_1 = count_1 +1;
                 end
            end
        end
    end
end
 
%解码
%若HDB3码为0，则源码为0
decode_HDB3 = []
count = 0; %计算连0的个数
for i = 1:100:length(HDB3)   
    if HDB3(i) == 0
        count = count+1;
        if(count == 3)%3连零前后非零脉冲同极性
            if HDB3(i+100) * HDB3(i-300) == 1
                HDB3( i+100:i+199 ) = zeros(1,100);
            end
        end
        
        if(count == 2)%2连零前后非零脉冲同极性
            if HDB3(i+100) * HDB3(i-200) == 1
                HDB3( i+100:i+199 ) = zeros(1,100);
                HDB3( i-200:i-101 ) = zeros(1,100);
            end
        end       
    else
        count = 0;
    end
end
 
%把-1极性修改为1 
for i=1:100:length(HDB3)
    if HDB3(i) == 0
        decode_HDB3 = [decode_HDB3 0]
    else
        decode_HDB3 = [decode_HDB3 1]
    end
end
                
 
disp("源码：")
disp(code)
disp("AMI解码：")
disp(decode_AMI)
disp("HDB3码解码：")
disp(decode_HDB3)
写的什么石山代码
关注
  • 7
    点赞
  • 踩
  • 10
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
am调制解调和相干检波matlab,AM信号的调制与解调仿真
weixin_36177953的博客
03-19 5478
AM信号的调制与解调仿真摘要：摘要：为了信号适合在信道中传输通常要将信号进行调制，本文介绍了一种AM信号的调制与解调的过程。并通过MATLAB仿真形象的看到各个过程中的图像，将抽象的公式及原理变化为可以直观观察的波形，是通信原理课程学习中非常有效的辅助手段，便于学生理解AM调制与解调的原理。关键词：关键词：AM；调制；解调；MATLAB作者：张婷作者单位：(山西大学商务学院，信息学院，山西 太原 ...
通信原理matlab版,通信原理matlab
weixin_36332394的博客
03-21 1174
帮助读者掌握通信的基本概念、基本理论和基本方法，理解通信系统的组成和设计过程。 在内容安排上循序渐进、由浅入深，以便教学内容组织和学习。 通过大量MATLAB、SystemView仿真实例，加深读者对通信原理的理解。 书中所有实例的仿真模型和程序源代码均可以提供下载。¥59.00定价：¥59.00教育部高等学校电工电子基础课程教学指导委员会推荐教材 电子信息学科基础课程系列教材——面向现代工程师培...
MATLAB中plot()函数的用法
热门推荐
qq_52544642的博客
05-02 2万+
文章目录前言一、plot()函数二、代码运行结果总结 前言 此为本人学习中的一些笔记总结，仅供参考。 提示：以下是本篇文章正文内容，下面案例可供参考 一、plot()函数 提示：关于plot()函数的部分叙述 1.plot(X,Y) 创建 Y 中数据对 X 中对应值的二维线图。plot 函数绘制 Y 对 X 的图。 2.plot(X,Y,LineSpec) 设置线型、标记符号和颜色。 3.plot(x,y1,x,y2,'--',x,y3,':') 表示绘制三条正弦曲线，第一条曲线使用默认的线型， 为第二
matlab制图
u013925378的博客
09-20 485
matlab制图，包括显示两组二维位置，标题，图例，x轴y轴标注等。 close all; clc;clear; filename1 = '1.txt'; filename2 = '2.txt'; data = load(filename1); datat = load(filename2); data1 = data(:,1); data2 = data(:,2); datat1 = da...
ASK、PSK、FSK的调制与解调
everything is possible，anything is OK.
11-13 5903
本文主要涉及数字信号的调制与解调，内容包括：2ASK、2PSK、2FSK的调制与解调以及频谱分析
matlab对信号进行DSB调制与解调（仿真）
moqisaonianqiong的博客
05-22 2万+
DSB （抑制载波双边带） 原理笔者不过多赘述，通信原理教材都有详细介绍，直接上干货 %AM调制信号的MATLAB实现 dt=0.001; %时间采样频谱 Fs=100; fm=1; %信源的最高频率 fc=10; %载波中心频率 T=4; %信号时长 N=T/dt; %采样点个数 t=[0:N-...
数字信号的调制与传输之二进制数字调制原理与matalab仿真
everything is possible，anything is OK.
10-12 619
数字信号的调制与传输中，二进制数字调制是一种常见的方式，它用来将数字信号转换为模拟信号或数字信号，以便在传输中传递信息。这种调制方式将二进制数据映射到不同的信号状态或符号上。二进制数字调制有多种类型，包括ASK（幅度移键控）、FSK（频率移键控）、PSK（相位移键控）等，下面本文将介绍它们的原理以及matalab仿真。
MASK、MPSK、MFSK信号的调制与解调+星座图
everything is possible，anything is OK.
11-13 3321
本文主要涉及多进制幅度键控（MASK）、多进制相移键控（MPSK）、多进频移键控（MFSK）的调制与解调，同时涉及到星座图的分析。
MATALAB信号处理——信号的变换（2）
qq_52379608的博客
02-04 1883
2-4 离散系统中的Z域描述 就像我们前面所说的线性时不变离散系统可以用线性常系数差分方程描述。 对两边同时使用Z变换，
2ASK调制与解调的MATALAB仿真实验.docx
07-07
2ASK调制与解调的MATALAB仿真实验.docx2ASK调制与解调的MATALAB仿真实验.docx2ASK调制与解调的MATALAB仿真实验.docx2ASK调制与解调的MATALAB仿真实验.docx2ASK调制与解调的MATALAB仿真实验.docx2ASK调制与解调的...
2ASK调制与解调的MATALAB仿真实验.pdf
07-06
2ASK调制与解调的MATALAB仿真实验.pdf2ASK调制与解调的MATALAB仿真实验.pdf2ASK调制与解调的MATALAB仿真实验.pdf2ASK调制与解调的MATALAB仿真实验.pdf2ASK调制与解调的MATALAB仿真实验.pdf2ASK调制与解调的MATALAB...
2ASK调制与解调的MATALAB仿真实验 (2).docx
07-06
2ASK调制与解调的MATALAB仿真实验 (2).docx2ASK调制与解调的MATALAB仿真实验 (2).docx2ASK调制与解调的MATALAB仿真实验 (2).docx2ASK调制与解调的MATALAB仿真实验 (2).docx2ASK调制与解调的MATALAB仿真实验 (2).docx...
2ASK调制与解调的MATALAB仿真实验 (2).pdf
07-06
2ASK调制与解调的MATALAB仿真实验 (2).pdf2ASK调制与解调的MATALAB仿真实验 (2).pdf2ASK调制与解调的MATALAB仿真实验 (2).pdf2ASK调制与解调的MATALAB仿真实验 (2).pdf2ASK调制与解调的MATALAB仿真实验 (2).pdf2ASK...
MATLAB编程实现2FSK信号的调制与解调（非相干解调）
09-08
MATLAB编程实现2FSK信号的调制与解调（非相干解调），本次资源主要针对非相干解调，下载之后解压，会有MATLAB源码，点开sydgy工程。第一次运行完后会导致数组内存装满，可以在命令行输入clear all回车。
Linux基础 - 环境变量.mp4
06-18
Linux基础 - 环境变量.mp4大厂内部培训资料，都是干货！
模拟退火巡检顺序判定 使用MATLAB实现的机器人导航模拟退火巡检顺序判定，课运行.zip
06-18
模拟退火巡检顺序判定 使用MATLAB实现的机器人导航模拟退火巡检顺序判定，课运行
mcuisp串口下载软件
最新发布
06-18
mcuisp串口下载软件
Java文本编辑器应用程序（含源码）+GUI界面+使用说明+可做毕设
06-18
项目简介 JavaFX文本编辑器应用程序是一个简单的类似记事本的应用程序。在这个应用程序中，您可以轻松创建文本文档。该系统在编辑文本文件或文档时为您提供了更大的灵活性。该系统内置了各种功能，提供了编辑和样式化文本的控制。 系统简介 顾名思义，该项目主要用于编辑和创建文本文档。该项目增加了许多新功能，例如选择字体和字体粗细，以及字体样式（如加粗、斜体或下划线）。不仅如此，您还可以使用颜色格式突出显示文本，或者直接选择文本并为其着色。您还可以进行单词搜索，这有助于控制段落中的拼写错误。您还可以使用菜单栏中的注释符号引用您喜欢的句子。 该项目为国外大神项目，可以作为毕业设计的项目，也可以作为大作业项目，不用担心代码重复，设计重复等，如果需要对项目进行修改，需要具备一定基础知识。 注意：如果装有360等杀毒软件，可能会出现误报的情况，源码本身并无病毒，使用源码时可以关闭360，或者添加信任。
HDB3编码c语言代码
10-29
HDB3编码是一种高密度双极性三级编码，可以将数字信号转换为双极性的信号。以下是HDB3编码的C语言代码实现： ``` #include <stdio.h> #include <string.h> void HDB3(char IN_ZIFU[]) { int i,j=0,k=0; int count=0; char OUT_ZIFU[100]={0}; for(i=0;i<strlen(IN_ZIFU);i++) { if(IN_ZIFU[i]=='1') { count++; OUT_ZIFU[j++]='1'; if(count%2==1) { OUT_ZIFU[j++]='+'; } else { OUT_ZIFU[j++]='-'; } } else { count=0; OUT_ZIFU[j++]='0'; } if(count==4) { OUT_ZIFU[j++]='-'; OUT_ZIFU[j++]='0'; OUT_ZIFU[j++]='-'; OUT_ZIFU[j++]='0'; count=0; } } printf("编码结果为：%s\n",OUT_ZIFU); } ``` 该代码实现了HDB3编码的功能，输入为一个字符串，输出为编码后的字符串。具体实现过程为：遍历输入字符串，如果当前字符为1，则计数器加1，并将1编码为+或-，如果计数器为奇数，则编码为+，否则编码为-；如果当前字符为0，则计数器清零，并将0编码为0；如果计数器连续出现4个1，则将前面的编码替换为000V，其中V为前面编码的极性。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

写的什么石山代码

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值