24小时热门版块排行榜    

查看: 1007  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

why小海

新虫 (初入文坛)


[交流] 求助各位大神,参数拟合结果不对,金币全送。谢谢

求助各位大神,用下面的程序做的参数拟合,拟合结果和图形都不对,希望大神们帮帮忙,金币全送。谢谢
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朵)

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by why小海 at 2019-03-26 19:51:16
谢谢,这个模型是方程的实数部分,如果加上虚数部分,就出错了。还有请问为什么改变初始值,结果也会变呢?谢谢
...

这个很正常呐,将每次算出的结果作为初始值,多试几次初始值,结果相对稳定就可以。

发自小木虫Android客户端
4楼2019-03-26 20:19:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

独孤神宇

版主 (知名作家)



小木虫: 金币+0.5, 给个红包,谢谢回帖
模型和数据不匹配,检查一下模型,或许书写哪里不对。
2楼2019-03-26 13:50:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

why小海

新虫 (初入文坛)


引用回帖:
2楼: Originally posted by 独孤神宇 at 2019-03-26 13:50:09
模型和数据不匹配,检查一下模型,或许书写哪里不对。

谢谢,这个模型是方程的实数部分,如果加上虚数部分,就出错了。还有请问为什么改变初始值,结果也会变呢?谢谢

发自小木虫IOS客户端
3楼2019-03-26 19:51:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

why小海

新虫 (初入文坛)


送红花一朵
谢谢,我试试^_^

发自小木虫IOS客户端
5楼2019-03-27 20:41:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料085600 303求调剂 +5 1bygone 2026-03-04 5/250 2026-03-04 16:34 by yuningshan
[考研] 085600材料与化工 298 调剂 +7 小西笑嘻嘻 2026-03-03 7/350 2026-03-04 16:31 by yuningshan
[考研] 0703 东华大学 理学化学均过 a 线,总分 281 求调剂 +3 香香 Lu 2026-03-03 5/250 2026-03-04 16:16 by zhukairuo
[考研] 301求调剂 +3 李LJR 2026-03-04 3/150 2026-03-04 16:15 by sslc1985
[考研] 295求调剂 +6 等春来, 2026-03-04 6/300 2026-03-04 16:11 by zhukairuo
[考研] 281求调剂 +3 我是小小葱葱 2026-03-03 5/250 2026-03-04 14:23 by kakakapanpan
[考研] 085600 英一数二272求调剂 5+6 vida_a 2026-03-01 47/2350 2026-03-04 13:35 by vida_a
[考研] 281求调剂 +3 Y?l?h 2026-03-04 3/150 2026-03-04 10:51 by 汪!?!
[考研] 291求调剂 +4 Afy123456 2026-03-03 7/350 2026-03-03 22:29 by barlinike
[考研] 276求调剂 +8 路lyh123 2026-02-28 10/500 2026-03-03 18:25 by xin吖
[考研] 成绩276,专业代码0856求调剂 +7 小陈朵 2026-03-03 7/350 2026-03-03 15:58 by tgxtgxtgx9
[考研] 计算机学硕 分数285求调剂 +4 glwshine 2026-03-02 5/250 2026-03-03 14:27 by king呀
[考研] 环境工程专硕307 一志愿211 四六级已过 求调剂 +3 ccc! 2026-03-03 3/150 2026-03-03 12:55 by 热情沙漠
[考研] 26考研报考西工大材料308分求调剂 +4 weizhong123 2026-03-01 5/250 2026-03-03 12:22 by weizhong123
[考研] 281电子信息求调剂 +3 jhtfeybgj 2026-03-02 6/300 2026-03-03 12:12 by 热情沙漠
[考研] 材料工程求调剂 +3 1431251 2026-03-03 3/150 2026-03-03 11:58 by EBSD
[考研] 284求调剂 +10 天下熯 2026-02-28 11/550 2026-03-02 11:03 by 无际的草原
[考研] 化工专硕342,一志愿大连理工大学,求调剂 +6 kyf化工 2026-02-28 7/350 2026-03-02 10:56 by 无际的草原
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
信息提示
请填处理意见