买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：长春理工大学

摘要：基于真三维显示系统的Roesser模型及其实现方法，属于多维系统控制理论技术领域，为了解决现有技术存在的问题，真三维显示系统的Roesser模型实现法将三维位置坐标与一维时间坐标相结合，对三维体空间中的体素点进行状态空间表示，根据系统的线性因果性，通过已知的传递函数，构造初始矩阵，对其进行初等变换和补充运算，从而得到Roesser实现矩阵；相较于已有方法，该实现方法得到的实现矩阵阶次更低，可以分析各变量对实现矩阵的影响，因此既能简化系统数学表达，又有利于系统分析与设计，具有极高的应用价值。

主权项：1.基于真三维显示系统的Roesser模型的实现方法，该模型为：在真三维显示系统中，将体素点的三维位置坐标与一维时间坐标结合，建立四维Roessor状态空间模型，如下所示： yn1,n2,n3,t＝Cxn1,n2,n3,t+Dun1,n2,n3,t式中： n1∈Z，n2∈Z，n3∈Z，t∈Z，向量xh∈Ra，xv∈Rb，xl∈Rc，xt∈Rd分别为x轴方向，y轴方向，z轴方向和时间t轴的向量；输入向量u∈Rp，输出向量y∈Rq；然后A1∈Ra×a，A2∈Ra×b，A3∈Ra×c，A4∈Ra×d，A5∈Rb×a，A6∈Rb×b，A7∈Rb×c，A8∈Rb×d，A9∈Rc×a，A10∈Rc×b，A11∈Rc×c，A12∈Rc×d，A13∈Rd×a，A14∈Rd×b，A15∈Rd×c，A16∈Rd×d，B1∈Ra×P，B2∈Rb×P，B3∈Rc×P，B4∈Rd×P，C∈Rq×a+b+c+d，D∈Rq×P；在时隙t时，对于体素点n1,n2,n3，分别从体素n1-1,n2,n3、n1,n2-1,n3和n1,n2,n3-1接收状态向量组xhn1,n2,n3,t，xvn1,n2,n3,t和xln1,n2,n3,t；用Roesser模型去计算向量xhn1+1,n2,n3,t，xvn1,n2+1,n3,t，xln1,n2,n3+1,t和xtn1,n2,n3,t+1；发送向量xhn1+1,n2,n3,t，xvn1,n2+1,n3,t和xln1,n2,n3+1,t；其特征是，该方法包括以下步骤：步骤一，令各系数矩阵：C＝[C1C2C3C4]，根据传递函数的定义，该四维Roesser模型的输出量yn1,n2,n3,t的z变换nz1,z2,z3,z4与输入量un1,n2,n3,t的z变换dz1,z2,z3,z4之比： 对系统进行z变换，则系统对应的传递函数为：Hz1,z2,z3,z4＝CZIr-AZ-1B+D3其中，对角阵Z＝diag{z1Ia,z2Ib,z3Ic,z4Id}，阶数r＝a+b+c+d；根据因果性D＝0；步骤二，定义四维多项式初始矩阵为 构造矩阵： 该矩阵应具有如下性质：a对角线上第一个元素只能是x；b对角线上其他元素只能是关于某一变量zk,k∈{1,2,3,4}的一维线性多项式，且常数项只能是1；c除第一行外，非对角线元素只能为关于某一变量zk,k∈{1,2,3,4}的线性单项式；d除x外，第一行中的元素均为常数项；e同一行的元素只能包含同一个元素zk,k∈{1,2,3,4}，且从第二行开始，所有的行都是按照z1,z2,z3,z4的顺序排列的；真三维显示系统的Roesser模型实现方法即为通过矩阵初等变换和补充运算，将初始矩阵M0变换为M；矩阵M0中对角线上第一个元素x在这里只是一个符号标志，而不是一个变量，在对初始矩阵M0变换的过程中，不能改变x的位置和表达式；这是因为，为了满足性质a和b，在对初始矩阵M0变换的过程中，不能对第一行进行任何变换运算，设置x就是为了阻止对第一行和第一列进行变换运算；步骤三，设任意一不含常数项的四维多项式为p'z1,z2,z3,z4，对其中某个变量zk,k∈{1,2,3,4}显然可以分解成如下形式：p'z1,z2,z3,z4＝p1zk+p2z1,...,zk-1,zk+1,...,z4+zkp3z1,z2,z3,z46其中，p1zk是只含有zk的一维线性多项式，p2z1,...,zk-1,zk+1,...,z4是不含有zk的四维线性多项式，且p1zk，p2z1,...,zk-1,zk+1,...,z4，p3z1,z2,z3,z4均不含常数项；令则初始矩阵M0为 令根据式6，将和分解为 式中，p1和q1是只含z1的一维线性单项式，p2和q2是不含z1的三维线性多项式，p3和q3是四维多项式；接着进行如下运算：M1＝augmentM0M1＝addrowM1,3,2,-z1M1＝addcolM1,3,2,p3 接着依次对p3、q3及M1中每次运算所产生的新的行进行运算，使得M1中含z1的项都变为关于z1的线性单项式；M1中的每一行依次对变量z1、z2、z3、z4进行同样的运算，使得M1中除了x之外的对角线元素均为常数项为1的4维线性多项式，而除第一行外的非对角元素均为不含常数项的4维线性多项式；步骤四，假设通过步骤一得到的矩阵M1为 式中，*和#都是线性多项式，ai,bi,ci,di,i＝{1,2,3,4}都是系数；将第一步中所得到的矩阵M1转化成M2，使得M2中除x外的对角线元素均为常数项为1的一维线性多项式，而除第一行外的非对角元素均为关于某个变量zk,k∈{1,2,3,4}的线性单项式；对式10中的M1进行如下运算M2＝augmentM1；M2＝addrowM2,4,2,-1；M2＝addcolM2,4,1,a2z2；M2＝addcolM2,4,2,b2z2； 在同一行中，对变量z3,z4进行相似的运算，最终得到矩阵M2，使得M2中第二行的前三个元素分别为a1z1，1+b1z1和c1z1，其余元素均为-1； 步骤五，通过适当的行变换和列变换，将M2中的每一行按照z1、z2、z3、z4的顺序排列，并将所有的一维线性多项式元素移到对角线位置，再通过列变换消去-1项，得到矩阵M3； 根据式5得出矩阵A,B,C。

全文数据：基于真三维显示系统的Roesser模型及其实现方法技术领域[0001]本发明涉及一种基于真三维显示系统的Roesser模型及其实现方法，属于多维系统控制理论技术领域。背景技术[0002]在真三维显示系统中，显示物体的每一个体素位于三维物理空间中的真实位置，且具有唯一的三维位置坐标和一维时间坐标。系统工作时，将原始三维图像数据转换成符合显示要求排列的体素点，送入显示单元，则每个体素点的亮度和色彩可控，体素点之间的相对空间位置关系被真实的体现，从而构成可直接观察的三维图像。[0003]Roesser模型是控制理论中最重要也是应用最广泛的模型之一，可以使用系统内部的状态变量来表征系统的动态特性。它由描述系统输出变量与状态变量间变换关系的输出方程组成。结合三维位置坐标与一维时间坐标，运用四维Roesser模型对真三维显示系统进行状态空间表示，能够简化系统的数学表达式，反映系统中四个独立变量的变化，进而确定系统的全部内部运动状态，能够反映系统动态行为的全部信息。因此，找到Roesser模型的实现方法，对真三维显示系统系统进行状态空间建模具有极高的应用价值。[0004]现有的RoeSser模型实现方法有徐粒-阎石EOA变换法与ψ矩阵法。徐粒-阎石EOA变换法采用ZllZ2,...，zn表示传递函数，引入EOA变换法给出了实现矩阵的求法，然而这种变换法中每次补充运算只能降低一阶。虽然可以通过一些分解来提高补充效率，但实际上只有极少数传递函数适用这种分解法。Ψ矩阵法所求得的GR模型实现矩阵阶次更低，但不能分析系数值对实现矩阵的影响，因此不宜采用。发明内容[0005]本发明为了解决现有技术存在的问题，提出一种基于真三维显示系统的Roesser模型及其实现方法。[0006]本发明的技术方案是：[0007]基于真三维显示系统的Roesser模型，其特征是，该模型为:在真三维显示系统中，将体素点的三维位置坐标与一维时间坐标结合，建立四维Roessor状态空间模型，如下所示：[0008][0009]ym，n2，n3，t=Cxm，n2，n3，t+Dum，n2，n3，t[0010]式中：[0011]mez，n2ez，n3ez，tez，向量xheR^xVeRbJ1eRc^xtGRd分别为X轴方向，y轴方向，z轴方向和时间t轴的向量;输入向量ueRp，输出向量yeRq;然后A1GRax^A2GRaxb,A3eRaXc，A4eRaxd，A5eRbxa，A6eRbxb，A?eRbxc，A8eRbxd，A9eRcXa，AioeRcxb，AneRcXc，Ai2eRcxd，Ai3eRdxa，A14eRdxb，Ai5eRdxc，Ai6eRdxd，BieRaxp，B2eRbxp，B3eRcxp，B4eRdxp，CeRqXa+b+c+d，DeRqXP。[0012]在时隙t时，对于体素点ηι，Π2，η3，分另Ij从体素ni_l，Π2，η3、（ηι，Π2-1，Π3和m,112，113-1接收状态向量组叉11111，112，113,1：，叉'?111，112，113,1：和叉1111，112，113,1：;用1?068861'模型去计算向量xhm+1，n2，n3,t，xvηι，Π2+1,Μ,t,X1ηι，Π2，n3+l，t和Xtni，n2，n3,t+1;发送向量xhm+1，Π2，η3,t，χνηι，η2+1，η3,t和X1ηι，Π2，η3+1，t。[0013]基于真三维显示系统的R0esser模型的实现方法，其特征是，包括以下步骤：[00Μ]步骤一，[0015]令各系数矩阵，C=[C1C2C3C4]，[0016]根据传递函数的定义，该四维Roesser模型的输出量yni，n2，n3，t的Z变换ηΖ1,Ζ2，Ζ3，Ζ4与输入量um，n2，n3，t的ζ变换dΖ1，Ζ2，Ζ3，Ζ4之比：[0017]2[0018]对系统进行ζ变换，则系统对应的传递函数为：[0019]Hzi,Z2,Z3,ζ4=CZIr-AZ_1B+D3[0020]其中，对角阵2=虹区{211£1，2211，231。，241]}，阶数1"=+匕+0+1;[0021]根据因果性D=〇;[0022]步骤二，定义四维多项式初始矩阵为[0023]4[0024]构造矩阵：[0025][0026]该矩阵应具有如下性质：[0027]a对角线上第一个元素只能是X;[0028]b对角线上其他元素只能是关于某一变量Zk，ke{1，2,3,4}的一维线性多项式，且常数项只能是I;[0029]c除第一行外，非对角线元素只能为关于某一变量zk，ke{1，2,3,4}的线性单项式；[0030]d除X外，第一行中的元素均为常数项；[0031]e同一行的元素只能包含同一个元素zk，ke{1，2,3,4}，且从第二行开始，所有的行都是按照Zl，Z2，Z3，Z4的顺序排列的；[0032]真三维显示系统的Roesser模型实现方法即为通过矩阵初等变换和补充运算，将初始矩阵Mo变换为M;[0033]矩阵Mo中对角线上第一个元素X在这里只是一个符号标志，而不是一个变量，在对初始矩阵Mo变换的过程中，不能改变X的位置和表达式;这是因为，为了满足性质a和b，在对初始矩阵Mo变换的过程中，不能对第一行进行任何变换运算，设置X就是为了阻止对第一行和第一列进行变换运算；[0034]步骤三，设任意一不含常数项的四维多项式为p’（Z1，Z2，Z3，Z4，对其中某个变量zk，ke{1，2,3,4}显然可以分解成如下形式：[0035]P，（Z1，Z2，Z3，Z4=PlZk+P2Zl,···，Zk-l，Zk+l,···，Z4+ZkP3Z1，Z2，Z3，Z46[0036]其中，PlZk是只含有Zk的一维线性多项式，P2Zl,…，Zk-l，Zk+l,…，Z4是不含有Zk的四维线性多项式，且PlZk，P2Zl,...，Zk-l，Zk+l,...，Z4，P3Z1，Z2，Z3，Z4均不含常数项；[0037]令：，则初始矩阵Mo为[0038]7[0039]令，根据式⑶，将》和#分解为[0040]p=+;：+Zi^,q=q+ty：+z-,q·.8[0041]式中，PdPq1是只含Z1的一维线性单项式，PdPq2是不含Z1的三维线性多项式，P3和q3是四维多项式；[0042]接着进行如下运算：[0043]Mi=augmentMo[0044]Mi=addrowΜι,3,2,-ζι[0045]Mi=addcoIΜι,3,2，p3[0046]9[0047]接着依次对p3、q3及施中每次运算所产生的新的行进行运算，使得M1中含Z1的项都变为关于Zi的线性单项式；[0048]Ml中的每一行依次对变量Zl、Z2、Z3、Z4进行同样的运算，使得Ml中除了X之外的对角线元素均为常数项为1的4维线性多项式，而除第一行外的非对角元素均为不含常数项的4维线性多项式；[0049]步骤四，假设通过步骤一得到的矩阵施为[0050][0051]式中，*和#都是线性多项式，1匕，^，^4={1，2,3,4}都是系数；[0052]将第一步中所得到的矩阵M1转化成M2,使得M2中除X外的对角线元素均为常数项为1的一维线性多项式，而除第一行外的非对角元素均为关于某个变量Zk，ke{1，2,3,4}的线性单项式；[0053]对式（10中的M1进行如下运算[0054]M2=augmentMi；[0055]M2=addrowM2，4，2，-1;[0056]M2=addcolM2,4,1，a2Z2;[0057]M2=addcolM2,4,2，b2Z2;[0058][0059]在同一行中，对变量z3，Z4进行相似的运算，最终得到矩阵跑，使得M2中第二行的前三个元素分别为aizi，Ι+bizi和Cizi，其余元素均为-1;[0060]12；[0061]步骤五，通过适当的行变换和列变换，将M2中的每一行按照Zl、Z2、Z3、Z4的顺序排列，并将所有的一维线性多项式元素移到对角线位置，再通过列变换消去-1项，得到矩阵M3;[0062]13[0063]根据式⑶得出矩阵A，B，C。[0064]本发明的有益效果：[0065]本发明运用四维Roesser模型对真三维显示系统进行状态空间表示，可以表征和计算体素变量，使得在时隙t计算的各变量只与相邻的前一体素或前一时刻有关，从而能够简化运算，便于微控制器的软件实现。[0066]与徐粒-阎石EOA变换法EOA变换法相比，本发明的四维Roesser模型实现方法得到的实现矩阵阶次更低，又克服了Ψ矩阵法不能分析系数值对实现矩阵的影响的缺点，因此既能简化系统数学表达，又有利于系统分析与设计，具有极高的应用价值。附图说明[0067]图1是本发明基于真三维显示系统的Roesser模型示意图。具体实施方式[0068]下面结合附图对本发明做进一步详细说明。[0069]如图1所示，基于真三维显示系统的Roesser模型，该模型为：在真三维显示系统中，显示物体的体素点按显示要求有序排列，且具有唯一的三维位置坐标和一维时间坐标，将体素点的三维位置坐标与一维时间坐标结合，建立四维Roessor状态空间模型，如下所示：[0070][0071]yni，m，M，t=Cxni，m，M，t+Duni，m，M，t[0072]式中：[0073]mez，n2ez，n3ez，tez，向量xheR^xVeRbJ1eRc^xtGRd分别为X轴方向，y轴方向，z轴方向和时间t轴的向量。输入向量ueRp，输出向量yeRq。然后A1GRax^A2GRaxb,A3eRaXc，A4eRaxd，A5eRbxa，A6eRbxb，A?eRbxc，A8eRbxd，A9eRcXa，AioeRcxb，AneRcXc，Ai2eRcxd，Ai3eRdxa，Ai4eRdxb，Ai5eRdxc，Ai6eRdxd，BieRaxp，B2eRbxp，B3eRcxp，B4eRdxp，CeRqXa+b+c+d，DeRqXP。[0074]在时隙t时，对于体素点m，n2，n3，分别从体素ni_l，Π2，η3、（ηι，Π2-1，n3和ni,112，113-1接收状态向量组叉11111，112，113,1：，叉'?111，112，113,1：和叉1111，112，113,1：;用1?068861'模型去计算向量xhm+1，Π2，η3,t，xvηι，Π2+1，η3,t,X1ηι，Π2，η3+1，t和Xtηι，Π2，η3,t+1;发送向量xhm+1，Π2，η3,t，χνηι，η2+1，η3,t和X1ηι，Π2，η3+1，t。[0075]基于真三维显示系统的Roesser模型的实现方法，包括以下步骤：[0076]步骤一，[0077]令各系数矩阵：£=[C1C2C3C4]，[0078]对系统进行z变换，则系统对应的传递函数为：[0079]Hzi，Z2，Z3，Z4=CZIr-AZ_1B+D2[0080]其中，对角阵2=虹区{211£1，2211，231。，241]}，阶数『=+匕+0+1〇[0081]根据传递函数的定义，该四维Roesser模型的输出量yni，n2，n3，t的z变换ηZ1,Ζ2，Ζ3，Ζ4与输入量um，n2，n3，t的ζ变换dΖ1，Ζ2，Ζ3，Ζ4之比：[0082]3[0083]每个体素点的显示是一个线性过程，由于物体的三维数学模型一定可实现可观测，即是一个因果过程，则整个体素空间的显示过程是一个第一卦限上的线性因果关系的过程，所以D=0。[0084]令M表示一个矩阵，定义以下符号表示的操作为：[0085]addrowM，i，j，k:将矩阵M的第i行所有元素乘以k后加到第j行对应的元素上；[0086]addcolM，i，j，k:将矩阵M的第i列所有元素乘以k后加到第j列对应的元素上；[0087]swapr〇WM，i，j:交换矩阵M的第i行和第j行对应的元素；[0088]swapcolM，i，j:交换矩阵M的第i列和第j列对应的元素；[0089][0090]步骤二，定义四维多项式初始矩阵为：[0091]⑷[0092]构造矩阵：[0095]该矩阵具有如下特点：[0096]a对角线上第一个元素只能是X;[0097]b对角线上其他元素只能是关于某一变量Zk，ke{1，2,3,4}的一维线性多项式，且常数项只能是1;[0098]c除第一行外，非对角线元素只能为关于某一变量Zk，keU，2,3,4}的线性单项式；[0099]d除X外，第一行中的元素均为常数项；[0100]e同一行的元素只能包含同一个元素Zk，ke{1，2,3,4}，且从第二行开始，所有的行都是按照Zl，Z2，Z3，Z4的顺序排列的。[0101]真三维显示系统的Roesser模型实现方法即为通过矩阵初等变换和补充运算，将初始矩阵Mo变换为M。[0102]矩阵Mo中对角线上第一个元素X在这里只是一个符号标志，而不是一个变量，在对初始矩阵Mo变换的过程中，不能改变X的位置和表达式。这是因为，为了满足性质a和b，在对初始矩阵Mo变换的过程中，不能对第一行进行任何变换运算，设置X就是为了阻止对第一行和第一列进行变换运算。[0103]步骤三，设任意一不含常数项的四维多项式为p’（Z1，Z2，Z3，Z4，对其中某个变量zk，ke{1，2,3,4}显然可以分解成如下形式：[0104]ρ’ζι，Z2，Z3，Z4=PlZk+P2Zl,···，Zk-l，Zk+l,···，Z4+ZkP3Z1，Z2，Z3，Z47[0105]其中，PlZk是只含有Zk的一维线性多项式，P2Zl,…，Zk-l，Zk+l,…，Z4是不含有Zk的四维线性多项式，且PlZk，P2Zl,...，Zk-l，Zk+l,...，Z4，P3Z1，Z2，Z3，Z4均不含常数项。[0106]令，23.,24.二與公,公，為,.24-1,则初始矩阵Mo为[0107]8[0108]令p.=czZl，c}=Zl，根据式7，将》和I分解为[0109]p=P'+P:+:、p'，q=q：+q：+z：q：Q[0110]式中，P^q1是只含Z1的一维线性单项式，P^q2是不含Z1的三维线性多项式，P3和Φ是四维多项式。[0111]接着进行如下运算：[0112]Mi=augmentMo[0113]Mi=addrowΜι,3,2,-ζι[0114]Mi=addcolΜι，3，2，p3[0115]C10^[0116]接着依次对P3、q3及施中每次运算所产生的新的行进行运算，使得M1中含Z1的项都变为关于Zi的线性单项式。[0117]Ml中的每一行依次对变量21、22、23、24进行同样的运算，使得施中除了1之外的对角线元素均为常数项为1的4维线性多项式，而除第一行外的非对角元素均为不含常数项的4维线性多项式。[0118]步骤四，假设通过步骤一得到的矩阵施为[0119][0120]式中，*和#都是线性多项式，ai，bi，Ci，cU，i={1，2，3，4}都是系数。[0121]将第一步中所得到的矩阵M1转化成M2,使得M2中除X外的对角线元素均为常数项为1的一维线性多项式，而除第一行外的非对角元素均为关于某个变量Zk，ke{1，2,3,4}的线性单项式。[0122]对式（11中的M1进行如下运算[0123]M2=augmentMi；[0124]M2=addrowM2,4,2,-l；[0125]M2=addcolM2,4,1，a2Z2;[0126]M2=addcolM2,4,2，b2Z2;[0127]V：=CiddcoliyM；,4,3,c'：z：[0128][0129]在同一行中，对变量Z3，Z4进行相似的运算，最终得到矩阵跑，使得M2中第二行的前三个元素分别为aizi，Ι+bizi和Cizi，其余元素均为-1。[0130]13[0131]步骤五，通过适当的行变换和列变换，将M2中的每一行按照Zl、Z2、Z3、Z4的顺序排列，并将所有的一维线性多项式元素移到对角线位置，再通过列变换消去-1项，得到矩阵M3〇[0132]14[0133]根据式⑶得出矩阵A，B，C。[0134]实施例：[0135]已知严密因果传递函数：[0136][0137]由严密因果性，得:D=H0,0,0,0=0。[0138]构造初始矩阵：[0139][0140]对其做如下变换：[0141]Mi=augmentMo[0142]Mi=addrowΜι,3,2,-ζι[0143]Mi=addcolMi，3,1，Z2+Z3+Z4[0144]M2=augmentMi[0145]M2=addrowM2,4,2,-Z2[0146]M2=addcolM2,4,2，Z3[0147]M3=augmentM2[0148]M3=addrowM3,5,2,-1[0149]M3=addcolM3,5,4，_Z2[0150]M4=augmentM3[0151]M4=addrowM4,6,3，_1[0152]M4=addcolM4,6,1，Z3+Z4[0153]M4=addcolM4,6,2，Z3+Z4[0154]Ms=augmentM4[0155]M5=addrowMs,7,6,-1[0156]]\fe=addcolMs，7,1，Z4[0157]]\fe=addcolMs，7,2，Z4[0158]Ms=swaprowMs,4,5[0159]Ms=swapcoIMs,4,5[0160]M5=addcolMs,2,4,1[0161]M5=addcolMs,3,6,1[0162]得:[0163][0164]根据式12，得:a=l氺=2,〇=2，1=1，即实现矩阵的阶数『=63[0165][0166]C=[II00I1],D=0〇

权利要求：1.基于真三维显示系统的Roesser模型，其特征是，该模型为:在真三维显示系统中，将体素点的三维位置坐标与一维时间坐标结合，建立四维Roessor状态空间模型，如下所示：式中：111£2，112£2，113£2八£2，向量叉11£俨4'^妒，叉1£1^，叉*£1?1分别为叉轴方向，7轴方向，2轴方向和时间t轴的向量;输入向量uGRp，输出向量A4eRaxd，A5eRbxa，A6eRbxb，A?eRbxc，A8eRbxd，A9eRcXa，AioeRcxb，AneRcXc，Ai2eRcxd，Ai3eRdxa，A14eRdxb，Ai5eRdxc，Ai6eRdxd，BieRaxp，B2eRbxp，B3eRcxp，B4eRdxp，CeRqxa+b+c+1，DGRqxp02.根据权利要求I所述的基于真三维显示系统的Roesser模型，其特征是，在时隙t时，对于体素点m，n2，n3，分别从体素ηι-1，Π2，η3、（ηι，Π2-1，n3和m，n2，n3-l接收状态向量组xhm，n2，n3,t，xvm，n2，n3,t和X1ηι，Π2，η3,t;用Roesser模型去计算向量xhm+1，]12，113,1：，叉'?111，112+1，113,1：，叉1111，112，113+1，1：和叉1：111，112，113，七+1;发送向量叉11111+1，112，113，t，xvni，n2+l，n3，t和X1ni，ri2，n3+l，t〇3.基于真三维显示系统的Roesser模型的实现方法，其特征是，包括以下步骤：步骤一，令各系数矩阵：根据传递函数的定义，该四维Roesser模型的输出量ym，Π2，m，t的z变换nZi，Z2，Z3，Ζ4与输入量um，n2，n3，t的ζ变换dζι，Ζ2，Ζ3，Ζ4之比：对系统进行ζ变换，则系统对应的传递函数为：HZ1,Z2,Z3,Z4=CZIr-AZ_1B+D3其中，对角阵7=118{211£1，2211，231。，241]}，阶数1=+匕+0+1;根据因果性D=O;步骤二，定义四维多项式初始矩阵为构造矩阵：该矩阵应具有如下性质：a对角线上第一个元素只能是X;b对角线上其他元素只能是关于某一变量Zk，ke{1，2,3,4}的一维线性多项式，且常数项只能是1;c除第一行外，非对角线元素只能为关于某一变量zk，ke{l，2,3,4}的线性单项式；d除X外，第一行中的元素均为常数项；e同一行的元素只能包含同一个元素zk，ke{1，2,3,4}，且从第二行开始，所有的行都是按照Zl，Z2，Z3，Z4的顺序排列的；真三维显示系统的Roesser模型实现方法即为通过矩阵初等变换和补充运算，将初始矩阵Mo变换为M;矩阵Mo中对角线上第一个元素X在这里只是一个符号标志，而不是一个变量，在对初始矩阵Mo变换的过程中，不能改变X的位置和表达式;这是因为，为了满足性质a和b，在对初始矩阵Mo变换的过程中，不能对第一行进行任何变换运算，设置X就是为了阻止对第一行和第一列进行变换运算；步骤三，设任意一不含常数项的四维多项式为P’（21，22,23,24，对其中某个变量zk，ke{1，2,3,4}显然可以分解成如下形式：P’（Z1，Z2，Z3，Z4=PlZk+P2Z1,···，Zk-l，Zk+l,···，Z4+ZkP3Z1，Z2，Z3，Z46其中，PlZk是只含有Zk的一维线性多项式，P2Zl,…，Zk-l，Zk+l,…，Z4是不含有Zk的四维线性多项式，且PlZk，P2Z1,...，Zk-l，Zk+l,...，Z4，P3Z1，Z2，Z3，Z4均不含常数项；令，则初始矩阵M0为令.根据式6，将分解为式中，P^q1是只含Z1的一维线性单项式，P^q2是不含Z1的三维线性多项式，P3和q3是四维多项式；接着进行如下运算：Mi=augmentMoMi=addrowΜι,3,2,-ζι接着依次对P3、q3及％中每次运算所产生的新的行进行运算，使得M1中含Z1的项都变为关于Zi的线性单项式；Ml中的每一行依次对变量Z1、Z2、Z3、Z4进行同样的运算，使得Ml中除了X之外的对角线元素均为常数项为1的4维线性多项式，而除第一行外的非对角元素均为不含常数项的4维线性多项式；步骤四，假设通过步骤一得到的矩阵施为式中，*和#都是线性多项式而，1^，^，1^={1，2,3,4}都是系数；将第一步中所得到的矩阵M1转化成M2,使得M冲除X外的对角线元素均为常数项为1的一维线性多项式，而除第一行外的非对角元素均为关于某个变量zk，ke{1，2,3,4}的线性单项式；对式（10中的M1进行如下运算在同一行中，对变量Z3，Z4进行相似的运算，最终得到矩阵M2，使得M2中第二行的前三个元素分别为aizi，Ι+bizi和CiZi，其余元素均为-1;步骤五，通过适当的行变换和列变换，将M2中的每一行按照Zl、Z2、Z3、Z4的顺序排列，并将所有的一维线性多项式元素移到对角线位置，再通过列变换消去-1项，得到矩阵M3;根据式⑶得出矩阵A，B，C。

百度查询： 长春理工大学 基于真三维显示系统的Roesser模型及其实现方法