| 查看: 2384 | 回复: 12 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
用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) 辛苦了  |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有123人回复
-
物理学I论文润色/翻译怎么收费?
已经有292人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 多元函数拟合
已经有9人回复
- 使用matlab最优化方法拟合获得多个动力学参数中的问题
已经有4人回复
- 求助：用matlab或其它软件进行如下计算模拟，多谢
已经有3人回复
- 拜求！matlab拟合动力学参数
已经有14人回复
- ]紧急求助复介电常数的模型拟合-Matlab or 1stOpt
已经有9人回复
- 求大神用matlab给我处理一组数据，有人没有啊？
已经有15人回复
- 求Origin或Matlab高手非线性拟合公式
已经有5人回复
- matlab 拟合求助
已经有8人回复
- 请帮忙拟合一组数据，求回归系数
已经有3人回复
- 三维图谱中，用matlab如何实现去除一部分数据并实现曲面拟合
已经有10人回复
- 用matlab里的nlinfit如何拟合分段函数
已经有8人回复
- matlab对数据进行最小二乘拟合，若果数据为负数怎么拟合？急急急！
已经有3人回复
- matlab数据拟合问题求助
已经有7人回复
- 用MATLAB或者ORIGIN怎么拟合复数数据
已经有10人回复
- 坐等----MAtlab中如何有三维数据拟合曲面，并得到曲面方程？？？？？？？
已经有10人回复
- matlab拟合求值
已经有17人回复
- 需要用matlab拟合数据的外行人急需帮助
已经有10人回复
- 【求助】matlab曲面拟合表达式
已经有10人回复
|
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
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68178.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
感谢参与，应助指数 +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 |
2楼2014-10-28 17:00:36
3楼2014-10-28 17:02:51
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68178.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
qj2006: 金币+20 2014-10-28 21:23:55
qj2006: 金币+20 2014-10-28 21:23:55
4楼2014-10-28 17:07:52
