开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 通讯编程 > 查看源码
power_allocate.m
资源名称：cooperative.rar　[点击查看]
上传用户：cjx922
上传日期：2017-12-13
资源大小：205k
文件大小：6k
源码类别：

通讯编程

开发平台：

Matlab

power_allocate.m：源码内容
  1. function [Ps,Pr]=power_allocate(P,asd,asr,ard,combining_type,al_method,transMode,sd_weight);
  2. %根据接收端不同的接收合并方式进行最优功率分配
  4. %MRC方式
  5. switch transMode
  6.     case 'AAF'
  7.         switch combining_type
  8.            case 'MRC'
  9.               for(i=1:1:size(asd,2))
  10.                  switch al_method
  11.                      case 'Best'
  12.                        if(((asd(i)*ard(i)+asr(i)*ard(i)-asd(i)*asr(i))>0)&&((asr(i)*ard(i)*P-asr(i)*asd(i)*P-asd(i))>0))
  13.                          Ps(i)=(asr(i)*ard(i)*P+asd(i)*ard(i)*P+asd(i))/(asd(i)*ard(i)+asr(i)*ard(i)-asd(i)*asr(i)+sqrt((asd(i)*ard(i)+asr(i)*ard(i)-asd(i)*asr(i))*asr(i)*ard(i)*(asr(i)*P+1)/(ard(i)*P+1)));
  14.                          Pr(i)=(asr(i)*ard(i)*P-asr(i)*asd(i)*P-asd(i))/(asd(i)*ard(i)+asr(i)*ard(i)-asd(i)*asr(i)+sqrt((asd(i)*ard(i)+asr(i)*ard(i)-asd(i)*asr(i))*asr(i)*ard(i)*(ard(i)*P+1)/(asr(i)*P+1)));
  15.                        else
  16.                         Ps(i)=P/2;
  17.                         Pr(i)=P/2;
  18.                        end
  19.                     case 'equal'
  20.                        Ps(i)=P/2;
  21.                        Pr(i)=P/2;
  22.                     case 'onlys'
  23.                        Ps(i)=P;
  24.                        Pr(i)=0;
  25.                     otherwise
  26.                        error('No such power allocate method'); 
  27.                  end
  28.               end
  29.             case 'ERC'
  30.                 for(i=1:1:size(asd,2))
  31.                     switch al_method
  32.                         case 'Best'
  33.                             if(((asd(i)-ard(i))*(asd(i)-2*asr(i))>0)&&(2*asr(i)*ard(i)*P-2*asd(i)-2*asd(i)*asr(i)*P-asd(i)*ard(i)*P>0))
  34.                                Ps(i)=(asr(i)*ard(i)*P+asd(i))/(asr(i)*ard(i)-asr(i)*asd(i)+sqrt((asr(i)*P+1)*(asd(i)-ard(i))*(asd(i)-2*asr(i))*asr(i)*ard(i)/(ard(i)*P+2)));
  35.                                Pr(i)=(2*asr(i)*ard(i)*P-2*asd(i)-2*asd(i)*asr(i)*P-asd(i)*ard(i)*P)/(2*asr(i)*ard(i)-2*asr(i)*asd(i)+sqrt((ard(i)*P+2)*(asd(i)-ard(i))*(asd(i)-2*asr(i))*ard(i)*asr(i)/(asr(i)*P+1)));
  36.                             else
  37.                                Ps(i)=P/2;
  38.                                Pr(i)=P/2;
  39.                             end
  40.                         case 'equal'
  41.                             Ps(i)=P/2;
  42.                             Pr(i)=P/2;
  43.                         case 'onlys'
  44.                             Ps(i)=P;
  45.                             Pr(i)=0;
  46.                         otherwise
  47.                             error('No such power allocate method'); 
  48.                     end                           
  49.                 end
  50.             case 'FRC'
  51.                 for(i=1:1:size(asd,2))
  52.                     switch al_method
  53.                         case 'Best'
  54.                             if(((asd(i)*sd_weight^2-ard(i))*(asd(i)*sd_weight^2-sd_weight^2*asr(i)-asr(i))>0)&&((P*asr(i)*ard(i)-asd(i)*sd_weight^2)*(sd_weight^2+1)-P*sd_weight^2*asd(i)*(asr(i)*sd_weight^2+asr(i)+ard(i))>0))
  55.                                Ps(i)=(asr(i)*ard(i)*P+sd_weight^2*asd(i))/(asr(i)*ard(i)-asr(i)*sd_weight^2*asd(i)+sqrt((asr(i)*P+1)*(asd(i)*sd_weight^2-ard(i))*(sd_weight^2*asd(i)-sd_weight^2*asr(i)-asr(i))*asr(i)*ard(i)/(ard(i)*P+sd_weight^2+1)));
  56.                                Pr(i)=((P*asr(i)*ard(i)-asd(i)*sd_weight^2)*(sd_weight^2+1)-P*sd_weight^2*asd(i)*(asr(i)*sd_weight^2+asr(i)+ard(i)))/((asr(i)*ard(i)-asd(i)*asr(i)*sd_weight^2)*(sd_weight^2+1)+sqrt((asd(i)*asr(i)*sd_weight^2-asr(i)*ard(i))*(ard(i)*P+sd_weight^2+1)*(sd_weight^2*asd(i)-sd_weight^2*asr(i)-asr(i))*ard(i)/(asr(i)*P+1)));
  57.                             else
  58.                                Ps(i)=P/2;
  59.                                Pr(i)=P/2;
  60.                             end
  61.                         case 'equal'
  62.                             Ps(i)=P/2;
  63.                             Pr(i)=P/2;
  64.                         case 'onlys'
  65.                             Ps(i)=P;
  66.                             Pr(i)=0;
  67.                         otherwise
  68.                             error('No such power allocate method'); 
  69.                     end                     
  70.                 end
  71.             case 'ESNRC'
  72.                 for(i=1:1:size(asd,2))
  73.                     switch al_method
  74.                       case 'Best'
  75.                            if(((asd(i)-ard(i))*(asd(i)-2*asr(i))>0)&&(2*asr(i)*ard(i)*P-2*asd(i)-2*asd(i)*asr(i)*P-asd(i)*ard(i)*P>0))
  76.                                basic_ps(i)=(asr(i)*ard(i)*P+asd(i))/(asr(i)*ard(i)-asr(i)*asd(i)+sqrt((asr(i)*P+1)*(asd(i)-ard(i))*(asd(i)-2*asr(i))*asr(i)*ard(i)/(ard(i)*P+2)));
  77.                                basic_pr(i)=(2*asr(i)*ard(i)*P-2*asd(i)-2*asd(i)*asr(i)*P-asd(i)*ard(i)*P)/(2*asr(i)*ard(i)-2*asr(i)*asd(i)+sqrt((ard(i)*P+2)*(asd(i)-ard(i))*(asd(i)-2*asr(i))*ard(i)*asr(i)/(asr(i)*P+1)));
  78.                                if(0.1<asd(i)*(asr(i)*basic_ps(i)+ard(i)*basic_pr(i)+1)/(basic_pr(i)*asr(i)*ard(i))&&(asd(i)*(asr(i)*basic_ps(i)+ard(i)*basic_pr(i)+1)/(asr(i)*ard(i)*basic_pr(i))<10))
  79.                                    Ps(i)=(asr(i)*ard(i)*P+asd(i))/(asr(i)*ard(i)-asr(i)*asd(i)+sqrt((asr(i)*P+1)*(asd(i)-ard(i))*(asd(i)-2*asr(i))*asr(i)*ard(i)/(ard(i)*P+2)));
  80.                                    Pr(i)=(2*asr(i)*ard(i)*P-2*asd(i)-2*asd(i)*asr(i)*P-asd(i)*ard(i)*P)/(2*asr(i)*ard(i)-2*asr(i)*asd(i)+sqrt((ard(i)*P+2)*(asd(i)-ard(i))*(asd(i)-2*asr(i))*ard(i)*asr(i)/(asr(i)*P+1)));
  81.                                else
  82.                                    Ps(i)=P/2;
  83.                                    Pr(i)=P/2;
  84.                                end
  85.                            else
  86.                                Ps(i)=P/2;
  87.                                Pr(i)=P/2;
  88.                            end  
  89.                         case 'equal'
  90.                            Ps(i)=P/2;
  91.                            Pr(i)=P/2;
  92.                         case 'onlys'
  93.                            Ps(i)=P;
  94.                            Pr(i)=0;
  95.                         otherwise
  96.                             error('No such power allocate method'); 
  97.                     end
  98.                 end                        
  99.             otherwise
  100.                error(['Unknown relay-type:',combining_type]);
  101.         end
  102.     case 'DAF'
  103.         for(i=1:1:size(asd,2))
  104.             Ps(i)=P/2;
  105.             Pr(i)=P/2;
  106.         end
  107.     otherwise
  108.         error('No such tran Mode'); 
  109. end
  110.             
  111.             
  112.             
  113.