mica2_bak.m

利用CC1000通信的仿真程序.这是适用于无线传感器网络Mica2平台的程序.

%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
% the program is written on path loss model for mica2 nodes by sunpg at  06-09-15.
%################################################################################
%**************************
%define constant variables.from chengdv.
%**************************
%format long;
%Pt=0.28183815;
%Gt = 1;
%Gr = 1;
%freq = 9.14e+8;
%L= 1.0;  
%d0 = 1.0; 
%std_var = 6.0;
%pathlossExp = 2.0; 
%speed_of_light = 3.0e+8;
%lambda = speed_of_light./freq;

%**************************
%define constant variables.from macro zunig..
%**************************
format long;
Pt=1.00;
Ptdb = 0.00;
Gt = 1;
Gr = 1;
freq = 9.14e+8;
%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
% the program is written on path loss model for mica2 nodes by sunpg at  06-09-18.
%################################################################################
format long;

%*************************************
%define wireless channel constants.
%*************************************
Gt = 1;
Gr = 1;
freq = 9.14e+8;
L= 1.0;  
speed_of_light = 3.0e+8;
lambda = speed_of_light./freq;
D = lambda/4;
df = 2*D.^2/lambda;%df = 0.041;
%**************************************
%parameter configuration.
%**************************************
%d0 = 0.8; 
std_var = 4.0;
pathlossExp = 3.5; 
f = 50;
pld0db = 55;
%**************************************
%transmission power configuration.
%**************************************
Pt=1;
Ptdb = 10*log10(Pt);
%**************************************
%the following find d0.(unit:m)
%**************************************
d0 = lambda/(4*pi*sqrt(power(10,-pld0db/10)));
%pld0db = -10*log(lambda.^2/(4*pi*d0).^2);
d = 15:0.01:100;
Xv = normrnd(0,std_var,size(d));
Prddb = 0 - pld0db - 10*pathlossExp*log(d./d0) + Xv;
Pr = power(10,Prddb./10);
plot(d,Prddb);
grid;

%*******************************************************
%prr is drawn.
%*******************************************************
%noise_db = -177;
%Pn = power(10,-177/10);
%snr_p = Pr/Pn;
%snr_db = exp(-Prddb/1.28);
%snr_g = exp(-snr_p/1.28);
%prr = power((1 - 0.5*snr_g),8*f);
%plot(d,Prddb);
%plot(d,prr,'o');
%axis([0 100 0 1]);
%grid;