½ºþ¾È¼±°¡£¡
Á½¸ö×Ó¿Õ¼äÖ®¼äµÄÖ÷½ÇÊÇÁ½¸ö×Ó¿Õ¼äÖ®¼ä²îÒì¶ÈµÄ¶ÈÁ¿¡£Ö÷½ÇԽС£¬Á½¸ö×Ó¿Õ¼äÖ®¼äµÄ²îÒì¶ÈԽС¡£µ±Á½ÕßÒ»ÖÂʱ£¬Ö÷½ÇΪÁã¡£×Ó¿Õ¼ä¸ú×ÙËã·¨µÄ×î´óÖ÷½ÇÖ¸µÄÊÇ×Ó¿Õ¼ä¸ú×ÙËã·¨¹À¼ÆµÄ×ÔÏà¹Ø¾ØÕó µÄÐźÅ×Ó¿Õ¼äÓëʵ¼Ê(»ùÓÚ×¼È·ÌØÕ÷Öµ·Ö½âµÃµ½£©µÄÐźÅ×Ó¿Õ¼äÖ®¼äµÄ×î´óÖ÷½Ç¡£
%% FSCAPI-FAPI×î´óÖ÷½Ç±È½Ï
clc,clear all
S=100;Ra=0;B=1000;L=0;SNR=20;
%Éú³ÉgoldseqÐòÁÐ
fbconnection=[0 1 0 0 1];
mseq=m_sequence(fbconnection);
fbconnection=[0 1 0 0 1];
mseq=m_sequence(fbconnection);
fbconnection1=[0 0 1 0 1];
fbconnection2=[0 1 1 1 1];
goldseq=gold_seq(fbconnection1,fbconnection2);
ind1=find(mseq==0);
mseq(ind1)=-1;
ind2=find(goldseq==0);
goldseq(ind2)=-1;
b=rand(6,10001); %Óû§Êý¾Ýb²úÉú
ind1=find(b<=0.5);
b(ind1)=-1;
ind2=find(b>0.5);
b(ind2)=1;
%GoldÂë¹éÒ»»¯
goldseq=1/sqrt(31)*goldseq;
%ÆÚÍûÓû§À©ÆµÐòÁÐ
s1=goldseq(1,: )';
s2=goldseq(2,: )';
s3=goldseq(3,: )';
s4=goldseq(4,: )';
s5=goldseq(5,: )';
s6=goldseq(6,: )';
A2=10;A3=10;A4=10;A5=10;A6=10;
N=length(s1);C=eye(N);
r1=goldseq(1,: )'*b(1,: );
r2=goldseq(2,: )'*b(2,: );
r3=goldseq(3,: )'*b(3,: );
r4=goldseq(4,: )'*b(4,: );
r5=goldseq(5,: )'*b(5,: );
r6=goldseq(6,: )'*b(6,: );
r=r1+r2*A2+r3*A3+r4*A4+r5*A5+r6*A6;
r=awgn(r,SNR);
% ÕâÀïµÄdsubspaceµ½µ×ÊÇʲôº¯Êý£¿ÇëÖªµÀµÄͬѧ°ïæ±àд£»ÓÐÎÊÌâÏÂÃæ»Ø¸´%ÎÒ£¬¸Ðл°¡£¡½ð±Ò²»ÊÇÎÊÌâ
[Us1,Un1,Ds1,Dn1,ds1,dn1]=dsubspace(r(:,1:B));
X=eye(31);Us2=X(:,1:6);Us3=X(:,1:6);[N,K]=size(Us2);beta=0.995;Z2=eye(K);Z3=eye(K);
for m=1:S
%[r1,r2,r3,r4,r5,r6,r7,r8,s1,s2,s3,s4,s5,s6,s7,s8,b1]=cgeneratingr(B,L);
%A1=1;A2=10;r=A1*r1+A2*(r2+r3+r4+r5+r6);
for t=L+1:L+B
x=r(:,t);
[Us2,Z2]=SFAPI_tracker(Us2,Z2,x);
[Us3,Z3]=FAPI_tracker(Us3,Z3,x);
[Q1,q1]=qr(Us1,0);
[Q2,q2]=qr(Us2,0);
[Q3,q3]=qr(Us3,0);
D1=Q1'*Q2;
[U1,E1,V1]=svd(D1,0);
D2=Q1'*Q3;
[U2,E2,V2]=svd(D2,0);
MA1(t-L)=acosd(E1(K,K));
MA2(t-L)=acosd(E2(K,K));
ratio(t-L)=20*log(MA1(t-L)/MA2(t-L));
end
Ra=Ra+ratio;
end
Ra=Ra/S;
t=1:B;
plot(t,Ra(t),'r');
xlabel('µü´ú´ÎÊý/´Î');ylabel('×î´óÖ÷½Ç±È/dB');
legend({'FSCAPI-FAPI'},'Location','Best');
![×Ó¿Õ¼ä¸ú×ÙËã·¨ÐÔÄܲâ¶È---×î´óÖ÷½Ç±È½Ï--Ïê¼ûÌû×Ó]()
11111.JPG |