请大家帮忙给看下个计算广义维数的程序，为什么q取不了负值啊

csuxiaofu · 发表于 2009-4-17 22:08

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

%本程序利用数每个盒子里的点数来计算多重分形
t=1:1:1000;
f=sin(0.03*t);
f=(f-min(f))/(max(f)-min(f));
   delta=0.05;%盒子边长
   numgrid=ceil(1/delta);%网格列数（此时行数与列数相等）
   buju=1000*delta;%为了编程的简便，在此取数据长度为1000，这样每段的数据点为1000*b，
   N_delta=zeros(numgrid,1);
      N_i=0;
   for m=1:1:numgrid
         if m==1
            y=f(1:buju);
         else
               y=f((m-1)*buju:m*buju);
         end
         maxy=max(y);
         miny=min(y);
         c=floor(maxy/delta);
         d=floor(miny/delta);
         if mod(maxy,delta)==0
               N_delta=abs(c-d);
         else
            N_delta=1+abs(c-d);%统计每一列网格里的盒子个数
         end
         N_m_i=zeros(1,N_delta);
            for s=1:1:N_delta
               for i=1:1:size(y,2)
                  if (y(i)>=(d+s-1)*delta) && (y(i)<=((d+s)*delta))%统计每列网格里的每个盒子内的点数
                     N_m_i(1,s)=N_m_i(1,s)+1;%
                  end
               end
            end
            N_i=[N_i N_m_i(1,:)];
      end
      lieshu=size(N_i,2);
      P=N_i(1,2:lieshu)./sum(sum(N_i(1,2:lieshu)));%计算每个盒子内的点数所占的概率

u=0;
for q=0:0.05:20%只是q为什么不能为负的呢？！
u=u+1;
N_q_delta=sum(sum(P.^q));
D(u)=(1/(q-1))*(log(N_q_delta)/log(delta));
end
q=0:0.05:20;
plot(q,D)

本程序的那个q取正值时，运行不出错，效果也还可以，但是为什么取不了负值呢？当取负值时P.^q为变得非常大，以致无穷大？那为什么有些文章却可以呢？感觉我的程序没有问题啊，希望大家能帮忙，谢谢

csuxiaofu · 发表于 2009-4-19 21:13

哈哈，发现问题了，原来矩阵P里有0元素，这样当q取负值时，相当于除0，这样就出错了。不过感觉这个论坛人气不旺啊，对于求助的问题都没有反应，哎……

kingsir · 发表于 2009-5-6 07:13

我以前也遇到过这类问题哦

苏舟 · 发表于 2014-6-11 17:33

我用这个程序进行仿真时，得到的结果只有q=0时有维数，其他的q值下维数都是1，不知道是我数据的原因还是程序中设置的网格大小不合适，请帮忙看一下qq：756047801

账号		自动登录	找回密码
密码			我要加入

[分形与混沌] 请大家帮忙给看下个计算广义维数的程序，为什么q取不了负值啊

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

评分