| 查看: 2648 | 回复: 4 | ||||
[交流]
多重广义分形维数计算已有3人参与
|
|
现已经有部分分形维数代码，但需要增加广义分形维数计算如果有小伙伴会等请赐教，现将已有分形代码进行分享 clc,clear; % function multifractal(A) format long g %输入一个列向量即可 DApath=['E:\MATLAB\fenxing code\geardata\zhengchang\']; % zc=[]; for I=1:2 DATAfile=['Data1-' int2str(I) '.txt']; % [DATAfile DApath]=uigetfile('*.txt','输入信号'); FILENAME=[DApath,DATAfile]; data=load(FILENAME); data=data(1:4096,1); A=data; %format long g L=length(A); i=1; modify=1; tmin=2;% 边框间距，“※” tmax=10; ttmin=-10; ttmax=10;% 自定义 q 的范围 for r=tmin:1:tmax c(i,1)=mod(L,r); i=i+1; end c'; % 计算不能被边长r整除的余数 a=L-c';% 计算并剔除掉不能被边长r整除的原始数据，得到有用数据的个数 n=length(a); % 求解格网化边长的个数，即为 n（网格边长可以有n=9种） TT=[]; j=1; r=tmin; % 自定义项 for i=1:1:n % 即n=25-10+1，自定义的结果 B=A(1:a(i),1); U=reshape(B,r,length(B)/r);%不明白？？？？？？？ T=mean(U);%默认的是求各列的均值 T=T'.*r^3; TT(1: length(T),i)=[T]; %概率密度分布函数 modifying(modify,1)=length(T); r=r+1; if r>= tmin+n % 自定义项 break; % 边长超过10+n，超过初始限制，则程序自动终止 end modify=modify+1; end modifying; % TT= nthroot(TT,1); %对TT开n次方根 开一次方根有意义吗？？？？ % 或者不缩小 % TT; %TT=nonzeros(TT); %矩阵TT中的非零元素 for cugb=1:1:n %cugb=9 modifying_modifying(cugb,1)=sum(modifying(1:cugb)); % 有影响的新加卷 % & * % ￥ # @ ！） …… end modifying_modifying; j=1; % q 为任意数，这里取1到n，为n，与 k取值保持一致，q过大，计算机无法识别，默认为无穷大，q过小，结果接近0，则意义不明确 for q=ttmin:ttmax %这里取 q=-10:1:10 for k=1:1:n X=TT(1:modifying_modifying(k,1),1).^q; if k>1 X=TT(modifying_modifying(k-1,1)+1:modifying_modifying(k,1),1).^q; end t=sum(X); XX(k,j)=[t]; % 这里用到两个循环，即考虑到了幂函数，又需考虑求和 end j=j+1; end XX; % 得到质量分配函数 X=log(tmin:1:tmax); % % X=log(tmax:-1:tmin); % 此系以前的自定义输入结果，“※—4” % Y=log(XX); % % figure(1) % plot(X',Y,'o-k') % 至此，计算多重分形谱的第一步，分配函数构建完毕 % %求质量指数τ(q) side_length= tmin:1:tmax; % 自定义网格边长，“※—5” side_length=side_length'; q=ttmin:ttmax; % q=-5:1:n-10，q=-10:1:10 m=1; [ha,hb]=size(XX); for i=1:1:hb % XX=XX’; % or not s=XX(:,i); % XX b=polyfit(log(side_length),log(s),1); % 在对数尺度下计算斜率 slope(m,1)=b(1,1); m=m+1; end slope; % 这里的Slope即为质量指数，τ(q) N=polyfit(q', slope,1); plot(q', slope) % 此步是考察τ(q)-q 之间的关系， % %求奇异指数 a=diff(slope)./diff(q'); % 第三步计算，diff函数求偏导确实少一列 % %求多重分形谱 q=q'; f_a=a.*q(1:end-1)-slope(1:end-1);%多重分形谱 % % figure(2) % a=sort(a,'ascend'); plot(a,f_a,'o-k') %xlabel('α','FontSize',12); %ylabel('f(α)','FontSize',12); % polyfit(a,f_a,3) a=sort(a,'ascend'); % a+1 % f_a+1 a % 奇异性指数 f_a % 不同层次的分维数 |
回复此楼» 收录本帖的淘贴专辑推荐
 资源收集 |
» 猜你喜欢
想问一下有没有需要文章但是没时间带学生的研究生导师。
已经有18人回复
-
颜宁:基础研究应顶天立地,立地即从需求端作为问题的出发点,立项背景就是面向真实需求
已经有3人回复
-
发现督导在听课，需要跟TA招呼示意吗？有同学看手机、课件有疏漏算大问题吗？
已经有4人回复
-
催稿后秒拒
已经有6人回复
-
前同事拿我以前未中的国自然标书申请今年的国自然，现在基金委发来调查函
已经有10人回复
-
298求调剂
已经有9人回复
-
申请省自然科学基金，研究区能否是省外区域
已经有13人回复
-
国内科研环境差导致国内科研水平差的根源
已经有3人回复
-
CSC在职申请求问
已经有7人回复
-
家乡二本高校/沿海传统私企，如何抉择？
已经有15人回复
2楼2018-06-12 01:42:08
3楼2018-06-23 15:16:03
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
4楼2018-09-26 12:10:47
5楼2018-10-11 09:41:45
15