RS 分形程序-经管之家官网！

坛友互助群

扫码加入各岗位、行业、专业交流群

请教高手 我这个程序运行出现点问题 ，哪位高手帮忙解决一下？谢谢

function[vmrs,VERS,H,C,R]=RS2(vector)
%this program makes use of the classical R/S method.
N=length(vector);
n=5;
vmrs=[];
%the vector vmrs is used to save the log of the mean value of R/S.
h=[];
%this vector is used to save the log value of length n.
VRS=[];
%it represents V-statistic
while n<N
a=floor(N/n);
s=0;
for i=1:a;
v1=zeros(n,1);
vl=vector(n*(i-1)+1:n*(i-1)+n);
m=mean(v1);
%find the mean value of the I-th sub-vector.
v2=v1-m;
v3=cumsum(v2);
%the cumulative deviations from mean of the I-th sub-vector.
r=max(v3)-min(v3)
%range calculation.;
sig=std(v1,1);
%sample standard deviations normalized by n.
rs=r/sig;
s=s+rs;
%the sum of rescaled range for‘a’sub-vectors correspond to the
%length n of subvector.
end
mrs=s/a;
vmrs=[vmrs;log10(mrs)];
h=[h;log10(n)];
VRS=[VRS;mrs/sqrt(n)];
n=n+1;
end
figure(1);
[H,C,R]=postreg(vmrs','h');
%to make use of BP network to make regression about‘vmrs’versus‘h’
figure(2)'
plot(h,vmrs,'*')
hold on;
y=C+H*h;
plot(h.y,'r-');
hold on;
polt(h,VRS,'r');
disp(VRS');
figure(3);
plot(h,vmrs,'*')
Cn=2^(2*H-1)-1;
disp('the correlation funtion Cm=');
disp(Cn);

扫码或添加微信号：坛友素质互助

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP