| 查看: 592 | 回复: 4 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
sgyin铁杆木虫 (正式写手)
|
[交流]
请教matlab问题
|
||
|
为什么我用matlab算出来的粒径数值最后一位都是7？ r 1.3857 1.4327 1.4727 1.5027 1.5237 1.5417 1.5627 1.5767 1.5807 1.5887 1.6027 1.6107 1.6177 1.6257 1.6357 1.6407 1.6447 1.6497 1.6587 1.6657 程序如下： [m,n]=size(aa); for t=1:8,%该for循环用于曲线平滑,以及将不必要的数据截掉.t为平滑次数. for s=(length(aa(:,2))-2):-1:2, sp=csapi([4,3,1],[aa(s+2,2),aa(s+1,2),aa(s-1,2)]);aa(s,2)=fnval(sp,2); end end wavelength=aa(1:m-40,1);absorptance=aa(1:m-40,2);E_eV=(h*c./(wavelength*1e-9))./e; particlesize=1/2./(50000*E_eV-157955).*(-5289+(-72372280279+22918000000*E_eV).^(1/2));%能量到粒径的转换 psize=particlesize(1):.001:particlesize(length(particlesize)); sp=csapi(particlesize,absorptance,psize); %%%%%%%%%plot(psize,sp,'--r'),xlabel('粒径'),ylabel('吸光度'),title('吸光度和粒径的关系(红色为插值结果)');%csapi()三次样条插值,fnplt()绘制插值结果。 spd=diff(sp)./diff(psize);psize1=psize(1:length(psize)-1);%spd即吸光度对粒径的导数。 V=4/3*3.1415926*psize1.^3; n_r=-spd./V; i=find(n_r<-1e-3&n_r>-1e+4);i=max(i);n_r(1:i)=NaN; j=find(n_r==max(n_r));psize1(j) [ Last edited by csfn on 2008-12-29 at 20:23 ] |
回复此楼» 猜你喜欢
刑法学论文投稿求助
已经有5人回复
-
大龄残疾硕士的一点执念
已经有8人回复
-
今年E04面上
已经有14人回复
-
Journal of Environmental Chemical Engineering
已经有3人回复
-
最后一年，祈求好运
已经有3人回复
-
风电环氧领域
已经有3人回复
-
2026年面上项目中了，2A+B, 会评顺利通过
已经有8人回复
-
昨日死，今日生
已经有8人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
sgyin
铁杆木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5519.4
- 散金: 100
- 红花: 8
- 帖子: 680
- 在线: 279.4小时
- 虫号: 405447
- 注册: 2007-06-17
- 性别: GG
- 专业: 材料物理化学
5楼2007-12-26 23:47:16
dnp
荣誉版主 (知名作家)
小木虫浪子
- 应助: 74 (初中生)
- 贵宾: 18.476
- 金币: 18827.9
- 散金: 2040
- 红花: 77
- 沙发: 6
- 帖子: 7319
- 在线: 598.4小时
- 虫号: 437184
- 注册: 2007-08-26
- 性别: GG
- 专业: 生物信息学
- 管辖: 有机交流
2楼2007-12-24 21:20:24
sgyin
铁杆木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5519.4
- 散金: 100
- 红花: 8
- 帖子: 680
- 在线: 279.4小时
- 虫号: 405447
- 注册: 2007-06-17
- 性别: GG
- 专业: 材料物理化学
|
下面是我用另外一组数据算出来的结果，第一列是long型，第二列是short，我觉得是那个样条平滑处理方法导致的，对它的平滑确定性没把握。明显看出数据的尾部都一样，这是不是说明这种算法有效数字只有前面的5位啊？不知道大家有什么更好的平滑算法？？？谢谢！ 1.31675586473640 1.3168 1.44455586473641 1.4446 1.54515586473641 1.5452 1.57215586473640 1.5722 1.58205586473640 1.5821 1.59265586473640 1.5927 1.60305586473641 1.6031 1.61725586473640 1.6173 1.62705586473641 1.6271 1.63455586473640 1.6346 1.64025586473640 1.6403 1.64855586473640 1.6486 1.66195586473641 1.662 1.66535586473641 1.6654 1.67405586473641 1.6741 1.68485586473640 1.6849 1.69125586473641 1.6913 1.70025586473641 1.7003 1.70755586473640 1.7076 1.71055586473641 1.7106 |
3楼2007-12-25 09:08:50
dnp
荣誉版主 (知名作家)
小木虫浪子
- 应助: 74 (初中生)
- 贵宾: 18.476
- 金币: 18827.9
- 散金: 2040
- 红花: 77
- 沙发: 6
- 帖子: 7319
- 在线: 598.4小时
- 虫号: 437184
- 注册: 2007-08-26
- 性别: GG
- 专业: 生物信息学
- 管辖: 有机交流
4楼2007-12-25 13:05:43