| 查看: 1194 | 回复: 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小海» 猜你喜欢
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有673人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
物理学I论文润色/翻译怎么收费?
已经有262人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有4人回复
» 抢金币啦！回帖就可以得到:
-
五大联赛看好哪支球队
+5/420
-
2026年财务管理、经济与信息技术国际会议 (FMEIT 2026)
+1/386
-
一点琐碎的家务事，却非常困扰我，让我内耗
+2/140
-
哈尔滨工业大学（校本部）电子与信息工程学院2027级博士招生
+1/96
-
坐标北京诚征男友
+1/81
-
5年了，散金祈愿面上得中，愿望成真再散十倍
+1/43
-
上海交通大学张航课题组招聘博士后（电化学能量存储与转换）
+1/34
-
出全新蔡司 Plan?APOCHROMAT 20×/0.8 平场复消色差物镜
+1/9
-
中国科学院生态环境研究中心环境工程相关方向招收保送硕士和联合培养博士和硕士研究生
+1/4
-
科研人员申请国际专利，到底该怎么走？（PCT vs 巴黎公约，一文讲透）
+1/3
-
河海大学海洋环境污染与生态毒理学研究组招收2027级硕士生 及博士研究生
+1/3
-
2026年生物学、大数据与信息工程国际会议 (BSBIE 2026)
+1/3
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/3
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/3
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/3
-
南京医科大学-冯璋课题组-招收2027级生化与分子生物学专业博士和硕士研究生
+1/3
-
27推免-广州医科大学封红亮教授招收精神医学应用心理学推免生
+1/3
-
中国科学院生物物理研究所田勇课题组招聘联合培养生
+1/2
-
北理工集成电路杰青团队 | 诚招科助理
+1/2
-
【产品测评分享】体外转录产量不理想、dsRNA 副产物高？T7 RNA 聚合酶实测体验
+1/1
3楼2019-03-26 19:51:16
2楼2019-03-26 13:50:09
4楼2019-03-26 20:19:44
5楼2019-03-27 20:41:08
