CN113660014B - 一种基于多谐振频率的信道容量优化方法及装置 - Google Patents

一种基于多谐振频率的信道容量优化方法及装置 Download PDF

Info

Publication number
CN113660014B
CN113660014B CN202110886466.6A CN202110886466A CN113660014B CN 113660014 B CN113660014 B CN 113660014B CN 202110886466 A CN202110886466 A CN 202110886466A CN 113660014 B CN113660014 B CN 113660014B
Authority
CN
China
Prior art keywords
matrix
magnetic induction
communication system
induction communication
channel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110886466.6A
Other languages
English (en)
Other versions
CN113660014A (zh
Inventor
陈洋
赵毓斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN202110886466.6A priority Critical patent/CN113660014B/zh
Publication of CN113660014A publication Critical patent/CN113660014A/zh
Application granted granted Critical
Publication of CN113660014B publication Critical patent/CN113660014B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/382—Monitoring; Testing of propagation channels for resource allocation, admission control or handover
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0426—Power distribution
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Power Engineering (AREA)
  • Near-Field Transmission Systems (AREA)

Abstract

本发明公开了一种基于多谐振频率的信道容量优化方法及装置，包括：计算磁感应通信系统中的互感矩阵；获取磁感应通信系统的第一等效电路模型和第二等效电路模型；根据基尔霍夫电压定律以及第一等效电路模型计算磁感应通信系统的通信信道；根据基尔霍夫电压定律以及第二等效电路模型计算磁感应通信系统的谐振频率；根据通信信道计算磁感应通信系统的接收信号向量以及谐振频率，并计算得到磁感应通信系统的信道容量；根据谐振频率对信道容量进行优化，得到最优信道容量。本实施例基于磁感应通信系统的等效电路模型计算通信信道和谐振频率，能够准确确定由于多线圈分频而产生的多共振频率，使得系统能够选择多个频率作为子载波信道，有效提高信道容量。

Description

一种基于多谐振频率的信道容量优化方法及装置
技术领域
本发明涉及通信技术领域，尤其是涉及一种基于多谐振频率的信道容量优化方法及装置。
背景技术
磁感应通信(MIC)是在射频信号无法穿透地下或水下环境等介质的具有挑战性或恶劣环境下，一种具有发展前景的新型无线通信模式。MIC系统中部署了几个耦合线圈，可以在相同的谐振频率下同时交换无线电源和信息。在这种情况下，MIC不仅提供了可预测、稳定的信号传播通道，而且传输延迟可以忽略不计。MIC由于其性价比高，被广泛应用于地下传感器网络、自主水下航行器、恶劣环境下的机器人以及一些军事武器。磁感应通信的主要挑战是有限的信道容量，主要受到接收功率的影响。具体来说，接收功率和信道容量是由线圈半径和传输距离决定的。在文献中首先对单载频单对收发线圈(SISO)的信道容量进行了分析。在此基础上，研究了单输入多输出(SIMO)、多输入单输入(MISO)和多输入多输出(MIMO)模型的信道容量。然而，MIMO信道容量的提高仍然是有限的。
现有的磁感应通信系统的信道容量优化方法主要有两类:一类是高质量的硬件和电路设计，另一类是开发优化方法。对于第一类，MIMO磁感应系统为原磁场或电流分布提供了更多的自由度，等效电路模型(ECM)是电路设计和性能分析的主要工具。在文献中，通过仿真和实际实验对双线圈和三线圈系统的性能进行了比较，证明了Tx和Rx线圈强耦合时存在分频现象。然后，Nguyen等人指出通过分裂频率的多通道通信可以用于数据传输。Dionigi等人设计了一种用于无线功率传输和近场通信的多波段谐振器。第二类是设计适当的信道容量优化方法，控制系统的电流或电压。这种方法不需要额外的电路。在文献中，为不同的MI继电器和收发器线圈分配不同的谐振频率，以计算系统的信道容量。
现有的多频率磁MIMO信道容量优化方法存在以下问题：
第一类方法需要额外的电路，例如控制电路和阻抗匹配网络，导致需要引入额外的功耗，而增加系统的复杂性。
第二类方法需要寄生电容和其他电路信息以获取共振频率，增加了计算复杂度。
发明内容
本发明提供了一种基于多谐振频率的信道容量优化方法及装置，以解决现有的空信道容量优化方法需要引入额外的功耗以及增加计算复杂度的问题。
本发明的第一实施例提供了一种基于多谐振频率的信道容量优化方法，包括：
计算磁感应通信系统中的发射端的阻抗矩阵和接收端的阻抗矩阵，根据所述发射端的阻抗矩阵和接收端的阻抗矩阵计算得到互感矩阵；
获取所述磁感应通信系统对应的等效电路模型，所述等效电路模型包括第一等效电路模型和第二等效电路模型；
根据基尔霍夫电压定律、所述互感矩阵以及所述第一等效电路模型计算所述磁感应通信系统的通信信道；
根据基尔霍夫电压定律以及所述第二等效电路模型构建所述磁感应通信系统的状态空间模型，并根据所述状态空间模型计算所述磁感应通信系统的谐振频率；
根据所述通信信道计算所述磁感应通信系统的接收信号向量，根据所述接收信号向量的协方差以及所述谐振频率，计算得到所述磁感应通信系统的信道容量；
在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率对所述信道容量进行优化，得到最优信道容量。
进一步的，所述根据基尔霍夫电压定律、所述互感矩阵以及所述第一等效电路模型计算所述磁感应通信系统的通信信道，具体为：
根据基尔霍夫电压定律计算所述第一等效电路模型的第一电路方程，根据所述第一电路方程计算得到电流矩阵；根据所述电流矩阵和所述互感矩阵计算得到接收电压，以所述接收电压定义接收信号和发送信号，根据所述接收电压计算所述接收信号与所述发射信号之间的系数作为通信信道。
进一步的，所述根据基尔霍夫电压定律以及所述第二等效电路模型构建所述磁感应通信系统的状态空间模型，并根据所述状态空间模型计算所述磁感应通信系统的谐振频率，具体为：
根据基尔霍夫电压定律计算所述第二等效电路模型的第二电路方程，根据所述第二电路方程计算所述磁感应通信系统的电容矩阵、电阻矩阵和电感矩阵，根据所述电容矩阵、电阻矩阵和电感矩阵构建所述状态空间模型，利用所述空间模型推导计算得到所述磁感应通信系统的谐振频率。
进一步的，所述在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率对所述信道容量进行优化，得到最优信道容量，具体为：
在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率在子信道上分配的Tx能量对所述信道容量的表达式进行优化，得到最优信道容量。
进一步的，所述根据所述谐振频率在子信道上分配的Tx能量对所述信道容量进行优化，得到最优信道容量，具体为：
根据所述通信信号将所有特征值按照大小降序排序，得到特征值排序集合；
通过对所述特征值排序结合进行迭代计算和子信道分配功率的非负性计算得到优化常数；
根据所述优化常数对子信道进行功率分配，当分配的子信道功率为负时将所述子信道丢弃，并重新运行迭代计算器加一；当分配给所述磁电感通信系统各个子型号的功率为负时，得到最优信道容量。
本发明的第二实施例提供了一种基于多谐振频率的信道容量优化装置，包括：
第一计算模块，用于计算磁感应通信系统中的发射端的阻抗矩阵和接收端的阻抗矩阵，根据所述发射端的阻抗矩阵和接收端的阻抗矩阵计算得到互感矩阵；
电路获取模块，用于获取所述磁感应通信系统对应的等效电路模型，所述等效电路模型包括第一等效电路模型和第二等效电路模型；
第二计算模块，用于根据基尔霍夫电压定律、所述互感矩阵以及所述第一等效电路模型计算所述磁感应通信系统的通信信道；
第三计算模块，用于根据基尔霍夫电压定律以及所述第二等效电路模型构建所述磁感应通信系统的状态空间模型，并根据所述状态空间模型计算所述磁感应通信系统的谐振频率；
第四计算模块，用于根据所述通信信道计算所述磁感应通信系统的接收信号向量，根据所述接收信号向量的协方差以及所述谐振频率，计算得到所述磁感应通信系统的信道容量；
信道容量优化模块，用于在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率对所述信道容量进行优化，得到最优信道容量。
进一步的，所述第二计算模块，具体用于：
根据基尔霍夫电压定律计算所述第一等效电路模型的第一电路方程，根据所述第一电路方程计算得到电流矩阵；根据所述电流矩阵和所述互感矩阵计算得到接收电压，以所述接收电压定义接收信号和发送信号，根据所述接收电压计算所述接收信号与所述发射信号之间的系数作为通信信道。
进一步的，所述第三计算模块，具体用于：
根据基尔霍夫电压定律计算所述第二等效电路模型的第二电路方程，根据所述第二电路方程计算所述磁感应通信系统的电容矩阵、电阻矩阵和电感矩阵，根据所述电容矩阵、电阻矩阵和电感矩阵构建所述状态空间模型，利用所述空间模型推导计算得到所述磁感应通信系统的谐振频率。
进一步的，所述信道容量优化模块，具体用于：
在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率在子信道上分配的Tx能量对所述信道容量的表达式进行优化，得到最优信道容量。
进一步的，所述根据所述谐振频率在子信道上分配的Tx能量对所述信道容量进行优化，得到最优信道容量，具体为：
根据所述通信信号将所有特征值按照大小降序排序，得到特征值排序集合；
通过对所述特征值排序结合进行迭代计算和子信道分配功率的非负性计算得到优化常数；
根据所述优化常数对子信道进行功率分配，当分配的子信道功率为负时将所述子信道丢弃，并重新运行迭代计算器加一；当分配给所述磁电感通信系统各个子型号的功率为负时，得到最优信道容量。
本发明实施例采用时频域分析模型方法，基于磁感应通信系统的等效电路模型计算通信信道和谐振频率，从而能够准确确定由于多线圈分频而产生的多共振频率，使得系统能够选择多个频率作为子载波信道，能够在不额外增加电路及功耗的情况下，有效扩展带宽，从而有效地提高了信道容量。
本发明实施例利用谐振频率作为子载波，推到出磁感应系统的信道容量，并且在给定的功率限制下，根据谐振频率对信道容量进行优化，得到最优信道容量，且本发明实施例能够使根据多个子载波来设计波形，使得系统的信道容量得到有效提高。
附图说明
图1是本发明实施例提供的一种基于多谐振频率的信道容量优化方法的流程示意图；
图2是本发明实施例提供的磁感应通信系统的结构式示意图；
图3是本发明实施例提供的第一等效电路模型的示意图；
图4是本发明实施例提供的第二等效电路模型的示意图；
图5是本发明实施例提供的一种基于多谐振频率的信道容量优化装置的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
在本申请的描述中，需要理解的是，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本申请的描述中，除非另有说明，“多个”的含义是两个或两个以上。
在本申请的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本申请中的具体含义。
请参阅图1-4，在本发明的第一实施例中，本发明的第一实施例提供了如1所示的一种基于多谐振频率的信道容量优化方法，包括：
S1、计算磁感应通信系统中的发射端的阻抗矩阵和接收端的阻抗矩阵，根据发射端的阻抗矩阵和接收端的阻抗矩阵计算得到互感矩阵；
S2、获取磁感应通信系统对应的等效电路模型，等效电路模型包括第一等效电路模型和第二等效电路模型；
S3、根据基尔霍夫电压定律、互感矩阵以及第一等效电路模型计算磁感应通信系统的通信信道；
S4、根据基尔霍夫电压定律以及第二等效电路模型构建磁感应通信系统的状态空间模型，并根据状态空间模型计算磁感应通信系统的谐振频率；
S5、根据通信信道计算磁感应通信系统的接收信号向量，根据接收信号向量的协方差以及谐振频率，计算得到磁感应通信系统的信道容量；
S6、在磁感应通信系统的总发射功率约束条件下，根据谐振频率对信道容量进行优化，得到最优信道容量。
可选地，本发明实施例可在如图2所示的磁感应通信系统中执行。请参阅图2，磁感应通信系统包括Tx端、Rx端、数字信号处理器(DSP)、数字模拟转换器(DAC)、放大器(AF)。在Tx端，使用数字信号处理器(DSP)产生发射的基带信号，数字模拟转换器(DAC)将离散的数字基带信号转换为连续的可变模拟基带信号，放大器(AF)用于调制和放大这种模拟信号，每个发射线圈与设置的一个特定谐振频率的电容器相连。检测模块用于测量每个发射线圈上的电压幅值和相位，并反馈给DSP。在接收侧，每个接收线圈与为了获得相同谐振频率而设置的接收侧电容器相连。模数转换器(ADC)对接收到的波形进行采样并发送给DSP。DSP用于对这些数据进行解码，以便进行通信和信号处理。
请参阅图3-4，为上述磁感应通信系统的两种等效电路模型。
请参阅表1，为本发明实施例提供的等效电路模型各参数定义。
表1 等效电路模型各参数定义
Figure GDA0003769637660000071
Figure GDA0003769637660000081
每个发射线圈的等效阻抗表示为
Figure GDA0003769637660000082
发射端的阻抗矩阵表示为
Figure GDA0003769637660000083
每个接收线圈的等效阻抗为
Figure GDA0003769637660000084
得到接收端的阻抗矩阵
Figure GDA0003769637660000091
互感矩阵定义为ZM，
Figure GDA0003769637660000092
表示Rx线圈对Tx线圈的影响。
Figure GDA0003769637660000093
作为本发明实施例的一种具体实施方式，根据基尔霍夫电压定律、互感矩阵以及第一等效电路模型计算磁感应通信系统的通信信道，具体为：
根据基尔霍夫电压定律计算第一等效电路模型的第一电路方程，根据第一电路方程计算得到电流矩阵；根据电流矩阵和互感矩阵计算得到接收电压，以接收电压定义接收信号和发送信号，根据接收电压计算接收信号与发射信号之间的系数作为通信信道。
示例性的，根据基尔霍夫电压定律(KVL)方程，得到图3所示MIMO磁感应通信系统的电路方程为
Figure GDA0003769637660000094
计算MIMO磁感应通信系统系统中所有电路的电流，将阻抗矩阵进行反求，得到电流矩阵：
Figure GDA0003769637660000095
Figure GDA0003769637660000096
Figure GDA0003769637660000097
根据电流矩阵和互感矩阵计算接收电压：
Figure GDA0003769637660000101
本发明实施例将电压类比为信号，则接收和发射的信号分别为Y和S，即Y和S之间的系数为通信信道H，定义为
Figure GDA0003769637660000102
通信信道高度依赖于线圈间的磁通密度，从电压角度看，MIMO磁感应通信系统类似于一般的无线通信系统，无线信号传播沿一定数量的衰落信道到达接收机，不同的传播机制对路径损耗和衰落模型有不同的影响。本发明实施例定义
Figure GDA0003769637660000103
作为H的元素,其中
Figure GDA0003769637660000104
是衰落因子。
作为本发明实施例的一种具体实施方式，根据基尔霍夫电压定律以及第二等效电路模型构建磁感应通信系统的状态空间模型，并根据状态空间模型计算磁感应通信系统的谐振频率，具体为：
根据基尔霍夫电压定律计算第二等效电路模型的第二电路方程，根据第二电路方程计算磁感应通信系统的电容矩阵、电阻矩阵和电感矩阵，根据电容矩阵、电阻矩阵和电感矩阵构建状态空间模型，利用空间模型推导计算得到磁感应通信系统的谐振频率。
示例性的，本发明实施例根据KVL，图4所示等效电路模型的电路方程为
Figure GDA0003769637660000105
其中Sn为
Figure GDA0003769637660000106
将上述的电路方程改写为
Figure GDA0003769637660000107
其中矩阵C为MIMO磁感应系统的电容矩阵，如下式所示：
Figure GDA0003769637660000111
矩阵R为MIMO磁感应通信系统的电阻矩阵，如下式所示
Figure GDA0003769637660000112
矩阵L为MIMO磁感应通信系统的电感矩阵，如下式所示
Figure GDA0003769637660000113
引入方程
Figure GDA0003769637660000114
构建系统的状态空间模型
Figure GDA0003769637660000115
其中
Figure GDA0003769637660000116
是状态矩阵，S′＝[ST,01×N]T为系统的输入矩阵，A和B具体为
Figure GDA0003769637660000117
其中IN×N表示一个N阶单位矩阵。
利用A和B矩阵构建该系统的系统矩阵D为
Figure GDA0003769637660000118
利用构建的状态空间模型和系统矩阵，进行推导得到如下引理：：MIMO MIC系统存在Nf个谐振频率
Figure GDA0003769637660000119
它们是系统矩阵D特征值的虚部的绝对值，如下式所示。
Figure GDA0003769637660000121
作为本发明实施例的一种具体实施方式，在磁感应通信系统的总发射功率约束条件下，根据谐振频率对信道容量进行优化，得到最优信道容量，具体为：
在磁感应通信系统的总发射功率约束条件下，根据谐振频率在子信道上分配的Tx能量对信道容量的表达式进行优化，得到最优信道容量。
示例性的，考虑MIMO磁感应通信系统信的道由NT个发射线圈和NR个接收线圈组成，本发明实施例假设通信信道已经做归一化处理的并且在这个频带内频率是平坦的，经过推导可知，MIMO磁感应通信系统的接收信号向量为：
Y＝HS+N-
其中H是通信信道，N-是高斯白噪声。
经过推导可知，多频MIMO磁感应通信系统的输入输出关系表示为
Figure GDA0003769637660000122
经过推导可知，多频MIMO磁感应通信系统的信道容量表示为
Figure GDA0003769637660000123
其中，Nf为分裂频率数，mf为在通过式(11)得到频率ωf的H(ωf)H(ωf)H正特征值的个数，λk(ωf)是这些特征值之一，
Figure GDA0003769637660000124
为信号向量在ωf频率处的协方差，N0为通带噪声谱功率。
在本发明实施例中，对于磁感应通信系统，多线圈发射机的主要目标是设计S在总发射功率约束下达到最大的信道容量，即
Figure GDA0003769637660000125
Figure GDA0003769637660000126
SHS≤PC
其中PC为最大允许的总传输功率，
Figure GDA0003769637660000131
为S的协方差.
Figure GDA0003769637660000132
为向量S的元素。
可选地，由于发射端可以接入MIMO磁感应通信系统的子信道，在总发射功率约束下，能够跨越子信道分配可变能量，使信道容量最大化。因此公式(12)的优化如下：
Figure GDA0003769637660000133
其中
Figure GDA0003769637660000134
表示表示在频率为ωf的第k个子信道上分配的Tx能量，满足
Figure GDA0003769637660000135
且最优能量分配方案
Figure GDA0003769637660000136
满足
Figure GDA0003769637660000137
其中η是一个常数，根据特征值的集合的迭代计算和子信道分配能量的非负性获得，x+是指
Figure GDA0003769637660000138
作为本发明实施例的一种具体实施方式，根据谐振频率在子信道上分配的Tx能量对信道容量进行优化，得到最优信道容量，具体为：
根据通信信号将所有特征值按照大小降序排序，得到特征值排序集合；
通过对特征值排序结合进行迭代计算和子信道分配功率的非负性计算得到优化常数；
根据优化常数对子信道进行功率分配，当分配的子信道功率为负时将子信道丢弃，并重新运行迭代计算器加一；当分配给磁电感通信系统各个子型号的功率为负时，得到最优信道容量。
在本发明实施例中，对自适应多频信道容量优化方案进行迭代设计，，在已知H的前提下，将所有的特征值λk(ωf)按大小降序排列，从而形成一个
Figure GDA0003769637660000139
的集合，并利用γl来表示相关的λl的功率分配。设置计数器p等于1，在每次迭代时，根据受约束的功率计算常数1/η*
Figure GDA0003769637660000141
每个模式在每次迭代中分配的功率按下式计算
Figure GDA0003769637660000142
若分配给增益最小信道的能量为负，如γi＜0，通过设置γi＝0丢弃这个信道，并重新运行迭代计数器p加1的方案。当分配给各磁感应通信系统子信道的功率为非负时，为最优的资源分配方案，即得到最优信道容量。
本发明实施例还提供了一种信道估计方法，具体为：
当Tx和Rx线圈固定时，通道几乎是固定的。在这种情况下，Rx线圈的电压通过互感与Tx线圈相互作用，信道估计主要是根据Rx侧的电压和负载计算互感，得到磁信道信息，接收端ADC模块对接收到的信号矢量进行采样。
在磁感应通信系统中，信道由接收机利用发射机发出的训练信号进行估计。接收机知道训练信号序列
Figure GDA0003769637660000143
定义
Figure GDA0003769637660000144
为估计的信道，则接收信号为
Figure GDA0003769637660000145
多线圈发射机的主要目标是在功率受限时控制
Figure GDA0003769637660000146
使得
Figure GDA0003769637660000147
误差达到最小值，即
Figure GDA0003769637660000148
Figure GDA0003769637660000149
Figure GDA00037696376600001410
其中
Figure GDA00037696376600001411
为为第Ti个Tx线圈上的信号，本发明实施例将训练的Tx信号矢量限制在总Tx功率PC范围内。
首先，将Tx训练信号序列改写为对角矩阵形式，即
Figure GDA00037696376600001412
然后对
Figure GDA00037696376600001413
进行奇异值分解：
Figure GDA00037696376600001414
其中∑＝diag(σ1，σ2，…，σz，…，σr)，且有σz≥σz+1≥0，z＝1，2，…，r，r为
Figure GDA00037696376600001415
的秩，σz为矩阵
Figure GDA0003769637660000151
第z个奇异值，E和F为酉矩阵。
因为E和F为酉矩阵，有
Figure GDA0003769637660000152
其中ez为E的列向量，fz为F的列向量。
令
Figure GDA0003769637660000153
其中矩阵
Figure GDA0003769637660000154
是
Figure GDA0003769637660000155
的Moorepenrose广义逆。其中
Figure GDA0003769637660000156
Figure GDA0003769637660000157
根据奇异值的性质，得到
Figure GDA0003769637660000158
由
Figure GDA0003769637660000159
联立(31)和(32)可得
Figure GDA00037696376600001510
本发明实施例只需要通过有限额电压反馈信息即可进行信道估计，基于此能够动态获取共振频率以进行资源分配，不仅能够有效降低信道容量优化的复杂度，还能够保证信道容量优化的可靠性。
实施本发明实施例，具有以下有益效果：
本发明实施例采用时频域分析模型方法，基于磁感应通信系统的等效电路模型计算通信信道和谐振频率，从而能够准确确定由于多线圈分频而产生的多共振频率，使得系统能够选择多个频率作为子载波信道，能够在不额外增加电路及功耗的情况下，有效扩展带宽，从而有效地提高了信道容量。
本发明实施例利用谐振频率作为子载波，推到出磁感应系统的信道容量，并且在给定的功率限制下，根据谐振频率对信道容量进行优化，得到最优信道容量，且本发明实施例能够使根据多个子载波来设计波形，使得系统的信道容量得到有效提高。
请参阅图5，本发明的第二实施例提供了一种基于多谐振频率的信道容量优化装置，包括：
第一计算模块10，用于计算磁感应通信系统中的发射端的阻抗矩阵和接收端的阻抗矩阵，根据发射端的阻抗矩阵和接收端的阻抗矩阵计算得到互感矩阵；
电路获取模块20，用于获取磁感应通信系统对应的等效电路模型，等效电路模型包括第一等效电路模型和第二等效电路模型；
第二计算模块30，用于根据基尔霍夫电压定律、互感矩阵以及第一等效电路模型计算磁感应通信系统的通信信道；
第三计算模块40，用于根据基尔霍夫电压定律以及第二等效电路模型构建磁感应通信系统的状态空间模型，并根据状态空间模型计算磁感应通信系统的谐振频率；
第四计算模块50，用于根据通信信道计算磁感应通信系统的接收信号向量，根据接收信号向量的协方差以及谐振频率，计算得到磁感应通信系统的信道容量；
信道容量优化模块60，用于在磁感应通信系统的总发射功率约束条件下，根据谐振频率对信道容量进行优化，得到最优信道容量。
可选地，本发明实施例可在如图2所示的磁感应通信系统中执行。请参阅图2，磁感应通信系统包括Tx端、Rx端、数字信号处理器(DSP)、数字模拟转换器(DAC)、放大器(AF)。在Tx端，使用数字信号处理器(DSP)产生发射的基带信号，数字模拟转换器(DAC)将离散的数字基带信号转换为连续的可变模拟基带信号，放大器(AF)用于调制和放大这种模拟信号，每个发射线圈与设置的一个特定谐振频率的电容器相连。检测模块用于测量每个发射线圈上的电压幅值和相位，并反馈给DSP。在接收侧，每个接收线圈与为了获得相同谐振频率而设置的接收侧电容器相连。模数转换器(ADC)对接收到的波形进行采样并发送给DSP。DSP用于对这些数据进行解码，以便进行通信和信号处理。
请参阅图3-4，为上述磁感应通信系统的两种等效电路模型。
请参阅表1，为本发明实施例提供的等效电路模型各参数定义。
表1 等效电路模型各参数定义
Figure GDA0003769637660000171
Figure GDA0003769637660000181
每个发射线圈的等效阻抗表示为
Figure GDA0003769637660000182
发射端的阻抗矩阵表示为
Figure GDA0003769637660000183
每个接收线圈的等效阻抗为
Figure GDA0003769637660000184
得到接收端的阻抗矩阵
Figure GDA0003769637660000185
互感矩阵定义为ZM，
Figure GDA0003769637660000186
表示Rx线圈对Tx线圈的影响。
Figure GDA0003769637660000187
作为本发明实施例的一种具体实施方式，第二计算模块30，具体用于：
根据基尔霍夫电压定律计算第一等效电路模型的第一电路方程，根据第一电路方程计算得到电流矩阵；根据电流矩阵和互感矩阵计算得到接收电压，以接收电压定义接收信号和发送信号，根据接收电压计算接收信号与发射信号之间的系数作为通信信道。
示例性的，根据基尔霍夫电压定律(KVL)方程，得到图3所示MIMO磁感应通信系统的电路方程为
Figure GDA0003769637660000191
计算MIMO磁感应通信系统系统中所有电路的电流，将阻抗矩阵进行反求，得到电流矩阵：
Figure GDA0003769637660000192
Figure GDA0003769637660000193
Figure GDA0003769637660000194
根据电流矩阵和互感矩阵计算接收电压：
Figure GDA0003769637660000195
本发明实施例将电压类比为信号，则接收和发射的信号分别为Y和S，即Y和S之间的系数为通信信道H，定义为
Figure GDA0003769637660000196
通信信道高度依赖于线圈间的磁通密度，从电压角度看，MIMO磁感应通信系统类似于一般的无线通信系统，无线信号传播沿一定数量的衰落信道到达接收机，不同的传播机制对路径损耗和衰落模型有不同的影响。本发明实施例定义
Figure GDA0003769637660000197
作为H的元素,其中
Figure GDA0003769637660000198
是衰落因子。
作为本发明实施例的一种具体实施方式，第三计算模块40，具体用于：
根据基尔霍夫电压定律计算第二等效电路模型的第二电路方程，根据第二电路方程计算磁感应通信系统的电容矩阵、电阻矩阵和电感矩阵，根据电容矩阵、电阻矩阵和电感矩阵构建状态空间模型，利用空间模型推导计算得到磁感应通信系统的谐振频率。
示例性的，本发明实施例根据KVL，图4所示等效电路模型的电路方程为
Figure GDA0003769637660000201
其中Sn为
Figure GDA0003769637660000202
将上述的电路方程改写为
Figure GDA0003769637660000203
其中矩阵C为MIMO磁感应系统的电容矩阵，如下式所示：
Figure GDA0003769637660000204
矩阵R为MIMO磁感应通信系统的电阻矩阵，如下式所示
Figure GDA0003769637660000205
矩阵L为MIMO磁感应通信系统的电感矩阵，如下式所示
Figure GDA0003769637660000206
引入方程
Figure GDA0003769637660000207
构建系统的状态空间模型
Figure GDA0003769637660000208
其中
Figure GDA0003769637660000211
是状态矩阵，S′＝[ST，01×N]T为系统的输入矩阵，A和B具体为
Figure GDA0003769637660000212
其中IN×N表示一个N阶单位矩阵。
利用A和B矩阵构建该系统的系统矩阵D为
Figure GDA0003769637660000213
利用构建的状态空间模型和系统矩阵，进行推导得到如下引理：：MIMO MIC系统存在Nf个谐振频率
Figure GDA0003769637660000214
它们是系统矩阵D特征值的虚部的绝对值，如下式所示。
Figure GDA0003769637660000215
作为本发明实施例的一种具体实施方式，信道容量优化模块60，具体用于：
在磁感应通信系统的总发射功率约束条件下，根据谐振频率在子信道上分配的Tx能量对信道容量的表达式进行优化，得到最优信道容量。
示例性的，考虑MIMO磁感应通信系统信的道由NT个发射线圈和NR个接收线圈组成，本发明实施例假设通信信道已经做归一化处理的并且在这个频带内频率是平坦的，经过推导可知，MIMO磁感应通信系统的接收信号向量为：
Y＝HS+N-
其中H是通信信道，N-是高斯白噪声。
经过推导可知，多频MIMO磁感应通信系统的输入输出关系表示为
Figure GDA0003769637660000216
经过推导可知，多频MIMO磁感应通信系统的信道容量表示为
Figure GDA0003769637660000217
其中，Nf为分裂频率数，mf为在通过式(11)得到频率ωf的H(ωf)H(ωf)H正特征值的个数，λk(ωf)是这些特征值之一，
Figure GDA0003769637660000221
为信号向量在ωf频率处的协方差，N0为通带噪声谱功率。
在本发明实施例中，对于磁感应通信系统，多线圈发射机的主要目标是设计S在总发射功率约束下达到最大的信道容量，即
Figure GDA0003769637660000222
Figure GDA0003769637660000223
SHS≤PC
其中PC为最大允许的总传输功率，
Figure GDA0003769637660000224
为S的协方差.
Figure GDA0003769637660000225
为向量S的元素。
可选地，由于发射端可以接入MIMO磁感应通信系统的子信道，在总发射功率约束下，能够跨越子信道分配可变能量，使信道容量最大化。因此公式(12)的优化如下：
Figure GDA0003769637660000226
其中
Figure GDA0003769637660000227
表示表示在频率为ωf的第k个子信道上分配的Tx能量，满足
Figure GDA0003769637660000228
且最优能量分配方案
Figure GDA0003769637660000229
满足
Figure GDA00037696376600002210
其中η*是一个常数，根据特征值的集合的迭代计算和子信道分配能量的非负性获得，x+是指
Figure GDA00037696376600002211
作为本发明实施例的一种具体实施方式，根据谐振频率在子信道上分配的Tx能量对信道容量进行优化，得到最优信道容量，具体为：
根据通信信号将所有特征值按照大小降序排序，得到特征值排序集合；
通过对特征值排序结合进行迭代计算和子信道分配功率的非负性计算得到优化常数；
根据优化常数对子信道进行功率分配，当分配的子信道功率为负时将子信道丢弃，并重新运行迭代计算器加一；当分配给磁电感通信系统各个子型号的功率为负时，得到最优信道容量。
在本发明实施例中，对自适应多频信道容量优化方案进行迭代设计，，在已知H的前提下，将所有的特征值λk(ωf)按大小降序排列，从而形成一个
Figure GDA0003769637660000231
的集合，并利用γl来表示相关的λl的功率分配。设置计数器p等于1，在每次迭代时，根据受约束的功率计算常数1/η*
Figure GDA0003769637660000232
每个模式在每次迭代中分配的功率按下式计算
Figure GDA0003769637660000233
若分配给增益最小信道的能量为负，如γi＜0，通过设置γi＝0丢弃这个信道，并重新运行迭代计数器p加1的方案。当分配给各磁感应通信系统子信道的功率为非负时，为最优的资源分配方案，即得到最优信道容量。
本发明实施例还提供了一种信道估计方法，具体为：
当Tx和Rx线圈固定时，通道几乎是固定的。在这种情况下，Rx线圈的电压通过互感与Tx线圈相互作用，信道估计主要是根据Rx侧的电压和负载计算互感，得到磁信道信息，接收端ADC模块对接收到的信号矢量进行采样。
在磁感应通信系统中，信道由接收机利用发射机发出的训练信号进行估计。接收机知道训练信号序列
Figure GDA0003769637660000234
定义
Figure GDA0003769637660000235
为估计的信道，则接收信号为
Figure GDA0003769637660000236
多线圈发射机的主要目标是在功率受限时控制
Figure GDA0003769637660000237
使得
Figure GDA0003769637660000238
误差达到最小值，即
Figure GDA0003769637660000239
Figure GDA00037696376600002310
Figure GDA0003769637660000241
其中
Figure GDA0003769637660000242
为为第Ti个Tx线圈上的信号，本发明实施例将训练的Tx信号矢量限制在总Tx功率PC范围内。
首先，将Tx训练信号序列改写为对角矩阵形式，即
Figure GDA0003769637660000243
然后对
Figure GDA0003769637660000244
进行奇异值分解：
Figure GDA0003769637660000245
其中∑＝diag(σ1，σ2，…，σz，…，σr)，且有σz≥σz+1≥0，z＝1，2，…，r，r为
Figure GDA00037696376600002417
的秩，σz为矩阵
Figure GDA0003769637660000247
第z个奇异值，E和F为酉矩阵。
因为E和F为酉矩阵，有
Figure GDA0003769637660000248
其中ez为E的列向量，fz为F的列向量。
令
Figure GDA0003769637660000249
其中矩阵
Figure GDA00037696376600002410
是
Figure GDA00037696376600002411
的Moorepenrose广义逆。其中
Figure GDA00037696376600002412
Figure GDA00037696376600002413
根据奇异值的性质，得到
Figure GDA00037696376600002414
由
Figure GDA00037696376600002415
联立(31)和(32)可得
Figure GDA00037696376600002416
本发明实施例只需要通过有限额电压反馈信息即可进行信道估计，基于此能够动态获取共振频率以进行资源分配，不仅能够有效降低信道容量优化的复杂度，还能够保证信道容量优化的可靠性。
实施本发明实施例，具有以下有益效果：
本发明实施例采用时频域分析模型方法，基于磁感应通信系统的等效电路模型计算通信信道和谐振频率，从而能够准确确定由于多线圈分频而产生的多共振频率，使得系统能够选择多个频率作为子载波信道，能够在不额外增加电路及功耗的情况下，有效扩展带宽，从而有效地提高了信道容量。
本发明实施例利用谐振频率作为子载波，推到出磁感应系统的信道容量，并且在给定的功率限制下，根据谐振频率对信道容量进行优化，得到最优信道容量，且本发明实施例能够使根据多个子载波来设计波形，使得系统的信道容量得到有效提高。
以上是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。

Claims (6)

1.一种基于多谐振频率的信道容量优化方法，其特征在于，包括：
计算磁感应通信系统中的发射端的阻抗矩阵和接收端的阻抗矩阵，根据所述发射端的阻抗矩阵和接收端的阻抗矩阵计算得到互感矩阵；
获取所述磁感应通信系统对应的等效电路模型，所述等效电路模型包括第一等效电路模型和第二等效电路模型；
根据基尔霍夫电压定律、所述互感矩阵以及所述第一等效电路模型计算所述磁感应通信系统的通信信道，包括：根据基尔霍夫电压定律方程，得到所述第一等效电路模型的第一电路方程为：
Figure 397105DEST_PATH_IMAGE001
，
其中，
Figure 102673DEST_PATH_IMAGE002
为发射端第
Figure 873052DEST_PATH_IMAGE003
线圈的发射信号，
Figure 686550DEST_PATH_IMAGE004
为发射端第
Figure 734140DEST_PATH_IMAGE005
个线圈的电阻，
Figure 106216DEST_PATH_IMAGE006
为接收 端第
Figure 176940DEST_PATH_IMAGE007
个线圈的电阻；
Figure 292663DEST_PATH_IMAGE008
为接收端第Ru个线圈的电流；
将所述第一电路方程改写为矩阵方程：
Figure 755612DEST_PATH_IMAGE009
根据上述的矩阵方程，得到发送电流向量：
Figure 236272DEST_PATH_IMAGE010
以及接收电流向量：
Figure 341763DEST_PATH_IMAGE011
根据所述接收电流矩阵和接收互感矩阵计算得到接收电压向量：
Figure 261177DEST_PATH_IMAGE012
将电压类比为信号，得到接收信号为Y和发射信号为S，将接收信号Y和发射信号S之间 的系数作为通信信道H，通信信道为：
Figure 533894DEST_PATH_IMAGE013
；
其中，通信信道H的元素为
Figure 247772DEST_PATH_IMAGE014
,
Figure 355405DEST_PATH_IMAGE015
表示第i个发射机与第u个接收机 之间的衰减衰落因子；
根据基尔霍夫电压定律以及所述第二等效电路模型构建所述磁感应通信系统的状态空间模型，并根据所述状态空间模型计算所述磁感应通信系统的谐振频率，包括：根据基尔霍夫电压定律，计算得到所述第二等效电路模型的第二电路方程为：
Figure 812931DEST_PATH_IMAGE016
其中
Figure 925506DEST_PATH_IMAGE017
为：
Figure 75865DEST_PATH_IMAGE018
将所述第二电路方程改写为：
Figure 405215DEST_PATH_IMAGE019
（1）
其中矩阵C为所述磁感应通信系统的电容矩阵，如下式所示：
Figure 400853DEST_PATH_IMAGE020
矩阵R为所述磁感应通信系统的电阻矩阵，如下式所示
Figure 69731DEST_PATH_IMAGE021
矩阵L为所述 磁感应通信系统的电感矩阵，如下式所示
Figure 889526DEST_PATH_IMAGE022
引入方程：
Figure 440594DEST_PATH_IMAGE023
（2）
构建所述磁感应通信系统的状态空间模型：
Figure 505502DEST_PATH_IMAGE024
（3）
其中
Figure 91204DEST_PATH_IMAGE025
是状态矩阵，
Figure 819250DEST_PATH_IMAGE026
为系统的输入矩阵，矩阵A和矩阵 B具体为：
Figure 60876DEST_PATH_IMAGE027
其中
Figure 726213DEST_PATH_IMAGE028
表示一个N阶单位矩阵；
利用矩阵A和矩阵B矩阵构建所述磁感应通信系统的系统矩阵D为
Figure 166421DEST_PATH_IMAGE029
利用所述状态空间模型和所述系统矩阵D，进行推导得到所述磁感应通信系统存在
Figure 767167DEST_PATH_IMAGE030
个谐振频率(
Figure 797221DEST_PATH_IMAGE031
)，所述谐振频率是系统矩阵D的虚部的绝对值，所述谐振频 率为：
Figure 203931DEST_PATH_IMAGE032
（11）；
根据所述通信信道计算所述磁感应通信系统的接收信号向量，根据所述接收信号向量的协方差以及所述谐振频率，计算得到所述磁感应通信系统的信道容量；
在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率对所述信道容量进行优化，得到最优信道容量，包括：推导得到所述磁感应通信系统的接收信号向量为：
Figure 764225DEST_PATH_IMAGE033
其中H是通信信道矩阵, N ̅是高斯白噪声向量；
所述磁感应通信系统的输入输出关系表示为：
Figure 332610DEST_PATH_IMAGE034
所述磁感应通信系统的信道容量表示为：
Figure 112609DEST_PATH_IMAGE035
（12）
其中，
Figure 57432DEST_PATH_IMAGE030
为分裂频率数，
Figure 409916DEST_PATH_IMAGE036
为在通过式(11)得到频率
Figure 149202DEST_PATH_IMAGE037
的
Figure 649453DEST_PATH_IMAGE038
特征值的 个数，
Figure 896501DEST_PATH_IMAGE039
为特征值，
Figure 228126DEST_PATH_IMAGE040
为
Figure 403892DEST_PATH_IMAGE037
频率对应的协方差，
Figure 955221DEST_PATH_IMAGE041
为通带噪声谱功率；
为求得总发射功率约束下达到最大的信道容量，即目标函数为：
Figure 241846DEST_PATH_IMAGE042
Figure 631239DEST_PATH_IMAGE043
；
Figure 470583DEST_PATH_IMAGE044
其中
Figure 679847DEST_PATH_IMAGE045
为最大允许的总发射功率，
Figure 832480DEST_PATH_IMAGE046
为发射信号S的协方差矩阵，
Figure 577844DEST_PATH_IMAGE047
为发射信号S向 量的元素；
将公式（12）优化为：
Figure 95413DEST_PATH_IMAGE048
；
其中
Figure 57553DEST_PATH_IMAGE049
表示在频率为
Figure 685980DEST_PATH_IMAGE037
的第k个子信道上分配的发送能量，满足
Figure 282921DEST_PATH_IMAGE050
，且最优能量分配方案
Figure 971392DEST_PATH_IMAGE051
满足
Figure 155248DEST_PATH_IMAGE052
；
其中，
Figure 587367DEST_PATH_IMAGE053
是常数，根据特征值的集合的迭代计算和子信道分配能量的非负性获得，
Figure 307323DEST_PATH_IMAGE054
是指
Figure 229012DEST_PATH_IMAGE055
。
2.如权利要求1所述的基于多谐振频率的信道容量优化方法，其特征在于，所述在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率对所述信道容量进行优化，得到最优信道容量，具体为：
在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率在所述磁感应通信系统的子信道上分配的可变能量对所述信道容量的表达式进行优化，得到最优信道容量。
3.如权利要求2所述的基于多谐振频率的信道容量优化方法，其特征在于，所述根据所述谐振频率在所述磁感应通信系统的子信道上分配的可变能量对所述信道容量进行优化，得到最优信道容量，具体为：
根据所述通信信道将所有特征值按照大小降序排序，得到特征值排序集合；
通过对所述特征值排序结合进行迭代计算和子信道分配功率的非负性计算得到优化常数；
根据所述优化常数对子信道进行功率分配，当分配的子信道功率为负时将所述子信道丢弃，并重新运行迭代计算器加一；当分配给所述磁感应通信系统各个子信道的功率为非负时，得到最优信道容量。
4.一种基于多谐振频率的信道容量优化装置，其特征在于，包括：
第一计算模块，用于计算磁感应通信系统中的发射端的阻抗矩阵和接收端的阻抗矩阵，根据所述发射端的阻抗矩阵和接收端的阻抗矩阵计算得到互感矩阵；
电路获取模块，用于获取所述磁感应通信系统对应的等效电路模型，所述等效电路模型包括第一等效电路模型和第二等效电路模型；
第二计算模块，用于根据基尔霍夫电压定律、所述互感矩阵以及所述第一等效电路模型计算所述磁感应通信系统的通信信道，具体用于：根据基尔霍夫电压定律方程，得到所述第一等效电路模型的第一电路方程为：
Figure 900165DEST_PATH_IMAGE001
，
其中，
Figure 374789DEST_PATH_IMAGE002
为发射端第
Figure 447787DEST_PATH_IMAGE003
线圈的发射信号，
Figure 478060DEST_PATH_IMAGE004
为发射端第
Figure 636509DEST_PATH_IMAGE005
个线圈的电阻，
Figure 348113DEST_PATH_IMAGE006
为接收 端第
Figure 42662DEST_PATH_IMAGE007
个线圈的电阻；
Figure 243836DEST_PATH_IMAGE008
为接收端第Ru个线圈的电流；
将所述第一电路方程改写为矩阵方程：
Figure 624002DEST_PATH_IMAGE009
根据上述的矩阵方程，得到发送电流向量：
Figure 139297DEST_PATH_IMAGE010
以及接收电流向量：
Figure 186887DEST_PATH_IMAGE011
根据所述接收电流矩阵和接收互感矩阵计算得到接收电压向量：
Figure 57498DEST_PATH_IMAGE012
将电压类比为信号，得到接收信号为Y和发射信号为S，将接收信号Y和发射信号S之间 的系数作为通信信道H，通信信道为：
Figure 987277DEST_PATH_IMAGE056
；
其中，通信信道H的元素为
Figure 103000DEST_PATH_IMAGE014
,
Figure 506562DEST_PATH_IMAGE015
表示第i个发射机与第u个接收机 之间的衰减衰落因子；
根据基尔霍夫电压定律以及所述第二等效电路模型构建所述磁感应通信系统的状态空间模型，并根据所述状态空间模型计算所述磁感应通信系统的谐振频率，包括：根据基尔霍夫电压定律，计算得到所述第二等效电路模型的第二电路方程为：
Figure 49539DEST_PATH_IMAGE016
其中
Figure 669876DEST_PATH_IMAGE017
为：
Figure 323711DEST_PATH_IMAGE018
将所述第二电路方程改写为：
Figure 924324DEST_PATH_IMAGE019
（1）
其中矩阵C为所述磁感应通信系统的电容矩阵，如下式所示：
Figure 638202DEST_PATH_IMAGE020
矩阵R为所述磁感应通信系统的电阻矩阵，如下式所示
Figure 480256DEST_PATH_IMAGE021
矩阵L为所述 磁感应通信系统的电感矩阵，如下式所示
Figure 501564DEST_PATH_IMAGE022
引入方程：
Figure 378253DEST_PATH_IMAGE023
（2）
构建所述磁感应通信系统的状态空间模型：
Figure 263033DEST_PATH_IMAGE024
（3）
其中
Figure 153235DEST_PATH_IMAGE025
是状态矩阵，
Figure 352135DEST_PATH_IMAGE026
为系统的输入矩阵，矩阵A和矩阵 B具体为：
Figure 83331DEST_PATH_IMAGE027
其中
Figure 139012DEST_PATH_IMAGE028
表示一个N阶单位矩阵；
利用矩阵A和矩阵B矩阵构建所述磁感应通信系统的系统矩阵D为
Figure 457123DEST_PATH_IMAGE029
利用所述状态空间模型和所述系统矩阵D，进行推导得到所述磁感应通信系统存在
Figure 584347DEST_PATH_IMAGE030
个谐振频率(
Figure 107733DEST_PATH_IMAGE031
)，所述谐振频率是系统矩阵D的虚部的绝对值，所述谐振频 率为：
Figure 334315DEST_PATH_IMAGE032
（11）；
根据所述通信信道计算所述磁感应通信系统的接收信号向量，根据所述接收信号向量的协方差以及所述谐振频率，计算得到所述磁感应通信系统的信道容量；
在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率对所述信道容量进行优化，得到最优信道容量，包括：推导得到所述磁感应通信系统的接收信号向量为：
Figure 142652DEST_PATH_IMAGE033
其中H是通信信道矩阵, N ̅是高斯白噪声向量；
所述磁感应通信系统的输入输出关系表示为：
Figure 745671DEST_PATH_IMAGE034
所述磁感应通信系统的信道容量表示为：
Figure 248197DEST_PATH_IMAGE035
（12）
其中，
Figure 599675DEST_PATH_IMAGE030
为分裂频率数，
Figure 623869DEST_PATH_IMAGE036
为在通过式(11)得到频率
Figure 30580DEST_PATH_IMAGE037
的
Figure 325295DEST_PATH_IMAGE038
特征值的 个数，
Figure 159259DEST_PATH_IMAGE039
为特征值，
Figure 673679DEST_PATH_IMAGE040
为
Figure 618501DEST_PATH_IMAGE037
频率对应的协方差，
Figure 970985DEST_PATH_IMAGE041
为通带噪声谱功率；
为求得总发射功率约束下达到最大的信道容量，即目标函数为：
Figure 975850DEST_PATH_IMAGE042
Figure 476102DEST_PATH_IMAGE043
；
Figure 717291DEST_PATH_IMAGE044
其中
Figure 986598DEST_PATH_IMAGE045
为最大允许的总发射功率，
Figure 100047DEST_PATH_IMAGE046
为发射信号S的协方差矩阵，
Figure 87595DEST_PATH_IMAGE047
为发射信号S向 量的元素；
将公式（12）优化为：
Figure 374220DEST_PATH_IMAGE048
；
其中
Figure 265078DEST_PATH_IMAGE049
表示在频率为
Figure 346166DEST_PATH_IMAGE037
的第k个子信道上分配的发送能量，满足
Figure 821010DEST_PATH_IMAGE050
，且最优能量分配方案
Figure 849009DEST_PATH_IMAGE051
满足
Figure 92908DEST_PATH_IMAGE052
；
其中，
Figure 109012DEST_PATH_IMAGE053
是常数，根据特征值的集合的迭代计算和子信道分配能量的非负性获得，
Figure 805573DEST_PATH_IMAGE054
是指
Figure 371684DEST_PATH_IMAGE055
。
5.如权利要求4所述的基于多谐振频率的信道容量优化装置，其特征在于，所述信道容量优化模块，具体用于：
在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率在所述磁感应通信系统的子信道上分配的可变能量对所述信道容量的表达式进行优化，得到最优信道容量。
6.如权利要求5所述的基于多谐振频率的信道容量优化装置，其特征在于，所述根据所述谐振频率在所述磁感应通信系统的子信道上分配的可变能量对所述信道容量进行优化，得到最优信道容量，具体为：
根据所述通信信道将所有特征值按照大小降序排序，得到特征值排序集合；
通过对所述特征值排序结合进行迭代计算和子信道分配功率的非负性计算得到优化常数；
根据所述优化常数对子信道进行功率分配，当分配的子信道功率为负时将所述子信道丢弃，并重新运行迭代计算器加一；当分配给所述磁感应通信系统各个子信道的功率为非负时，得到最优信道容量。
CN202110886466.6A 2021-08-03 2021-08-03 一种基于多谐振频率的信道容量优化方法及装置 Active CN113660014B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110886466.6A CN113660014B (zh) 2021-08-03 2021-08-03 一种基于多谐振频率的信道容量优化方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110886466.6A CN113660014B (zh) 2021-08-03 2021-08-03 一种基于多谐振频率的信道容量优化方法及装置

Publications (2)

Publication Number Publication Date
CN113660014A CN113660014A (zh) 2021-11-16
CN113660014B true CN113660014B (zh) 2023-03-24

Family

ID=78478297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110886466.6A Active CN113660014B (zh) 2021-08-03 2021-08-03 一种基于多谐振频率的信道容量优化方法及装置

Country Status (1)

Country Link
CN (1) CN113660014B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115801063B (zh) * 2022-10-21 2024-05-14 华中科技大学 一种地下磁感应无线通信方法、装置及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106849221A (zh) * 2016-12-30 2017-06-13 深圳市尧元科技有限公司 基于磁共振多设备无线充电功率最大化传输方法及系统
CN112235027A (zh) * 2020-09-14 2021-01-15 德清阿尔法创新研究院 一种基于磁谐振的mimo充电并行通信方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10249950B1 (en) * 2017-09-16 2019-04-02 Searete Llc Systems and methods for reduced control inputs in tunable meta-devices
CN209088633U (zh) * 2018-03-23 2019-07-09 天津工业大学 一种新型的高密度无线电能传输线圈结构
US11374618B2 (en) * 2019-10-19 2022-06-28 Meta Platforms, Inc. Surface waveguide with a two-dimensional conductive surface surrounded by a conductive wall
CN112202476B (zh) * 2020-09-29 2021-09-10 西安电子科技大学 基于线圈圆阵列的涡旋电磁波近场无线通信方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106849221A (zh) * 2016-12-30 2017-06-13 深圳市尧元科技有限公司 基于磁共振多设备无线充电功率最大化传输方法及系统
CN112235027A (zh) * 2020-09-14 2021-01-15 德清阿尔法创新研究院 一种基于磁谐振的mimo充电并行通信方法

Also Published As

Publication number Publication date
CN113660014A (zh) 2021-11-16

Similar Documents

Publication Publication Date Title
US9473223B2 (en) Methods for optimal collaborative MIMO-SDMA
CN103069759B (zh) Mimo信道状态信息估计的方法和接收器
US9236922B2 (en) Method of communicating between base station and terminal based on interference alignment in multi-cell multi-user multiple-input multiple-output (MIMO) interference channel and method and apparatus of communication using interference alignment and block successive interference pre-cancellation in multi-user multiple-input multiple-output interference channel
KR101312015B1 (ko) 무선 통신 시스템 내의 무선 디바이스에 전송될 프리코딩 데이터에 대한 프리코딩 벡터를 결정하는 방법 및 장치
US8649457B2 (en) Precoding process for a transmitter of a MU-MIMO communication system
CN101383633A (zh) 无线通信系统内对接收信号预测的系统和方法
KR101546177B1 (ko) 잡음 평가 필터
CN113162663B (zh) 非理想信道信息下的鲁棒多载波mimo无线携能系统设计方法
Couillet et al. Random beamforming over quasi-static and fading channels: A deterministic equivalent approach
US20100246715A1 (en) Wireless communication method and apparatus
CN103475401A (zh) 一种下行波束赋形方法与装置
CN102571296A (zh) 一种预编码的方法及装置
KR20120126572A (ko) 다중 사용자 다중 입출력 인지 무선 통신 시스템에서 인지 무선 기지국 및 인지 무선 기지국의 통신 방법 및 장치
EP2429087B1 (en) Interference cancellation method and data transmission device
CN103249080A (zh) 一种确定基站的天线校准系数的方法、系统以及装置
Guthy et al. Large system analysis of projection based algorithms for the MIMO broadcast channel
CN113660014B (zh) 一种基于多谐振频率的信道容量优化方法及装置
JP4227648B2 (ja) 無線通信システムにおけるトレーニングシーケンスの前置フィルタリング方法
CN108012272B (zh) 基于认知网络中动态功率分配的干扰对齐方法
CN107248878A (zh) 一种预编码的方法及装置
CN105099530B (zh) 认知无线电mimo-ofdm系统中基于认知用户漏功率的干扰抑制预编码方法
CN107979398A (zh) 结合天线选择和干扰对齐的联合优化算法
CN102215095A (zh) 一种实现数据重传控制的方法和系统
Kaleva et al. Decentralized sum MSE minimization for coordinated multi-point transmission
Dikmen et al. A new method in pilot reuse factor selection in spectrum efficient massive MIMO systems

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant