curv.m

来自「无线移动信道的matlab源码 用于对无线通信信道的编程实现」· M 代码 · 共 31 行

function out=curv(thita)
%      p=0.8;
%      sigma_1_2=1;
%      sigma_2_2=p*sigma_1_2;
sigma_1_2=0.10391;
sigma_2_2=0.030488;

     BEITA1=1103.4298;
     BEITA2=1091.5206;
     
     r=10;
     N_r=zeros(length(thita),length(r));
     for i=1:length(r)
         ri=10^(r(i)/20);
         a1=sigma_2_2*(cos(thita).*cos(thita));
         a2=sigma_1_2*(sin(thita).*sin(thita));
         c=a1+a2;
         b=sqrt(sigma_1_2*sigma_2_2);
         d=(2*pi)^(3/2);
         b1=ri/(d*b);
         b2=-(ri.*ri)/(2*sigma_1_2*sigma_2_2);
         b3=exp(b2.*c);
         a3=BEITA1*(cos(thita).*cos(thita));
         a4=BEITA2*(sin(thita).*sin(thita));
         c2=sqrt(a3+a4);
         d1=b3.*c2;
         N_r(:,i)=transpose(b1).*transpose(d1);
     end
     out=N_r;