dcd.m

大尺度衰落的仿镇

%dcd.m                                                          
% 大尺度损耗模型
% 参数设置
% distance 		: T-R距离
% frequency 	: 载波频率
% n         	: 路径损耗指数
% sigma     	: (0均值)正态随机变量（即高斯变量）的标准偏差
% d0        	: 近地参考距离
% m         	: 预测样本数目
% gt        	: 发射天线的增益
% gr        	: 接收天线的增益
hold off;
clear;
c=3*10^8;
disp('** Please enter para **');
disp('** Attention: d0 must <= distance **');
d0=input('d0=');
distance=input('distance=');
frequency=input('frequency=');
n=input('n=');
sigma=input('sigma=');
m=input('m=');
gt=input('gt=');
gr=input('gr=');
%这里按照课本要求仅仅对distance和d0之间的关系作检查，当然我们还可以添加检查程序对其它数据作检查。
while distance<d0|d0<=0
    disp('错误：近地参考距离必须小于等于T-R距离!');
    disp('      d0 must <= distance, 请重新输入参数!');
    d0=input('d0=');
    distance=input('distance=');
end;

%自取附加测试点distancex，步长100m，从d0到distance
%每100米步长内每隔10m取1个测试点distancexx,每点作平均损耗计算(不考虑对数正态阴影)，连成线
xstep=100;
xxstep=10;
for distancex=d0:xstep:distance;
    %不考虑对数正态阴影
    pathlossv = -10*log10(gt*gr*(c/frequency/4/pi)^2/d0^n) + 10*n*log10(distancex/d0);
    %产生1×m阶0均值sigma标准偏差的随机正态分布矩阵,m为预测样本数目
    R=normrnd(0,sigma,1,m);
    for i=1:m;
        %考虑对数正态阴影
        pathloss = pathlossv + R(i);
        %从d0到T-R距离上每100米取为附加测试点，画出其含对数阴影的损耗
        plot(distancex,pathloss,'b.');
        hold on;
    end;
    %每100米步长内每隔10m取1个测试点distancexx,每点作平均损耗计算(不考虑对数正态阴影)，连线
    if distancex<distance
        for distancexx=distancex:xxstep:100+distancex;
        pathlossvv = -10*log10(gt*gr*(c/frequency/4/pi)^2/d0^n) + 10*n*log10(distancexx/d0);
        %防止打印文档看不出此线，特设为星形，显示为一条很粗的红线
        plot(distancexx,pathlossvv,'r*');
        hold on;
        end;
    end;
end;
title('大尺度路径损耗(对数正态阴影)与T-R距离的关系');
xlabel('T-R距离(m)');
ylabel('路径损耗(dB)');