版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

1135725495

鐵桿木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 6868.7
散金: 18
紅花: 1
帖子: 1492
在線: 209.6小時
蟲號: 2756771
注冊: 2013-10-27
性別: MM
專業(yè): 有機化工

[求助] matlab 已有1人參與

用matlab軟件計算置信區(qū)間，算出來的結(jié)果過大怎么解決？？？

@月只藍

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有6人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有230人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學工程學碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有65人回復(fù)
華東師范大學芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

1樓 2017-04-11 09:38:07

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
1135725495: 金幣+5, ★★★很有幫助 2017-04-11 12:31:43

做參數(shù)擬合，首先先看決定系數(shù)（或相關(guān)系數(shù)）的高低，再看其他的。
置信區(qū)間跟設(shè)置的置信度也有關(guān)系，95%的置信度是最常用的，置信度越高，置信區(qū)間也越大。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2017-04-11 10:02:44

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

1135725495

鐵桿木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 6868.7
散金: 18
紅花: 1
帖子: 1492
在線: 209.6小時
蟲號: 2756771
注冊: 2013-10-27
性別: MM
專業(yè): 有機化工

引用回帖:

2樓: Originally posted by 月只藍 at 2017-04-11 10:02:44
做參數(shù)擬合，首先先看決定系數(shù)（或相關(guān)系數(shù)）的高低，再看其他的。
置信區(qū)間跟設(shè)置的置信度也有關(guān)系，95%的置信度是最常用的，置信度越高，置信區(qū)間也越大。

你好！以下是用matlab計算的結(jié)果,使用的ci = nlparci(k,residual,jacobian)函數(shù)，置信度為95%。算出來的決定系數(shù)在0.9以上呢。算出來的置信區(qū)間好大啊，不知問題出在哪？
ci =

1.0e+03 *

-0.1666 0.1667
-0.0000 0.0000
-0.1562 0.1562
-6.0471 6.0552

以fmincon()的結(jié)果為初值，使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:
k1 = 0.05881980024
k2 = 0.00000000000
k3 = 0.00396318725
k4 = 4.05955359944
  The sum of the squares is: 6.0e-06

output =
firstorderopt: 6.1534e-07
   iterations: 1
      funcCount: 10
   cgiterations: 0
      algorithm: 'trust-region-reflective'
      message: [1x425 char]

贊一下

回復(fù)此樓

3樓2017-04-11 11:16:48

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

引用回帖:

3樓: Originally posted by 1135725495 at 2017-04-11 11:16:48
你好！以下是用matlab計算的結(jié)果,使用的ci = nlparci(k,residual,jacobian)函數(shù)，置信度為95%。算出來的決定系數(shù)在0.9以上呢。算出來的置信區(qū)間好大啊，不知問題出在哪？
ci =

1.0e+03 *

-0.1666 ...

給出具體數(shù)據(jù)和方程，并指明因變量、自變量、待估計參數(shù)

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

4樓2017-04-11 11:22:39

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

1135725495

鐵桿木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 6868.7
散金: 18
紅花: 1
帖子: 1492
在線: 209.6小時
蟲號: 2756771
注冊: 2013-10-27
性別: MM
專業(yè): 有機化工

引用回帖:

4樓: Originally posted by 月只藍 at 2017-04-11 11:22:39
給出具體數(shù)據(jù)和方程，并指明因變量、自變量、待估計參數(shù)...

下面是程序，方程 r1 =-dc/dt=(k(1)*c^2-0.25*k(2)*(9.404-c)^2)/(1+k(3)*c+0.5*k(4)*(9.404-c))^2，數(shù)據(jù)是c關(guān)于t的方程，因變化量c,自變量是t,待估計參數(shù)k(1),k(2),k(3),k(4)
function parafit
%
% r1 =(k(1)*c^2-0.25*k(2)*(9.404-c)^2)/(1+k(3)*c+0.5*k(4)*(9.404-c))^2;
%
% dCAdt = - r1;
clear all
clc
%       t/min CA    / mol/L
  Kinetics=[0 8.8599786
3 7.81716904
6 7.3661532
9 7.062404
12 6.82862056
15 6.63696704
18 6.47380764
21 6.33133704
31 5.95301412
41 5.6663802
51 5.43447756
61 5.23953264];
k0 = [0.04539  0.008031  0.07300  3.0230];       % 參數(shù)初值
lb = [0 0  0  0];                % 參數(shù)下限
ub =[1  1 1  10];    % 參數(shù)上限
x0 = [8.8599786];
yexp = Kinetics;
warning off
% 使用函數(shù)fmincon()進行參數(shù)估計
[k,fval,flag] = fmincon(@ObjFunc7Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp);
fprintf('\n使用函數(shù)fmincon()估計得到的參數(shù)值為:\n')
fprintf('\tk1 = %.11f\n',k(1))
fprintf('\tk2 = %.11f\n',k(2))
fprintf('\tk3 = %.11f\n',k(3))
fprintf('\tk4 = %.11f\n',k(4))
fprintf('  The sum of the squares is: %.1e\n\n',fval)
k_fm= k;
warning off
% 使用函數(shù)lsqnonlin()進行參數(shù)估計
[k,resnorm,residual,exitflag,output,lambda, jacobian,] = ...
lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp);
ci = nlparci(k,residual,jacobian)
fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:\n')
fprintf('\tk1 = %.11f\n',k(1))
fprintf('\tk2 = %.11f\n',k(2))
fprintf('\tk3 = %.11f\n',k(3))
fprintf('\tk4 = %.11f\n',k(4))
fprintf('  The sum of the squares is: %.1e\n\n',resnorm)
k_ls = k;
output
warning off
% 以函數(shù)fmincon()估計得到的結(jié)果為初值，使用函數(shù)lsqnonlin()進行參數(shù)估計
k0 = k_fm;
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp);
ci = nlparci(k,residual,jacobian)
fprintf('\n\n以fmincon()的結(jié)果為初值，使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:\n')
fprintf('\tk1 = %.11f\n',k(1))
fprintf('\tk2 = %.11f\n',k(2))
fprintf('\tk3 = %.11f\n',k(3))
fprintf('\tk4 = %.11f\n',k(4))
fprintf('  The sum of the squares is: %.1e\n\n',resnorm)
k_fmls = k;
output
tspan = [0  3 6 9 12 15 18 21 31 41 51 61];
[t x] = ode45(@KineticEqs,tspan,x0,[],k_fmls);
figure;
plot(t,x,'o-',t,yexp(:,2),'*');legend('Predictive value','Experimental data')
x

% ------------------------------------------------------------------
function f = ObjFunc7Fmincon(k,x0,yexp)
tspan = [0  3 6 9 12 15 18 21 31 41 51 61];
[t x] = ode45(@KineticEqs,tspan,x0,[],k);
f =  sum((x-yexp(:,2)).^2) ;
% ------------------------------------------------------------------
function f = ObjFunc7LNL(k,x0,yexp)
tspan = [0  3 6 9 12 15 18 21 31 41 51 61];
[t x] = ode45(@KineticEqs,tspan,x0,[],k);
f1 = x - yexp(:,2);
f = [f1];
% ------------------------------------------------------------------
function dxdt = KineticEqs(t,x,k)
dxdt=-(k(1)*x^2-0.25*k(2)*(9.404-x)^2)/(1+k(3)*x+0.5*k(4)*(9.404-x))^2;

贊一下

回復(fù)此樓

5樓2017-04-11 12:31:25

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

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

★ ★ ★ ★ ★
1135725495: 金幣+5, ★★★很有幫助 2017-04-11 14:46:28

內(nèi)容已刪除

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

6樓2017-04-11 13:48:24

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

1135725495

鐵桿木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 6868.7
散金: 18
紅花: 1
帖子: 1492
在線: 209.6小時
蟲號: 2756771
注冊: 2013-10-27
性別: MM
專業(yè): 有機化工

引用回帖:

6樓: Originally posted by 月只藍 at 2017-04-11 13:48:24
換一下初值試試：

function parafit
%
% r1 =(k(1)*c^2-0.25*k(2)*(9.404-c)^2)/(1+k(3)*c+0.5*k(4)*(9.404-c))^2;
%
% dCAdt = - r1;
clear all
clc
%       t/min CA    / mol/L
  Kinetics= ...

你好，謝謝您的幫助，請問這個初值的選擇有什么規(guī)律嗎？我發(fā)現(xiàn)初值差一點，都會造成置信區(qū)間的范圍很大

贊一下

回復(fù)此樓

7樓2017-04-11 15:17:47

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

引用回帖:

7樓: Originally posted by 1135725495 at 2017-04-11 15:17:47
你好，謝謝您的幫助，請問這個初值的選擇有什么規(guī)律嗎？我發(fā)現(xiàn)初值差一點，都會造成置信區(qū)間的范圍很大...

多試

回復(fù)此樓

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

1135725495

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

8樓2017-04-11 15:19:29

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

1135725495

鐵桿木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 6868.7
散金: 18
紅花: 1
帖子: 1492
在線: 209.6小時
蟲號: 2756771
注冊: 2013-10-27
性別: MM
專業(yè): 有機化工

送紅花一朵

引用回帖:

8樓: Originally posted by 月只藍 at 2017-04-11 15:19:29
多試...

你好！按照您說的，我試了一下午到現(xiàn)在還是不行，置信區(qū)間的數(shù)量級沒有變化。還有好幾個方程需要計算，也是這種問題，就是換一下c的值。你能不能幫我試一下，或者你有沒有什么方法可以找到比較好的初值。謝謝！

贊一下

回復(fù)此樓

9樓2017-04-11 20:49:19

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

引用回帖:

9樓: Originally posted by 1135725495 at 2017-04-11 20:49:19
你好！按照您說的，我試了一下午到現(xiàn)在還是不行，置信區(qū)間的數(shù)量級沒有變化。還有好幾個方程需要計算，也是這種問題，就是換一下c的值。你能不能幫我試一下，或者你有沒有什么方法可以找到比較好的初值。謝謝！...

待估計參數(shù)不超過4個，都可以用5.5版本的1stopt軟件計算的結(jié)果當初值。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

10樓2017-04-11 21:05:54

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

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

返回列表

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

[考研] 一志愿西安交通大學學碩 354求調(diào)劑 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 317求調(diào)劑 +6	申子申申 2026-03-19	11/550	2026-03-20 20:07 by zhukairuo
[考研] 085700資源與環(huán)境308求調(diào)劑 +11	墨墨漠 2026-03-18	12/600	2026-03-20 19:43 by 丁丁*
[考研] 290求調(diào)劑 +4	^O^乜 2026-03-19	4/200	2026-03-20 11:51 by 學員8dgXkO
[考研] 材料與化工求調(diào)劑 +7	為學666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 一志愿福大288有機化學，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 一志愿985，本科211，0817化學工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 332求調(diào)劑 +3	ydfyh 2026-03-17	3/150	2026-03-19 10:14 by 功夫瘋狂
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 0703化學調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 考研化學學碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 302求調(diào)劑 +4	小賈同學123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學123
[考研] 藥學383 求調(diào)劑 +3	藥學chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike