zhangxiaodao

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1929.4
散金: 31
紅花: 1
帖子: 228
在線: 286小時
蟲號: 871428
注冊: 2009-10-13

[求助] 根據(jù)數(shù)據(jù)和方程的形式，求方程中的參數(shù)

RT，數(shù)據(jù)如下列，然后方程的形式為y=a*exp[-p/(x+q)],求方程中的p和q，謝謝各位大俠，急求啊~~
x=1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
2
2.1
2.2
2.3
2.4
y=10363972.36
5697599.795
3269017.372
1948866.775
1202604.284
765578.5464
501320.0508
336809.4864
231640.3371
162754.7914
116618.904
85080.3473
63109.81473
47535.43773

回復(fù)此樓

» 猜你喜歡

085600材料與化工求調(diào)劑已經(jīng)有16人回復(fù)
317求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有7人回復(fù)
265求調(diào)劑已經(jīng)有5人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有3人回復(fù)
求調(diào)劑已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)

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

【500金幣本版活動】“經(jīng)驗力場分子模擬+第一性原理計算問題討論互助及資源共享” 已經(jīng)有82人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
求非線性方程中的參數(shù)問題已經(jīng)有13人回復(fù)

1樓 2011-12-19 22:50:41

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

715211229

木蟲 (正式寫手)

應(yīng)助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時
蟲號: 1357002
注冊: 2011-07-27
性別: GG
專業(yè): 傳動機(jī)械學(xué)

★ ★
臭水溝(金幣+2): 謝謝交流~~ 2011-12-20 11:03:15

matlab ：ctool工具箱
%%%擬合效果很差

Fit computation did not converge:
Maximum number of function evaluations exceeded. Increasing
MaxFunEvals (in fit options) may allow for a better fit, or
the current equation may not be a good model for the data.

Fit found when optimization terminated:

General model:
   f(x) = a*exp(-p/(x+q))
Coefficients (with 95% confidence bounds):
   a =    11.28  (-2.505e+004, 2.507e+004)
   p =    -7.497  (-3191, 3176)
   q =    -0.4499  (-151.6, 150.7)

Goodness of fit:
  SSE: 1.309e+014
  R-square: -0.159
  Adjusted R-square: -0.3697
  RMSE: 3.45e+006

贊一下(1人)

回復(fù)此樓

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

zhangxiaodao

我是蝸牛

2樓2011-12-20 05:28:41

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

715211229

木蟲 (正式寫手)

應(yīng)助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時
蟲號: 1357002
注冊: 2011-07-27
性別: GG
專業(yè): 傳動機(jī)械學(xué)

cftool工具箱

回復(fù)此樓

我是蝸牛

3樓2011-12-20 05:29:19

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

zhangxiaodao

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1929.4
散金: 31
紅花: 1
帖子: 228
在線: 286小時
蟲號: 871428
注冊: 2009-10-13

送鮮花一朵

引用回帖:

樓: Originally posted by 715211229 at 2011-12-20 05:28:41:
matlab ：ctool工具箱
%%%擬合效果很差

Fit computation did not converge:
Maximum number of function evaluations exceeded. Increasing
MaxFunEvals (in fit options) may allow ...

謝謝

回復(fù)此樓

4樓2011-12-20 09:00:53

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

715211229

木蟲 (正式寫手)

應(yīng)助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時
蟲號: 1357002
注冊: 2011-07-27
性別: GG
專業(yè): 傳動機(jī)械學(xué)

★
臭水溝(金幣+1): 謝謝交流~~ 2011-12-20 11:03:26

你看下這個帖子，我這個可能是錯誤的
http://www.gaoyang168.com/bbs/viewthread.php?tid=3866180

贊一下(1人)

回復(fù)此樓

我是蝸牛

5樓2011-12-20 10:48:24

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

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

★
感謝參與，應(yīng)助指數(shù) +1
zhangxiaodao(金幣+20): 2011-12-20 13:52:42
xiegangmai(金幣+1): 謝謝應(yīng)助！ 2011-12-20 23:02:49

General model:
   f(x) = a*exp(-p/(x+q))
Coefficients (with 95% confidence bounds):
   a =          1  (1, 1)
   p =       -42  (-42, -42)
   q =       1.5  (1.5, 1.5)

Goodness of fit:
  SSE: 1.767e-007
  R-square: 1
  Adjusted R-square: 1
  RMSE: 0.0001267

初值a=10000 p=-2  q=-1

贊一下(1人)

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

6樓2011-12-20 11:22:06

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

★
xiegangmai(金幣+1): 辛苦了！ 2011-12-20 23:03:02

引用回帖:

5樓: Originally posted by 715211229 at 2011-12-20 10:48:24:
你看下這個帖子，我這個可能是錯誤的
http://www.gaoyang168.com/bbs/viewthread.php?tid=3866180

圖如下

贊一下(1人)

回復(fù)此樓

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

zhangxiaodao

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

7樓2011-12-20 11:23:32

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

zhangxiaodao

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1929.4
散金: 31
紅花: 1
帖子: 228
在線: 286小時
蟲號: 871428
注冊: 2009-10-13

送鮮花一朵

引用回帖:

樓: Originally posted by dbb627 at 2011-12-20 11:23:32:
圖如下

非常感謝，您擬合的太棒了~~
ps: 怎么送金幣呢？然后能講程序發(fā)給我嗎，我是MATLAB門外漢。。。。

贊一下

回復(fù)此樓

8樓2011-12-20 13:48:57

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

CODE:

 x=[1.1

1.2

1.3

1.4

1.5

1.6

1.7

1.8

1.9

2

2.1

2.2

2.3

2.4];

y=[10363972.36

5697599.795

3269017.372

1948866.775

1202604.284

765578.5464

501320.0508

336809.4864

231640.3371

162754.7914

116618.904

85080.3473

63109.81473

47535.43773];

st_ = [1000 -2 -1];

ft_ = fittype('a*exp(-p/(x+q))','dependent',{'y'},'independent',{'x'},'coefficients',{'a', 'p','q'});

[cf_,good]= fit(x,y,ft_ ,'Startpoint',st_)

h_ = plot(cf_,'fit',0.95);

legend off;  % turn off legend from plot method call

set(h_(1),'Color',[1 0 0],...

     'LineStyle','-', 'LineWidth',2,...

     'Marker','none', 'MarkerSize',6);

hold on,plot(x,y,'*')

cf_ =

   General model:
   cf_(x) = a*exp(-p/(x+q))
   Coefficients (with 95% confidence bounds):
   a =          1  (1, 1)
   p =       -42  (-42, -42)
   q =       1.5  (1.5, 1.5)

good =

         sse: 1.7672e-007
   rsquare: 1
         dfe: 11
adjrsquare: 1
      rmse: 1.2675e-004

贊一下

回復(fù)此樓

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

zhangxiaodao

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

9樓2011-12-20 14:38:54

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

zhangxiaodao

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1929.4
散金: 31
紅花: 1
帖子: 228
在線: 286小時
蟲號: 871428
注冊: 2009-10-13

送鮮花一朵

引用回帖:

樓: Originally posted by dbb627 at 2011-12-20 14:38:54:
[code]
x=[1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
2
2.1
2.2
2.3
2.4];
y=[10363972.36
5697599.795
3269017.372
1948866.775
1202604.284
765578.5464
501320.0508
336809.4864
231640.3 ...

非常感謝dbb627版主，可能您也知道我問您要源程序的目的了，我手里還有一些數(shù)據(jù)，是跟上面一樣的方程，希望得到a，p，q值，請問能直接套用這個程序嗎，或者需要做哪些調(diào)整，謝謝您。

贊一下

回復(fù)此樓

10樓2011-12-20 15:41:40

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

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

[考研] 279分求調(diào)劑一志愿211 +3	chaojifeixia 2026-03-19	3/150	2026-03-20 05:47 by laoshidan
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 296求調(diào)劑 +3	www_q 2026-03-18	6/300	2026-03-19 22:28 by zhq0425
[考研] 321求調(diào)劑 +8	何潤采123 2026-03-18	10/500	2026-03-19 16:46 by 何潤采123
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +5	吃吃吃才有意義 2026-03-19	5/250	2026-03-19 16:18 by 30660438
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	10/500	2026-03-19 16:05 by 余麟余
[考研] 本人考085602 化學(xué)工程專碩 +17	不知道叫什么！ 2026-03-15	19/950	2026-03-19 15:06 by 盡舜堯1
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 求調(diào)劑 +3	Mqqqqqq 2026-03-19	3/150	2026-03-19 14:11 by peike
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[基金申請] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky