CN113343453B - 一种基于小步长离散化的电力电子级联变换器的建模方法 - Google Patents

一种基于小步长离散化的电力电子级联变换器的建模方法 Download PDF

Info

Publication number
CN113343453B
CN113343453B CN202110595176.6A CN202110595176A CN113343453B CN 113343453 B CN113343453 B CN 113343453B CN 202110595176 A CN202110595176 A CN 202110595176A CN 113343453 B CN113343453 B CN 113343453B
Authority
CN
China
Prior art keywords
small
iteration
power electronic
state
modeling method
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
CN202110595176.6A
Other languages
English (en)
Other versions
CN113343453A (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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN202110595176.6A priority Critical patent/CN113343453B/zh
Publication of CN113343453A publication Critical patent/CN113343453A/zh
Application granted granted Critical
Publication of CN113343453B publication Critical patent/CN113343453B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M3/00—Conversion of dc power input into dc power output
    • H02M3/02—Conversion of dc power input into dc power output without intermediate conversion into ac
    • H02M3/04—Conversion of dc power input into dc power output without intermediate conversion into ac by static converters
    • H02M3/10—Conversion of dc power input into dc power output without intermediate conversion into ac by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
    • H02M3/145—Conversion of dc power input into dc power output without intermediate conversion into ac by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
    • H02M3/155—Conversion of dc power input into dc power output without intermediate conversion into ac by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Power Engineering (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Dc-Dc Converters (AREA)

Abstract

本发明公开了一种基于小步长离散化的电力电子级联变换器的建模方法，包括以下步骤：S1、根据电力电子变换器级联系统的具体拓扑及控制方式，获得状态方程及各状态检测向量的表达式；S2、选取合适的步长，在每个小步长区间中，使用欧拉法解出各状态方程，将各状态方程变为线性方程；S3、设置标志位、迭代初值进行迭代运算，根据每次迭代标志位及状态检测向量的值，选取合适工作状态方程计算各状态变量的值；S4、输出满足条件的迭代结果。本发明考虑了级联系统中所有变换器完整的动态过程，所得到的数值仿真结果将更加贴合实际情况。

Description

一种基于小步长离散化的电力电子级联变换器的建模方法
技术领域
本发明涉及电力电子变换器的技术领域，尤其涉及一种基于小步长离散化的电力电子级联变换器的建模方法。
背景技术
传统的火力发电造成了大量污染，已经不适合当今时代的发展。为了解决能源问题，各国都在积极发展分布式发电系统。而在分布式电源中，级联的电力电子变换器是必不可少的。要对级联的电力电子变换器进行详尽的分析，其模型是必不可少的。
针对级联的电力电子变换器系统，传统的离散映射方法过于复杂，实际可操作性不高；状态空间平均模型只能针对同频率的级联系统使用，且无法分析系统的快时标失稳现象；降阶的恒功率源法在级联系统每级的频率相差不大时，准确度不高。因此，必须提出一种模型能够完整的考虑级联系统中的每一级变换器且不进行平均化或降阶处理，得到的结果才会接近实际。
Fan Xie,Bo Zhang,Dongyuan Qiu等人对与电力电子级联变换器系统的建模已有过研究，他们利用降阶思想，使用了恒功率源法对级联系统进行了建模，但这种方法没有考虑级联系统中每一级变换器完整的动态行为，无法对级联系统进行准确的分析。而本发明所提出的小步长离散化建模法克服了Fan Xie等人所使用的恒功率源法的缺陷，完整的考虑了电力电子级联变换器系统中的每一级变换器，没有进行平均化或降阶处理，可以准确捕捉整个级联系统的动态行为及稳定性。(Fan Xie,Bo Zhang,Dongyuan Qiu and YanweiJiang.“Non-linear dynamic behaviours of DC cascaded converters system withmulti-load converters.”IET Power Electronic,Vol.9,pp.1093-1102,2016.)
电力电子级联变换器系统建模的难点在于实际系统的工作状态太多，且其状态排布没有规律，若不进行平均化或降阶处理，很难对其工作规律进行捕捉。随着计算机计算能力的提高和数值算法的完善，小步长离散化的理论在建模问题中一直有广泛的应用。
发明内容
本发明的目的在于克服上述现有模型的不足，使求得的模型更加精确、更加接近实际情况，提出了一种基于小步长离散化的电力电子级联变换器的建模方法，其完整的考虑了级联系统中的每一级变换器，基本没有做近似处理，模型数值仿真的结果将会更加精确。
本发明至少通过如下技术方案之一实现。
一种基于小步长离散化的电力电子级联变换器的建模方法，包括以下步骤：
S1、根据电力电子变换器级联系统的具体拓扑及控制方式，获得状态方程及各状态检测向量的表达式；
S2、设置步长，在每个小步长区间中，使用欧拉法解出各状态方程，将各状态方程变为线性方程；
S3、设置标志位、迭代初值进行迭代运算，根据每次迭代标志位及状态检测向量的值选取工作状态方程，计算各状态变量的值；
S4、输出满足条件的迭代结果。
优选的，所述状态方程为：
Figure GDA0003881546040000031
其中X为电力电子级联变换器系统的状态变量；Ai、Bi均为常量矩阵，与级联变换器的具体拓扑有关，i为级联变换器系统的工作模态个数；
Figure GDA0003881546040000032
Vin为输入电压。
优选的，将各状态方程变为线性方程分为以下两种情况：
第一种情况，在一个小步长内系统的工作状态没有发生改变
Figure GDA0003881546040000033
其中，Xa是离散区间起始时状态变量的值；Xc是该离散区间结束时状态变量的值；i为系统当前的工作状态；定义Hi,i(Xa)为第一种情况的离散迭代方程；Te为小步长的长度；
第二种情况，在一个步长内系统的工作状态由i变为j：
Figure GDA0003881546040000034
其中d为虚拟占空比，根据峰值电流控制的原理获得；定义Hi,j(Xa,d)为第二种情况的离散迭代方程。
优选的，标志位设置为：
P＝[b1 b0 a1 a0]
其中b1、b0代表一个小步长内前部分系统的工作状态开关器件通断情况，为0则通，为1则断；a1、a0则代表一个小步长内后部分系统的开关的通断情况，设置两个十进制数p1、p2以指示系统开关的通断的情况，其值根据标志位P的值设置：
Figure GDA0003881546040000041
优选的，所述状态检测向量的设置如下：
Q3＝[q3,0 q3,1 q3,2 q3,3]
其中
Figure GDA0003881546040000042
q3,1＝Te-mod((k1+1)Te,T1)；q3,2＝Te-mod((k1+1)Te,T2)；q3,3＝-1；qi,j(i,j∈{0,1,2,3})为状态检测向量Qi中的各元素，其下标i表示电力电子级联变换器系统此时所处的工作模态，j表示若系统满足qi,j所对应的切换条件所跳转至的工作模态；T1为前级周期；T2为后级周期；Te为小步长的长度；k1为常数。
优选的，在迭代中，p1＝p2，则该次迭代完成后，须检查系统所处的工作模态所对应的检测向量Qi中的各个元素，如果Qi中的所有元素均小于0，则迭代继续，如果Qi中有元数大于0，则系统在该次迭代中应进行工作模态的切换，由模态i切换到模态s，令p2＝s重新进行迭代，如果在某次迭代中，p1≠p2，则该次迭代完成后，令p1＝p2，再进行接下来的迭代。
优选的，在重新进行的迭代中，根据开关切换条件gi,s获取出虚拟占空比d。
优选的，所述各状态变量的值保存在记忆体Xr中。
优选的，所述小步长基于前级周期T1、后级周期T2的最大公约数进行设置。
优选的，所述小步长取≤5e-6。
与现有技术相比，本发明具有如下优点：
1、通过小步长离散化的方法，没有进行近似的对电力电子级联变换器系统进行了建模，能够精确的描述系统的工作情况。
2、在每个小步长区间中，利用欧拉法将系统的状态方程写成了线性方程，这大大减少了计算量。
3、在基于小步长离散化的级联变换器模型的基础上，还可对系统进行稳定性分析，这在实际工程设计中具有重要指导意义。
附图说明
图1为本发明实施例一种基于小步长离散化的电力电子级联变换器的建模方法的流程图；
图2为本发明实施例基于小步长离散化的不同频率级联Boost变换器建模的等效电路图；
图3为实施例一种基于小步长离散化的电力电子级联变换器的建模方法原理图；
图4为实施例中由该建模方法得到各状态变量的时域波形与Simulink仿真得到的时域波形的对比图；
图5为实施例中前级电感电流随前级参考电流变化的分岔图。
具体实施方式
本部分将详细描述本发明的具体实施例，本发明之较佳实施例在附图中示出，附图的作用在于用图形补充说明书文字部分的描述，使人能够直观地、形象地理解本发明的每个技术特征和整体技术方案，但其不能理解为对本发明保护范围的限制。
如图1所示，本实施例一种基于小步长离散化的电力电子级联变换器的建模方法，包括以下步骤：
首先，根据电力电子变换器级联系统的具体拓扑及控制方式，得到电力电子变换器级联系统的状态方程及各状态检测向量的表达式；
然后选取合适的步长(一般可取≤5e-6)，在每个小步长区间中，使用欧拉法解出各状态方程，将各状态方程变为简单的线性方程。
接着设置迭代次数u和迭代初值进行迭代运算，根据每次迭代标志位及状态检测向量的值选择表达式计算各状态变量的值X0、X1、Xn，并将结果保存在记忆体Xr中。
如果在某次迭代中，p1＝p2，则该次迭代完成后，须检查此时电力电子变换器系统所处的工作模态所对应的检测向量Qi中的各个元素，如果Qi中的所有元素均小于0，则迭代继续，如果Qi中有元数大于0，比如qi,s，则系统在该次迭代中应进行工作模态的切换，由模态i切换到模态s，所以令p2＝s重新进行迭代，在重新进行的迭代中，可根据电力电子级联变换器系统的开关切换条件gi,s计算出虚拟占空比d。
如果在某次迭代中，p1≠p2，则该次迭代完成后，令p1＝p2，再进行接下来的迭代。
最后，当达到满足工程技术人员要求的迭代次数u后，输出满足条件的迭代结果。
下面以峰值电流控制不同频率级联的Boost变换器为例，对实施上述方法进行具体说明，峰值电流控制不同频率级联的Boost变换器的等效电路如图2所示。为了验证方便，电路结构中的所有器件均视为理想器件，电路共有四种状态，其状态方程为：
Figure GDA0003881546040000071
X＝[iL1 vC1 iL2 vC2]T为电力电子级联变换器系统的状态变量；Ai、Bi为常量矩阵，i为级联变换器系统的工作模态个数；
Figure GDA0003881546040000072
Figure GDA0003881546040000073
B3＝B2＝B1＝B0＝[1/L1 0 0 0]T；Vin为输入电压；L1为图2所对应的系统的前级电感，L2为后级电感，C1为前级电容，C2为后级电容，R为负载电阻，iL1为前级电感电流，iL2为后级电感电流，vC1为前级电容电压，vC2为后级电容电压。
根据欧拉法，在每一个离散的小步长中，可分为两种情况将其状态方程写成线性形式。
第一种情况，在一个小步长内系统的工作状态没有发生改变
Figure GDA0003881546040000074
其中，Xa是离散区间起始时状态变量的值；Xc是该离散区间结束时状态变量的值；i为系统当前的工作状态。定义Hi,i(Xa)为此种情况下的离散迭代方程。
第二种情况，在一个小步长内系统的工作状态由i变为j：
Figure GDA0003881546040000081
其中d为虚拟占空比，可根据峰值电流控制的原理算得。定义Hi,j(Xa,d)为此种情况下的离散迭代方程。
选择标志位：
P＝[b1 b0 a1 a0]
其中b1、b0代表一个小步长内前部分峰值电流控制Boost变换器级联系统的工作状态开关器件通断情况(即图2中的S1、S2)，为0则通，为1则断；a1、a0则代表一个小步长内后部分峰值电流控制Boost变换器级联系统的工作状态开关器件的通断情况，可设置两个十进制数p1、p2来方便状态的切换，其计算方式如下：
Figure GDA0003881546040000082
每次迭代都要根据p1、p2来选择合适的离散迭代方程
Figure GDA0003881546040000083
以来计算状态变量的值。
而标志位P的设置与状态检测向量Qi息息相关，以图3所示的级联系统先处于两级均断的模态3后跳转到两级均通的模态0为例说明如何计算状态检测向量Qi。假设变换器运行到时间t为nT1＝mT2(n、m均为常数)前的一个Te(第k1-1个Te，k1为常数)内，S1、S2均关断，系统工作在两级均断的模态3，此时标志位P＝[1111]、所以指示开关状态的两个十进制数p1＝3、p2＝3。因此离散迭代方程H3,3被选择来迭代计算此个小步长内状态变量的值。依次检测状态检测向量Q3中的元素是否都小0。所述依次是要先检测多个开关同时动作的情况，所以，Q3可按下式设置：
Q3＝[q3,0 q3,1 q3,2 q3,3]
其中
Figure GDA0003881546040000091
q3,1＝Te-mod((k1+1)Te,T1)；q3,2＝Te-mod((k1+1)Te,T2)；q3,3＝-1，qi,j(i,j∈{0,1,2,3})为Qi中的各元素，其下标i表示电力电子级联变换器系统此时所处的工作模态，j表示若系统满足qi,j所对应的切换条件所跳转至的工作模态；T1为前级周期；T2为后级周期；Te为小步长的长度。
例如当电力电子级联变换器系统处于工作模态3时，此时标志位P所对应的两个十进制数为p1＝p2＝3，所以对应的离散迭代方程H3,3先被选取来计算状态变量的值，根据H3,3进行迭代运算完成后，依次状态检测向量Q3中的值，此时q3,0大于0，则本次迭代舍去，并令p2＝0，重新进行迭代。
在图3中，当t＝k1Te时，因为p1＝p2＝3，所以H3,3先被选取以来计算状态变量的值，计算完成后，依次检测Q3中的值，此时q3,0大于0了，则本次迭代舍去，并令p2＝0(q3,0下标与这里的一致性，即p2的值应为q3,0下标后面的0)，重新进行迭代。其余情况与此类似。由上述流程建模并进行数值仿真得到的各状态变量的时域波形与Simulink仿真的结果对比如图4所示，可以看出，两者基本一致，这也验证了该模型的正确性。
另外，小步长可特殊地基于T1、T2的最大公约数来选取，这样可以获得一个周期表达式，通过该周期表达式可对系统进行非线性动力学稳定性分析，图5所示的分岔图由此得到。
上述实施案例为本发明较佳的实施方式，但本发明的实施方式并不受所述实施案例的限制，其它的任何未背离本发明的精神实质与原理下所做的改动、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (9)

1.一种基于小步长离散化的电力电子级联变换器的建模方法，其特征在于：包括以下步骤：
S1、根据电力电子变换器级联系统的具体拓扑及控制方式，获得状态方程及各状态检测向量的表达式；
S2、设置步长，在每个小步长区间中，使用欧拉法解出各状态方程，将各状态方程变为线性方程；
S3、设置标志位、迭代初值进行迭代运算，根据每次迭代标志位及状态检测向量的值选取工作状态方程，计算各状态变量的值；
将各状态方程变为线性方程分为以下两种情况：
第一种情况，在一个小步长内系统的工作状态没有发生改变
Figure FDA0003888825340000011
其中，Xa是离散区间起始时状态变量的值；Xc是该离散区间结束时状态变量的值；Ai、Bi均为常量矩阵，与级联变换器的具体拓扑有关，其下标i表示电力电子级联变换器系统此时所处的工作模态；Vin为输入电压；定义Hi,i(Xa)为第一种情况的离散迭代方程；Te为小步长的长度；
第二种情况，在一个步长内系统的工作模态由i变为j：
Figure FDA0003888825340000012
其中d为虚拟占空比，根据峰值电流控制的原理获得；定义Hi,j(Xa,d)为第二种情况的离散迭代方程；
S4、输出满足条件的迭代结果。
2.根据权利要求1所述的一种基于小步长离散化的电力电子级联变换器的建模方法，其特征在于，所述状态方程为：
Figure FDA0003888825340000021
其中X为电力电子级联变换器系统的状态变量；
Figure FDA0003888825340000022
3.根据权利要求2所述的一种基于小步长离散化的电力电子级联变换器的建模方法，其特征在于，标志位设置为：
P＝[b1 b0 a1 a0]
其中b1、b0代表一个小步长内前部分系统的工作状态开关器件通断情况，为0则通，为1则断；a1、a0则代表一个小步长内后部分系统的开关的通断情况，设置两个十进制数p1、p2以指示系统开关的通断的情况，其值根据标志位P的值设置：
Figure FDA0003888825340000023
4.根据权利要求3所述的一种基于小步长离散化的电力电子级联变换器的建模方法，其特征在于，所述状态检测向量的设置如下：
Q3＝[q3,0 q3,1 q3,2 q3,3]
其中
Figure FDA0003888825340000024
q3,1＝Te-mod((k1+1)Te,T1)；q3,2＝Te-mod((k1+1)Te,T2)；q3,3＝-1；qi,j为状态检测向量Qi中的各元素，i,j∈{0,1,2,3}，j表示若系统满足qi,j所对应的切换条件所跳转至的工作模态；T1为前级周期；T2为后级周期；Te为小步长的长度；k1为常数。
5.根据权利要求4所述的一种基于小步长离散化的电力电子级联变换器的建模方法，其特征在于，在迭代中，如果p1＝p2，则该次迭代完成后，须检查系统所处的工作模态所对应的检测向量Qi中的各个元素，如果Qi中的所有元素均小于0，则迭代继续，如果Qi中有元数大于0，则系统在该次迭代中应进行工作模态的切换，由模态i切换到模态s，令p2＝s重新进行迭代，如果在某次迭代中，p1≠p2，则该次迭代完成后，令p1＝p2，再进行接下来的迭代。
6.根据权利要求5所述的一种基于小步长离散化的电力电子级联变换器的建模方法，其特征在于，在重新进行的迭代中，根据开关切换条件gi,s获取出虚拟占空比d。
7.根据权利要求1～6任一项所述的一种基于小步长离散化的电力电子级联变换器的建模方法，其特征在于，所述各状态变量的值保存在记忆体Xr中。
8.根据权利要求7所述的一种基于小步长离散化的电力电子级联变换器的建模方法，其特征在于，所述小步长基于前级周期T1、后级周期T2的最大公约数进行设置。
9.根据权利要求8所述的一种基于小步长离散化的电力电子级联变换器的建模方法，其特征在于，所述小步长取≤5e-6。
CN202110595176.6A 2021-05-28 2021-05-28 一种基于小步长离散化的电力电子级联变换器的建模方法 Active CN113343453B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110595176.6A CN113343453B (zh) 2021-05-28 2021-05-28 一种基于小步长离散化的电力电子级联变换器的建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110595176.6A CN113343453B (zh) 2021-05-28 2021-05-28 一种基于小步长离散化的电力电子级联变换器的建模方法

Publications (2)

Publication Number Publication Date
CN113343453A CN113343453A (zh) 2021-09-03
CN113343453B true CN113343453B (zh) 2023-02-10

Family

ID=77472101

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110595176.6A Active CN113343453B (zh) 2021-05-28 2021-05-28 一种基于小步长离散化的电力电子级联变换器的建模方法

Country Status (1)

Country Link
CN (1) CN113343453B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013114512A (ja) * 2011-11-29 2013-06-10 Chuo Univ 回路シミュレーション方法、回路シミュレーション装置、および回路シミュレーションプログラム

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108347180B (zh) * 2017-01-24 2019-11-05 台达电子企业管理（上海）有限公司 级联变换器系统及其变换器模块投入运行的方法
CN107290977B (zh) * 2017-06-07 2019-09-27 清华大学 后向离散状态事件驱动电力电子仿真方法、设备和介质
EP3640822A4 (en) * 2017-06-14 2020-07-08 Mitsubishi Electric Corporation DEVICE FOR ESTIMATING THE CONDITION
CN110348057B (zh) * 2019-06-12 2020-12-22 华南理工大学 快速获得电力电子dc/dc变换器动态和稳态特性的方法
CN112398402A (zh) * 2020-11-30 2021-02-23 潍柴动力股份有限公司 基于快速有效集的电机模型预测控制方法、装置及控制器

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013114512A (ja) * 2011-11-29 2013-06-10 Chuo Univ 回路シミュレーション方法、回路シミュレーション装置、および回路シミュレーションプログラム

Also Published As

Publication number Publication date
CN113343453A (zh) 2021-09-03

Similar Documents

Publication Publication Date Title
US20200125686A1 (en) Discrete state event-driven simulation method for simulation of power electronic system
US6789237B1 (en) Efficient model order reduction via multi-point moment matching
CN110765584B (zh) 一种含多开关元件的电磁暂态仿真方法、系统、介质及设备
CN113343453B (zh) 一种基于小步长离散化的电力电子级联变换器的建模方法
Lantsov A new algorithm for solving of harmonic balance equations by using the model order reduction method
Iordache et al. ACAP-Analog Circuit Analysis Program
Wichmann et al. On the simplification of nonlinear DAE systems in analog circuit design
US9582622B1 (en) Evaluating on-chip voltage regulation
CN110516276B (zh) 基于fpga的高频开关功率变换器实时仿真方法
Pejovic et al. An algorithm for solving piecewise-linear networks that include elements with discontinuous characteristics
Maffezzoni et al. Event-driven time-domain simulation of closed-loop switched circuits
Li et al. A comparison of steady state methods for power electronic circuits
Tant et al. Accurate second-order interpolation for power electronic circuit simulation
Gil et al. SystemC AMS power electronic modeling with ideal instantaneous switches
CN108846247B (zh) 用于电力变流器系统电磁暂态仿真加速的坐标变换方法
Moroz et al. Why High-Order Integrators Not Rational on Electrical Systems' Computer Calculation
Chung Simulation of PWM switching regulators using linear output predictions and corrections
Sagan et al. Transient adjoint dae sensitivities: a complete, rigorous, and numerically accurate formulation
Nahvi et al. A Quasi-linearisation approach to trajectory based methods for nonlinear MOR
CN115859661A (zh) 一种励磁机的动态仿真方法、装置、设备和介质
CN116011372B (zh) 一种基于分数阶电感的双有源桥变换器分数阶建模方法
Schmidt et al. Fast converter simulation method including parasitic nonlinear capacitances
Durbaba et al. Fast prototyping of a photovoltaic system by using DSP in MATLAB simulation loop
Li et al. PECS. An efficacious solution for simulating switched networks with nonlinear elements
Chiarantoni et al. Applying fixed point homotopy to nonlinear DAEs deriving from switching circuits

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