| 查看: 961 | 回复: 5 | |||
[交流]
求助各位大神，参数拟合结果不对，金币全送。谢谢
|
 求助各位大神，用下面的程序做的参数拟合，拟合结果和图形都不对，希望大神们帮帮忙，金币全送。谢谢function fit_nonl12 clear all;clc format long data=[ 2561.0000 10.0000 2560.0000 15.0000 2560.0000 25.0000 2558.0000 40.0000 2557.0000 60.0000 2556.0000 100.0000 2555.0000 160.0000 2554.0000 250.0000 2552.0000 400.0000 2551.0000 630.0000 2550.0000 1000.0000 ]; xspan=data(:,2); Texp=data(:,1); k0=[2.0 30.0 6.0 0.5]; lb=[0.1 0.1 0.1 0.0]; ub=[10.0 100.0 20.0 1.0]; % 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)) fprintf('\n\t参数 E2 = %.16f',k(2)) fprintf('\n\t参数 E3 = %.16f',k(3)) fprintf('\n\t参数 E4 = %.16f',k(4)) 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 x1 =kineticsEqs(x,k) E1=k(1);E2=k(2);E3=k(3);E4=k(4); x1=E1+(((2.*pi.*E3).^(-E4)+x.^(E4).*cos(pi.*E4./2)).*E2.*(2.*pi.*E3).^(-E4))./(((2.*pi.*E3).^(-E4)+x.^(E4).*cos(pi.*E4./2)).^2+(x.*(E4).*sin(pi.*E4./2)).^2); |
回复此楼» 本帖已获得的红花（最新10朵）
why小海» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有195人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
坐标济南，来碰碰运气
+1/450
-
中国科学院大学功能多孔组装材料实验室招聘启事
+2/340
-
供应德国EXAKT艾卡特半导体导热散热材料三辊研磨机50 PLUS
+1/85
-
科瑞赛生物内皮细胞培养基试用装限时大放送，助力你的实验高效进阶！
+1/84
-
南京理工大学曾海波/李伟金 招聘博士后（电磁响应：介电调控等方向）
+1/79
-
坐标济南，山东农科院招 有机合成 or 药物化学 联培硕士研究生
+1/37
-
国重点实验室双一流A类长江学者团队招2026年全日制博士1-2名/博后1-2名
+2/34
-
SCI，计算机相关可以写
+1/23
-
SCI，计算机相关可以写
+1/19
-
杨老师招收联合培养硕士、博士生或客座学生
+1/17
-
吉林大学材料物理本科生求问调剂信息
+1/11
-
2026年黄河科技学院纳米功能材料研究所招聘
+2/8
-
【博士后/科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+1/8
-
2026 博士自荐-机器人机构学方向
+1/5
-
中国科学院大学-杨晗课题组-诚聘-博士后、副研究员
+1/5
-
北京工业大学材料学院吴玉锋教授、王长龙研究员招收博士研究生
+1/3
-
海南大学化学院—功能分子器件团队博士后招聘
+1/2
-
香港中文大学（深圳）下一代半导体激光器概念验证中心（筹）招聘工程师
+1/2
-
广东医科大学干细胞与再生组织工程重点实验室联合博士后招聘
+1/2
-
中国矿业大学黄赳课题组联合中国科学院南京土壤研究所朱晓芳研究员诚聘博士后
+1/1
2楼2019-03-26 13:50:09
3楼2019-03-26 19:51:16
4楼2019-03-26 20:19:44
5楼2019-03-27 20:41:08
6楼2019-03-28 11:22:35