 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1320 | 回復(fù): 12 | |||
1135725495鐵桿木蟲 (著名寫手)
|
[求助]
matlab 已有1人參與
|
||
用matlab軟件計(jì)算置信區(qū)間，算出來的結(jié)果過大怎么解決？？？ @月只藍(lán) |
回復(fù)此樓» 猜你喜歡
2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有294人回復(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ù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有30人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有63人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
6樓2017-04-11 13:48:24
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
2樓2017-04-11 10:02:44
1135725495
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時(shí)
- 蟲號: 2756771
- 注冊: 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
|
你好！以下是用matlab計(jì)算的結(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()估計(jì)得到的參數(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] |
3樓2017-04-11 11:16:48
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
4樓2017-04-11 11:22:39
1135725495
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時(shí)
- 蟲號: 2756771
- 注冊: 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
|
下面是程序，方程 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,待估計(jì)參數(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()進(jìn)行參數(shù)估計(jì) [k,fval,flag] = fmincon(@ObjFunc7Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp); fprintf('\n使用函數(shù)fmincon()估計(jì)得到的參數(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()進(jìn)行參數(shù)估計(jì) [k,resnorm,residual,exitflag,output,lambda, jacobian,] = ... lsqnonlin(@ObjFunc7LNL,k0,lb,ub,[],x0,yexp); ci = nlparci(k,residual,jacobian) fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(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()估計(jì)得到的結(jié)果為初值，使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計(jì) 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()估計(jì)得到的參數(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; |
5樓2017-04-11 12:31:25
1135725495
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時(shí)
- 蟲號: 2756771
- 注冊: 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
7樓2017-04-11 15:17:47
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
11357254958樓2017-04-11 15:19:29
1135725495
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 6868.7
- 散金: 18
- 紅花: 1
- 帖子: 1492
- 在線: 209.6小時(shí)
- 蟲號: 2756771
- 注冊: 2013-10-27
- 性別: MM
- 專業(yè): 有機(jī)化工
9樓2017-04-11 20:49:19
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
10樓2017-04-11 21:05:54
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 1135725495 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +3 | lemonzzn 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|