% function [z]=MPointM(t,H);
%中點(diǎn)隨機(jī)位移法生成具有確定維數(shù)的曲線
%H:Hurst指數(shù)；維數(shù)D=2-H；改變此值得到不同維數(shù)的曲線
%t:迭代次數(shù),此值表明最后所需要數(shù)據(jù)量，比如t為10則最后生成2^10=1024個(gè)數(shù)據(jù)
%
%
%
% clear
% t=10 %迭代次數(shù)
% 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%計(jì)算總點(diǎn)數(shù)        
end
tt=[1:length(z)]
plot(tt,z)
grid on