%%K从1到15的辅助信号都可以计算
load('Data0.mat');
k=1+Data0;
Data0=k;
save('Data0.mat','Data0');
x=fminbnd('F',0,1);
H=Data{k};
P0=Data5{k};
P1=Data6{k};
A=H'*(1/(1-x)*P0+1/x*P1)^-1*H;
a=eig(A);%特征值矩阵
[m,p]=max(a);%最大特征值
Data10{k}=m;
Data9{k}=x;
b(:,k)=Data10{k};
[n,y]=eig(A);%求矩阵的特征值和特征向量,n为特征向量矩阵,y为特征值矩阵
eigenvalue=diag(y);%取矩阵y的对角线元素
lamda=max(eigenvalue);%求最大特征值
for i=1:length(A)%求最大特征值对应的序数
if lamda==eigenvalue(i)
break;
end
end
y_lamda=n(:,i);%求矩阵最大特征值对应的特征向量
w=sqrt(1/b(:,k))*y_lamda/norm(y_lamda);
Data12{k}=w;
W=w'*w;
Data11{k}=W;