有高斯噪声的RSSI值仿真代码

   close all;
   clear all;
   clc;


    PtW = 10e-6; %单位是W
    Pt = 10*log10（PtW）; %单位是dB
    f = 9e8; %载频单位是Hz
    n = 2; %路径损耗指数
    d0 = 20; %近地参考距离单位是m
    %d = 100*sqrt（13）; %选定信标节点与未知节点之间的精确T-R距离单位是m
    c = 3*10^8; %光速单位是m/s
    lamida = c/f; %波长单位是m
    Gt = 1;Gr = 1;L = 1; %Gt为发射天线增益;Gr为接收天线增益;L为与传播无关的系统损耗因子（不小于1）

    %PL0为近地参考距离的路径损耗
    %PrW = PtW*Gt*Gr*lamida^2/（（4*pi）^2*d0^2*L） %单位是W
    %PL0 = 10*log10（Pt/Pr） %单位是dB
    PL0 = -10*log10（Gt*Gr*lamida^2/（（4*pi）^2*d0^2*L））; %单位是dB
    Pr0 = Pt-PL0; %单位是dB
    
    %引入坐标矩阵，计算距离
    anchor=[1 1 1 0 0;
        2 2 1 0 0;
        3 3 1 0 0;
        4 1 2 0 0;
        5 2 2 0 0;
        6 3 2 0 0;
        7 1 3 0 0;
        8 2 3 0 0;
        9 3 3 0 0];
    polastx=2.2;
    polasty=2.35;
    ponowx=2.5;
    ponowy=2.5;
    n=2;
    count=0;
    %t-1时刻接收到的RSSI值
   
        for m=1:9
            d1=sqrt（（polastx-anchor（mn））.^2+（polasty-anchor（mn+1））.^2）;
            
            %PL为精确T-R距离的路径损耗;Pr为信标节点的接收功率
            PL = PL0+10*n*log10（d1/d0）; %单位是dB