%广义最小方差控制(显示控制)考虑如下系统：式中x(k)为方差为0.1的白噪声。取，期望输出yr(k)为幅值为10的方波信号。clear all;close all;a=1 -1.7 0.7;b=1 2;c=1 0.2;d=4;na=length(a)-1;nb=length(b)-1;nc=length(c)-1;nf=nb+d-1;ng=na-1;P=1;R=1;Q=2; %加权多项式np=length(P)-1;nr=length(R)-1;nq=length(Q)-1;L=400;uk=zeros(d+nb,1);yk=zeros(na,1);yrk=zeros(nc,1);xik=zeros(nc,1);yr=10*ones(L/4,1);-ones(L/4,1);ones(L/4,1);-ones(L/4+d,1);xi=sqrt(0.1)*randn(L,1);e,f,g=singlediophantine(a,b,c,d);CQ=conv(c,Q);FP
