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

1.本发明涉及通信技术领域，尤其是涉及一种基于多谐振频率的信道容量优化方法及装置。


背景技术：

2.磁感应通信(mic)是在射频信号无法穿透地下或水下环境等介质的具有挑战性或恶劣环境下，一种具有发展前景的新型无线通信模式。mic系统中部署了几个耦合线圈，可以在相同的谐振频率下同时交换无线电源和信息。在这种情况下，mic不仅提供了可预测、稳定的信号传播通道，而且传输延迟可以忽略不计。mic由于其性价比高，被广泛应用于地下传感器网络、自主水下航行器、恶劣环境下的机器人以及一些军事武器。磁感应通信的主要挑战是有限的信道容量，主要受到接收功率的影响。具体来说，接收功率和信道容量是由线圈半径和传输距离决定的。在文献中首先对单载频单对收发线圈(siso)的信道容量进行了分析。在此基础上，研究了单输入多输出(simo)、多输入单输入(miso)和多输入多输出(mimo)模型的信道容量。然而，mimo信道容量的提高仍然是有限的。
3.现有的磁感应通信系统的信道容量优化方法主要有两类:一类是高质量的硬件和电路设计，另一类是开发优化方法。对于第一类，mimo磁感应系统为原磁场或电流分布提供了更多的自由度，等效电路模型(ecm)是电路设计和性能分析的主要工具。在文献中，通过仿真和实际实验对双线圈和三线圈系统的性能进行了比较，证明了tx和rx线圈强耦合时存在分频现象。然后，nguyen等人指出通过分裂频率的多通道通信可以用于数据传输。dionigi等人设计了一种用于无线功率传输和近场通信的多波段谐振器。第二类是设计适当的信道容量优化方法，控制系统的电流或电压。这种方法不需要额外的电路。在文献中，为不同的mi继电器和收发器线圈分配不同的谐振频率，以计算系统的信道容量。
4.现有的多频率磁mimo信道容量优化方法存在以下问题：
5.第一类方法需要额外的电路，例如控制电路和阻抗匹配网络，导致需要引入额外的功耗，而增加系统的复杂性。
6.第二类方法需要寄生电容和其他电路信息以获取共振频率，增加了计算复杂度。


技术实现要素：

7.本发明提供了一种基于多谐振频率的信道容量优化方法及装置，以解决现有的空信道容量优化方法需要引入额外的功耗以及增加计算复杂度的问题。
8.本发明的第一实施例提供了一种基于多谐振频率的信道容量优化方法，包括：
9.计算磁感应通信系统中的发射端的阻抗矩阵和接收端的阻抗矩阵，根据所述发射端的阻抗矩阵和接收端的阻抗矩阵计算得到互感矩阵；
10.获取所述磁感应通信系统对应的等效电路模型，所述等效电路模型包括第一等效电路模型和第二等效电路模型；
11.根据基尔霍夫电压定律、所述互感矩阵以及所述第一等效电路模型计算所述磁感
应通信系统的通信信道；
12.根据基尔霍夫电压定律以及所述第二等效电路模型构建所述磁感应通信系统的状态空间模型，并根据所述状态空间模型计算所述磁感应通信系统的谐振频率；
13.根据所述通信信道计算所述磁感应通信系统的接收信号向量，根据所述接收信号向量的协方差以及所述谐振频率，计算得到所述磁感应通信系统的信道容量；
14.在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率对所述信道容量进行优化，得到最优信道容量。
15.进一步的，所述根据基尔霍夫电压定律、所述互感矩阵以及所述第一等效电路模型计算所述磁感应通信系统的通信信道，具体为：
16.根据基尔霍夫电压定律计算所述第一等效电路模型的第一电路方程，根据所述第一电路方程计算得到电流矩阵；根据所述电流矩阵和所述互感矩阵计算得到接收电压，以所述接收电压定义接收信号和发送信号，根据所述接收电压计算所述接收信号与所述发射信号之间的系数作为通信信道。
17.进一步的，所述根据基尔霍夫电压定律以及所述第二等效电路模型构建所述磁感应通信系统的状态空间模型，并根据所述状态空间模型计算所述磁感应通信系统的谐振频率，具体为：
18.根据基尔霍夫电压定律计算所述第二等效电路模型的第二电路方程，根据所述第二电路方程计算所述磁感应通信系统的电容矩阵、电阻矩阵和电感矩阵，根据所述电容矩阵、电阻矩阵和电感矩阵构建所述状态空间模型，利用所述空间模型推导计算得到所述磁感应通信系统的谐振频率。
19.进一步的，所述在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率对所述信道容量进行优化，得到最优信道容量，具体为：
20.在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率在子信道上分配的tx能量对所述信道容量的表达式进行优化，得到最优信道容量。
21.进一步的，所述根据所述谐振频率在子信道上分配的tx能量对所述信道容量进行优化，得到最优信道容量，具体为：
22.根据所述通信信号将所有特征值按照大小降序排序，得到特征值排序集合；
23.通过对所述特征值排序结合进行迭代计算和子信道分配功率的非负性计算得到优化常数；
24.根据所述优化常数对子信道进行功率分配，当分配的子信道功率为负时将所述子信道丢弃，并重新运行迭代计算器加一；当分配给所述磁电感通信系统各个子型号的功率为负时，得到最优信道容量。
25.本发明的第二实施例提供了一种基于多谐振频率的信道容量优化装置，包括：
26.第一计算模块，用于计算磁感应通信系统中的发射端的阻抗矩阵和接收端的阻抗矩阵，根据所述发射端的阻抗矩阵和接收端的阻抗矩阵计算得到互感矩阵；
27.电路获取模块，用于获取所述磁感应通信系统对应的等效电路模型，所述等效电路模型包括第一等效电路模型和第二等效电路模型；
28.第二计算模块，用于根据基尔霍夫电压定律、所述互感矩阵以及所述第一等效电路模型计算所述磁感应通信系统的通信信道；
29.第三计算模块，用于根据基尔霍夫电压定律以及所述第二等效电路模型构建所述磁感应通信系统的状态空间模型，并根据所述状态空间模型计算所述磁感应通信系统的谐振频率；
30.第四计算模块，用于根据所述通信信道计算所述磁感应通信系统的接收信号向量，根据所述接收信号向量的协方差以及所述谐振频率，计算得到所述磁感应通信系统的信道容量；
31.信道容量优化模块，用于在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率对所述信道容量进行优化，得到最优信道容量。
32.进一步的，所述第二计算模块，具体用于：
33.根据基尔霍夫电压定律计算所述第一等效电路模型的第一电路方程，根据所述第一电路方程计算得到电流矩阵；根据所述电流矩阵和所述互感矩阵计算得到接收电压，以所述接收电压定义接收信号和发送信号，根据所述接收电压计算所述接收信号与所述发射信号之间的系数作为通信信道。
34.进一步的，所述第三计算模块，具体用于：
35.根据基尔霍夫电压定律计算所述第二等效电路模型的第二电路方程，根据所述第二电路方程计算所述磁感应通信系统的电容矩阵、电阻矩阵和电感矩阵，根据所述电容矩阵、电阻矩阵和电感矩阵构建所述状态空间模型，利用所述空间模型推导计算得到所述磁感应通信系统的谐振频率。
36.进一步的，所述信道容量优化模块，具体用于：
37.在所述磁感应通信系统的总发射功率约束条件下，根据所述谐振频率在子信道上分配的tx能量对所述信道容量的表达式进行优化，得到最优信道容量。
38.进一步的，所述根据所述谐振频率在子信道上分配的tx能量对所述信道容量进行优化，得到最优信道容量，具体为：
39.根据所述通信信号将所有特征值按照大小降序排序，得到特征值排序集合；
40.通过对所述特征值排序结合进行迭代计算和子信道分配功率的非负性计算得到优化常数；
41.根据所述优化常数对子信道进行功率分配，当分配的子信道功率为负时将所述子信道丢弃，并重新运行迭代计算器加一；当分配给所述磁电感通信系统各个子型号的功率为负时，得到最优信道容量。
42.本发明实施例采用时频域分析模型方法，基于磁感应通信系统的等效电路模型计算通信信道和谐振频率，从而能够准确确定由于多线圈分频而产生的多共振频率，使得系统能够选择多个频率作为子载波信道，能够在不额外增加电路及功耗的情况下，有效扩展带宽，从而有效地提高了信道容量。
43.本发明实施例利用谐振频率作为子载波，推到出磁感应系统的信道容量，并且在给定的功率限制下，根据谐振频率对信道容量进行优化，得到最优信道容量，且本发明实施例能够使根据多个子载波来设计波形，使得系统的信道容量得到有效提高。
附图说明
44.图1是本发明实施例提供的一种基于多谐振频率的信道容量优化方法的流程示意
图；
45.图2是本发明实施例提供的磁感应通信系统的结构式示意图；
46.图3是本发明实施例提供的第一等效电路模型的示意图；
47.图4是本发明实施例提供的第二等效电路模型的示意图；
48.图5是本发明实施例提供的一种基于多谐振频率的信道容量优化装置的结构示意图。
具体实施方式
49.下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
50.在本技术的描述中，需要理解的是，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本技术的描述中，除非另有说明，“多个”的含义是两个或两个以上。
51.在本技术的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本技术中的具体含义。
52.请参阅图1
‑
4，在本发明的第一实施例中，本发明的第一实施例提供了如1所示的一种基于多谐振频率的信道容量优化方法，包括：
53.s1、计算磁感应通信系统中的发射端的阻抗矩阵和接收端的阻抗矩阵，根据发射端的阻抗矩阵和接收端的阻抗矩阵计算得到互感矩阵；
54.s2、获取磁感应通信系统对应的等效电路模型，等效电路模型包括第一等效电路模型和第二等效电路模型；
55.s3、根据基尔霍夫电压定律、互感矩阵以及第一等效电路模型计算磁感应通信系统的通信信道；
56.s4、根据基尔霍夫电压定律以及第二等效电路模型构建磁感应通信系统的状态空间模型，并根据状态空间模型计算磁感应通信系统的谐振频率；
57.s5、根据通信信道计算磁感应通信系统的接收信号向量，根据接收信号向量的协方差以及谐振频率，计算得到磁感应通信系统的信道容量；
58.s6、在磁感应通信系统的总发射功率约束条件下，根据谐振频率对信道容量进行优化，得到最优信道容量。
59.可选地，本发明实施例可在如图2所示的磁感应通信系统中执行。请参阅图2，磁感应通信系统包括tx端、rx端、数字信号处理器(dsp)、数字模拟转换器(dac)、放大器(af)。在tx端，使用数字信号处理器(dsp)产生发射的基带信号，数字模拟转换器(dac)将离散的数字基带信号转换为连续的可变模拟基带信号，放大器(af)用于调制和放大这种模拟信号，
每个发射线圈与设置的一个特定谐振频率的电容器相连。检测模块用于测量每个发射线圈上的电压幅值和相位，并反馈给dsp。在接收侧，每个接收线圈与为了获得相同谐振频率而设置的接收侧电容器相连。模数转换器(adc)对接收到的波形进行采样并发送给dsp。dsp用于对这些数据进行解码，以便进行通信和信号处理。
60.请参阅图3
‑
4，为上述磁感应通信系统的两种等效电路模型。
61.请参阅表1，为本发明实施例提供的等效电路模型各参数定义。
62.表1等效电路模型各参数定义
63.[0064][0065]
每个发射线圈的等效阻抗表示为
[0066][0067]
发射端的阻抗矩阵表示为
[0068][0069]
每个接收线圈的等效阻抗为
[0070][0071]
得到接收端的阻抗矩阵
[0072][0073]
互感矩阵定义为z
m
，表示rx线圈对tx线圈的影响。
[0074][0075]
作为本发明实施例的一种具体实施方式，根据基尔霍夫电压定律、互感矩阵以及第一等效电路模型计算磁感应通信系统的通信信道，具体为：
[0076]
根据基尔霍夫电压定律计算第一等效电路模型的第一电路方程，根据第一电路方程计算得到电流矩阵；根据电流矩阵和互感矩阵计算得到接收电压，以接收电压定义接收信号和发送信号，根据接收电压计算接收信号与发射信号之间的系数作为通信信道。
[0077]
示例性的，根据基尔霍夫电压定律(kvl)方程，得到图3所示m i mo磁感应通信系统的电路方程为
[0078][0079]
计算m i mo磁感应通信系统系统中所有电路的电流，将阻抗矩阵进行反求，得到电流矩阵：
[0080][0081][0082][0083]
根据电流矩阵和护肝矩阵计算接收电压：
[0084][0085]
本发明实施例将电压类比为信号，则接收和发射的信号分别为y和s，即y和s之间的系数为通信信道h，定义为
[0086]
通信信道高度依赖于线圈间的磁通密度，从电压角度看，m i mo磁感应通信系统类似于一般的无线通信系统，无线信号传播沿一定数量的衰落信道到达接收机，不同的传播机制对路径损耗和衰落模型有不同的影响。本发明实施例定义作为h的元素,其中是衰落因子。
[0087]
作为本发明实施例的一种具体实施方式，根据基尔霍夫电压定律以及第二等效电路模型构建磁感应通信系统的状态空间模型，并根据状态空间模型计算磁感应通信系统的谐振频率，具体为：
[0088]
根据基尔霍夫电压定律计算第二等效电路模型的第二电路方程，根据第二电路方程计算磁感应通信系统的电容矩阵、电阻矩阵和电感矩阵，根据电容矩阵、电阻矩阵和电感矩阵构建状态空间模型，利用空间模型推导计算得到磁感应通信系统的谐振频率。
[0089]
示例性的，本发明实施例根据kvl，图4所示等效电路模型的电路方程为
[0090][0091]
其中s
n
为
[0092][0093]
将上述的电路方程改写为
[0094][0095]
其中矩阵c为m i mo磁感应系统的电容矩阵，如下式所示：
[0096][0097]
矩阵r为m i mo磁感应通信系统的电阻矩阵，如下式所示
[0098][0099]
矩阵l为m i mo磁感应通信系统的电感矩阵，如下式所示
[0100][0101]
引入方程
[0102][0103]
构建系统的状态空间模型
[0104][0105]
其中是状态矩阵，s
′
＝[s
t
,01×
n
]
t
为系统的输入矩阵，a和b具体为
[0106][0107]
其中i
n
×
n
表示一个n阶单位矩阵。
[0108]
利用a和b矩阵构建该系统的系统矩阵d为
[0109][0110]
利用构建的状态空间模型和系统矩阵，进行推导得到如下引理：：m i mo m i c系统存在n
f
个谐振频率它们是系统矩阵d特征值的虚部的绝对值，如下式所示。
[0111]
ω1,ω2,
…
,ω
n
＝|im(eigenvalue(d))|
ꢀꢀꢀꢀ
(11)
[0112]
作为本发明实施例的一种具体实施方式，在磁感应通信系统的总发射功率约束条件下，根据谐振频率对信道容量进行优化，得到最优信道容量，具体为：
[0113]
在磁感应通信系统的总发射功率约束条件下，根据谐振频率在子信道上分配的tx能量对信道容量的表达式进行优化，得到最优信道容量。
[0114]
示例性的，考虑mi mo磁感应通信系统信的道由n
t
个发射线圈和n
r
个接收线圈组成，本发明实施例假设通信信道已经做归一化处理的并且在这个频带内频率是平坦的，经过推导可知，mi mo磁感应通信系统的接收信号向量为：
[0115]
y＝hs n
[0116]
其中h是通信信道,n是高斯白噪声。
[0117]
经过推导可知，多频m i mo磁感应通信系统的输入输出关系表示为
[0118][0119]
经过推导可知，多频m i mo磁感应通信系统的信道容量表示为
[0120][0121]
其中，n
f
为分裂频率数，m
f
为在通过式(11)得到频率ω
f
的h(ω
f
)h(ω
f
)
h
正特征值的个数，λ
k
(ω
f
)是这些特征值之一，为信号向量在ω
f
频率处的协方差，n0为通带噪声谱功率。
[0122]
在本发明实施例中，对于磁感应通信系统，多线圈发射机的主要目标是设计s在总发射功率约束下达到最大的信道容量，即
[0123][0124][0125]
s
h
s≤p
c
[0126]
其中p
c
为最大允许的总传输功率，为s的协方差.为向量s的元素。
[0127]
可选地，由于发射端可以接入mimo磁感应通信系统的子信道，在总发射功率约束下，能够跨越子信道分配可变能量，使信道容量最大化。因此公式(12)的优化如下：
[0128][0129]
其中表示表示在频率为ω
f
的第k个子信道上分配的tx能量，满足且最优能量分配方案满足
[0130][0131]
其中η
*
是一个常数，根据特征值的集合的迭代计算和子信道分配能量的非负性获得，x

是指
[0132][0133]
作为本发明实施例的一种具体实施方式，根据谐振频率在子信道上分配的tx能量对信道容量进行优化，得到最优信道容量，具体为：
[0134]
根据通信信号将所有特征值按照大小降序排序，得到特征值排序集合；
[0135]
通过对特征值排序结合进行迭代计算和子信道分配功率的非负性计算得到优化常数；
[0136]
根据优化常数对子信道进行功率分配，当分配的子信道功率为负时将子信道丢弃，并重新运行迭代计算器加一；当分配给磁电感通信系统各个子型号的功率为负时，得到最优信道容量。
[0137]
在本发明实施例中，对自适应多频信道容量优化方案进行迭代设计，，在已知h的前提下，将所有的特征值λ
k
(ω
f
)按大小降序排列，从而形成一个的集合，并利用γ
l
来表示相关的λ
l
的功率分配。设置计数器p等于1，在每次迭代时，根据受约束的功率计算常数1/η
*
[0138][0139]
每个模式在每次迭代中分配的功率按下式计算
[0140][0141]
若分配给增益最小信道的能量为负，如γ
i
＜0，通过设置γ
i
＝0丢弃这个信道，并重新运行迭代计数器p加1的方案。当分配给各磁感应通信系统子信道的功率为非负时，为最优的资源分配方案，即得到最优信道容量。
[0142]
本发明实施例还提供了一种信道估计方法，具体为：
[0143]
当tx和rx线圈固定时，通道几乎是固定的。在这种情况下，rx线圈的电压通过互感与tx线圈相互作用，信道估计主要是根据rx侧的电压和负载计算互感，得到磁信道信息，接
收端adc模块对接收到的信号矢量进行采样。
[0144]
在磁感应通信系统中，信道由接收机利用发射机发出的训练信号进行估计。接收机知道训练信号序列定义为估计的信道，则接收信号为多线圈发射机的主要目标是在功率受限时控制使得误差达到最小值，即
[0145][0146][0147][0148]
其中为为第t
i
个tx线圈上的信号，本发明实施例将训练的tx信号矢量限制在总tx功率p
c
范围内。
[0149]
首先，将tx训练信号序列改写为对角矩阵形式，即然后对进行奇异值分解：
[0150][0151]
其中∑＝diag(σ1，σ2，
…
，σ
z
，
…
，σ
r
)，且有σ
z
≥σ
z 1
≥0，z＝1，2，
…
，r，r为的秩，σ
z
为矩阵第z个奇异值，e和f为酉矩阵。
[0152]
因为e和f为酉矩阵，有
[0153][0154]
其中e
z
为e的列向量，f
z
为f的列向量。
[0155]
令其中矩阵是的moorepenrose广义逆。其中的moorepenrose广义逆。其中
[0156]
根据奇异值的性质，得到
[0157][0158]
由
[0159][0160]
联立(31)和(32)可得
[0161][0162]
本发明实施例只需要通过有限额电压反馈信息即可进行信道估计，基于此能够动态获取共振频率以进行资源分配，不仅能够有效降低信道容量优化的复杂度，还能够保证信道容量优化的可靠性。
[0163]
实施本发明实施例，具有以下有益效果：
[0164]
本发明实施例采用时频域分析模型方法，基于磁感应通信系统的等效电路模型计算通信信道和谐振频率，从而能够准确确定由于多线圈分频而产生的多共振频率，使得系统能够选择多个频率作为子载波信道，能够在不额外增加电路及功耗的情况下，有效扩展带宽，从而有效地提高了信道容量。
[0165]
本发明实施例利用谐振频率作为子载波，推到出磁感应系统的信道容量，并且在给定的功率限制下，根据谐振频率对信道容量进行优化，得到最优信道容量，且本发明实施例能够使根据多个子载波来设计波形，使得系统的信道容量得到有效提高。
[0166]
请参阅图5，本发明的第二实施例提供了一种基于多谐振频率的信道容量优化装置，包括：
[0167]
第一计算模块10，用于计算磁感应通信系统中的发射端的阻抗矩阵和接收端的阻抗矩阵，根据发射端的阻抗矩阵和接收端的阻抗矩阵计算得到互感矩阵；
[0168]
电路获取模块20，用于获取磁感应通信系统对应的等效电路模型，等效电路模型包括第一等效电路模型和第二等效电路模型；
[0169]
第二计算模块30，用于根据基尔霍夫电压定律、互感矩阵以及第一等效电路模型计算磁感应通信系统的通信信道；
[0170]
第三计算模块40，用于根据基尔霍夫电压定律以及第二等效电路模型构建磁感应通信系统的状态空间模型，并根据状态空间模型计算磁感应通信系统的谐振频率；
[0171]
第四计算模块50，用于根据通信信道计算磁感应通信系统的接收信号向量，根据接收信号向量的协方差以及谐振频率，计算得到磁感应通信系统的信道容量；
[0172]
信道容量优化模块60，用于在磁感应通信系统的总发射功率约束条件下，根据谐振频率对信道容量进行优化，得到最优信道容量。
[0173]
可选地，本发明实施例可在如图2所示的磁感应通信系统中执行。请参阅图2，磁感应通信系统包括tx端、rx端、数字信号处理器(dsp)、数字模拟转换器(dac)、放大器(af)。在tx端，使用数字信号处理器(dsp)产生发射的基带信号，数字模拟转换器(dac)将离散的数字基带信号转换为连续的可变模拟基带信号，放大器(af)用于调制和放大这种模拟信号，每个发射线圈与设置的一个特定谐振频率的电容器相连。检测模块用于测量每个发射线圈上的电压幅值和相位，并反馈给dsp。在接收侧，每个接收线圈与为了获得相同谐振频率而设置的接收侧电容器相连。模数转换器(adc)对接收到的波形进行采样并发送给dsp。dsp用于对这些数据进行解码，以便进行通信和信号处理。
[0174]
请参阅图3
‑
4，为上述磁感应通信系统的两种等效电路模型。
[0175]
请参阅表1，为本发明实施例提供的等效电路模型各参数定义。
[0176]
表1等效电路模型各参数定义
[0177][0178][0179]
每个发射线圈的等效阻抗表示为
[0180][0181]
发射端的阻抗矩阵表示为
[0182][0183]
每个接收线圈的等效阻抗为
[0184][0185]
得到接收端的阻抗矩阵
[0186][0187]
互感矩阵定义为zm，表示rx线圈对tx线圈的影响。
[0188][0189]
作为本发明实施例的一种具体实施方式，第二计算模块30，具体用于：
[0190]
根据基尔霍夫电压定律计算第一等效电路模型的第一电路方程，根据第一电路方程计算得到电流矩阵；根据电流矩阵和互感矩阵计算得到接收电压，以接收电压定义接收信号和发送信号，根据接收电压计算接收信号与发射信号之间的系数作为通信信道。
[0191]
示例性的，根据基尔霍夫电压定律(kvl)方程，得到图3所示mimo磁感应通信系统的电路方程为
[0192][0193]
计算mimo磁感应通信系统系统中所有电路的电流，将阻抗矩阵进行反求，得到电流矩阵：
[0194][0195][0196][0197]
根据电流矩阵和护肝矩阵计算接收电压：
[0198][0199]
本发明实施例将电压类比为信号，则接收和发射的信号分别为y和s，即y和s之间的系数为通信信道h，定义为
[0200]
通信信道高度依赖于线圈间的磁通密度，从电压角度看，mimo磁感应通信系统类似于一般的无线通信系统，无线信号传播沿一定数量的衰落信道到达接收机，不同的传播机制对路径损耗和衰落模型有不同的影响。本发明实施例定义作为h的元素，其中是衰落因子。
[0201]
作为本发明实施例的一种具体实施方式，第三计算模块40，具体用于：
[0202]
根据基尔霍夫电压定律计算第二等效电路模型的第二电路方程，根据第二电路方程计算磁感应通信系统的电容矩阵、电阻矩阵和电感矩阵，根据电容矩阵、电阻矩阵和电感矩阵构建状态空间模型，利用空间模型推导计算得到磁感应通信系统的谐振频率。
[0203]
示例性的，本发明实施例根据kvl，图4所示等效电路模型的电路方程为
[0204][0205]
其中sn为
[0206][0207]
将上述的电路方程改写为
[0208][0209]
其中矩阵c为mimo磁感应系统的电容矩阵，如下式所示：
[0210][0211]
矩阵r为mimo磁感应通信系统的电阻矩阵，如下式所示
[0212][0213]
矩阵l为mimo磁感应通信系统的电感矩阵，如下式所示
[0214][0215]
引入方程
[0216][0217]
构建系统的状态空间模型
[0218][0219]
其中是状态矩阵，s
′
＝[s
t
，01×
n
]
t
为系统的输入矩阵，a和b具体为
[0220][0221]
其中i
n
×
n
表示一个n阶单位矩阵。
[0222]
利用a和b矩阵构建该系统的系统矩阵d为
[0223][0224]
利用构建的状态空间模型和系统矩阵，进行推导得到如下引理：：mim0mic系统存在nf个谐振频率它们是系统矩阵d特征值的虚部的绝对值，如下式所示。
[0225]
ω1，ω2，
…
，ω
n
＝|im(eigenvalue(d))|
ꢀꢀꢀꢀ
(11)。
[0226]
作为本发明实施例的一种具体实施方式，信道容量优化模块60，具体用于：
[0227]
在磁感应通信系统的总发射功率约束条件下，根据谐振频率在子信道上分配的tx能量对信道容量的表达式进行优化，得到最优信道容量。
[0228]
示例性的，考虑mimo磁感应通信系统信的道由nt个发射线圈和nr个接收线圈组成，本发明实施例假设通信信道已经做归一化处理的并且在这个频带内频率是平坦的，经过推导可知，mimo磁感应通信系统的接收信号向量为：
[0229]
y＝hs n
[0230]
其中h是通信信道，n是高斯白噪声。
[0231]
经过推导可知，多频mimo磁感应通信系统的输入输出关系表示为
[0232][0233]
经过推导可知，多频mimo磁感应通信系统的信道容量表示为
[0234][0235]
其中，n
f
为分裂频率数，m
f
为在通过式(11)得到频率ω
f
的h(ω
f
)h(ω
f
)
h
正特征值的个数，λ
k
(ω
f
)是这些特征值之一，为信号向量在ω
f
频率处的协方差，n0为通
带噪声谱功率。
[0236]
在本发明实施例中，对于磁感应通信系统，多线圈发射机的主要目标是设计s在总发射功率约束下达到最大的信道容量，即
[0237][0238][0239]
s
h
s≤p
c
[0240]
其中pc为最大允许的总传输功率，为s的协方差，为向量s的元素。
[0241]
可选地，由于发射端可以接入mimo磁感应通信系统的子信道，在总发射功率约束下，能够跨越子信道分配可变能量，使信道容量最大化。因此公式(12)的优化如下：
[0242][0243]
其中表示表示在频率为ω
f
的第k个子信道上分配的tx能量，满足且最优能量分配方案满足
[0244][0245]
其中η
*
是一个常数，根据特征值的集合的迭代计算和子信道分配能量的非负性获得，x

是指
[0246][0247]
作为本发明实施例的一种具体实施方式，根据谐振频率在子信道上分配的tx能量对信道容量进行优化，得到最优信道容量，具体为：
[0248]
根据通信信号将所有特征值按照大小降序排序，得到特征值排序集合；
[0249]
通过对特征值排序结合进行迭代计算和子信道分配功率的非负性计算得到优化常数；
[0250]
根据优化常数对子信道进行功率分配，当分配的子信道功率为负时将子信道丢弃，并重新运行迭代计算器加一；当分配给磁电感通信系统各个子型号的功率为负时，得到最优信道容量。
[0251]
在本发明实施例中，对自适应多频信道容量优化方案进行迭代设计，，在已知h的前提下，将所有的特征值λ
k
(ω
f
)按大小降序排列，从而形成一个的集合，并利用γ
l
来表示相关的λ
l
的功率分配。设置计数器p等于1，在每次迭代时，根据受约束的功率计算常数1/η
*
[0252][0253]
每个模式在每次迭代中分配的功率按下式计算
[0254][0255]
若分配给增益最小信道的能量为负，如γ
i
＜0，通过设置γ
i
＝0丢弃这个信道，并重新运行迭代计数器p加1的方案。当分配给各磁感应通信系统子信道的功率为非负时，为最优的资源分配方案，即得到最优信道容量。
[0256]
本发明实施例还提供了一种信道估计方法，具体为：
[0257]
当tx和rx线圈固定时，通道几乎是固定的。在这种情况下，rx线圈的电压通过互感与tx线圈相互作用，信道估计主要是根据rx侧的电压和负载计算互感，得到磁信道信息，接收端adc模块对接收到的信号矢量进行采样。
[0258]
在磁感应通信系统中，信道由接收机利用发射机发出的训练信号进行估计。接收机知道训练信号序列定义为估计的信道，则接收信号为多线圈发射机的主要目标是在功率受限时控制使得误差达到最小值，即
[0259][0260][0261][0262]
其中为为第t
i
个tx线圈上的信号，本发明实施例将训练的tx信号矢量限制在总tx功率p
c
范围内。
[0263]
首先，将tx训练信号序列改写为对角矩阵形式，即然后对进行奇异值分解：
[0264][0265]
其中∑＝diag(σ1，σ2，
…
，σ
z
，
…
，σ
r
)，且有σ
z
≥σ
z 1
≥0，z＝1，2，
…
，r，r为的秩，σ
z
为矩阵第z个奇异值，e和f为酉矩阵。
[0266]
因为e和f为酉矩阵，有
[0267][0268]
其中e
z
为e的列向量，f
z
为f的列向量。
[0269]
令其中矩阵是的moorepenrose广义逆。其中的moorepenrose广义逆。其中
[0270]
根据奇异值的性质，得到
[0271][0272]
由
[0273][0274]
联立(31)和(32)可得
[0275][0276]
本发明实施例只需要通过有限额电压反馈信息即可进行信道估计，基于此能够动态获取共振频率以进行资源分配，不仅能够有效降低信道容量优化的复杂度，还能够保证信道容量优化的可靠性。
[0277]
实施本发明实施例，具有以下有益效果：
[0278]
本发明实施例采用时频域分析模型方法，基于磁感应通信系统的等效电路模型计算通信信道和谐振频率，从而能够准确确定由于多线圈分频而产生的多共振频率，使得系统能够选择多个频率作为子载波信道，能够在不额外增加电路及功耗的情况下，有效扩展带宽，从而有效地提高了信道容量。
[0279]
本发明实施例利用谐振频率作为子载波，推到出磁感应系统的信道容量，并且在给定的功率限制下，根据谐振频率对信道容量进行优化，得到最优信道容量，且本发明实施例能够使根据多个子载波来设计波形，使得系统的信道容量得到有效提高。
[0280]
以上是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。