objectif.m

来自「Tensor MIMO system simulation using MATL」· M 代码 · 共 21 行

function s=objectif(X)
% function s=objectif(X)
% P.Comon, June 16, 1998

global VECTEURSNOYAU
global POLYNOMECIRCULAIRE

PC=POLYNOMECIRCULAIRE;
v1=VECTEURSNOYAU(:,1);v2=VECTEURSNOYAU(:,2);
theta=X(1);phi=X(2);ii=sqrt(-1);n=length(v1); % n=r+1=4 generiquement

  u=v1*cos(theta)+v2*sin(theta)*exp(ii*phi);
  % construction des 3 formes lineaires
  if abs(u(1))<2*eps, q=roots(u(2:n));L=[1 q(:).';0 ones(1,n-2)];
  else q=roots(u); L=[q(:).';ones(1,n-1)]; end;
  P1=carr2PR2(mod2LC2PR(L(:,1)));
  P2=carr2PR2(mod2LC2PR(L(:,2)));
  P3=carr2PR2(mod2LC2PR(L(:,3)));
  pppPC=[P1 P2 P3 PC];
  s=min(svd(pppPC));