下面是一段产生log-normal分布的代码,以此进行说明。
clear all;
clc;
for t=1:100
Traffic(t) =curve(t);
end
MaxTraffic = max(Traffic);
w = 0.2;
Wmax = 2*pi*w/3000;
x=[0:10:300];
y=[0:10:300];
Nx=length(x);
Ny=length(y);
Sigma = 0.53;
t = 0;
M = 10*curve(t)/MaxTraffic;
sum = 0;
for i=1:Nx
forj=1:Ny
Mu = log(M)-0.5*Sigma^2;
Rho(i,j) = RhoFromCoordination(x(i),y(j),Wmax,Sigma,Mu);
Lognrnd(i,j) = round(exp(Sigma*Rho(i,j)+Mu));
sum = Lognrnd(i,j)+sum;
end
end
sum
[xi,yi]=meshgrid(0:2:300,0:2:300);
z1=interp2(x,y,Lognrnd,xi,yi,'spline');%三次样条插值
surf(xi,yi,z1)
(1)首先理解meshgrid的原理和用法。简单地说,就是产生Oxy平面的网格坐标。