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

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

Ronny_chou

木蟲 (正式寫手)

應(yīng)助: 43 (小學(xué)生)
金幣: 3031.5
散金: 16
紅花: 14
帖子: 485
在線: 313小時(shí)
蟲號(hào): 2537599
注冊(cè): 2013-07-08
專業(yè): 分離過程

[求助] Origin自定義模擬與matlab中cftool工具箱擬合結(jié)果比較已有1人參與

P=[0.5070 1.318 2.772 4.338 6.241 7.839 9.849 14.66 19.35 24.56 29.93  34.13 39.04  46.20] ;
n=[0.4429 1.085 1.901 2.533 3.162 3.697 4.209 5.052 6.043 6.754 7.489 8.279 9.184 10.16];
想要擬合成曲線f(x) = nmax*(b*x)^(1/k)/(1+(b*x)^(1/k)) 與f(x) = nmax*x/(1+(a*x)^k)^(1/k) 的形式，然后我
通過Origin中的自定義曲線擬合與matlab中cftool工具箱模擬結(jié)果，比較如下：
（1）  matlab中cftool工具箱模擬結(jié)果
1.1. General model:
   f(x) = nmax*(b*x)^(1/k)/(1+(b*x)^(1/k))
Coefficients (with 95% confidence bounds):
   b = 0.0001122  (-0.001508, 0.001732)
   k =    1.658  (1.348, 1.968)
   nmax =    248.9  (-1699, 2197)

Goodness of fit:
  SSE: 0.3045
  R-square: 0.9975
  Adjusted R-square: 0.9971
  RMSE: 0.1664

1.2 General model:
General model:
   f(x) = nmax*a*x/(1+(a*x)^k)^(1/k)
Coefficients (with 95% confidence bounds):
   a = 0.001947  (-0.01345, 0.01735)
   k =    0.1019  (-0.1258, 0.3297)
   nmax = 3.203e+04  (-5.908e+05, 6.549e+05)

Goodness of fit:
  SSE: 0.3006
  R-square: 0.9976
  Adjusted R-square: 0.9971
  RMSE: 0.1653

但是用origin模擬結(jié)果如下：
（1）    f(x) = nmax*(b*x)^(1/k)/(1+(b*x)^(1/k)) 模擬結(jié)果

（2） f(x) = nmax*x/(1+(a*x)^k)^(1/k) 模擬結(jié)果

Origin自定義模擬與matlab中cftool工具箱擬合結(jié)果比較-1

用兩者所模擬出來的nmax是有區(qū)別的，而nmax在我這里是有物理意義的，因此在收斂時(shí)不管用什么方法所模擬出來的結(jié)果應(yīng)該是一樣的。朋友們幫忙看看？我發(fā)現(xiàn)如果曲線擬合收斂比較困難時(shí)，調(diào)節(jié)初始化參數(shù)會(huì)對(duì)結(jié)果影響比較大。

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : simulation2.jpg

2015-09-23 11:32:52, 666.02 K

附件 2 : simulation2.jpg

2015-09-23 11:34:19, 666.02 K

» 猜你喜歡

計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有161人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有29人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有26人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名已經(jīng)有0人回復(fù)

1樓 2015-09-23 11:49:39

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

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

引用回帖:

3樓: Originally posted by Ronny_chou at 2015-09-24 10:32:53
謝謝，我想請(qǐng)問下用1stOpt求解出來的是沒有設(shè)初始值的，然后求解出來的結(jié)果穩(wěn)定唯一怎么理解？因?yàn)橛胢atlab用您簡(jiǎn)化的那個(gè)函數(shù)一開始是收斂的，當(dāng)我改變初始值時(shí)就會(huì)出現(xiàn)不一樣的結(jié)果，即使最后通過固定某一個(gè)值達(dá) ...

論優(yōu)化擬合，Matlab和Origin都無法與1stOpt相比。在1stOpt中，參數(shù)范圍限制也是很容易的。

贊一下

回復(fù)此樓

4樓2015-09-24 14:02:21

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

查看全部 5 個(gè)回答

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
Ronny_chou: 金幣+20, ★★★★★最佳答案 2015-10-29 10:03:09

1：
第一個(gè)公式：y = nmax*(b*x)^(1/k)/(1+(b*x)^(1/k))，不能簡(jiǎn)化下嗎，完全等同于：

y = nmax/(1+(b*x)^k);，看著都清爽些。直接擬合這個(gè)公式的話，目標(biāo)函數(shù)（SSE）基本一樣但nmax和b不穩(wěn)定：
1）
均方差(RMSE):0.146688964176323
殘差平方和(SSE):0.301247130955714
相關(guān)系數(shù)(R): 0.998801066757962
相關(guān)系數(shù)之平方(R^2): 0.997603570956843
確定系數(shù)(DC): 0.997545478102313
卡方系數(shù)(Chi-Square): 0.085846619465631
F統(tǒng)計(jì)(F-Statistic): 2240.76224587093

參數(shù) 最佳估算
---------- -------------
b 6.38790682352766E-17
nmax 3715236827.67985
k -0.589760207223845

2）
均方差(RMSE):0.146688964157574
殘差平方和(SSE):0.301247130878707
相關(guān)系數(shù)(R): 0.998801066740756
相關(guān)系數(shù)之平方(R^2): 0.997603570922471
確定系數(shù)(DC): 0.997545478102941
卡方系數(shù)(Chi-Square): 0.0858466176890851
F統(tǒng)計(jì)(F-Statistic): 2240.76224644373

參數(shù) 最佳估算
---------- -------------
b 2.37230022879993E-20
nmax 391719737466.736
k -0.589760208502794

應(yīng)該是方程的形式造成這種結(jié)局，如果不考慮參數(shù)物理意義，其實(shí)公式：y = nmax/(1+(b*x)^k); --》y = nmax/(1+b^k*x^k); -->y = nmax/(1+b1*x^k);

此時(shí)穩(wěn)定唯一：
均方差(RMSE):0.146648461643796
殘差平方和(SSE):0.301080798234885
相關(guān)系數(shù)(R): 0.998804382533961
相關(guān)系數(shù)之平方(R^2): 0.997610194569047
確定系數(shù)(DC): 0.997546833359388
卡方系數(shù)(Chi-Square): 0.0881571807051252
F統(tǒng)計(jì)(F-Statistic): 2242.00015969098

參數(shù) 最佳估算
---------- -------------
nmax -863.166782962704
b -823.26708610708
k -0.586003926097037

2：第二個(gè)式子： y = nmax*a*x/(1+(a*x)^k)^(1/k)

1stOpt求解穩(wěn)定唯一，Matlab誤差太大，Origin也有差距：

均方差(RMSE):0.144604300178649
殘差平方和(SSE):0.292745650822194
相關(guān)系數(shù)(R): 0.998819337279252
相關(guān)系數(shù)之平方(R^2): 0.997640068522964
確定系數(shù)(DC): 0.997614747041354
卡方系數(shù)(Chi-Square): 0.06651177594336
F統(tǒng)計(jì)(F-Statistic): 2305.83510233767

參數(shù) 最佳估算
---------- -------------
nmax 185924116.646796
a 5.71405015369486E-5
k 0.0511661433683439

贊一下

回復(fù)此樓

2樓2015-09-23 22:21:18

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

Ronny_chou

木蟲 (正式寫手)

應(yīng)助: 43 (小學(xué)生)
金幣: 3031.5
散金: 16
紅花: 14
帖子: 485
在線: 313小時(shí)
蟲號(hào): 2537599
注冊(cè): 2013-07-08
專業(yè): 分離過程

引用回帖:

2樓: Originally posted by dingd at 2015-09-23 22:21:18
1：
第一個(gè)公式：y = nmax*(b*x)^(1/k)/(1+(b*x)^(1/k))，不能簡(jiǎn)化下嗎，完全等同于：

y = nmax/(1+(b*x)^k);，看著都清爽些。直接擬合這個(gè)公式的話，目標(biāo)函數(shù)（SSE）基本一樣但nmax和b不穩(wěn)定：
1）
均方差(R ...

謝謝，我想請(qǐng)問下用1stOpt求解出來的是沒有設(shè)初始值的，然后求解出來的結(jié)果穩(wěn)定唯一怎么理解？因?yàn)橛胢atlab用您簡(jiǎn)化的那個(gè)函數(shù)一開始是收斂的，當(dāng)我改變初始值時(shí)就會(huì)出現(xiàn)不一樣的結(jié)果，即使最后通過固定某一個(gè)值達(dá)到收斂，但跟1stOpt結(jié)果相差較大，那究竟哪一個(gè)具有說服力？此外1stOpt中能夠?qū)讉€(gè)系數(shù)都限制為大于零么？謝謝啦~

贊一下

回復(fù)此樓

3樓2015-09-24 10:32:53

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

Ronny_chou

木蟲 (正式寫手)

應(yīng)助: 43 (小學(xué)生)
金幣: 3031.5
散金: 16
紅花: 14
帖子: 485
在線: 313小時(shí)
蟲號(hào): 2537599
注冊(cè): 2013-07-08
專業(yè): 分離過程

引用回帖:

4樓: Originally posted by dingd at 2015-09-24 14:02:21
論優(yōu)化擬合，Matlab和Origin都無法與1stOpt相比。在1stOpt中，參數(shù)范圍限制也是很容易的。...

好的，謝謝，我再找1stOpt相關(guān)文檔看看~

贊一下

回復(fù)此樓

5樓2015-09-24 15:26:27

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

查看全部 5 個(gè)回答

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

[考研] 08工科求調(diào)劑290分 +5	1314捧花 2026-04-02	6/300	2026-04-02 12:43 by 見光死
[考研] 一志愿鄭大材料工程290求調(diào)劑 +19	Youth_ 2026-03-30	19/950	2026-04-02 11:21 by Sammy2
[考研] 07生物學(xué)求調(diào)劑一志愿同濟(jì)大學(xué)359分 +3	LAMC. 2026-03-30	3/150	2026-04-02 10:26 by 18828373951
[考研] 348環(huán)境工程調(diào)劑 +3	吳彥祖24k 2026-04-01	3/150	2026-04-02 09:14 by nanaliuyun
[考研] 805600專碩材料與化工348分求調(diào)劑 +6	上學(xué)啦！ 2026-04-01	6/300	2026-04-02 07:49 by 2026材料調(diào)劑
[考博] 26年申博 +3	staryer 2026-03-30	4/200	2026-04-01 23:21 by ai4pharm
[考研] 材料調(diào)劑 +14	一樣YWY 2026-04-01	14/700	2026-04-01 21:07 by lijunpoly
[考研] 285求調(diào)劑 +7	AZMK 2026-03-30	13/650	2026-04-01 17:00 by 七度不信任
[考研] 材料專業(yè)調(diào)劑 +5	啦啦啦哭 2026-03-31	6/300	2026-04-01 16:48 by JourneyLucky
[考研] 350求調(diào)劑 +7	阿佳～ 2026-03-31	7/350	2026-04-01 16:12 by yanflower7133
[考研] 求調(diào)劑 +5	零八# 2026-03-27	5/250	2026-04-01 14:40 by yulian1987
[考研] 263求調(diào)劑 +3	DDDDuu 2026-03-27	3/150	2026-03-31 16:21 by 土木碩士招生
[考研] 求調(diào)劑 +8	11ggg 2026-03-30	8/400	2026-03-31 13:56 by nanaliuyun
[考研] 一志愿浙江大學(xué)工科動(dòng)力工程370,數(shù)一121,專業(yè)課135，現(xiàn)在能去哪里 +3	080700調(diào)劑 2026-03-30	4/200	2026-03-31 12:00 by KLMY666
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-29	3/150	2026-03-31 10:06 by cal0306
[考研] 328求調(diào)劑 +8	嗯滴的基本都 2026-03-27	8/400	2026-03-30 17:20 by Wang200018
[考研] 071010 323 分求調(diào)劑 +3	Baekzhy 2026-03-27	3/150	2026-03-30 14:24 by andresqi
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4	qianbai11 2026-03-29	4/200	2026-03-30 11:34 by 唐沐兒
[考研] 086000生物與醫(yī)藥調(diào)劑 +5	Feisty。 2026-03-28	9/450	2026-03-29 12:02 by longlotian
[考研] 08開頭275求調(diào)劑 +4	拉誰不重要 2026-03-26	4/200	2026-03-27 14:12 by Delta2012

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

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

Ronny_chou

[求助] Origin自定義模擬與matlab中cftool工具箱擬合結(jié)果比較 已有1人參與

» 本帖附件資源列表

» 猜你喜歡

dingd

dingd

【答案】應(yīng)助回帖

Ronny_chou

Ronny_chou

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

[求助] Origin自定義模擬與matlab中cftool工具箱擬合結(jié)果比較已有1人參與