24小时热门版块排行榜    

查看: 3275  |  回复: 29
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

飞也如铃

新虫 (小有名气)

[求助] 求助大侠帮忙,有关matlab模拟数据

x = [0.768;
0.711;
0.685;
0.671;
0.659;
0.657;
0.655;
0.654;
0.653;
0.654];

y = [4.584;
4.375;
4.234;
4.128;
4.041;
3.972;
3.913;
3.861;
3.815;
3,774];

上面是x,y的数据,公式是
请问哪位大侠帮忙给编辑一下公式代码啊,我自己用matlab模拟不出来,请哪位好心人帮帮忙!
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

dbb627

荣誉版主 (著名写手)

引用回帖:
23楼: Originally posted by guojunhang at 2011-12-03 17:41:42:
程序运行完毕后
        a =       4.167
        b =      0.1139
我想检验这个拟合效果的时候,使用这样一段程序,
但是执行效果显示,很差,说明这个拟合结果严重不靠谱啊。
期待解释。
其实上面也有 ...

这个x的数据做过归一化处理后再回归的
CODE:
a =  4.167; b = 0.1139;
x = [0.768;
0.711;
0.685;
0.671;
0.659;
0.657;
0.655;
0.654;
0.653;
0.654];
y = [4.584;
4.375;
4.234;
4.128;
4.041;
3.972;
3.913;
3.861;
3.815;
3.774];
figure;
plot(x,y,'k-')
x1=x-mean(x);
x1=x1/std(x1);
hold on;
y_fit = a+b*log10(1-b*x1)+(1-b)*log10(1-b*x1+x1);
plot(x,y_fit,'ro-')

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
24楼2011-12-03 20:10:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 30 个回答

juloong

木虫 (职业作家)

X射线衍/散射测试

【答案】应助回帖

飞也如铃(金币+2): 非常感谢您提供的这个信息,但是我找不到我需要的公式,请问我的公式适合fitting里面的哪个类型的公式呢,谢谢回复! 2011-11-22 21:58:19
matlab里面有个cftool工具,拟合起来非常方便的
2楼2011-11-22 21:15:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juloong

木虫 (职业作家)

X射线衍/散射测试

【答案】应助回帖

拟合出来的结果
General model:
       f(x) = a+b*log(1-b*x)+(1-b)*log(1-b*x+x)
Coefficients (with 95% confidence bounds):
       a =      -1.309  (-9942, 9939)
       b =      -117.5  (-1.173e+006, 1.173e+006)

Goodness of fit:
  SSE: 0.4103
  R-square: 0.337
  Adjusted R-square: 0.2541
  RMSE: 0.2265
3楼2011-11-22 21:25:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juloong

木虫 (职业作家)

X射线衍/散射测试

【答案】应助回帖

用这个式子拟合貌似不好
4楼2011-11-22 21:27:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:40 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 7/350 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 两块石头 +6 阿美_Lml888 2026-09-03 8/400 2026-09-05 13:52 by 阿美_Lml888
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +4 SCSIOyxguo 2026-09-03 7/350 2026-09-05 10:39 by Honglianlian
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见