24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 540  |  回复: 2

purplelotus

新虫 (初入文坛)

[交流] 相关对数阴影下合作频谱感知已有1人参与

最近读了“Asymptotic performance of collaborative spectrum sensing under correlated log-normal shadowing”这篇文章,并进行了仿真,但仿真结果与理论结果不符,不知道错在哪里,请大家指教。



理论值:

uR=1;
delta=-5.19;
sigma=2.3;
alpha=0.1;
D1=100;
beta=0.05;
for n=1:100
    rho=exp(-alpha*D1/(n-1));
    pm(n)=1-qfunc(delta/sigma*(((1-rho)*n+2*rho)/(1+rho))^0.5+qfuncinv(beta));
    rho=exp(-alpha*D1*2/(n-1));   
    pm1(n)=1-qfunc(delta/sigma*(((1-rho)*n+2*rho)/(1+rho))^0.5+qfuncinv(beta));
end
hold on;
semilogy([1:100],pm);
semilogy([1:100],pm1);
b=1-qfunc(delta/sigma*(alpha*100/2+1)^0.5+qfuncinv(beta));
plot([1:100],b);
b1=1-qfunc(delta/sigma*(alpha*200/2+1)^0.5+qfuncinv(beta));
plot([1:100],b1);



仿真值:

uR=1;
delta=-5.19;
sigma=2.3;
alpha=0.1;
D=100;
beta=0.01;
for n=2:20
    pm(n)=0;
    for k=1:500
        rho=exp(-alpha*D/(n-1));
        y=sigma*randn(n-1,10000)+uR+delta; %
        r=[0:1:n-2];
        r1=rho.^r;
        sigma_b=toeplitz(r1);% correlation matrix
        oneT=ones(1,n-1);
        one=oneT';
        temp=oneT*sigma_b^(-1)*one;
        lemda=uR-sigma/temp^0.5*qfuncinv(beta);%decision threshold
        C = chol(sigma_b,'lower');
        y1=C*y; % correlated Log-normal SNRs
        lemda_b=oneT*sigma_b^(-1)*y1/temp;
        lemda_a=ones(1,10000)*lemda;
        p=ge(lemda_b,lemda_a);
        num=sum(p);
        p=num/10000;
        pm(n)=pm(n)+p;
    end
    pm(n)=pm(n)/500;
end
semilogy([1:100],pm);
hold on;





仿真结果相关对数阴影下合作频谱感知
1.png
回复此楼

» 本帖附件资源列表

» 猜你喜欢

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

sunyuanxin

禁虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

2楼2018-09-18 15:51:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

purplelotus

新虫 (初入文坛)

不是,红线是D=100d的仿真值,仿真值与理论值不一致
3楼2018-09-18 18:49:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 purplelotus 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 普通院校药学硕士,做合成的,感觉找不到工作 +8 pom戴墨镜 2024-04-24 15/750 2024-04-25 20:50 by pom戴墨镜
[教师之家] 刚博士毕业,入职普通一本高校 +18 保持肌肤 2024-04-20 19/950 2024-04-25 18:30 by lyfbangong
[考博] 25年博士申请 +5 Changzixuan 2024-04-25 5/250 2024-04-25 17:37 by 那年的旧城
[基金申请] 基金开始函评了吗? +11 wych1103 2024-04-25 11/550 2024-04-25 16:54 by 未岸亡灵
[论文投稿] with editor 两个月了,什么原因? +6 yiersan9 2024-04-24 8/400 2024-04-25 15:45 by topedit
[硕博家园] 聊天 +11 暮色恋伊人 2024-04-22 12/600 2024-04-25 13:53 by UCTS
[论文投稿] 研二光催化6月底四篇二区什么水平 5+4 wjtab 2024-04-22 11/550 2024-04-25 13:20 by 椰蛋树的曾经
[基金申请] 两类问题算是白选了~ +7 jurkat.1640 2024-04-23 11/550 2024-04-25 13:06 by 立诛杀曹无伤
[基金申请] 国社科项目,你们学校都限额申报吗? +7 屡战屡败 2024-04-21 10/500 2024-04-25 12:10 by 屡战屡败
[博后之家] 南京大学-广州大学联合招聘博士后 欢迎广大优秀人才!!! +4 黑魔变身啾 2024-04-20 12/600 2024-04-25 11:18 by dodonaomi
[论文投稿] Chemical Engineering Journal投稿3周了,一直显示With editor状态。这是送审了吗? 10+4 yifeng11 2024-04-20 10/500 2024-04-25 02:14 by yifeng11
[基金申请] 企业博后是否能申请CSC博士后项目?谢谢 66+3 龙包子 2024-04-19 6/300 2024-04-24 21:59 by 3115321
[访问学者] CSC的访问学者申请,没有个评审意见,也不知道怎么改,还有必要申请吗 +4 flyingship 2024-04-20 4/200 2024-04-24 21:09 by 59038mute
[硕博家园] 考研,求职还是考编? +13 xizj 2024-04-21 21/1050 2024-04-24 17:23 by jxdic
[教师之家] 大家访学都是怎么找的啊? +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
[论文投稿] 期刊推荐 20+4 木颜尘ip 2024-04-22 7/350 2024-04-24 10:06 by bobvan
[教师之家] 找工作过程中碰到好多人问南京某新换校长211的政策,规定如下 +15 ZHONGWU_U 2024-04-21 17/850 2024-04-24 05:12 by juta256
[公派出国] 東京大学 应用化学专业柳田研究室 招收修士/博士生 +5 蝶梦123 2024-04-19 5/250 2024-04-23 22:40 by 90472jvvh
[论文投稿] 编辑是选国外的好还是国内的好。 +8 lizhengke06 2024-04-20 8/400 2024-04-22 08:58 by cuiyunjian
[考研] 312求调剂 +3 Leroic000 2024-04-19 3/150 2024-04-20 18:19 by xinaishan
信息提示
请填处理意见