clear all;
clc;
pfa = 0.05; %虚警概率
PH1=0.7;
PH0=1-PH1;
SNRp=-15;
SNRs=20;%信噪比
T=100e-3;
SNRP=10.^(SNRp/10);
SNRS=10.^(SNRs/10);
C0=log2(1+SNRS);
fs=100e+3;
ts=1e-3:0.5e-3:100e-3;
M=3;
N=ts*fs/M
PathGains=(randn(1,M)+j*randn(1,M))
for iii=1:length(ts)
for kk=1:M
D=(N(iii)+SNRs*(abs(PathGains(kk))).^2).^0.5.*eye(kk,kk)
SNR=(ones(1,kk)*SNRs.*(abs(PathGains(1:kk))).^2)'
Dinv=inv(D)
cooper=Dinv*SNR*SNR'*Dinv
[v ,d]=eig(cooper)
[a,b]=max(max(d))
eigen=Dinv*v(:,b,:)
Q=norm(abs(eigen))