| 查看: 2275 | 回复: 9 | ||
lvubaby木虫 (正式写手)
|
[求助]
二维粒子分布的分形维数 已有1人参与
|
|
大家好，我得到了一个二维的粒子分布，想计算他的分形维数。采用回转半径法，输入粒子坐标，编了一个小程序，可是计算得到的分形维数为2.262大于2。理论上讲二维图像下的分形维数应该小于2啊，是哪里出错了，请各位大神指点。 程序为： clc clear dim=128; load final_paticle_distribution_2D.txt particle_num=size(final_paticle_distribution_2D,1); c1=final_paticle_distribution_2D(:,1); c2=final_paticle_distribution_2D(:,2); figure(1) plot(c1,c2,'.') %%求得重心坐标 mc1=round(sum(c1)/particle_num); mc2=round(sum(c2)/particle_num); num=zeros(dim,1);%粒子数目 r=zeros(dim,1);%半径 s=zeros(particle_num,1);%粒子与重心距离 for i=1:particle_num %%计算粒子与质心的距离 s(i)=sqrt((c1(i)-mc1)^2+(c2(i)-mc2)^2); for j=1:dim r(j)=j-0.5; if s(i)<=r(j) num(j)=num(j)+1; end end end r_num=[r,num]; save('radius_particlenum.txt', 'r_num','-ASCII') %%将不含有粒子的空间删去 %%外侧 ii=1; for jj=1:dim-1 if num(jj)==num(jj+1) num_zero_ind_high(ii)=jj+1; ii=ii+1; end end num(num_zero_ind_high)=[]; r(num_zero_ind_high)=[]; %%内侧 particle_ind=find(num); particle_num=num(particle_ind); particle_r=r(particle_ind); log_num=log(particle_num); log_r=log(particle_r); fitobject=fit(log_r,log_num,'poly1') figure(2); plot(log_r,log_num,'o');  2D particle distribution.png |
回复此楼» 猜你喜欢
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 粒子之间的距离怎么求解
已经有6人回复
- 有半个自由度一说吗？
已经有13人回复
- 拍摄到粒子的照片，怎样得到样品的粒径分布图？可以应用什么软件？
已经有8人回复
- 纳米粒子尺寸分布
已经有20人回复
- XRD和TEM能统计纳米粒子尺寸的分布吗？
已经有10人回复
- origin8制作粒径分布的视频（自做的），旧版和新版视频
已经有19人回复
- 用什么仪器表征纳米粒子的粒径分布？
已经有10人回复
- fluent中颗粒速度求助
已经有19人回复
- 【求助】FITC荧光检测的文献，要看到纳米粒子在细胞内的分布
已经有18人回复
- 【转帖】动力学蒙卡方法（KMC）模拟薄膜生长的基本知识
已经有11人回复
- 【求助】Materials explorer可以模拟预测晶体生长的形貌吗？
已经有9人回复
- 【求助】有关无限深势阱及一维谐振子的一些新的简介及分析
已经有7人回复
- 【求助】控制聚合物基纳米复合材料中无机纳米粒子空间分布的方法
已经有6人回复
alorsmail
铁杆木虫 (著名写手)
- 应助: 52 (初中生)
- 金币: 7398.3
- 红花: 41
- 帖子: 2104
- 在线: 237小时
- 虫号: 2537936
- 注册: 2013-07-09
- 专业: 系统生物学
2楼2014-04-20 04:06:51
lvubaby
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 2417.5
- 散金: 20
- 红花: 2
- 帖子: 557
- 在线: 179小时
- 虫号: 1172683
- 注册: 2010-12-18
- 性别: MM
- 专业: 传热传质学
3楼2014-04-20 16:39:55
alorsmail
铁杆木虫 (著名写手)
- 应助: 52 (初中生)
- 金币: 7398.3
- 红花: 41
- 帖子: 2104
- 在线: 237小时
- 虫号: 2537936
- 注册: 2013-07-09
- 专业: 系统生物学
4楼2014-04-22 17:50:23
alorsmail
铁杆木虫 (著名写手)
- 应助: 52 (初中生)
- 金币: 7398.3
- 红花: 41
- 帖子: 2104
- 在线: 237小时
- 虫号: 2537936
- 注册: 2013-07-09
- 专业: 系统生物学
5楼2014-04-22 21:33:37
alorsmail
铁杆木虫 (著名写手)
- 应助: 52 (初中生)
- 金币: 7398.3
- 红花: 41
- 帖子: 2104
- 在线: 237小时
- 虫号: 2537936
- 注册: 2013-07-09
- 专业: 系统生物学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lvubaby: 金币+20, ★★★★★最佳答案 2014-04-25 07:52:38
lvubaby: 金币+20, ★★★★★最佳答案 2014-04-25 07:52:38
|
你好，我找到问题所在了， 你的程序没有问题。之所以出错，是因为一个细节： 重心坐标 你这里用的是平均值，但是这个并非真正的坐标。 所以建议用 “核密度估计 ” Kernel density estimation 来找更加合理的重心 比如 ================================= [F,XI]=ksdensity(c1) ; mc1 = floor( XI(find(F==max(F))) ) ; %#ok<*FNDSB> [F,XI]=ksdensity(c2) ; mc2 = floor( XI(find(F==max(F))) ) ; ================================= 然后运行程序，得到的分形维数是 1.313 （大致是这个数值，每次运行可能有变化，这个是由于MATLAB的线性逼近引起的） 这个结果比较合理， 用Box counting得到的fractal dimension 是 1.46左右。考虑到用盒子的误差，你的结果应该是对的。 所以，恭喜你！   |
6楼2014-04-25 05:37:27
lvubaby
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 2417.5
- 散金: 20
- 红花: 2
- 帖子: 557
- 在线: 179小时
- 虫号: 1172683
- 注册: 2010-12-18
- 性别: MM
- 专业: 传热传质学
7楼2014-04-25 07:54:25
alorsmail
铁杆木虫 (著名写手)
- 应助: 52 (初中生)
- 金币: 7398.3
- 红花: 41
- 帖子: 2104
- 在线: 237小时
- 虫号: 2537936
- 注册: 2013-07-09
- 专业: 系统生物学
|
不客气，实际上，这里还有一个更小的细节 ： 重心坐标并不一定要取整值，这里用整值只是为了保持和原数据一直。 如果直接用计算ksdensity的计算结果的话， ================================= [F,XI]=ksdensity(c1); mc1 = ( XI(find(F==max(F))) ); %#ok<*FNDSB> [F,XI]=ksdensity(c2); mc2 = ( XI(find(F==max(F))) ); ================================= 得到的将是 1.353 (1.32, 1.385)， |
8楼2014-04-25 12:24:15
9楼2014-08-08 16:04:24
10楼2014-08-08 16:05:36