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

24小时热门版块排行榜    

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

铁杆木虫 (正式写手)
下面是我用另外一组数据算出来的结果，第一列是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的回帖
查看全部 5 个回答

dnp

荣誉版主 (知名作家)

小木虫浪子

优秀版主
★ ★
csfn(金币+2,VIP+0):平安夜，圣诞，元旦节日发主题帖或应助，额外加金币一个祝福~~
你将所有的r中的数据类型在数据显示窗口中改为longG试试看，我看不懂你写的程序，如果不是数据显示格式的问题的话，就建议你看一下你的计算过程是不是有问题。
你在那里面既用平滑又用微分的，计算数据很难把握，而且还平滑很多次，因此，什么情况都能发生。
赞一下(10人) 回复此楼
What would Jesus do?
2楼2007-12-24 21:20:24
已阅   回复此楼   关注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的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +16 Lanmanbaby 2026-08-09 25/1250 2026-08-11 15:51 by 医学老男孩
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +7 Tide man 2026-08-10 8/400 2026-08-11 15:49 by 铁帽子农民
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 帮忙看看fileCode +6 wwncly 2026-08-10 12/600 2026-08-11 13:01 by wwncly
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
关闭
请填处理意见
关闭