AF协议DF协议对通信性能的影响分析

最新推荐文章于 2024-05-11 20:40:56 发布
最新推荐文章于 2024-05-11 20:40:56 发布
阅读量2.2k 收藏 9
点赞数 2
分类专栏： 板块1:通信与信号处理 文章标签： AF DF
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ccsss22/article/details/118666711
版权

clear all;
SINRO=0.1,T=10^21,k=1.38*10^(-23),refer=0.4,B=10^6,noise=10^(-6),N=50,Pp=10,rr=1.5*10^6;%基本元素的设定，此处的refer值得是第二时隙次用户发送自身信号和主用户信号的功率比，如ps1=（1-refer）ps，ps2=refer*ps
%设定各种信道的信道衰落
K=1,M=15;%规定信道矩阵的大小
V1=10^(-1)/sqrt(3.14/2);%瑞利信道的参数，注意这里不是平均值
V2=10^(-2)/sqrt(3.14/2);
H=15;%Hb表示整个系统程序循环的次数
W=ones(H,N);%每一横排表示一次循环的速率值，不同的纵坐标表示不同的信噪比下系统的速率值
W1=ones(H,N);%case 3 r df模式
W2=ones(H,N);%af模式
Pout=ones(H,N);%case 2 pout
Pout1=ones(H,N);%case 3 pout df模式
Pout2=ones(H,N);%af模式
Gss=random('rayleigh',10^(-1)/sqrt(3.14/2),K,M);%次用户发送端到次用户接收端的信道增益
Gpss=random('rayleigh',10^(-1)/sqrt(3.14/2),K,M);%主用户发送端到次用户发送端的信道增益
Gpp=random('rayleigh',10^(-2)/sqrt(3.14/2),K,M);%主用户发送端到主用户接收端的信道增益
Gsp=random('rayleigh',10^(-1)/sqrt(3.14/2),K,M);%次用户发送端到主用户接收端的信道增益
Gpssd=random('rayleigh',10^(-2)/sqrt(3.14/2),K,M);%主用户发送端到次用户接收端的信道增益
for j=1:1:H

%设定横坐标的范围
SNR=linspace(-10,30,N);                                          %横坐标为信噪比[-10,30]db,SNR=Gss*ps2/(B*noise)
snr=ones(size(SNR));                                                   %snr为信噪比的瓦特表示
for i=1:1:N                                                      %SNR到snr的转换
    snr(i)=10^(SNR(i)/10);
end%横坐标是次用户接收端的信噪比
%带入公式仿真
%第二种情形的系统速率和中断概率
a=snr.*B*noise;%a=Gss*ps2
 ps2=a./Gss(j);
for i=1:1:N
    if Gsp(j)*ps2(i)<=T*k*B 
        ps1(i)=ps2(i)/refer*(1-refer);
rs(i)=B*log2(1+a(i)/(B*noise+Gss(j)*ps1(i)));%rs=B*log2(1+a/b),a=Gss*ps2,b=B*noise+Gss*ps1
temp(i)=ps2(i)/(2^(rr/B)-1);
temp1(i)=-((B^2)*(noise^2))/(2*((temp(i)-ps1(i))^2)*(V1^2));
pout(i)=1-exp(temp1(i));
    else
     p2o(i)=T*k*B/Gsp(j);%次用户发送自身信号的边界功率
     p1o(i)=p2o(i)/refer*(1-refer);
 rs(i)=B*log2(1+Gss(j)*p2o(i)/(B*noise+Gss(j)*p1o(i)));   
 temp(i)=p2o(i)/(2^(rr/B)-1);
temp1(i)=-((B^2)*(noise^2))/(2*((temp(i)-p1o(i))^2)*(V1^2));
pout(i)=1-exp(temp1(i));

    end
end
for i=1:1:N
W(j,i)=rs(i);
Pout(j,i)=pout(i);
end
%rs2=ones(size(rs));
%rs2=rs;%把情形二下的速率值付给rs2

 %第三种情形的系统速率
 a=snr.*B*noise;%a=Gss*ps2
 ps2=a./Gss(j);
 for i=1:1:N
     if Gsp(j)*ps2(i)<=T*k*B %满足对主用户的影响小于特定值
          ps1(i)=ps2(i)/refer*(1-refer);
         if a(i)/(Pp*Gpssd(j)+B*noise+Gss(j)*ps1(i))>=SINRO%满足次用户的最小通信要求
             rs3(i)=B*log2(1+a(i)/(Pp*Gpssd(j)+B*noise+Gss(j)*ps1(i)));%只有同时满足两个条件才能用这个表达式表达。
             c=B*log2(1+Gpss(j)*Pp/(B*noise));
             d1=1+Gpp(j)*Pp/(B*noise);
             d2(i)=Gsp(j)*ps1(i)/(B*noise+Gsp(j)*ps2(i));
             e= ps1(i)*Pp*Gpss(j)*Gsp(j)/(B*noise+Pp*Gpss(j));
             e1=B*noise*(1+ps1(i)*Gsp(j)/(B*noise+Pp*Gpss(j)))+ps2(i)*Gsp(j);
             d3(i)=e/e1;
           
             d(i)=B*log2(d1+d2(i));%df模式
             rp(i)=min(c,d(i));
             r(i)=rs3(i)+rp(i);
              da(i)=B*log2(d1+d3(i));%af模式
               rpa(i)=min(c,da(i));
             ra(i)=rs3(i)+rpa(i);
             
          %   首先得是df模式主用户的中断概率
             if c<d(i)%如果较小的是前者，则主用户的中断概率是一个常数如下
           temp4=1;
           temp2=(2^(rr/B)-temp4)^2;
           temp3=-(temp2*(B^2)*(noise^2))/(2*(Pp^2)*(V1^2));
             else
          temp4=1+Gsp(j)*ps1(i)/(B*noise+Gsp(j)*ps2(i));
          temp2=(2^(rr/B)-temp4)^2;
           temp3=-(temp2*(B^2)*(noise^2))/(2*(Pp^2)*(V2^2));
             end
           %df主用户的中断概率
        
           poutp(i)=1-exp(temp3);
           %然后是af模式主用户的中断概率
            if c<da(i)%如果较小的是前者，则主用户的中断概率是一个常数如下
           temp41=1;
            temp21=(2^(rr/B)-temp41)^2;
           temp31=-(temp21*(B^2)*(noise^2))/(2*(Pp^2)*(V1^2));
             else
          temp41=1+d3(i);
           temp21=(2^(rr/B)-temp41)^2;
           temp31=-(temp21*(B^2)*(noise^2))/(2*(Pp^2)*(V2^2));
            end
               %af主用户的中断概率
        
           poutpa(i)=1-exp(temp31);
            
            %接下来是次用户的中断概率
            temp(i)=ps2(i)/(2^(rr/B)-1);
            temp1(i)=-((Pp*Gpssd(j)+B*noise)^2)/(2*((temp(i)-ps1(i))^2)*(V1^2));
            pout(i)=1-exp(temp1(i));

            %系统的中断概率
           poutsys(i)=pout(i)*poutp(i);
           poutsysa(i)=pout(i)*poutpa(i);
         else%不满足次用户的最小通信要求，即ps2(i)的值太小了
             rs3(i)=B*log2(1+SINRO);
             refer1=1/refer-1;
             ps22=B*noise*SINRO/(Gss(j)*(1-refer1*SINRO));%边界功率值
             ps11=ps22/refer*(1-refer);
             c=B*log2(1+Gpss(j)*Pp/(B*noise));
             d1=1+Gpp(j)*Pp/(B*noise);
             d2=Gsp(j)*ps11/(B*noise+Gsp(j)*ps22);
              
              
              e= ps11*Pp*Gpss(j)*Gsp(j)/(B*noise+Pp*Gpss(j));
             e1=B*noise*(1+ps11*Gsp(j)/(B*noise+Pp*Gpss(j)))+ps22*Gsp(j);
             d3(i)=e/e1;
             %df模式
             d=B*log2(d1+d2);
             rp(i)=min(c,d);
             r(i)=rs3(i)+rp(i);
             %af模式
             da(i)=B*log2(d1+d3(i));
               rpa(i)=min(c,da(i));
             ra(i)=rs3(i)+rpa(i);
             
             
             %df模式主用户的中断概率
             if c<d%如果较小的是前者，则主用户的中断概率是一个常数如下
            temp4=1;
             temp2=(2^(rr/B)-temp4)^2;
           temp3=-(temp2*(B^2)*(noise^2))/(2*(Pp^2)*(V1^2));
             else
            temp4=1+Gsp(j)*ps11/(B*noise+Gsp(j)*ps22);
             temp2=(2^(rr/B)-temp4)^2;
           temp3=-(temp2*(B^2)*(noise^2))/(2*(Pp^2)*(V2^2));
             end
           %主用户的中断概率
       
           poutp(i)=1-exp(temp3);
           
            %然后是af模式主用户的中断概率
            if c<da(i)%如果较小的是前者，则主用户的中断概率是一个常数如下
           temp41=1;
            temp21=(2^(rr/B)-temp41)^2;
           temp31=-(temp21*(B^2)*(noise^2))/(2*(Pp^2)*(V1^2));
             else
          temp41=1+d3(i);
           temp21=(2^(rr/B)-temp41)^2;
           temp31=-(temp21*(B^2)*(noise^2))/(2*(Pp^2)*(V2^2));
            end
               %af主用户的中断概率
       
           poutpa(i)=1-exp(temp31);
   
            %接下来是次用户的中断概率
            temp(i)=ps22/(2^(rr/B)-1);
            temp1(i)=-((Pp*Gpssd(j)+B*noise)^2)/(2*((temp(i)-ps11)^2)*(V1^2));
            pout(i)=1-exp(temp1(i));

            %系统的中断概率
            poutsys(i)=pout(i)*poutp(i);
            poutsysa(i)=pout(i)*poutpa(i); 
             
         end
     else%不满足对主用户的影响小于特定值，即ps2（i）的值太大了
         p2oo(i)=T*k*B/Gsp(j);%次用户发送自身信号的边界功率
         p1oo(i)=p2oo(i)/refer*(1-refer);
         rs3(i)=B*log2(1+Gss(j)*p2oo(i)/(Pp*Gpssd(j)+B*noise+Gss(j)*p1oo(i)));   
         c=B*log2(1+Gpss(j)*Pp/(B*noise));
         d1=1+Gpp(j)*Pp/(B*noise);
             d2(i)=Gsp(j)*p1oo(i)/(B*noise+Gsp(j)*p2oo(i));
             e= p1oo(i)*Pp*Gpss(j)*Gsp(j)/(B*noise+Pp*Gpss(j));
             e1=B*noise*(1+p1oo(i)*Gsp(j)/(B*noise+Pp*Gpss(j)))+p2oo(i)*Gsp(j);
             d3(i)=e/e1;
             %df模式
             d(i)=B*log2(d1+d2(i));
             rp(i)=min(c,d(i));
             r(i)=rs3(i)+rp(i);
              %af模式
             da(i)=B*log2(d1+d3(i));
             rpa(i)=min(c,da(i));
             ra(i)=rs3(i)+rpa(i);
             
             %df模式主用户中断概率
           
             if c<d(i)%如果较小的是前者，则主用户的中断概率是一个常数如下
            temp4=1;
            temp2=(2^(rr/B)-temp4)^2;
            temp3=-(temp2*(B^2)*(noise^2))/(2*(Pp^2)*(V1^2));
             else
             temp4=1+Gsp(j)*p1oo(i)/(B*noise+Gsp(j)*p2oo(i));
             temp2=(2^(rr/B)-temp4)^2;
            temp3=-(temp2*(B^2)*(noise^2))/(2*(Pp^2)*(V2^2));
             end
            %主用户的中断概率 
           
            poutp(i)=1-exp(temp3);
             %af模式主用户中断概率
           
             if c<da(i)%如果较小的是前者，则主用户的中断概率是一个常数如下
            temp41=1;
             temp21=(2^(rr/B)-temp41)^2;
            temp31=-(temp21*(B^2)*(noise^2))/(2*(Pp^2)*(V1^2));
             else
             temp41=1+d3(i);
              temp21=(2^(rr/B)-temp41)^2;
            temp31=-(temp21*(B^2)*(noise^2))/(2*(Pp^2)*(V2^2));
             end
            %主用户的中断概率 
           
            poutpa(i)=1-exp(temp31);
            
            %接下来是次用户的中断概率
            temp(i)=p2oo(i)/(2^(rr/B)-1);
            temp1(i)=-((Pp*Gpssd(j)+B*noise)^2)/(2*((temp(i)-p1oo(i))^2)*(V1^2));
            pout(i)=1-exp(temp1(i));

           %系统的中断概率
            poutsys(i)=pout(i)*poutp(i);
             poutsysa(i)=pout(i)*poutpa(i); 
     end
 end
 for i=1:1:N
W1(j,i)=r(i);%df模式下系统的容量
W2(j,i)=ra(i);%af模式下系统的容量
Pout1(j,i)=poutsys(i);%df模式下系统的中断概率
Pout2(j,i)=poutsysa(i);%af模式下系统的中断概率
 end

end%大循环的结束

for i=1:1:N
    r2(i)=mean(W(1:1:H,i));%H次循环的平均值
    r3(i)=mean(W1(1:1:H,i));
     r3a(i)=mean(W2(1:1:H,i));
    Pou(i)=mean(Pout(1:1:H,i));
    Pou1(i)=mean(Pout1(1:1:H,i));
    Pou1a(i)=mean(Pout2(1:1:H,i));
end
plot(SNR,r2,'o',SNR,r3,'.',SNR,r3a,':');


xlabel('SNR/db');
ylabel('r/(bit/second)');
legend('case2 存/不存','case3 存/存,DF模式','case3 存/存,AF模式');
figure(2)
plot(SNR,Pou,'o',SNR,Pou1,'.',SNR,Pou1a,':');
xlabel('SNR/db');
ylabel('Pout');
legend('case2 存/不存','case3 存/存,DF模式','case3 存/存,AF模式');

D100

fpga和matlab
关注
  • 2
    点赞
  • 踩
  • 9
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
专栏目录
基于MATLAB的中继系统在瑞利信道下的误码率仿真分析，matlab2021a仿真
04-29
基于MATLAB的中继系统在瑞利信道下的误码率仿真分析，matlab2021a仿真 %设定各种信道的信道衰落 K=1,M=15;%规定信道矩阵的大小 V1=10^(-1)/sqrt(3.14/2);%瑞利信道的参数，注意这里不是平均值 V2=10^(-2)/sqrt(3.14/2); H=15;%Hb表示整个系统程序循环的次数 W=ones(H,N);%每一横排表示一次循环的速率值，不同的纵坐标表示不同的信噪比下系统的速率值 W1=ones(H,N);%case 3 r df模式 W2=ones(H,N);%af模式 Pout=ones(H,N);%case 2 pout Pout1=ones(H,N);%case 3 pout df模式 Pout2=ones(H,N);%af模式 Gss=random('rayleigh',10^(-1)/sqrt(3.14/2),K,M);%次用户发送端到次用户接收端的信道增益 Gpss=random('rayleigh',10^(-1)/sqrt(3.14/2),K,M);%主用户发送端到次用户发送端的信道增益
瑞利信道误码率
03-01
瑞利信道误码率 的matlab仿真，可实现BPSK条件下的仿真。
中继误码率 matlab,关于误码率的问题 急！！！！！
weixin_42367233的博客
03-23 133
该楼层疑似违规已被系统折叠隐藏此楼查看此楼求助各位大佬要求是2ask调制 通过滚降系数为0,0.5,1的升余弦滚降滤波器，信道加入1-15db的高斯白噪声，相干解调，匹配滤波后计算出误码率与信噪比的关系图与理论值拟合。现在问题是计算出的误码率过高无法拟合，请教各位大佬应该怎么修改。明天就要交论文了 曲线一直拟合不上好急%%---baseband signal parameters---%%D_R...
瑞利衰落信道的误码率分析
06-30
了解瑞利衰落信道的基本特征，使用 Matlab 进行 QPSK 调制与解调，使用 Matlab 分析信道的误码率，信道编解码的基本原理及其 Matlab 实现。
瑞利信道中 QPSK 的误码率：瑞利信道中 QPSK 的误码率-matlab开发
06-01
瑞利信道中 QPSK 的 BER
协作通信-AF、DF的MATLAB仿真_MATLAB仿真AF/DF协议_协议仿真_
09-29
DF协议提供了错误检测和纠正能力，因此在信噪比较低的环境下，其性能通常优于AF协议，但其缺点是增加了处理复杂度和延迟。 **MATLAB仿真：** MATLAB是一种强大的数学计算和数据分析工具，广泛应用于通信系统的建模...
协作通信AF与DF协议性能分析
05-16
协作通信AF与DF协议性能分析，通过matlab仿真分别绘出误码率与误比特率与信噪比的关系，得出AF与DF协作的性能。
协作通信-AF、DF的MATLAB仿真_MATLAB仿真AFDF协议_协议仿真_源码.zip
10-18
总的来说，这个MATLAB仿真项目为研究和教学提供了实用的工具，帮助通信专业人员和学生更好地理解和评估协作通信中的AF和DF协议。通过实际操作和调试代码，学习者能够深入理解无线通信理论，并提高问题解决和编程能力...
04016535_龙政兴_AF-DF协作中功率分配对协作网络性能的影响1
08-08
实验总结强调了功率分配对协作通信性能的重要性，以及DF协议在降低误码率方面的优势。通过这次实验，学生不仅深化了对AF-DF协作技术原理的理解，还通过编程实践增强了对协作通信实现的认识。实验结果为优化功率分配...
16QAM在瑞利信道下的误码率
05-03
16QAM在瑞利信道下的误码率
BPSK瑞利和高斯信道衰落误码率
07-16
本资源有一个matlab程序段，是仿真BPSK分别在高斯噪声和瑞利衰落下的误码率，产生图形对仿真值和理论值进行比较
bpsk在瑞利信道下的误码率与信噪比的关系
12-15
bpsk在瑞利信道下的误码率与信噪比的关系，通信相关专业用处较大。
matlab仿真qpsk在高斯信道和瑞利衰落信道下的误码率
01-12
本资源有两个matlab程序段，都是仿真QPSK分别在高斯噪声和瑞利衰落下的误码率，产生图形对仿真值和理论值进行比较。
DF1协议详解(PDF)
01-05
ab PLC的DF1协议详解。对于希望了解DF1协议或相关编程解析等工作来说十分具有参考价值。
matlab信号经过瑞利,Matlab关于 BPSK信号通过瑞利信道的误码性能分析的仿真实验....
weixin_29056781的博客
03-20 1902
主函数。SNRindB1=0:2:12;%定义信噪比序列SNRindB2=0:2:12;%采用dB作为单位for i=1:length(SNRindB1)s_e(i)=ray(SNRindB1(i));%计算误码率endfor i=1:length(SNRindB2)SNR=exp(SNRindB2(i)*log(10)/10);%转化信噪比为数值表示。。theo_err_prb(i)= 1/2*...
2BPSK调制下瑞利信道三种分集合并技术的误码率分析
wangwttwtt的博客
11-18 311
【代码】2BPSK调制下瑞利信道三种分集合并技术的误码率分析。
瑞利信道下BPSK的误码率
cz20010101的博客
04-09 2276
瑞利信道条件下BPSK的误码率 通信系统 sss为随机产生的信号，为0或者1 调制后的信号为xxx,0映射为−Eb-\sqrt{E_b}−Eb​​,1映射为+Eb+\sqrt{E_b}+Eb​​ 高斯白噪声为nnn,满足均值为0，方差为σ2=N02\sigma^2=\frac{N0}{2}σ2=2N0​的高斯分布 收到的信号y=hx+ny=hx+ny=hx+n,h用于描述信道，这里的信道是瑞利信道，固h有实部和虚部，且实部和虚部是独立同分布的高斯变量，均值为0，方差为12\frac{1}{2}21​(这
QPSK信号在高斯信道、瑞利信道、Ricean信道下的误码率仿真（matlab）
如切如磋，如琢如磨
06-23 5190
这里只展示最终的结果，需要源码的联系上面的二维码免费获取
砸饭碗系列一、BPSK高斯、瑞利、莱斯信道误码率对比（手把手推导瑞利信道理论误码率公式）
最新发布
王翊的博客
05-11 2290
做这个的原因是，想学习瑞利跟莱斯信道，结果发现csdn 上拥有这个曲线的博主，居然要收80多块钱！但是美中不足的是缺少莱斯信道的理论误码率曲线。先说结论，由于莱斯信道是在瑞利信道基础上加入一条直射路径，因此功率分为两个部分，直射路径功率P1以及瑞利多径的功率P2，用莱斯因子K进行衡量,且K趋于无穷代表没有衰落，K=0等与瑞利衰落。本文是做了BPSK高斯信道理论以及仿真的误码率曲线、瑞利信道理论以及仿真的误码率曲线、莱斯信道仿真的误码率曲线。（注意没有莱斯信道理论的误码率曲线公式，为什么呢，因为我没做出来！
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

fpga和matlab

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值