通信原理与MATLAB(十三):AMI的编解码

已于 2023-06-17 12:26:02 修改
阅读量9.6k 收藏 84
点赞数 17
分类专栏： MATLAB 通信原理 文章标签： matlab 学习
于 2023-01-14 11:30:34 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_47598782/article/details/128683271
版权

目录

1.AMI的的编解码原理

1.1 AMI编码原理

如下图所示，AMI的编码原理:将原始码元的1交替转换成1和-1,0保持不变。

1.2 AMI解码原理

AMI解码原理:将AMI码元中的1和-1转换成1，0转换成0。
如果是经过信道后的AMI码元，那么解码之前先要经过低通滤波器，然后抽样判决后，恢复出AMI码，然后再进行解码。
在这里插入图片描述

2.AMI编解码的代码

clear all;                  % 清除所有变量
close all;                  % 关闭所有窗口
clc;                        % 清屏
%% 基本参数
M=10;                       % 产生码元数    
L=100;                      % 每码元复制L次,每个码元采样次数
Ts=0.001;                   % 每个码元的宽度,即码元的持续时间
Rb=1/Ts;                    % 码元速率1K
dt=Ts/L;                    % 采样间隔
TotalT=M*Ts;                % 总时间
t=0:dt:TotalT-dt;           % 时间
Fs=1/dt;                    % 采样间隔的倒数即采样频率

%% 产生单极性波形
wave=randi([0,1],1,M);      % 产生二进制随机码,M为码元个数
%% AMI编码，将原码中的1交替转换成-1和1,0不变
ami=wave;
a_sign=-1;
for i=1:M
    if ami(i)==0
        continue;
    elseif ami(i)==1&&a_sign==-1
        ami(i)=-ami(i);
        a_sign=1;
    elseif ami(i)==1&&a_sign==1
        ami(i)=ami(i);
        a_sign=-1;
    end
end

fz=ones(1,L);               % 定义复制的次数L,L为每码元的采样点数
x1=wave(fz,:);              % 将原来wave的第一行复制L次，称为L*M的矩阵
jidai=reshape(x1,1,L*M);    % 产生单极性不归零矩形脉冲波形，将刚得到的L*M矩阵，按列重新排列形成1*(L*M)的矩阵

x2=ami(fz,:);               % 将原来ami的第一行复制L次，称为L*M的矩阵
ami_tz=reshape(x2,1,L*M);   % 将刚得到的L*M矩阵，按列重新排列形成1*(L*M)的矩阵

%% 绘制波形
figure(1);                  % 绘制第1幅图
subplot(311);               % 窗口分割成3*1的，当前是第1个子图 
plot(t,jidai,'LineWidth',2);% 绘制基带码元波形，线宽为2
title('基带信号波形');      % 标题
xlabel('时间/s');           % x轴标签
ylabel('幅度');             % y轴标签
axis([0,TotalT,-0.1,1.1])   % 坐标范围限制

subplot(312)                % 窗口分割成3*1的，当前是第2个子图 
plot(t,ami_tz,'LineWidth',2);% 绘制ami的波形 
title('AMI信号波形')        % 标题
axis([0,TotalT,-1.1,1.1]);  % 坐标范围限制
xlabel('时间/s');           % x轴标签
ylabel('幅度');             % y轴标签
%% 信号经过高斯白噪声信道
tz=awgn(ami_tz,20);         % 信号ami_tz中加入白噪声，信噪比为SNR=20dB
subplot(313);               % 窗口分割成3*1的，当前是第1个子图 
plot(t,tz,'LineWidth',2);   % 绘制2ASK信号加入白噪声的波形
axis([0,TotalT,-1.5,1.5]);  % 坐标范围设置
title('通过高斯白噪声信道后的信号');% 标题
xlabel('时间/s');           % x轴标签
ylabel('幅度');             % y轴标签
%% 解码部分
%% 加噪信号经过滤波器
% 低通滤波器设计
fp=2*Rb;                    % 低通滤波器截止频率，乘以2是因为下面要将模拟频率转换成数字频率wp=Rb/(Fs/2)
b=fir1(30, fp/Fs, boxcar(31));% 生成fir滤波器系统函数中分子多项式的系数
% fir1函数三个参数分别是阶数，数字截止频率，滤波器类型
% 这里是生成了30阶(31个抽头系数)的矩形窗滤波器
[h,w]=freqz(b, 1,512);      % 生成fir滤波器的频率响应
% freqz函数的三个参数分别是滤波器系统函数的分子多项式的系数，分母多项式的系数(fir滤波器分母系数为1)和采样点数(默认)512
lvbo=fftfilt(b,tz);         % 对信号进行滤波，tz是等待滤波的信号，b是fir滤波器的系统函数的分子多项式系数
figure(2);                  % 绘制第2幅图  
subplot(411);               % 窗口分割成4*1的，当前是第1个子图 
plot(w/pi*Fs/2,20*log(abs(h)),'LineWidth',2); % 绘制滤波器的幅频响应
title('低通滤波器的频谱');  % 标题
xlabel('频率/Hz');          % x轴标签
ylabel('幅度/dB');          % y轴标签

subplot(412)                % 窗口分割成4*1的，当前是第2个子图 
plot(t,lvbo,'LineWidth',2); % 绘制经过低通滤波器后的信号
axis([0,TotalT,-1.5,1.5]);  % 设置坐标范围
title("经过低通滤波器后的信号");% 标题
xlabel('时间/s');           % x轴标签
ylabel('幅度');             % y轴标签

%% 抽样判决
pdst=1*(lvbo>0);            % 滤波后的向量的每个元素大于0.5为1，小于-0.5为-1，否则为0
for ii=1:length(lvbo)
    if lvbo(ii)>=0.5
        pdst(i)=1;
    elseif lvbo(ii)<=-0.5
        pdst(ii)=-1;
    else
        pdst(ii)=0;
    end
end
% 取码元中间时刻值为判决值
panjue=[];

for j=(L/2):L:(L*M)
    if pdst(j)==1
        panjue=[panjue,1];
    elseif pdst(j)==-1
        panjue=[panjue,-1];
    else
        panjue=[panjue,0];
    end
end

x3=panjue(fz,:);            % 将原来panjue的第一行复制L次，称为L*M的矩阵
huifu=reshape(x3,1,L*M);    % 将刚得到的L*M矩阵，按列重新排列形成1*(L*M)的矩阵

subplot(413)                % 窗口分割成4*1的，当前是第3个子图 
plot(t,huifu,'LineWidth',2) % 画出经过抽样判决后的信号
axis([0,TotalT,-1.1,1.1]);  % 设置坐标范用
title("经过抽样判决后的信号")% 标题
xlabel('时间/s');           % x轴标签
ylabel('幅度');             % y轴标签

%% AMI解码，1和-1对应转换为1，0转换成0
for i=1:M
    if panjue(i)~=0
        dout(i)=1;
    else 
        dout(i)=0;
    end
end

x4=dout(fz,:);              % 将原来dout的第一行复制L次，称为L*M的矩阵
dout_wave=reshape(x4,1,L*M);% 将刚得到的L*M矩阵，按列重新排列形成1*(L*M)的矩阵

subplot(414)                % 窗口分割成4*1的，当前是第4个子图 
plot(t,dout_wave,'LineWidth',2)% 画出经过抽样判决后的信号
axis([0,TotalT,-0.1,1.1]);  % 设置坐标范用
title("AMI解码后波形")      % 标题
xlabel('时间/s');           % x轴标签
ylabel('幅度');             % y轴标签


%% 绘制频谱
%% 基带信号频谱
T=t(end);                   % 时间
df=1/T;                     % 频谱分辨率
N=length(jidai);            % 采样长度
f=(-N/2:N/2-1)*df;          % 频率范围
mf=fftshift(abs(fft(jidai)));%对信源信号采用快速傅里叶变换并移到矩阵中心
figure(3)                   % 绘制第3幅图
subplot(211);               % 窗口分割成2*1的，当前是第1个子图 
plot(f,mf,'LineWidth',2);   % 绘制信源频谱波形
title("基带信号频谱");      % 标题
xlabel('频率/Hz');          % x轴标签
ylabel('幅度');             % y轴标签

%% AMI信号频谱
sf=fftshift(abs(fft(ami_tz))); % 对AMI信号采用快速傅里叶变换并将0-fs频谱移动到-fs/2-fs/2
subplot(212)                % 窗口分割成2*1的，当前是第1个子图 
plot(f,sf,'LineWidth',2)    % 绘制调制信号频谱
title("AMI信号频谱")        % 标题
xlabel('频率/Hz');          % x轴标签
ylabel('幅度');             % y轴标签


%% 加入噪声后的频谱
mmf=fftshift(abs(fft(tz))); % 对信号采用快速傅里叶变换并移到矩阵中心
figure(4)                   % 绘制第4幅图
subplot(211);               % 窗口分割成2*1的，当前是第1个子图 
plot(f,mmf,'LineWidth',2)   % 画出乘以相干载波后的频谱
title("加入噪声后的频谱")   % 标题
xlabel('频率/Hz');          % x轴标签
ylabel('幅度');             % y轴标签

%% 经过低通滤波后的频谱
dmf=fftshift(abs(fft(lvbo)));%对低通滤波信号采用快速傅里叶变换并移到矩阵中心
subplot(212);               % 窗口分割成2*1的，当前是第2个子图 
plot(f,dmf,'LineWidth',2)   % 画出经过低通滤波后的频谱
title("经过低通滤波后的频谱");% 标题
xlabel('频率/Hz');          % x轴标签
ylabel('幅度');             % y轴标签

3.AMI编解码结果图

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

4.AMI的误码率曲线

4.1 原理

计算某信噪比下的错误的码元个数，然后除以总码元个数，即该信噪比下的误码率。然后更改信噪比，计算新的误码率，一个信噪比对应一个误码率，连点成线。

4.2 AMI的误码率曲线代码

clear all;                  % 清除所有变量
close all;                  % 关闭所有窗口
clc;                        % 清屏
%% 基本参数
M=100000;                   % 产生码元数    
L=100;                      % 每码元复制L次,每个码元采样次数
Ts=0.001;                   % 每个码元的宽度,即码元的持续时间
Rb=1/Ts;                    % 码元速率1K
dt=Ts/L;                    % 采样间隔
TotalT=M*Ts;                % 总时间
t=0:dt:TotalT-dt;           % 时间
Fs=1/dt;                    % 采样间隔的倒数即采样频率

%% 产生单极性波形
wave=randi([0,1],1,M);      % 产生二进制随机码,M为码元个数
%% AMI编码，将原码中的1交替转换成-1和1,0不变
ami=wave;
a_sign=-1;
for i=1:M
    if ami(i)==0
        continue;
    elseif ami(i)==1&&a_sign==-1
        ami(i)=-ami(i);
        a_sign=1;
    elseif ami(i)==1&&a_sign==1
        ami(i)=ami(i);
        a_sign=-1;
    end
end

fz=ones(1,L);               % 定义复制的次数L,L为每码元的采样点数
x1=wave(fz,:);              % 将原来wave的第一行复制L次，称为L*M的矩阵
jidai=reshape(x1,1,L*M);    % 产生单极性不归零矩形脉冲波形，将刚得到的L*M矩阵，按列重新排列形成1*(L*M)的矩阵

x2=ami(fz,:);               % 将原来ami的第一行复制L次，称为L*M的矩阵
ami_tz=reshape(x2,1,L*M);   % 将刚得到的L*M矩阵，按列重新排列形成1*(L*M)的矩阵

%% 信号经过高斯白噪声信道
EbN0_dB = -15:1:5;
for i=1:21
tz=awgn(ami_tz,EbN0_dB(i)); % 信号ami_tz中加入白噪声

%% 解码部分
%% 加噪信号经过滤波器
% 低通滤波器设计
fp=2*Rb;                    % 低通滤波器截止频率，乘以2是因为下面要将模拟频率转换成数字频率wp=Rb/(Fs/2)
b=fir1(30, fp/Fs, boxcar(31));% 生成fir滤波器系统函数中分子多项式的系数
% fir1函数三个参数分别是阶数，数字截止频率，滤波器类型
% 这里是生成了30阶(31个抽头系数)的矩形窗滤波器
[h,w]=freqz(b, 1,512);      % 生成fir滤波器的频率响应
% freqz函数的三个参数分别是滤波器系统函数的分子多项式的系数，分母多项式的系数(fir滤波器分母系数为1)和采样点数(默认)512
lvbo=fftfilt(b,tz);         % 对信号进行滤波，tz是等待滤波的信号，b是fir滤波器的系统函数的分子多项式系数

%% 抽样判决
pdst=1*(lvbo>0);            % 滤波后的向量的每个元素大于0.5为1，小于-0.5为-1，否则为0
for ii=1:length(lvbo)
    if lvbo(ii)>=0.5
        pdst(i)=1;
    elseif lvbo(ii)<=-0.5
        pdst(ii)=-1;
    else
        pdst(ii)=0;
    end
end
% 取码元中间时刻值为判决值
panjue=[];

for j=(L/2):L:(L*M)
    if pdst(j)==1
        panjue=[panjue,1];
    elseif pdst(j)==-1
        panjue=[panjue,-1];
    else
        panjue=[panjue,0];
    end
end

%% AMI解码，1和-1对应转换为1，0转换成0
for k=1:M
    if panjue(k)~=0
        dout(k)=1;
    else 
        dout(k)=0;
    end
end

% 计算错误码元数，然后除以总码元数得到误码率
error_num=length(find(dout~=wave));

error_rate(i)=error_num/M;

end

figure(1)
semilogy(EbN0_dB,error_rate,'-ob');% 绘制误码率曲线
xlabel('Eb/n0');                   % x轴标签
ylabel('BRE');                     % y轴标签
title('AMI的误码率曲线')           % 标题

4.3 误码率曲线图

在这里插入图片描述

通信工程实验四 AMI/HDB3编译码实验
yibolsx的博客
02-25 955
正确译码之后21TP01与20P01的波形应一致，但由于HDB3码的编译码规则较复杂，当前的输出HDB3码字与前4个码字有关，因而HDB3码的编译码时延较大。HDB3码保留了AMI码所有的优点（如前所述），还可将连“0”码限制在3个以内，克服了AMI码出现长连“0”过多，对提取定时钟不利的缺点。1．二进制序列中的“0”码在HDB3码中仍编为“0”码，但当出现四个连“0”码时，用取代节000V或B00V代替四个连“0”码。取代节中的V码、B码均代表“1”码，它们可正可负（即V。以恢复成原二进制码序列）。
AMI编码演示
02-23
AMI编码中的1码通常称为传号，0码则叫空号，这是沿用了早期电报通信中的叫法。从形态上看，它已是三状态信号，所以AMI码是伪三进制码。 一、编码规则: 消息代码中的0 传输码中的0 消息代码中的1 传输码中的+1、-1交替 例如: 消息代码:1 0 1 0 1 0 0 0 1 0 1 1 1 AMI码: +1 0 -1 0 +1 0 0 0 -1 0 +1 -1 +1 此例程演示AMI码的生成过程：message中的字符先先被转换成ASC编码，并以二进制形式翻译成消息码，再按AMI规则转换成AMI码
Matlab技术的AMI编码和解码
04-12
基于Matlab技术的AMI编码和解码的实现
AMI编码matlab
05-18
在matlab上的AMI编码实现及频谱分析
用matlab，simulink实现的ami编码和解码
11-30
用matlab编程实现的ami编码解码，用simulink实现的显示波形，没想到其他的办法，这样勉强实现了，欢迎高手和我讨论，邮箱chen0510566@gmail.com
AMI2编码 matlab实现
11-19
AMI2编码的matlab实现 本代码是老师上课要求写的，自己编码实现
通信原理ami码c语言实现,通信原理AMI码型变换实验
weixin_35343148的博客
05-23 5011
《通信原理AMI码型变换实验》由会员分享，可在线阅读，更多相关《通信原理AMI码型变换实验(16页珍藏版)》请在人人文库网上搜索。1、实验二 AMI码型变换实验一、实验目的1、 了解几种常用的数字基带信号的特征和作用。2、 掌握AMI码的编译规则。3、 了解滤波法位同步在的码变换过程中的作用。二、实验器材1、 主控&信号源、2号、8号、13号模块 各一块2、 双踪示波器 一台3、 连接线 ...
MATLAB教程：HDB3编解码仿真与AMI码转换详解
用户将学会如何使用MATLAB强大的矩阵运算和图形显示功能来设计HDB3编码器和解码器，以及如何将AMI码转换为HDB3码，并能够进行相应的解码过程。 教程中的编程部分不仅包括函数的编写，还会详细讲解如何处理边界条件...
通信原理实验必看：AMI码与HDB3码型变换的关键步骤与要点
![通信原理实验必看：AMI码与HDB3码型变换的关键步骤与要点]... # 摘要 本论文旨在回顾数字通信基础，并深入探讨AMI码与HDB3码型的基本概念、变换机制及实践技巧。文章首先介绍了AMI码与HDB3码的定义、原理、编码规则...
Matlab技术的AMI编码和解码的实现
04-09
在本文中，我们将深入探讨如何使用Matlab技术实现AMI（ Alternate Mark Inversion，交替标志反转）编码和解码。AMI编码是一种广泛应用于数字通信中的线性脉冲编码调制（PCM）方法，它通过改变脉冲的位置来表示二进制...
matlab用于HDB3的编解码，用matlab仿真，有完整的程序.zip
04-20
标题中的“matlab用于HDB3的编解码，用matlab仿真，有完整的程序”指的是利用MATLAB语言实现HDB3编码和解码的过程，并且提供了一整套的仿真程序。下面将详细解释HDB3编码的基本原理以及MATLAB如何进行仿真。 HDB3...
matlabami代码-matlab:MATLAB
06-07
matlab ami代码MatLab 备忘单 :memo: 目录 MatLab 的东西 ans变量存储上次操作的结果。 clear命令清除内存。 clc命令清除控制台历史记录。 help 命令显示一个缩写的文档。 （这不是帮助中的内容） 当除以0而不是抛出错误时，计算除法的“限制”（ 1/0 == Inf ， -1/0 == -Inf ） label 和 title 函数也可以在 latex 模式下工作： f(..., 'interpreter', 'latex') 控制台编写 disp(x)函数可用于将预定义值写入控制台： disp('Hello world!')或变量的值： disp(a) 。 如果你不结束一行; -it，它的值将显示在控制台上。 使用这个技巧和sprintf('format string', param1, param2, ...)函数，您可以sprintf('format string', param1, param2, ...)包含变量的文本 ( = template ) sprintf('format string', param1, p
利用MATLAB实现AMI、HDB3码-画图
07-14
数字通信基础实验一，利用MATLAB实现AMI、HDB3码-画图，有详细代码注释！好用！
matlabami代码-line-coding:模拟流行线路编码技术的Matlab函数
06-07
matlab ami代码关于 这个存储库包含模拟一些流行的线路编码技术的 matlab 函数。 我把它们写成大学作业的一部分。 该代码的灵感来自 . 职能 每个函数都接受一个包含“1”和“0”值的数组。 p_nrzl()极性不归零 u_nrzl()单极不归零 nrzi()不归零倒置 manchester()曼彻斯特 diff_manchester()差分曼彻斯特 b_ami()双极交替标记反转 pseudoternary()三元 b8zs()具有 8 个零替换的双极 hdb3()高密度双极 3 零 截图 查看项目的 . 执照 请参阅复制文件。 作者 :copyright:2012 Tasos Latsas
通信原理AMIHDB3码型变换实验
06-07
AMI码的全称是传号交替反转码。这是一种将消息代码0（空号）和1（传号）按如下规则进行编码的码：代码的0仍变换为传输码的0，而把代码中的1交替地变换为传输码的+1、–1、+1、–1…
HDB3码、AMI码的MATLAB实现
01-10
HDB3码、AMI码编码和译码的MATLAB实现
hdb3.rar_AMI HDB3_HDB3_HDB3 matlab_ami matlab_hdb3 ami matlab
09-24
6. AMI解码：将AMI码还原为二进制数据。 标签"ami_hdb3 hdb3 hdb3_matlab ami_matlab hdb3_ami_matlab"进一步确认了文件内容的重点，即MATLAB实现的HDB3码和AMI码的相关工作。 综上所述，这个MATLAB项目提供了对...
ami码matlab,matlab编程实现ami码的相关问题
weixin_27531501的博客
03-17 1653
错误是数组索引必须为正整数或逻辑值，但我不知道咋改，哪位大佬指点一下完整代码：%%%%% 参数初始化 %%%%%%%%%%%%%Ts=1; %码元周期N_sample=8; %每个码元的抽样点数9dt=Ts/N_sample; %抽样时间间隔N=30; %码元数t=0:dt:(N*N_sample-1)*dt;T=t(end);Again=100; %迭代次数%%%数据初始化 %%%%%%%%%%...
【MATLAB源码-第291期】基于matlab的AMI编码解码系统仿真，输出各个节点波形。
最新发布
Matlab程序猿i的博客
11-12 568
AMI编码的特点是在传输过程中，对于0信号使用相同的电平，而对于1信号则交替使用不同的电平。它通过避免长时间的直流分量来提高信号的抗干扰能力，减少了直流漂移的影响，从而使信号的频谱特性更加有利于传输。具体来说，对于二进制流中的每一个“1”，代码生成了一个+1或-1的极性，而“0”则映射为0。在通信系统中，AMI编码后的信号通常是通过脉冲信号的方式进行传输的，这些脉冲信号是经过基带处理的。在AMI编码中，对于输入信号中的每一个“1”，都会对应一个不同的电平，通常是交替的+1和-1；例如，如果输入的比特流为。
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

菜菜菜~~

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

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

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

打赏作者

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

抵扣说明：

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

余额充值