pointlmi.m

线性时变系统控制器设计的工具包

function [vgam,vclp,vctrl] = pointlmi(vlpv,nmeas,nctrl)

% [vgam,vX,vY] = pointlmi(vlpv,nmeas,nctrl)

viv = getiv(vlpv);
npts = length(viv);
vgam = zeros(npts,1);
nx = xnum(xtracti(vlpv,1,1));
vX = zeros(npts*nx,nx);
vY = zeros(npts*nx,nx);

for k = 1:npts
 disp(['Point ' int2str(k) ' of ' int2str(npts) ':']);
 [gam,xmat,ymat] = vlpvsol(xtracti(vlpv,k),nmeas,nctrl);
 vgam(k) = gam;
 vX((k-1)*nx+1:k*nx,1:nx) = xmat;
 vY((k-1)*nx+1:k*nx,1:nx) = ymat;
end

vgam = vpck(vgam,viv);
vX = vpck(vX,viv);
vY = vpck(vY,viv);