24小时热门版块排行榜    

查看: 2365  |  回复: 8

smileyjo

新虫 (初入文坛)

[求助] MATLAB作图拟合求曲线 已有6人参与

x=[0.066667 0.133333 0.2 0.266667 0.333333 0.4 0.466667 0.533333 0.6 0.666667]'; y=[0.073 0.0145 0.023 0.0017 0.0013 0.0013 0.0023 0.0005 0.0003 0.0009]'; 已知方程满足y= a*x^-11+b*x^-31,要怎么拟合求解系数呢,跪求代码T T~~~~ 如果要y= a*x^-11+b*x^-31+c 加上常数项,又应该怎么写代码呢? 求助啊,太感谢了!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
1stOpt求解:
CODE:
Function  y= a*x^(-11)+b*x^(-31);
Data;
x=[0.066667,0.133333,0.2,0.266667,0.333333,0.4,0.466667,0.533333,0.6,0.666667];
y=[0.073,0.0145,0.023,0.0017,0.0013,0.0013,0.0023,0.0005,0.0003,0.0009];

均方差(RMSE):0.00730705877927612
残差平方和(SSE):0.000533931080037963
相关系数(R): 0.953599387787668
相关系数之平方(R^2): 0.909351792389015
确定系数(DC): 0.885658460987378
卡方系数(Chi-Square): 0.0307716782978097
F统计(F-Statistic): 69.9658240485718

参数        最佳估算
----------        -------------
a        3.49576759701372E-12
b        -1.04884423089604E-35
CODE:
Function  y= a*x^(-11)+b*x^(-31)+c;
Data;
x=[0.066667,0.133333,0.2,0.266667,0.333333,0.4,0.466667,0.533333,0.6,0.666667];
y=[0.073,0.0145,0.023,0.0017,0.0013,0.0013,0.0023,0.0005,0.0003,0.0009];

均方差(RMSE):0.00643858647334295
残差平方和(SSE):0.000414553957747148
相关系数(R): 0.954580073144648
相关系数之平方(R^2): 0.911223116044842
确定系数(DC): 0.911223116044842
卡方系数(Chi-Square): 0.0246845977988711
F统计(F-Statistic): 39.4246772816209

参数        最佳估算
----------        -------------
a        2.56882567356543E-12
b        -7.70192879141E-36
c        0.00386877580554191
2楼2014-07-11 16:40:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smileyjo

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2014-07-11 16:40:01
1stOpt求解:

Function  y= a*x^(-11)+b*x^(-31);
Data;
x=;
y=;

均方差(RMSE):0.00730705877927612
残差平方和(SSE):0.000533931080037963
相关系数(R): 0.953599387787668
相关系数之平方(R^2): 0.90 ...

谢谢,请问是用什么方法计算的?可以提供代码么?
3楼2014-07-11 17:26:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
2楼的专家已经给出代码,而且他的计算结果应该是最优了。
不过就方法而言,包括Origin在内的4款软件操作步骤参见:
http://muchong.com/bbs/viewthread.php?tid=6572963&authorid=1122189
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2014-07-11 19:31:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smileyjo

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by 月只蓝 at 2014-07-11 19:31:26
2楼的专家已经给出代码,而且他的计算结果应该是最优了。
不过就方法而言,包括Origin在内的4款软件操作步骤参见:
http://muchong.com/bbs/viewthread.php?tid=6572963&authorid=1122189

呃。。刚刚接触MATLAB,比较小白,为什么把代码复制进去运行不起来呢?。。。因为Y的数据有好几组,如果y分别是 y=[0.296 0.156 0.101 0.064 0.0453 0.0288 0.0252 0.0183 0.012 0.01]’; y=[0.029 0.012 0.0083 0.0067 0.0053 0.0049 0.0038 0.0032 0.003 0.0024]’; y=[0.066 0.02 0.011 0.0092 0.0064 0.0051 0.0038 0.0032 0.0034 0.0024]'; 的话结果是多少呢? 太感谢了T-T
5楼2014-07-11 22:52:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

waterque

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
用matlab拟合,直接用其拟合工具,cftool (在命令窗口输入该命令即可),在弹出窗口中选择自变量x,因变量y, 以及拟合函数(可以自定义)
行自己的路:waterque@hotmail.com
6楼2014-07-12 14:47:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kabaisun

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
推荐用matlab的拟合工具箱,不必自己写代码。
严实诚勇
7楼2014-07-12 15:17:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shikang999

新虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
你这个是线性方程组,完全没有必要采用现成拟合算法进行拟合。采用最小二乘求解一个超定方程组即可得到最优系数!
PS:
y=a*x^-11+b*x^-31
y=A*X
A=[a    b]
X=[x^-11    x^-31]
我还是那个我,过去如此,未来亦如此!
8楼2014-07-13 00:54:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyong9572

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
用1stopt,非常好用的一个小软件,不用安装,网上搜搜看,拟合曲线和matlab完全一样,亲测。

[ 发自小木虫客户端 ]
9楼2014-07-13 07:18:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 smileyjo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 h4CP7TrQR8Lg 2026-08-22 4/200 2026-08-23 20:55 by 0p9y1isJp2jM
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 19:34 by lXgU0RvNTDC9
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 19:34 by lXgU0RvNTDC9
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 19:22 by lXgU0RvNTDC9
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 19:10 by lXgU0RvNTDC9
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 18:55 by lXgU0RvNTDC9
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 18:46 by lXgU0RvNTDC9
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 15:30 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:21 by KM5EcsNQRBPn
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +3 炎黄贵胄 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
信息提示
请填处理意见