| 查看: 2729 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
感觉是下周放榜了
已经有8人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有81人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 粒子之间的距离怎么求解
已经有6人回复
- 有半个自由度一说吗？
已经有13人回复
- 拍摄到粒子的照片，怎样得到样品的粒径分布图？可以应用什么软件？
已经有8人回复
- 纳米粒子尺寸分布
已经有20人回复
- XRD和TEM能统计纳米粒子尺寸的分布吗？
已经有10人回复
- origin8制作粒径分布的视频（自做的），旧版和新版视频
已经有19人回复
- 用什么仪器表征纳米粒子的粒径分布？
已经有10人回复
- fluent中颗粒速度求助
已经有19人回复
- 【求助】FITC荧光检测的文献，要看到纳米粒子在细胞内的分布
已经有18人回复
- 【转帖】动力学蒙卡方法（KMC）模拟薄膜生长的基本知识
已经有11人回复
- 【求助】Materials explorer可以模拟预测晶体生长的形貌吗？
已经有9人回复
- 【求助】有关无限深势阱及一维谐振子的一些新的简介及分析
已经有7人回复
- 【求助】控制聚合物基纳米复合材料中无机纳米粒子空间分布的方法
已经有6人回复
9楼2014-08-08 16:04:24
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
25