动态系统建模理论与应用习题

《动态系统建模理论与应用》课程习题一、选择题：答案唯一，在（）内填入正确答案的编号。1.对于批量最小二乘格式LLLEY，其最小二乘无偏估计的必要条件是（）。A.输入序列}{ku为“持续激励”信号B.LE与TLLTL1)(正交C.LE为非白噪声向量D.0}{LEE2.对象模型为Tkkkye时，采用递推最小二乘估计后的残差序列的计算式为（）。A.1ˆTkkkkyB.1ˆTkkkkyC.ˆTkkkkyD.11ˆTkkkky3.在上题的条件下，递推最小二乘算法中的增益矩阵kK可以写成（）。A.11kkPB.1kkPC.1kkPD.kkP4.可以同时得到对象参数和干扰噪声模型参数的估计算法是（）。A.辅助变量法B.广义最小二乘法C.最小二乘限定记忆法D.相关最小二乘两步法5.增广最小二乘估计的关键是（）。A.将控制项增广进k中，并用残差项取代进行估计B.将输出项增广进k中，并用残差项取代进行估计C.将噪声项增广进k中，并用残差项取代进行估计D.将噪声项增广进k中，并用输出项取代进行估计答案：1.B2.C3.D4.B5.C■二、判断题：以○表示正确或×表示错误。1．估计残差平方和最小是确定辨识过程对象结构的唯一标准。（）2．最小二乘估计的批量算法和递推算法在数学上是等价的。（）3．广义最小二乘法就是辅助变量法和增广最小二乘法交替试用。（）4．在递推最小二乘算法中，若置0TkPPP，则该算法也能克服“数据饱和”现象，进而可适用于时变系统。（）5．用神经网络对SISO非线性系统辨识，采用的是输入层和输出层均为一个神经元的三层前馈神经元网络结构。（）答案：1.×2.○3.×4.○5.×■三、设y和n21x,x,x之间满足关系)xaxaxa(expynn2211，试图利用y和n21x,x,x的观测值来估计参数n21a,a,a，请将该模型化成最小二乘格式。答案：Tnn2211xaxaxaln(y)z其中，n21Tn21Tx,,x,xa,,a,a■四、对于多输入单输出（MISO）系统可由下面的模型描述kkkeuzByzA111)()(其中，ku为系统的m×1维输入向量；ky为系统的标量输出；ke为标量i.i.d随机噪声；1z为延迟算子，即11kkyyz；)(1zA为标量参数多项式，)(1zB为1×m的参数多项式向量：aannzazazA.1)(111bbnnzBzBBzB.)(1101请写出：最小二乘递推算法公式和计算步骤或流程。答案：根据题意，可写出最小二乘格式为：kTkkey其中，TnkTkTknkkkTkbauuuyyy12121,,;,,1201,,,;,,,abTnnaaaBBB因此，采用批量最小二乘法估计时，设采集数据时刻为k=1,2,…,L，则有批量最小二乘格式为：LLLEY其中，LLyyyY21，TLTTL21，LLeeeE21从而，批量最小二乘估计公式为：LTLLTLY1)(ˆ递推最小二乘估计公式为：)ˆ(ˆˆ11kTkkkkkyKkkTkkkkPPK111，kkTkkTkkkkkPPPPP11111初始估计：0ˆ0，IP20，2是一个充分大的正数。计算流程为：（0）给定0,,ˆ00kP；（1）量测1ky，组成Tk1；（2）计算1kK；（3）计算1ˆk；（4）输出估计结果，并由误差限或数据长度L来确定是否停止估计。若条件满足，则停止估计；否则，继续进行。（5）计算1kP；（6）1kk，返回到（1）。■五、对于SISO系统的数学模型kkkvuzByzA111)()(其中，ku和ky分别为系统的输入输出量，kv为干扰噪声，)(1zA和)(1zB为参数多项式：aannzazazA1111)(bbnnzbzbbzB1101)(且bann，1z为延迟算子，即11kkyyz。1．对于量测ku、ky，Nk,2,1，写出估计系统参数的最小二乘批量算法详细公式。2．给出最小二乘法无偏估计的条件并加以证明。3．简述辨识动态系统数学模型的一般步骤。答案：1．由题意可知，采用L次测量的批量最小二乘格式可写为：NNNVY其中，12121,,;,,abTkkkknkkknyyyuuu1201,,,;,,,abTnnaaabbbNNyyyY21，TkTTN21，NNvvvV21因此，最小二乘批量算法公式为：NTNNTNY1)(ˆ■2．证明：))()()()()()ˆ(111NTNNTNNNTNNTNNTNNTNVEEVEYEE当1)(NTN和NTNV不相关时，上式第二项为零，最小二乘估计为无偏估计，}{kv为零均值独立随机序列时，此条件自然满足。此时，)()ˆ(EE。■3．辨识动态系统数学模型的一般步骤为：Step1:确定建模目的，并由工艺和物理/化学过程初步确定模型形式和结构；Step2:试验设计：包括试验信号设计、采样周期选择、实验数据长度选定、试验方式（离线/在线）等；Step3:实际系统试验，采集输入输出数据，并进行数据的预处理；Step4:模型结构假设，选定阶次范围；Step5:选供适用算法进行参数估计，得到一组数学模型；Step6:模型结构的确定，得到一个数学模型；Step7:模型检验；根据检验结果，可能要从Step2到Step6中的任何一步重新做起。Step8:若模型检验合格，则得到最终模型。■六、某系统的动态模型为11kkkkyaybue，假设：系统是稳定的，且ke和ku都为零均值广义平稳随机序列。采用辅助变量法进行参数估计，进行L次量测，且L充分大,试证明：0111LkkLLuuZuuuu是一个合适的辅助变量矩阵。答案：证明：辅助变量法的计算公式为1111ˆ()()TTTTIVLLLLLLLLZZYZZYLL根据题义有0001111111kLTLLkkkLLLyuuuuZyuuuuyu1111111111LLiiiiiiLLiiiiiiyuuuyuuu10111kLTLLkkLLyuuuZYyuuuy111LiiiLiiiyuyu当对象的输入输出序列都是具有遍历性的平稳随机序列时，则输入信号的自相关函数和输入输出的互相关函数为111lim1limLuuk-lkLkLuyk-lkLkRluuLRluyL因{vk}、{uk}和{yk}均为是零均值广义平稳随机序列，所以，(0)(0)1lim(1)(1)yuuuTLLLyuuuRRZRRL(1)1lim(0)yuTLLLyuRZYRL1ˆ(0)(0)(1)(1)(0)(1)1ˆlimlimˆ(1)(1)(0)(1)(0)(0)yuuuyuuuuuyuIVLLyuuuyuyuyuyuaRRRRRRppRRRRRRb2(1)(1)(0)(0)1(1)(1)(0)uuyuuuyuyuyuyuRRRRRRR式中，(0)(1)(1)(0)yuuuyuuuRRRR又相关函数kk-1(1)(0)(0)yuyuuuREyuaRbRkk(0)(1)(1)yuyuuuREyuaRbR可得(1)(0)(0)(0)(1)(1)ˆˆ1limˆˆ(1)(0)(0)(0)(1)(1)uuyuuuuuyuuuNyuyuuuyuyuuuRaRbRRaRbRaapbbRaRbRRaRbR由此可知，ZN矩阵是一个合适的辅助变量矩阵。■七、在递推最小二乘估计中，新息的表达式为11,ˆ~kTkkkkyy。1.请写出残差的表达式kky,~2.证明：)1/(~~11,,kkTkkkkkPyy答案：1.kTkkkkyyˆ~,2.证明：)1/(~))1/(1(~)1(~~~))ˆ(ˆ(ˆ~11,111,1,1,1,11,kkTkkkkkTkkkTkkkkTkkkkkkTkkkkTkkkkTkkkTkkkkPyPPyKyyKyyKyyy■八、请证明：在递推最小二乘估计中kkkKP。证明：在递推最小二乘估计中kkTkkkkPPK111，kkTkkTkkkkkPPPPP11111111111111(1)111TTkkkkkkkkkkkkkkTTkkkkkkkkkTkkkPPPPPPPPPKP■九、在数据预处理环节，去除数据中的一价趋势项。设系统输出的测量值为1,2,...,kykL;，其一价趋势项可用式01kyhhk表示。试用最小二乘法确定0h和1h。解：令，01Tkyhhk，其中，01[1],[]TTkhh对于1,2,...,kL，将上式写成批量最小二乘格式LLY，其中11121LL，12LLyyYy应用LS法1ˆ()TTLLLLY122111(1)/2(1)/2(1)/2(1)/2(1)(21)/6LkTLLLLLkkkLkLLLLLLLLkLLLLLkk1(1)(21)/6(1)/21(1)/21(1)(1)/12TLLLLLLLLL11LkkTLLLkkyYky所以1011ˆ(1)(21)/6(1)/21ˆˆ(1)/21(1)(1)/12LkkLkkyhLLLLLLLhky1102(21)6ˆ(1)LLkkkkLykyhLL111126(1)ˆ(1)(1)LLkkkkkyLyhLLL则去除数据中的一价趋势项的操作：*01ˆˆ,1,2,...kkkkyyyyhhkkL■十、考虑一个SISO闭环系统如图所示，其中ku和ky分别为前向通道过程输入和输出量，}{ke为白噪声扰动序列，过程参数多项式)(1zA、)(1zB、)(1zC和已知的调节器参数多项式)(1zP、)(1zQ分别表示为：aannzazazA.1)(111bannnnzbzbbzBbb,.)(11011,.1)(111cnnnzczczCccppnnzpzpzP.1)(111qqnnzqzqqzQ.)(1101