xiaoqiia

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 730.9
散金: 60
紅花: 3
帖子: 79
在線: 106.8小時
蟲號: 2603013
注冊: 2013-08-18
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

[交流] 求大神幫忙matlab擬合函數(shù)求參數(shù)

如題，x= 1 7 14 21 35
y=17.24 16.31 14.21 13.16 11.87
方程為y=18+a[1-exp(-bx)]-cln(1+dx)，求參數(shù)abcd，要擬合的函數(shù)過程啊，在下第一次接觸matlab，完全新手~求大神幫忙，金幣不多，愿相送50個。

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有64人回復(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)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有73人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

拜求！matlab擬合動力學(xué)參數(shù) 已經(jīng)有14人回復(fù)
求助用matlab擬合動力學(xué)方程已經(jīng)有9人回復(fù)
如何用已知函數(shù)去擬合實驗數(shù)據(jù)，得出參數(shù) 已經(jīng)有10人回復(fù)
matlab 擬合求助已經(jīng)有8人回復(fù)
關(guān)于matlab擬合已經(jīng)有4人回復(fù)
用matlab里的nlinfit如何擬合分段函數(shù) 已經(jīng)有8人回復(fù)
用Matlab做多變量非線性擬合求參數(shù)值，求各位大俠鼎力相助！已經(jīng)有3人回復(fù)
Matlab6.5精簡版可擬合數(shù)據(jù)繪圖已經(jīng)有7人回復(fù)
MATLAB微分方程參數(shù)擬合問題，求大神已經(jīng)有7人回復(fù)
matlab數(shù)據(jù)擬合希望高手幫忙寫一下matlab的程序跪謝已經(jīng)有12人回復(fù)
matlab數(shù)據(jù)擬合問題求助已經(jīng)有7人回復(fù)
Matlab多項式擬合已經(jīng)有6人回復(fù)
求matlab畫一個函數(shù)的圖已經(jīng)有6人回復(fù)
matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

matlab擬合方程參數(shù)時初值的選擇已經(jīng)有15人回復(fù)
matlab擬合求值已經(jīng)有17人回復(fù)
需要用matlab擬合數(shù)據(jù)的外行人急需幫助已經(jīng)有10人回復(fù)
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復(fù)
求助matlab隱函數(shù)參數(shù)擬合已經(jīng)有16人回復(fù)
【求助】matlab曲面擬合表達(dá)式已經(jīng)有10人回復(fù)
【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過一個點已經(jīng)有8人回復(fù)

奔跑吧少年

1樓 2014-04-19 22:01:18

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

qq2537742567

新蟲 (小有名氣)

應(yīng)助: 46 (小學(xué)生)
金幣: 164.5
帖子: 93
在線: 43.5小時
蟲號: 3059261
注冊: 2014-03-16

★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
xiaoqiia: 金幣+2, 能來個過程嗎？小弟剛接觸，不大懂~ 2014-04-20 10:02:55
fegg7502: 金幣+1, 應(yīng)助指數(shù)+1, 鼓勵交流 2014-04-21 07:50:36

lsqcurvefit你參考這個函數(shù)來做

贊一下

回復(fù)此樓

2樓2014-04-20 09:45:09

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

NicoleLas

木蟲 (正式寫手)

應(yīng)助: 66 (初中生)
金幣: 4288.7
紅花: 26
帖子: 514
在線: 177.2小時
蟲號: 2872321
注冊: 2013-12-15
性別: GG
專業(yè): 計算機(jī)軟件

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
xiaoqiia: 金幣+10, 用了為什么參數(shù)不能顯示d呢？只有abc跟non。。求大神繼續(xù)指導(dǎo)下 2014-04-20 20:53:33
fegg7502: 金幣+1, 應(yīng)助指數(shù)+1, 鼓勵交流 2014-04-21 07:50:45

用1stOpt計算挺好的，結(jié)果好像不止一種：
代碼為：
Title "Curve";
Parameters a,b,c,d;
Variable x,y;
Function y=18+a*(1-exp(-b*x))-c*log(1+d*x);
Data;
1 17.24
7 16.31
14 14.21
21 13.16
35 11.87

贊一下

回復(fù)此樓

Be water, my friend.

3樓2014-04-20 12:56:42

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

xiaoqiia

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 730.9
散金: 60
紅花: 3
帖子: 79
在線: 106.8小時
蟲號: 2603013
注冊: 2013-08-18
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

引用回帖:

3樓: Originally posted by NicoleLas at 2014-04-20 12:56:42
用1stOpt計算挺好的，結(jié)果好像不止一種：
代碼為：
Title "Curve";
Parameters a,b,c,d;
Variable x,y;
Function y=18+a*(1-exp(-b*x))-c*log(1+d*x);
Data;
1 17.24
7 16.31
14 14. ...

是LN不是LOG~

回復(fù)此樓

奔跑吧少年

4樓2014-04-20 20:54:35

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

mzhyan

至尊木蟲 (文壇精英)

應(yīng)助: 712 (博后)
金幣: 37461.7
紅花: 24
帖子: 45756
在線: 287.3小時
蟲號: 1191620
注冊: 2011-01-18
專業(yè): 環(huán)境污染化學(xué)

★
xiaoqiia: 金幣+1 2014-04-20 21:53:25

blessing

回復(fù)此樓

5樓2014-04-20 21:45:12

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

亦岳9863

6樓2014-04-20 22:17:26

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

NicoleLas

木蟲 (正式寫手)

應(yīng)助: 66 (初中生)
金幣: 4288.7
紅花: 26
帖子: 514
在線: 177.2小時
蟲號: 2872321
注冊: 2013-12-15
性別: GG
專業(yè): 計算機(jī)軟件

★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
fegg7502: 金幣+1, 鼓勵交流 2014-04-21 07:50:53
fegg7502: 應(yīng)助指數(shù)+1, 鼓勵交流 2014-04-21 07:51:03

比如這是一組解。
"Curve"

迭代數(shù): 66
計算用時(時:分:秒:微秒): 00:00:01:60
優(yōu)化算法: 準(zhǔn)牛頓法(BFGS) + 通用全局優(yōu)化法
計算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn)
均方差(RMSE): 0.0720037079046967
殘差平方和(SSE): 0.0259226697601244
相關(guān)系數(shù)(R): 0.999337253425239
相關(guān)系數(shù)之平方(R^2): 0.9986749460835
決定系數(shù)(DC): 0.998674945880788
卡方系數(shù)(Chi-Square): 0.000930938251442266
F統(tǒng)計(F-Statistic): 2261.05881499843

參數(shù) 最佳估算
---------- -------------
a       31.1225668620202
b       1.96505984749032
c       2.75886155053898
d       21490.6088747969

====== 結(jié)果輸出 =====

No 實測值y 計算值y
1 17.24 17.2400202
2 16.31 16.2333504
3 14.21 14.3210955
4 13.16 13.2024765
5 11.87 11.7931817

不好意思，換成ln了，以前在使用其它軟件時，LOG就是LN，而LG為LOG10()。
我沒有1stOpt正版的，你最好求助一下dingd版主。

贊一下

回復(fù)此樓

Be water, my friend.

7樓2014-04-20 23:00:35

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

waterque

金蟲 (正式寫手)

應(yīng)助: 59 (初中生)
金幣: 1989.1
散金: 96
紅花: 10
帖子: 759
在線: 191.7小時
蟲號: 2095273
注冊: 2012-10-30
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
xiaoqiia: 金幣+5, 嗯，謝謝 2014-04-29 10:01:05

matlab 擬合還是很強(qiáng)大的，有一個curve fitting tool很強(qiáng)大。建議使用
輸入命令cftool即彈出擬合窗口，設(shè)置擬合的參數(shù)。
你這個數(shù)據(jù)點太少，擬合有一定的不確定性，需要根據(jù)實際進(jìn)行分析對比。

贊一下(1人)

回復(fù)此樓

行自己的路：waterque@hotmail.com

8樓2014-04-21 11:28:40

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

spiderone

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 3474.1
散金: 12
紅花: 1
帖子: 157
在線: 108.4小時
蟲號: 2520847
注冊: 2013-06-25
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
xiaoqiia: 金幣+10, 非常感謝 2014-04-29 09:59:08

function jisuan
clear all
clc
x=[1 7 14 21 35]';  % x數(shù)據(jù)
yexp=[17.24 16.31 14.21 13.16 11.87]';  % y數(shù)據(jù)
lb=[];
ub=[];
% 非線性擬合
k0=[26.1 1.8 2.8 3429.0];  %  參數(shù)初始值
[k,resnorm,residual,exitflag] = ...
lsqnonlin(@f,k0,[],[],[],x,yexp)

y=18+k(1).*(1-exp(-k(2).*x))-k(3).*log(1+k(4).*x)  % 計算結(jié)果

% ------------------------------------------------------------------
function y=f(k,x,yexp)
y=18+k(1).*(1-exp(-k(2).*x))-k(3).*log(1+k(4).*x)-yexp;

%--------------------------------------------------------------------
計算結(jié)果如下：
k =

  1.0e+003 *

0.0261 0.0018 0.0028 3.4290

resnorm =

0.0259

residual =

0.0000
-0.0767
0.1111
0.0425
-0.0769
計算值y =

17.2400
16.2333
14.3211
13.2025
11.7931

贊一下(2人)

回復(fù)此樓

9樓2014-04-22 19:21:30

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

sjmr1221

10樓2014-04-22 19:51:56

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

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

[考研] 307求調(diào)劑 +3	余意卿 2026-03-21	3/150	2026-03-23 10:32 by Iveryant
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚.. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 08工科 320總分求調(diào)劑 +11	梨花珞晚風(fēng) 2026-03-17	11/550	2026-03-22 17:42 by luoyongfeng
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 305分求調(diào)劑（食品工程） +4	Sxy112 2026-03-21	6/300	2026-03-22 15:26 by 無懈可擊111
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	4/200	2026-03-22 10:25 by QGZDSYS
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 求調(diào)劑 +4	要好好無聊 2026-03-21	4/200	2026-03-21 18:57 by 學(xué)員8dgXkO
[基金申請] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 268求調(diào)劑 +9	簡單點0 2026-03-17	9/450	2026-03-21 15:37 by lature00
[考研] 299求調(diào)劑 +6	△小透明* 2026-03-17	6/300	2026-03-21 02:42 by JourneyLucky
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 274求調(diào)劑 +10	S.H1 2026-03-18	10/500	2026-03-20 23:51 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無懈可擊111
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero