load Nominalmatrices.mat;
i=4; %number of measurments for combination 
[B,sset,ops,ctime,flag]=pb3wc(G1,F,Wd,Wn,Juu,Jud,i,tlimit,1);
% meas{i,1} = sset
% Gy = G1(meas{i},:);
Gy=G1(sset,:);
Y1=[F*Wd Wn];
% Y = Y1(meas{i},:);
Y=Y1(sset,:);
H1=inv(Y*Y')*Gy*inv(Gy'*inv(Y*Y')*Gy)*sqrtm(Juu);
B*1e8/(365*24)
H = H1'




