您现在的位置:首页 >> ★免费资源 >> 源码下载 >> 内容

MATLAB代做|FPGA代做|Python代做--没有协作分支的频谱效率数值计算

时间:2019-1-5 15:10:18 点击:

  核心提示:MATLAB代做|FPGA代做|Python代做--没有协作分支的频谱效率数值计算...
%%%%%%%%%%%%%%%%%%%%%%%%%没有协作分支的频谱效率数值计算%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%       
%%%%%%%%%%%%%%%%%%%%%%%%%没有协作分支的频谱效率数值计算%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

    sen0=[];
    Rn=[];
    an=[67.7328,73.8279,58.7332,55.9137,50.0552];    %%拟合函数参数
    gn=[0.9819,0.4945,0.1641,0.0989,0.0381];         %%拟合函数参数
    rn=[ 4.2935, 8.6986,24.8142,40.6799,102.5723,1000];    %%预设信噪比阈值rn=rp,线性比值%
    p0=1080*10^(-6);                                %%%%%%目标分组错误概率%%%%%%%  
    for snrdB=0:30;
        syms f t; 
        f=cooperation0pdfint(snrdB);
        I=int(f,t,0,rn(6));
        e=eval(I);
        while e<0.9999          %%%%求信噪比转换阈值rn(6),令截断误差小于0.0001
             rn(6)=rn(6)+100;
             I=int(f,t,0,rn(6));
             e=eval(I);
        end
         Rn(snrdB+1,:)=rn
   
          
    for n=5:-1:1           %%%%求信噪比转换阈值rn(1)-rn(5),令截断误差小于0.0001
        syms f h t;  
        f=cooperation0pdf(snrdB);
        digits(4)
        g=vpa(f);
        h=an(n)*exp(-gn(n)*t)*g; 
        h=vpa(h);
        g=char(g);                 % 将 g 转换成字符串表达式
        h=char(h);
        g=strrep(g, '*', '.*');
        g=strrep(g, '/', './');    % 替换运算符
        g=strrep(g, '^', '.^');
        h=strrep(h, '*', '.*');
        h=strrep(h, '/', './');    % 替换运算符
        h=strrep(h, '^', '.^');    % 最终满足要求的 f
        PERn=(1/quad(g,rn(n),rn(n+1)))*quad(h,rn(n),rn(n+1)); %用 quad函数计算数值积分!
        while PERn>p0
              rn(n)=rn(n)+0.001*(rn(n+1)-rn(n));
              PERn=(1/quad(g,rn(n),rn(n+1)))*quad(h,rn(n),rn(n+1));
        end
        rn(n)=rn(n);
    end
    Rn((snrdB+1),:)=rn;
    u=0,v=0;
        for i=1:5             % 求一个分支的协作信道在平均信噪比为snr dB时的频谱效率
            syms f t; 
            f=cooperation0pdfint(snrdB);
            I=int(f,t,Rn(snrdB+1,i),Rn(snrdB+1,i+1));
            u=i*eval(I);
            v=u+v;
        end
        sen0(snrdB+1,:)=v;
        rn=[ 4.2935, 8.6986,24.8142,40.6799,102.5723,1000];
  end

  Rn

  sen0=sort(sen0')
 
  
  
  
  plot(sen0)
  xlabel('SNR(dB)');
  ylabel('Average Spectral Efficiency(bits/symbol)');
  title('没有协作的频谱效率');
 
  

作者:频谱效率数值 来源:频谱效率数值
  • 您是如何找到本站的?
  • 百度搜索
  • Google搜索
  • 查阅资料过程中
  • 论坛发现
  • 百度贴吧发现
  • 朋友介绍
本站最新成功开发工程项目案例
相关文章
  • 没有相关文章
相关评论
发表我的评论
  • 大名:
  • 内容:
  • matlab代做|matlab专业代做|matlab淘宝代做(www.hslogic.com) © 2019 版权所有 All Rights Reserved.
  • Email:highspeed_logic@163.com 站长QQ: 1224848052

    专业代做/代写/承接、MATLAB、SIMULINK、FPGA项目、博士/硕士/本科毕业设计、课题设计、论文,毕业论文,Coursework、Eassy、Assignment