hlayumi的个人空间 https://blog.eetop.cn/hslogic [收藏] [复制] [分享] [RSS]

空间首页 动态 记录 日志 相册 主题 分享 留言板 个人资料

日志

ML接收器

已有 1569 次阅读| 2019-4-27 23:33

%ML detector
j =sqrt(-1);

    s1 = 1/sqrt(2)*(1+j);
    s2 = 1/sqrt(2)*(1-j);
    s3 = 1/sqrt(2)*(-1-j);
    s4 = 1/sqrt(2)*(-1+j);
    s = [s1 s2 s3 s4];
    b = [0 0 0 1 1 1 1 0];
    bhat = zeros(1,2*length(r_SD));
    %N0 = 1;
    
    
    b_tr = RALEY_SETTINGS.b_tr_I;
    cc1 = ahat_SD* qpsk(b_tr);
    cc2 = ahat_RD* qpsk(b_tr);
    N0_SD = norm( r_SD(1:length(qpsk(b_tr)))-cc1 )^2/length(cc1);
    N0_RD = norm( r_RD(1:length(qpsk(b_tr)))-cc2 )^2/length(cc2);
%     N0_SD = 0.5;
%     N0_RD = 0.5;
    for i =1:length(r_SD)
                d1 = (1/(pi*N0_SD))*exp(-norm(r_SD(i)-ahat_SD*s(1))^2/ N0_SD)*...
                    ( (1-eps)^2*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(1))^2/N0_RD) +...
                    eps*(1-eps)*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(2))^2/N0_RD) +...
                    eps*(1-eps)*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(4))^2/N0_RD) +...
                    eps^2*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(3))^2/N0_RD) );
        
                d2 = (1/(pi*N0_SD))*exp(-norm(r_SD(i)-ahat_SD*s(2))^2/ N0_SD)*...
                    ( (1-eps)^2*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(2))^2/N0_RD) +...
                    eps*(1-eps)*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(1))^2/N0_RD) +...
                    eps*(1-eps)*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(3))^2/N0_RD) +...
                    eps^2*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(4))^2/N0_RD) );
        
                d3 = (1/(pi*N0_SD))*exp(-norm(r_SD(i)-ahat_SD*s(3))^2/ N0_SD)*...
                    ( (1-eps)^2*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(3))^2/N0_RD) +...
                    eps*(1-eps)*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(2))^2/N0_RD) +...
                    eps*(1-eps)*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(4))^2/N0_RD) +...
                    eps^2*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(1))^2/N0_RD) );
        
                d4 = (1/(pi*N0_SD))*exp(-norm(r_SD(i)-ahat_SD*s(4))^2/ N0_SD)*...
                    ( (1-eps)^2*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(4))^2/N0_RD) +...
                    eps*(1-eps)*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(1))^2/N0_RD) +...
                    eps*(1-eps)*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(3))^2/N0_RD) +...
                    eps^2*(1/(pi*N0_RD))*exp(-norm(r_RD(i)-ahat_RD*s(2))^2/N0_RD) );
        
        d = [d1 d2 d3 d4];
        
        bhat(2*i-1) = b(1);
        bhat(2*i)   = b(2);
        d_max= d(1);
        for j =2:4
            if  (d_max <d(j))
                d_max = d(j);
                bhat(2*i-1) = b(2*j-1);
                bhat(2*i)   = b(2*j);
            end
        end
        
    end

点赞

评论 (0 个评论)

facelist

您需要登录后才可以评论 登录 | 注册

  • 关注TA
  • 加好友
  • 联系TA
  • 0

    周排名
  • 0

    月排名
  • 0

    总排名
  • 0

    关注
  • 5

    粉丝
  • 0

    好友
  • 1

    获赞
  • 16

    评论
  • 5260

    访问数
关闭

站长推荐 上一条 /1 下一条

小黑屋| 关于我们| 联系我们| 在线咨询| 隐私声明| EETOP 创芯网
( 京ICP备:10050787号 京公网安备:11010502037710 )

GMT+8, 2024-3-29 05:50 , Processed in 0.015396 second(s), 6 queries , Gzip On, Redis On.

eetop公众号 创芯大讲堂 创芯人才网
返回顶部