mpointm.m

来自「终点位移法生成具有确定维数的程序」· M 代码 · 共 33 行

% function [z]=MPointM(t,H);
%中点随机位移法生成具有确定维数的曲线
%H:Hurst指数；维数D=2-H；改变此值得到不同维数的曲线
%t:迭代次数,此值表明最后所需要数据量，比如t为10则最后生成2^10=1024个数据
%
%
%
% clear
% t=10 %迭代次数
% H=0.3 
Number=2
B=normrnd(0,(0.5^(2*H))*4)

z(1)=0
z(2)=B
for i=1:t        
    zz=z
    z=0
    for m=1:Number
        z(2*m-1)=zz(m)
    end
    for n=1:Number-1
        B=normrnd(0,0.5^((i+1)*H))
        z(2*n)=(z(2*n-1)+z(2*n+1))/2+B
    end
    Number=Number*2-1%计算总点数        
end
tt=[1:length(z)]
plot(tt,z)
grid on