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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 用MATLAB拟合如下数据（求程序代码和拟合结果）
131/2返回列表12下一页
查看: 2368  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

qj2006

金虫 (著名写手)

[求助] 用MATLAB拟合如下数据（求程序代码和拟合结果） 已有1人参与

现有如下数据
x=[2.07931 ,1.80644,1.61,1.48184,1.4116];
y=[1.74024,1.63,1.49,1.4,1.36];
根据如下关系式拟合 y=atan((3.14^2*459*1.449e-29*x*A/6.626E-34*3753^2)^0.5)*x/(3.14^2*459*1.449e-29*x*A/6.626E-34*3753^2)^0.5)
辛苦了
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-10-28 16:44:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
感觉你给出的公式有点问题，一下结果仅供参考：
均方差(RMSE): 0.0590588265039104
残差平方和(SSE): 0.0174397249400949
相关系数(R): 0.993602844977662
相关系数之平方(R^2): 0.987246613547703
决定系数(DC): 0.828021511100056
卡方系数(Chi-Square): 0.00560692501155114
F统计(F-Statistic): 232.231639158851

参数        最佳估算
----------        -------------
a                  1.52138772416843E-16

====== 结果输出 =====

No        实测值y        计算值y
1        1.7402        1.8354974
2        1.63        1.6175395
3        1.49        1.4570489
4        1.4        1.3505577
5        1.36        1.2916440
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2014-10-28 17:00:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qj2006

金虫 (著名写手)
引用回帖:
2楼: Originally posted by 月只蓝 at 2014-10-28 17:00:36
感觉你给出的公式有点问题，一下结果仅供参考：
均方差(RMSE): 0.0590588265039104
残差平方和(SSE): 0.0174397249400949
相关系数(R): 0.993602844977662
相关系数之平方(R^2): 0.987246613547703
决定系数(D ...

谢谢你，你给我发个拟合的代码，其他的我想自己拟合
赞一下 回复此楼
3楼2014-10-28 17:02:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
qj2006: 金币+20 2014-10-28 21:23:55
引用回帖:
3楼: Originally posted by qj2006 at 2014-10-28 17:02:51
谢谢你，你给我发个拟合的代码，其他的我想自己拟合...

http://muchong.com/bbs/viewthread.php?tid=6572963&authorid=1122189
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2014-10-28 17:07:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qj2006

金虫 (著名写手)
引用回帖:
2楼: Originally posted by 月只蓝 at 2014-10-28 17:00:36
感觉你给出的公式有点问题，一下结果仅供参考：
均方差(RMSE): 0.0590588265039104
残差平方和(SSE): 0.0174397249400949
相关系数(R): 0.993602844977662
相关系数之平方(R^2): 0.987246613547703
决定系数(D ...

function fit_nonl
clear all;clc
format long

data=[    2.07931   1.74024
   1.80644   1.63
   1.61   1.49
   1.48184   1.4
   1.4116   1.36
   ];
xspan=data(:,2);    %x的数据，在此输入
Texp=data(:,1);    %T的数据，在此输入

k0=[0.1 0.1 1 100];
lb=-[1 1 1 1]*1e9;
ub=[1 1 1 1]*1e9;


%-------------------------------------------------------------------------

% 使用函数lsqnonlin()进行参数估计

OPTIONS=optimset('MaxFunEvals',1000);
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
    lsqnonlin(@ObjFunc,k0,lb,ub,OPTIONS,xspan,Texp);

ci = nlparci(k,residual,jacobian);
%residual;
fprintf('\n\n拟合结果:\n')
fprintf('\n\t参数 E1 = %.16f',k(1))
y=KineticsEqs(xspan,k);
R2=1-sum((Texp-y).^2)./sum((Texp-mean(y)).^2);
fprintf('\n\tR^2 = %.16f',R2);
figure
plot(xspan,KineticsEqs(xspan,k),'b',xspan,Texp,'or'),legend('计算值','实验值','Location','Best')


%-------------------------------------------------------------------------

function f = ObjFunc(k,xspan,Texp)
f=KineticsEqs(xspan,k)-Texp;

%------------------------------------------------------------------------
function xt = KineticsEqs(x,k)
E1=k(1);
xt=atan((3.14^2*459*1.449e-29*x*E1/6.626E-34*3753^2)^0.5)*x/(3.14^2*459*1.449e-29*x*E1/6.626E-34*3753^2)^0.5);
我这个为啥拟合不出结果，系统老提示有问题，谢谢你的帮忙
赞一下 回复此楼
5楼2014-10-28 17:27:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qj2006

金虫 (著名写手)
引用回帖:
4楼: Originally posted by 月只蓝 at 2014-10-28 17:07:52
http://muchong.com/bbs/viewthread.php?tid=6572963&authorid=1122189...

你好，我的还是拟合不出来，帮我改改那个代码吧，麻烦你了，谢谢
赞一下 回复此楼
6楼2014-10-28 18:24:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
6楼: Originally posted by qj2006 at 2014-10-28 18:24:58
你好，我的还是拟合不出来，帮我改改那个代码吧，麻烦你了，谢谢...

用1stopt软件试试。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
7楼2014-10-28 19:15:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qj2006

金虫 (著名写手)
引用回帖:
7楼: Originally posted by 月只蓝 at 2014-10-28 19:15:53
用1stopt软件试试。...

我没有这个软件
回复此楼
8楼2014-10-28 19:53:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
8楼: Originally posted by qj2006 at 2014-10-28 19:53:04
我没有这个软件...

http://muchong.com/bbs/viewthread.php?tid=4586018&fpage=1&target=blank
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
9楼2014-10-29 08:28:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qj2006

金虫 (著名写手)
引用回帖:
9楼: Originally posted by 月只蓝 at 2014-10-29 08:28:13
http://muchong.com/bbs/viewthread.php?tid=4586018&fpage=1&target=blank...

你好帮我把剩下几组拟合一下，我实在不会，谢谢你
第一组
x=[2.07931 ,1.80644,1.61,1.48184,1.4116];
y1=[1.74024,1.63,1.49,1.4,1.36];
根据如下关系式拟合 y=atan((3.14^2*459*1.449e-29*x*A/6.626E-34*3753^2)^0.5)*x/(3.14^2*459*1.449e-29*x*A/6.626E-34*3753^2)^0.5)；

第二组
x=[2.07931 ,1.80644,1.61,1.48184,1.4116];
y2=[1.73395,1.58,1.45,1.37,1.34];
根据如下关系式拟合 y=atan((3.14^2*460*1.444e-29*x*A/6.626E-34*3756^2)^0.5)*x/(3.14^2*460*1.444e-29*x*A/6.626E-34*3756^2)^0.5)；

第三组
x=[2.07931 ,1.80644,1.61,1.48184,1.4116];
y3=[1.55275,1.51,1.41,1.34,1.3];
根据如下关系式拟合 y=atan((3.14^2*460*1.431-29*x*A/6.626E-34*3750^2)^0.5)*x/(3.14^2*460*1.431e-29*x*A/6.626E-34*3750^2)^0.5)；
第四组
x=[2.07931 ,1.80644,1.61,1.48184,1.4116];
y4=[1.55959,1.52,1.42,1.35,1.33];
根据如下关系式拟合 y=atan((3.14^2*461*1.418-29*x*A/6.626E-34*3742^2)^0.5)*x/(3.14^2*461*1.418e-29*x*A/6.626E-34*3742^2)^0.5)；
赞一下 回复此楼
10楼2014-10-29 09:30:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qj2006 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 C79jjtjAKjEn 2026-09-15 3/150 2026-09-18 03:40 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:16 by cITaGg3p5edr
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:06 by cITaGg3p5edr
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-18 01:54 by cITaGg3p5edr
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-18 00:28 by cITaGg3p5edr
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-17 22:39 by cuPZTDXS3VOd
[有机交流] 各位大神，目前国内有哪些比较好用的逆合成软件？ 20+3 zapen 2026-09-15 7/350 2026-09-17 16:40 by 北京莱茵编辑
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-17 16:03 by 0ngsZVSCTS7z
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-17 13:27 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 yaCh4X3wd045 2026-09-14 4/200 2026-09-17 13:17 by ObVzqQQh1rrE
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-17 13:17 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-17 12:39 by ObVzqQQh1rrE
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 04:52 by k06BKNblrGRH
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-17 02:38 by Tql5LQhh5rLK
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 02:02 by Tql5LQhh5rLK
[考博] 上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:36 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭