2021年3期
科技创新与应用
Technology Innovation and Application
研究视界
丁同禹，王孟于
（集美大学信息工程学院，福建厦门361021）
引言
随着电子与远程通信技术的飞速发展，电子电气设备与系统的集成度不断提高，其在制程及使用过程中受工艺或外界电磁干扰影响而产生参数不确定性（parame - ter uncertainty ）的现象日益严重[1-3]，比如，设备、元件及线路在长期使用后发生老化；各种外界因素（温度，空气湿
度，静电放电等）对元件、线束造成性能改变或损坏等。一旦这些内部或外界的干扰因素造成了设备或元件电参数（如微带线的介电常数）的变化，将会直接影响整个系统的性能，造成电路或系统可靠性的下降和安全系数的大幅降低。
不仅如此，不确定参数问题还会对电磁信号传播特性的研究造成危害，尤其是室内环境的无线信道建模。在利用射线追踪（ray-tracing ）、时域有限差分（FDTD ）等方法建立指定室内场景无线信道模型时，目标环境中的不确定参数（如空气湿度导致的建筑物墙体电导率及介电常数变化，施工过程产生的建筑物几何尺寸误差等），会在信道模型中引入大量的随机变量，导致该模型的输出（如路径损耗等）变成不确定的范围区间，从而降低模型可信度，严重影响其在无线系统布局及优化中的应用价值。
因此，当电路、元件或电磁环境中某一媒质的电磁参数发生变化时，亟需对系统响应的随机特性进行预测和评估，有效地降低不确定参数对电路与系统可靠性造成的危害；在信道建模研究等领域则可获取更为完整的电波传播特性数据和信息，提高无线信道模型的可信度。
1最坏分析算法
最坏分析算法作为一种高效的数值仿真技术能够准确分析不确定参数系统的随机特性并预测其响应（输出）的变化范围，不仅可为电路系统的可靠性评估提供依据，还能缩短研发周期，提高设备与系
统的安全系数[4]。用于不确定参数系统分析的传统数值技术大都以概率类算法为基础，如以Monte Carlo 算法为代表的基于抽样理论的系列方法，通过大量的随机试验来覆盖不确定参数可能出现的所有取值，从而不断接近所要预测的真实边界，其缺点是当不确定参数变化范围较大时，所需的样本数量异常巨大，导致效率急剧下降（Monte Carlo 算法的收敛速度为1/M √，M 为样本数量），即便如此，其预测结果仍然可能出现遗漏真实值的情况。
相比之下，基于最坏分析算法的数值技术以绝对保守的方式对系统输出进行预测，其所提供的预测结果能
摘要：电子线路与系统的不确定参数问题对系统可靠性与稳定性造成严重危害。当电路、元件，甚至电磁环境中某一媒质的电磁参数发生变化时，亟需对系统响应的随机特性进行预测和评估，从而降低不确定参数可能导致的安全性问题。文章提出基于仿射算法的最坏分析架构，并模拟了集总参数电路在受到不确定因素影响时，其系统响应的变化范围。相比于传统的随机试验方法，本模型具有更高的准确度和计算效率。
关键词：仿射算法；不确定参数；可靠性中图分类号:O231
文献标志码:A
文章编号:2095-2945(2021)03-0076-03
Abstract :The uncertain parameters of electronic circuits and systems are harmful to the reliability and stability of the
system.When the electromagnetic parameters of circuits,components and even a medium in electromagnetic environment change,it is necessary to predict and evaluate the random characteristics of the system response,so as to reduce the security problems caused by uncertain parameters.In this paper,a worst-case analysis framework based on affine algorithm is proposed,and the variation range of the system response of lumped parameter circuit under the influence of uncertain factors is simulated.Compared with the traditional random test methods,this model achieves higher accuracy and computational
efficiency.
Keywords :affine algorithm;uncertain parameter;reliability
*基金项目：福建省自然科学基金青年基金(编号:2018J05109)
作者简介：丁同禹（1987-），男，博士，副教授，研究方向：电路与系统可靠性分析；王孟于（1996-），男，硕士研究生，研究方向：
电路与系统可靠性分析。
76--
2021年3期
科技创新与应用
Technology Innovation and Application
研究视界
够覆盖所有的真实情况，即“始终包含全部可能出现的真实结果”，不会产生任何遗漏，从而提供比经典的概率类算法更加可靠的预测结果。同时，最坏分析类算法的收敛速度在变量数相对较少时远大于概率类算法，在效率上也具有明显的优势。
音乐枕头2仿射算法的基本原理
作为最坏分析算法的典型代表，仿射算法
（affine arithmetic ，AA ）的概念在20世纪90年代由Comba 和
Stolfi 首次提出[5]。仿射算法能够实现对预测对象浮动范围的保守估计，即直接预测电路与系统可能出现的最坏情况。该算法从区间算法（interval analysis ）变形而来，因
此也能够自动记录浮点的截尾和舍入误差。此外，仿射算法增加了一些附加的噪声系数，若多个仿射形式有相同的噪声系数，那么证明两者具有关联性，仿射形式之间关联性和依赖性可以根据变量之间相同噪声系数的多少决定。这种明确变量之间关联性的方法，有助于辨别计算过程中所记录的各个不确定参数的依赖关系，如果变量之间的噪声系数呈负相关时，可以通过计算变量相互抵消达到减小区间扩张值的目的，得到更小的区间结果，提高计算结果的准确性，同时缩短计算时间。
在区间算法中，变量是使用两个临界点的极值来表示，例如，x ⎺=[a ，b]代表一个变量x ，其取值范围落在区间
[a ，b]之内。该方法的主要缺陷是无法捕捉到各个变量之
间的关系，从而导致对真实区间边界的“高估”。仿射算法
对不确定参数的表达形式进行了修订，同样对于一个区间变量x ∈[a ，b]，其对应的AA 变量表达为：x ~=x 0+x 1ϵ，其中x 0=（a+b ）/2，而x 1=（b-a ）/2称为其“半径”，噪声系数ϵ
是一个取值范围始终在±1之间的独立区间变量。可见，
相比于传统的区间算法，仿射算法在运算过程中考虑了变量之间的关联，并将其引入到计算链中，从而有效消除了变量裕量叠加所引起的过度估计，对于任意一个给定的含有n 个不确定因素的随机变量，其仿射表达形式为：
x ~=x 0+∑i=1n
x i ϵi
（1）
在此定义下给出的变量x ~的最坏边界为（上下边界）：
x 0±∑i |x i |。若已知两个独立的区间变量：
x ~和y ~，其仿射形式的乘法运算如公式（2）所示。（2）
显然地，在此过程中产生了一个新的高阶项。
为了消除该项所引入的计算复杂度，引入一个新的变量：ζ。令：
则新产生的z ~仍然满足AA 形式，同时引入了一项新
的变量ζ，其半径幅值为R （x ~y ~）。在使用仿射算法进行函数界限的计算时，不同形式的运算所引起的运算量、运算精度以及二者之间的关系具有很大的差异，清楚在什么情况下使用什么形式的区间方法可以获得何种精度的结果是十分重要的，因此本研究建立了包括乘法运算、指数、对数、倒数以及除法运算等在内的各种基本函数的仿射运算模型，形成了仿射函数库，并在此基础上实现了参数矩阵的指数和求逆运算，其实现方法在第3节简述。
3基于仿射算法的矩阵运算
3.1参数矩阵求逆
矩阵求逆是电路频域求解过程中十分关键的一项内容。对于仿射形式的变量矩阵，其求逆运算需要结合其他数值算法，为了叙述简便，本文只讨论一个不确定因素的情况，例如求解：
（4）
其中，X 0是中心矩阵，X 1代表不确定干扰因素ϵ1的
“半径”矩阵。当X 1的秩为1时，使用Sherman Morrison 公式对其求逆，过程如式（5）所示。
（5）其中，g=trace （X 0-1X 1）为矩阵的迹。以上运算将求得
无线抄表的Y ~矩阵转换成标准的仿射形式。当X ~矩阵的秩不为1，或其同时受到多个不确定因素干扰时，例如：
（6）
则需要将其拆分成若干个秩为1的子矩阵的加和形式，然后逐项使用公式（5）完成求逆运算。
3.2参数矩阵的指数运算
集总和分布式电路求解过程中经常使用到的另一项运算就是求矩阵指数。本文使用了bilinear 近似公式作为辅助技术，完成仿射形式的矩阵指数运算，如公式（7）所示：
（7）
其中，I 为单位矩阵，Y 0=exp （X 0）为新的中心矩阵，
ΔX=X 1ϵ1+…+X n ϵn 为原矩阵的全部干扰项。进一步化简
n n n
n
i i i i i i
i 1i 1i 1
i 1
x y x y R( xy)
%%òò（3）
°° 11011
Y X X X ò° 101111100110011
Y X X 1X (X X X )
1g
Y Y òòò°° 11011n n
Y X X X X L òò°0
0011n n 01X X 0Y exp(X X X )Y Y 11
e
I (X X)X(e I)
22
L òòn n 0i i 0i i i 1i 1
n
矩阵干扰
n n 000i 0i i i i i i i 1
i 1i 1x y x x y y x y (x y y x )x y %
%òòòòòz %77--
2021年3期
科技创新与应用
Technology Innovation and Application
研究视界
ΔY 可得：
（8）
4不确定参数电路的时域分析
本节应用仿射算法完成一个元件参数为不确定量的
RLC 电路的时域仿真和分析，电路图如图1所示。其中，L=1mH ，R L =0.5Ω，C=10μF ，R C =0.5Ω，R=20Ω。设电容C 存在±50%的容差，应用仿射算法预测输出（电感电流i L （t
防爆鼠标））的响应范围。
图1含有不确定参数元件的RLC 电路
首先对该电路的状态方程：x ̇（t ）=Ax （t ）+bu （t ）（其中
状态向量x （t ）=[i L （t ），v C （t ）]T 为电感电流和电容电压，且x （0）=[0，0]T ）进行时间离散：
x （k+1）=A d x （k ）+B d u （k ）（9）其中，A d =e A ·Δt ，B d =A -1（e A ·Δt -I ）b
；
[A]=
-（[R L ]+[R]//[R C ]）/[L]
-[R]//[R C ]（[L][R C ]）
蚝排[R]//[R C ]/（[C][R C ]）-[R]//[R C ]/（[C][R][R C ]）]
[
]
[b]=[1/[L]，0]T
以上公式中的a//b 表示并联：a ·b/（a+b ）。
应用上文所述的仿射算法数值结构，依次完成本电路求解运算中所需要的对参数矩阵A 的求逆和指数运算，并解得仿射形式的电感电流i L （t ）在时域的输出波形如图2所示。
如图2所示，本文所使用的基于仿射算法的数值架构能够有效预测在RLC 电路的电容值出现变化时，电路的时域输出波形可能产生的变化，并成功估计了电感电流在时域可能出现的最大值和最小值两种极端情况。为了验证算法的有效性，与经典的Monte Carlo 随机试验方法（图2灰实线）进行了对比验证，结果表明仿射算法（图2黑实线）的预测范围准确，同时在Matlab 环境下运行时仿真速度达到Monte Carlo 算法
（10000次随机试验）的3倍以上。
5结束语
以电路与系统的可靠性预测为目标，本文对一种典型的RLC 电路进行了基于仿射算法的最坏分析。该研究通过将仿射算法与其他数值技术，如Sherman Morrison 公式，bilinear 数值近似等技术结合，解决了仿射形式矩
阵的运算问题，进而求解了目标电路的电容元件参数发
生变化时，系统响应的变化范围。该预测结果与经典的
Monte Carlo 方法进行了对比，不仅达到了较高的精度，且在运算效率上得到了显著提高。
参考文献院
[1]谢永强，陈建军.不确定系统稳定性的仿射不等式分析[J].高技术通讯，2009，19（8）：101-105.
[2]谢永强，陈建军，朱增青.不确定系统响应上下界分析的改进仿射算法[J].电子科技大学学报，2011，40（4）：634-640.[3]乐健，廖小兵，李奔，等.基于扩展仿射模型的不确定性静态电压稳定性全局灵敏度分析[J].电工技术学报，2020-10-30.[4]T.Ding ，R.Trinchero ，P.Manfredi ，I.S.Stievano ，F.G.Canavero.How affine arithmetic helps beat uncertainties in electrical systems [J].IEEE Circuits Syst.Mag.，2015，
15（4）：70-79.
[5]STOLFI J.Stolfi ，L.H.Figueiredo.Self-validated reduced numerical methods and applications [C].Monograph Brazilian Mathematics Colloquium ，IMPA ，Rio De Janeiro ，Brazil ，1997.
1
0001122n n
1
j 00j 01Y (Y I)X X(Y I)
2
Y Y Y 1Y (Y I)X X (Y I),j 1,2,,n.
2
L òò
ò图2电感电流i L （t
）的预测结果78--1.08b
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议