版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

why小海

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 131.5
帖子: 33
在線: 8.1小時
蟲號: 14670528

[交流] 求助各位大神，參數(shù)擬合結(jié)果不對，金幣全送。謝謝

求助各位大神，用下面的程序做的參數(shù)擬合，擬合結(jié)果和圖形都不對，希望大神們幫幫忙，金幣全送。謝謝

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擬合結(jié)果;\n')
fprintf('\n\t參數(shù) E1 = %.16f',k(1))
fprintf('\n\t參數(shù) E2 = %.16f',k(2))
fprintf('\n\t參數(shù) E3 = %.16f',k(3))
fprintf('\n\t參數(shù) 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);

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

why小海

» 猜你喜歡

2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有299人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有26人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有34人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2019-03-25 15:50:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

獨孤神宇

版主 (知名作家)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

模型和數(shù)據(jù)不匹配，檢查一下模型，或許書寫哪里不對。

贊一下

回復(fù)此樓

2樓2019-03-26 13:50:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

why小海

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 131.5
帖子: 33
在線: 8.1小時
蟲號: 14670528

引用回帖:

2樓: Originally posted by 獨孤神宇 at 2019-03-26 13:50:09
模型和數(shù)據(jù)不匹配，檢查一下模型，或許書寫哪里不對。

謝謝，這個模型是方程的實數(shù)部分，如果加上虛數(shù)部分，就出錯了。還有請問為什么改變初始值，結(jié)果也會變呢？謝謝

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

3樓2019-03-26 19:51:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

獨孤神宇

版主 (知名作家)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

3樓: Originally posted by why小海 at 2019-03-26 19:51:16
謝謝，這個模型是方程的實數(shù)部分，如果加上虛數(shù)部分，就出錯了。還有請問為什么改變初始值，結(jié)果也會變呢？謝謝
...

這個很正常吶，將每次算出的結(jié)果作為初始值，多試幾次初始值，結(jié)果相對穩(wěn)定就可以。

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

4樓2019-03-26 20:19:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

why小海

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 131.5
帖子: 33
在線: 8.1小時
蟲號: 14670528

送紅花一朵

謝謝，我試試^_^

發(fā)自小木蟲IOS客戶端

回復(fù)此樓

5樓2019-03-27 20:41:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

why小海

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 131.5
帖子: 33
在線: 8.1小時
蟲號: 14670528

請問介電常數(shù)擬合參數(shù)是Cole-Cole模型的實部還是整個函數(shù)？謝謝

贊一下

回復(fù)此樓

6樓2019-03-28 11:22:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 why小海的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3	脫穎而出 2026-03-16	6/300	2026-03-17 21:19 by z1z2z3879
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 312求調(diào)劑 +4	陌宸希 2026-03-16	5/250	2026-03-17 17:09 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 工科材料085601 279求調(diào)劑 +3	困于星晨 2026-03-17	3/150	2026-03-17 14:08 by ms629
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8	秋有木北 2026-03-14	8/400	2026-03-16 15:21 by 哦哦123
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 308求調(diào)劑 +5	是Lupa啊 2026-03-11	5/250	2026-03-13 22:13 by JourneyLucky
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 307求調(diào)劑 +5	超級伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[論文投稿] 投稿問題 5+4	星光燦爛xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

why小海

[交流] 求助各位大神，參數(shù)擬合結(jié)果不對，金幣全送。謝謝

» 本帖已獲得的紅花（最新10朵）

» 猜你喜歡

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

獨孤神宇

why小海

獨孤神宇

why小海

why小海

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[交流] 求助各位大神，參數(shù)擬合結(jié)果不對，金幣全送。謝謝

» 搶金幣啦！回帖就可以得到:

查看全部散金貼