| 查看: 1014 | 回复: 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小海» 猜你喜欢
西安电子科学大学杭州研究院刘丽香教授招收智能多模态传感器和微型储能器件方向博士
已经有13人回复
-
重庆交通大学光子学微结构与器件课题组招收2026年硕士研究生信息
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有248人回复
-
一志愿郑大材料学硕298分，求调剂
已经有6人回复
-
寻合作：应力腐蚀多尺度模拟
已经有4人回复
-
考研交流
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
» 抢金币啦！回帖就可以得到:
-
Ph.D. Opportunity in Mechanical Engineering – Fall 2026
+2/154
-
江西科技师范大学微生物天然药物2026研究生调剂招生
+1/87
-
北京信息科技大学智能机器人技术方向，联合清华培养硕生、博生，并长期招博士后
+1/87
-
北京信息科技大学联合清华大学某国家重点实验室，联合培养硕生、博生，并长期招博士后
+1/84
-
化生学院蒋老师课题组招收研究生
+1/81
-
中国计量科学研究院长度量子基准实验室2026年调剂招生
+1/32
-
南京林业大学-国家级青年人才团队 招2026级申请考核制博士（化学方向）
+1/17
-
26届计算机、电子信息类、电科、控制、通信考研T j信息pp骏：74+08+76+6+20/
+1/14
-
电子科技大学材料学院SFT创新中心招收2026年博士生及科研助理
+1/13
-
求DNAstar 破解版
+1/9
-
中南林业科技大学万才超教授课题组现招收2026年秋季入学博士研究生
+1/9
-
地理学数二英一316求调剂
+1/9
-
济南大学耿兵氟化学课题组现招收2026级化学工程学硕专硕共4名
+1/8
-
温州大学王娟课题组招收2026年硕博生
+1/8
-
中科院微生物所-肿瘤免疫或微生物组方向的博士后（特别研究助理）
+1/7
-
085704测绘工程 一志愿华东师大 本科双非 总分340求调剂
+1/4
-
食品考研交流
+1/3
-
双一流天津工业大学电信学院李鸿强教授招收2026年申请审核制博士
+1/2
-
招调剂
+1/1
-
140300设计，一志愿985，354求调剂
+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
