24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lvubaby

木虫 (正式写手)

应助: 4 (幼儿园)
金币: 2417.5
散金: 20
红花: 2
帖子: 557
在线: 179小时
虫号: 1172683
注册: 2010-12-18
性别: MM
专业: 传热传质学

[求助] 二维粒子分布的分形维数已有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人回复

海风吹

1楼 2014-04-16 19:32:51

已阅关注TA 给TA发消息送TA红花 TA的回帖

alorsmail

铁杆木虫 (著名写手)

应助: 52 (初中生)
金币: 7398.3
红花: 41
帖子: 2104
在线: 237小时
虫号: 2537936
注册: 2013-07-09
专业: 系统生物学

引用回帖:

7楼: Originally posted by lvubaby at 2014-04-25 07:54:25
你真是大牛啊，按照你说的方法改后，结果是1.313.文献中的结果是1.38+-0.06. 非常温和. 多谢指教！！！...

不客气，实际上，这里还有一个更小的细节：
重心坐标并不一定要取整值，这里用整值只是为了保持和原数据一直。
如果直接用计算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)，