版塊導(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í)熱門(mén)版塊排行榜

返回列表

why小海

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 131.5
帖子: 33
在線: 8.1小時(shí)
蟲(chóng)號(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)有73人回復(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寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有24人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有6人回復(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í)
蟲(chóng)號(hào): 3522474

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

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

贊一下

回復(fù)此樓

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

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

why小海

新蟲(chóng) (初入文壇)

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

引用回帖:

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

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

發(fā)自小木蟲(chóng)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í)
蟲(chóng)號(hào): 3522474

★
小木蟲(chóng): 金幣+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ā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

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

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

why小海

新蟲(chóng) (初入文壇)

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

送紅花一朵

謝謝，我試試^_^

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

回復(fù)此樓

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

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

why小海

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 131.5
帖子: 33
在線: 8.1小時(shí)
蟲(chóng)號(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ā)表

[考研] 301求調(diào)劑 +3	yy要上岸呀 2026-03-17	3/150	2026-03-17 14:04 by laoshidan
[考研] 268求調(diào)劑 +7	一定有學(xué)上- 2026-03-14	8/400	2026-03-17 13:10 by zz820
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 材料專碩326求調(diào)劑 +5	墨煜姒莘 2026-03-15	5/250	2026-03-16 21:30 by 木瓜膏
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475 +5	Becho359 2026-03-11	5/250	2026-03-14 11:35 by 哦哦123
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 材料工程,326分，求調(diào)劑 +6	KRSLSR 2026-03-10	6/300	2026-03-13 23:47 by JourneyLucky
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 求調(diào)劑（材料與化工327） +4	愛(ài)吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過(guò) 本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考研] 420求調(diào)劑 +4	莫向外求11 2026-03-10	6/300	2026-03-12 14:41 by ruiyingmiao
[考研] 大連大學(xué)化學(xué)專業(yè)研究生調(diào)劑 +3	琪久. 2026-03-10	8/400	2026-03-11 10:02 by 琪久.