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

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

24小時(shí)熱門版塊排行榜

返回列表

why小海

新蟲 (初入文壇)

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

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

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

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('計(jì)算值','實(shí)驗(yàn)值','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)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有228人回復(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ù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有23人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)

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

查看全部散金貼

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

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474

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

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

贊一下

回復(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小時(shí)
蟲號(hào): 14670528

引用回帖:

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

謝謝，這個(gè)模型是方程的實(shí)數(shù)部分，如果加上虛數(shù)部分，就出錯(cuò)了。還有請(qǐng)問(wèn)為什么改變初始值，結(jié)果也會(huì)變呢？謝謝

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

贊一下

回復(fù)此樓

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

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474

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

引用回帖:

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

這個(gè)很正常吶，將每次算出的結(jié)果作為初始值，多試幾次初始值，結(jié)果相對(duì)穩(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小時(shí)
蟲號(hào): 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小時(shí)
蟲號(hào): 14670528

請(qǐng)問(wèn)介電常數(shù)擬合參數(shù)是Cole-Cole模型的實(shí)部還是整個(gè)函數(shù)？謝謝

贊一下

回復(fù)此樓

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

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

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

返回列表

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

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

[考研] 271求調(diào)劑 +12	生如夏花… 2026-03-11	14/700	2026-03-17 10:56 by lovewei0727
[考研] 290求調(diào)劑 +6	@將就將就看 2026-03-10	11/550	2026-03-17 10:22 by @將就將就看
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 283求調(diào)劑 +3	聽(tīng)風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 085601求調(diào)劑 +3	Du.11 2026-03-16	3/150	2026-03-16 20:42 by 無(wú)際的草原
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過(guò)，有科研經(jīng)歷 +7	曦熙兮 2026-03-15	7/350	2026-03-16 16:34 by houyaoxu
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 311求調(diào)劑 +5	牛乳糖的卡卡 2026-03-10	5/250	2026-03-14 00:05 by JourneyLucky
[考研] 311求調(diào)劑 +8	zchqwer 2026-03-10	8/400	2026-03-14 00:01 by JourneyLucky
[考研] 材料工程,326分，求調(diào)劑 +6	KRSLSR 2026-03-10	6/300	2026-03-13 23:47 by JourneyLucky
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12	11啦啦啦 2026-03-11	35/1750	2026-03-13 21:25 by JourneyLucky
[考研] 【考研調(diào)劑求收留】 +3	Ceciilia 2026-03-11	3/150	2026-03-13 20:18 by JourneyLucky
[考研] 301求調(diào)劑 +6	Liyouyumairs 2026-03-11	6/300	2026-03-13 20:11 by JourneyLucky