24小时热门版块排行榜    

Znn3bq.jpeg
查看: 835  |  回复: 1

美沙黄路夜

新虫 (正式写手)

[求助] 子空间跟踪算法性能测度---最大主角比较--详见帖子

江湖救急啊!
两个子空间之间的主角是两个子空间之间差异度的度量。主角越小,两个子空间之间的差异度越小。当两者一致时,主角为零。子空间跟踪算法的最大主角指的是子空间跟踪算法估计的自相关矩阵 的信号子空间与实际(基于准确特征值分解得到)的信号子空间之间的最大主角。

%% 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
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyuanxin

禁虫 (著名写手)

本帖内容被屏蔽

2楼2015-05-30 00:38:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 美沙黄路夜 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 河北省自然科学基金 +5 Peterchao 2026-05-18 7/350 2026-05-23 16:57 by articlefan
[基金申请] 揭秘青基评审内幕:几个A才能顺利中标 +3 国自然国社科中 2026-05-23 4/200 2026-05-23 15:37 by 2000zf36392
[论文投稿] 投MDPI的三个期刊都被桌拒 +4 木虫研究牲 2026-05-17 5/250 2026-05-23 14:44 by 期刊投稿求助
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 6/300 2026-05-23 14:32 by tfang
[基金申请] 青B发送上会通知了吗 +5 chemBioBro 2026-05-22 7/350 2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +3 zju2000 2026-05-21 4/200 2026-05-23 08:33 by xli1984
[基金申请] 今年审到国自然15份,谈谈感受 +27 国自然国社科中 2026-05-17 28/1400 2026-05-22 16:25 by hunter无悔
[论文投稿] 投稿求助,期刊 +4 希冀,有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀,有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +14 1234567wang 2026-05-17 16/800 2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[基金申请] 评审有感 +15 popular289 2026-05-18 26/1300 2026-05-21 10:35 by 西葫芦炒鸡蛋
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +4 光伏爱好者 2026-05-17 6/300 2026-05-21 01:37 by 17520222914
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
请填处理意见