一种基于微分求积法的时变状态方程解法制造技术
【技术实现步骤摘要】
一种基于微分求积法的时变状态方程解法
本专利技术属于时间域控制系统数值计算领域，特别涉及一种基于微分求积法的时变状态方程解法。
技术介绍
状态向量所满足的向量常微分方程称为控制系统的状态方程，状态方程是控制系统数学模型的重要组成部分，状态方程常常被用来研究复杂的控制系统的动态特性。实际复杂的控制系统大多都是时间域上的控制系统，所建立状态方程的系统矩阵和输入矩阵的元素都是时间的函数，难以求解状态方程的解析解，因此亟需提出一种能够求解复杂时间域控制系统所建立的时变状态方程的求。
技术实现思路
本专利技术为了克服上述现有技术的不足，提供了一种基于微分求积法的时变状态方程解法，本求解方法具备求解精度高、不受时变状态方程矩阵对于时间变量复杂度的限制、易于编程的特点。为实现上述目的，本专利技术采用了以下技术措施：一种基于微分求积法的时变状态方程解法包括以下步骤：S1、建立待求解电路系统的时变状态方程；其中，fi(i＝1,2,…k)为第i个电路系统的时变状态变量，且fi＝fi(t)，fi′(i＝1,2,…k)为第i个电路系统时变状态变量关于时间t的一阶导函数，系统矩阵和输入矩阵均为...
【技术保护点】
一种基于微分求积法的时变状态方程解法，其特征在于，包括以下步骤：S1、建立待求解电路系统的时变状态方程；
【技术特征摘要】
1.一种基于微分求积法的时变状态方程解法，其特征在于，包括以下步骤：S1、建立待求解电路系统的时变状态方程；其中，fi(i＝1,2,…k)为第i个电路系统的时变状态变量，且fi＝fi(t)，fi′(i＝1,2,…k)为第i个电路系统时变状态变量关于时间t的一阶导函数，系统矩阵和输入矩阵均为关于时间t的函数矩阵，其中，aij(i＝1,2,…,k；j＝1,2,…k)为fi′和fj之间的相关系数，mi(i＝1,2,…,k)为第i个输入参数，us(t)为关于时间t的激励函数；S2、在时间域[a,b]上取N个非均匀节点，通过取点公式确定非均匀节点，得到时间域上N个非均匀节点，即a＝t1＜t2＜…＜tN＝b，其中，ti为所取的第i个节点；S3、将待求解的时变状态方程用微分求积的形式代替，即为将时变状态方程中的函数及函数导数用全域上的全部非均匀节点的函数值进行加权求和来表示；S4、将微分求积形式代替后的线性方程组进行移项整理，提取方程组的系数矩阵TEM和BEM；S5、确定初值条件，即待求解电路系统中时变状态变量的初值向量[f1(1)f2(1)…fk(1)]，其中，fi(1)＝fi(t1),(i＝1,2,…,k)，对线性方程组的系数矩阵进行修正，去掉由每一个初始微分方程形成线性方程组的第一个方程，将与初值有关的项移到等式的右端，得到修正后线性方程组求解的系数矩阵TEM和BEM；S7、对修正后线性方程组进行求解，得到各个时变状态变量在节点处的值，即为f1(1),f1(2),…,f1(N)；f2(1),f2(2),…,f2(N)；………；fk(1),fk(2)，…fk(N)，其中，fi(j)＝fi(tj),(i＝1,2,…,k；j＝1,2,…,N)，n为状态变量的个数，N为节点的个数。2.如权利要求1所述的一种基于微分求积法的时变状态方程解法，其特征在于，步骤S3的具体操作步骤包括：S31、一维函数f(t)为时变状态变量fi(i＝1,2,…,k)中的任意一个，它在时间域[a,b]上连续可微，则其中，L为线性微分算子，Wj(t)为插值基函数，tj为N个互异节点a＝t1＜t2＜…＜tN＝b中第j个节点的坐标值，令Aij＝Wj(ti)、fj＝f(tj)，则其中，Aij为函数f(t)一阶导数的加权系数，则[Aij]为一阶导数的加权系数矩阵；S32、通过拉格朗日插值方法得到微分求积法加权系数的显式表达式即为Aij，Aij由下式计算得到：S33、将步骤S1中的待求解电路系统的时变状态方程中的微分项用微分求积的形式代替，即fij＝fi(tj)(i＝1,2,…,k；j＝1,2,…,N)，得到如下的方程组：3.如权利要求2所述的一种基于微分求积法的时变状态方程解法，其特征在于，步骤S4的具体操作步骤包括：S41、将微分求积形式代替后的线性方程组进行移项整理，可得如下标准形...
【专利技术属性】
技术研发人员：陈波，张斌斌，储昭碧，胡永新，王有康，李华，
申请(专利权)人：合肥工业大学，
类型：发明
国别省市：安徽,34
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。