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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » 请教matlab问题
51/1返回列表
查看: 610  |  回复: 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的回帖

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的回帖
查看全部 5 个回答

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的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 yPv6ndDtiRxQ 2026-08-10 8/400 2026-08-11 19:56 by thgx1vqMqrLw
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 GpvwGNVshCOp 2026-08-10 7/350 2026-08-11 19:39 by thgx1vqMqrLw
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +16 Lanmanbaby 2026-08-09 26/1300 2026-08-11 17:16 by 天神眷顾
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] FileCode能看出啥？ +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
信息提示
关闭
请填处理意见
关闭