应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » 请教matlab问题
51/1返回列表
查看: 561  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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 ]
回复此楼

» 猜你喜欢
&amp;amp;lt;ahref=&amp;amp;quot;http://www.TickerFactory.com/&amp;amp;quot;&amp;amp;gt;&amp;amp;lt;imgborder=&amp;amp;quot;0&amp;amp;quot;src=&amp;am
1楼 2007-12-23 22:13:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dnp

荣誉版主 (知名作家)

小木虫浪子

优秀版主
★ ★
csfn(金币+2,VIP+0):平安夜，圣诞，元旦节日发主题帖或应助，额外加金币一个祝福~~
你将所有的r中的数据类型在数据显示窗口中改为longG试试看，我看不懂你写的程序，如果不是数据显示格式的问题的话，就建议你看一下你的计算过程是不是有问题。
你在那里面既用平滑又用微分的，计算数据很难把握，而且还平滑很多次，因此，什么情况都能发生。
赞一下(10人) 回复此楼
What would Jesus do?
2楼2007-12-24 21:20:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sgyin

铁杆木虫 (正式写手)
下面是我用另外一组数据算出来的结果，第一列是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
赞一下 回复此楼
&amp;amp;lt;ahref=&amp;amp;quot;http://www.TickerFactory.com/&amp;amp;quot;&amp;amp;gt;&amp;amp;lt;imgborder=&amp;amp;quot;0&amp;amp;quot;src=&amp;am
3楼2007-12-25 09:08:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dnp

荣誉版主 (知名作家)

小木虫浪子

优秀版主
★
csfn(金币+1,VIP+0):Thanks~
平滑算法非常多，这个可能是跟你的数据有点关系。
你先分析一下你的数据，看一下有什么特点然后再选择适当的平滑方法，最简单的就是移动窗口平均法。
赞一下(10人) 回复此楼
What would Jesus do?
4楼2007-12-25 13:05:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sgyin

铁杆木虫 (正式写手)
那请问虫虫dnp，移动窗口平均法的具体算法是怎样的呢？谢谢了啊，百度，google都没找到
还有我现在知道了为什么最后粒径数据后几位都一样了，因为前面的插值精度只有0.001
赞一下 回复此楼
&amp;amp;lt;ahref=&amp;amp;quot;http://www.TickerFactory.com/&amp;amp;quot;&amp;amp;gt;&amp;amp;lt;imgborder=&amp;amp;quot;0&amp;amp;quot;src=&amp;am
5楼2007-12-26 23:47:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sgyin 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 评审有感 +15 popular289 2026-05-18 24/1200 2026-05-21 06:04 by dandongli
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +4 光伏爱好者 2026-05-17 6/300 2026-05-21 01:37 by 17520222914
[有机交流] 反应很差，大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[基金申请] 面上本子正文33页，违规吗？会被低分嘛？ +11 1234567wang 2026-05-17 13/650 2026-05-20 18:55 by tanfancy
[基金申请] 提交了我也来说说感想 +6 fummck 2026-05-20 6/300 2026-05-20 17:29 by 姜太公烤鱼
[论文投稿] Sci. Bull. 悲剧经验 +6 jyang1999 2026-05-16 7/350 2026-05-20 16:39 by feng6531
[基金申请] 今年审到国自然15份，谈谈感受 +19 国自然国社科中 2026-05-17 20/1000 2026-05-20 14:14 by 仲夏夜的星星
[基金申请] 国自然评分 +3 无名者登山 2026-05-20 4/200 2026-05-20 13:29 by vito刘
[基金申请] 河北省自然科学基金 +3 Peterchao 2026-05-18 3/150 2026-05-20 08:57 by 霸_霸
[教师之家] 上海大学实验技术岗位非升即走 +12 嘻嘻哈哈乐呵呵 2026-05-15 13/650 2026-05-20 08:34 by xli1984
[文学芳草园] 献血感触 +6 呀呀好傻 2026-05-19 11/550 2026-05-19 22:26 by 而立得乐
[基金申请] 2026山东省优青 +3 luckinging 2026-05-14 4/200 2026-05-19 16:15 by 52wls
[考博] 找博士生导师 +6 小代想上岸 2026-05-15 7/350 2026-05-19 10:22 by free_fisher
[考博] 26/27申博自荐-锂/钠电池方向 5+3 狗头军师. 2026-05-15 4/200 2026-05-19 09:10 by moonboat
[基金申请] 别被青基扩招骗了！26年科研内卷才刚刚开始 +3 国自然国社科中 2026-05-14 4/200 2026-05-19 08:48 by archvillain
[考博] 2026博士还有哪些学校有名额 +7 小王求读研 2026-05-15 8/400 2026-05-19 08:27 by zhyzzh
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[考博] 光量子物理方向 博士招生 1人（2026.09） +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
信息提示
关闭
请填处理意见
关闭