 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 819 | 回復(fù): 4 | |||
你的CEO鐵蟲(chóng) (小有名氣)
|
[求助]
動(dòng)力學(xué)參數(shù)擬合matlab程序有問(wèn)題 已有1人參與
|
| 求大神幫看下程序問(wèn)題在哪？ |
回復(fù)此樓» 猜你喜歡
265求調(diào)劑
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有12人回復(fù)
-
286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù)
已經(jīng)有4人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
081700化工學(xué)碩調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一）
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 向大神求教一個(gè)用matlab擬合動(dòng)力學(xué)參數(shù)的問(wèn)題，求指導(dǎo)
已經(jīng)有5人回復(fù)
- 給我大神求教一個(gè)用matlab擬合動(dòng)力學(xué)參數(shù)的問(wèn)題，求指導(dǎo)
已經(jīng)有0人回復(fù)
- Matlab擬合動(dòng)力學(xué)參數(shù)問(wèn)題，求指導(dǎo)！
已經(jīng)有9人回復(fù)
- 各位大神能不能幫小弟看看動(dòng)力學(xué)編程的問(wèn)題
已經(jīng)有8人回復(fù)
- matlab 非線(xiàn)性擬合 數(shù)據(jù)擬合
已經(jīng)有3人回復(fù)
- 最小二乘 平方和初值小于函數(shù)公差的默認(rèn)值（動(dòng)力學(xué)方程擬合）
已經(jīng)有3人回復(fù)
- 使用matlab最優(yōu)化方法擬合獲得多個(gè)動(dòng)力學(xué)參數(shù)中的問(wèn)題
已經(jīng)有4人回復(fù)
- 求助1stopt擬合動(dòng)力學(xué)參數(shù)
已經(jīng)有4人回復(fù)
- Matlab動(dòng)力學(xué)擬合問(wèn)題！求助�。。。�
已經(jīng)有0人回復(fù)
- 使用lsqnonlin函數(shù)優(yōu)化動(dòng)力學(xué)參數(shù)，總是得不到合理的結(jié)果
已經(jīng)有6人回復(fù)
- matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- matlab擬合擬均相動(dòng)力學(xué)參數(shù)
已經(jīng)有6人回復(fù)
- 請(qǐng)那位幫我看看這個(gè)MATLAB程序那里有問(wèn)題！先謝謝
已經(jīng)有2人回復(fù)
- 請(qǐng)那位幫我看看這個(gè)MATLAB程序那里有問(wèn)題！先謝謝
已經(jīng)有0人回復(fù)
- 【求助】使用Matlab預(yù)估動(dòng)力學(xué)方程問(wèn)題
已經(jīng)有13人回復(fù)
你的CEO
鐵蟲(chóng) (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 175.6
- 帖子: 71
- 在線(xiàn): 18小時(shí)
- 蟲(chóng)號(hào): 3071920
- 注冊(cè): 2014-03-19
- 性別: GG
- 專(zhuān)業(yè): 能源化工
|
程序如下 ： function PX4800602 clear all clc global x0 % k0 = [1.8887 0.0036 0.0070 0.7557 0.1455]; % k0 = [4.4163008451 0.0033686022 0.0065669052 1.6201160059 16.2647816]; k0= [ 5.2766233638 0.0030934069 0.0060459974 1.8656945060 20.6241874]; % 參數(shù)初值 lb = [0 0 0 0 0]; % 參數(shù)下限 ub = [+inf +inf +inf +inf +inf ]; % 參數(shù)上限 x0 = [0 0 0 0 0 0]; tspn=[0 43.2 72]; % yexp =[0 0 0 0 0;0.2282 0.2067 0.0057 0.0025 0.0132;0.2394 0.2170 0.0063 0.0037 0.0123; 0.2515 0.2278 0.0069 0.0040 0.0128;0.2638 0.2378 0.0081 0.0043 0.0135;0.2738109 0.2461347 0.0049171 0.011883 0.01087607]; % yexp: 實(shí)驗(yàn)數(shù)據(jù) % yexp =[0 0 0 0 0;22.81824 20.67075 0.571198 0.253676 1.322614; 23.94047 21.7016 0.633422 0.372055 1.233396;25.14798 22.77878 0.692198 0.39614 1.280859;26.37669 23.78016 0.811828 0.432807 1.351893;27.38109 24.61347 0.8883 0.49171 1.387607]*0.01; % yexp: 實(shí)驗(yàn)數(shù)據(jù) yexp=[0 0 0 0 0 0; 22.253 20.067 0.340 0.633 1.213 53.06753634; 23.532 21.175 0.385 0.662 1.310 50.31636493]; % yexp=[0 0 0 0 0 0;23.679 21.655 0.421 0.715 0.888 50.9; 24.415 22.131 0.463 0.724 1.096 48.978;25.148 22.580 0.481 0.751 1.336 47.073;26.377 23.580 0.564 0.869 1.364 44.518;27.110 23.948 0.594 1.177 1.392 43; % 15.51 14.1 0.2476 0.545 0.614 67.75;15.78 14.24 0.265 0.6352 0.632 67.18;16.4 14.51 0.292 0.672 0.937 65.32;]*0.01; % yexp: 兩個(gè)配比 % yexp=[0 0 0 0 0 0;23.679 21.655 0.421 0.715 0.888 50.9; 24.415 22.131 0.463 0.724 1.096 48.978;25.148 22.580 0.481 0.751 1.336 47.073;26.377 23.580 0.564 0.869 1.364 44.518;27.110 23.948 0.594 1.177 1.392 43;] % yexp: 一個(gè)配比 % 使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計(jì) [k,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@ObjFunc4LNL,k0,lb,ub,optimset('TolFun',1.0000e-12),x0,yexp); ci = nlparci(k,residual,jacobian) fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n') fprintf('\tk1 = %.10f ± %.7f\n',k(1),ci(1,2)-k(1)) fprintf('\tk2 = %.10f ± %.7f\n',k(2),ci(2,2)-k(2)) fprintf('\tk3 = %.10f ± %.7f\n',k(3),ci(3,2)-k(3)) fprintf('\tk4 = %.10f ± %.7f\n',k(4),ci(4,2)-k(4)) fprintf('\tk5 = %.7f ± %.7f\n',k(5),ci(5,2)-k(5)) % fprintf('\tk6 = %.7f ± %.7f\n',k(4),ci(6,2)-k(6)) fprintf(' The sum of the squares is: %.1e\n\n',resnorm) % ------------------------------------------------------------------ function f = ObjFunc4LNL(k,x0,yexp) tspn=[0 43.2 72]; [t x] = ode45(@KineticEqs1,tspn,x0,[],k); y = x f1 = y(2:3,1) - yexp(2:3,1); f2 = y(2:3,2) - yexp(2:3,2); f3 = y(2:3,3) - yexp(2:3,3); f4 = y(2:3,4) - yexp(2:3,4); f5 = y(2:3,5) - yexp(2:3,5); f6 = y(2:3,6) - yexp(2:3,6); z1=x % -------------------------------------------不同配比的方程組 function dxdt = KineticEqs1(t,x,k) p0=1.3*1; dxdt = ... [(k(1)*(2*p0*(1-x(1))/(19+0.5*x(6)))*(1-2*x(1)-2*x(5)-x(6))*p0/(19+0.5*x(6))) (k(1)*(2*p0*(1-x(1))/(19+0.5*x(6)))*((1-2*x(1)-2*x(5)-x(6))*p0/(19+0.5*x(6)))-2*(k(2)+k(3))*p0*x(2)/(19+0.5*x(6))-2*k(4)*p0*x(2)/(19+0.5*x(6))*((1-2*x(1)-2*x(5)-x(6))*p0/(19+0.5*x(6)))) (2*k(2)*p0*x(2)/(19+0.5*x(6))-2*k(4)*p0*x(3)/(19+0.5*x(6))*((1-2*x(1)-2*x(5)-x(6))*p0/(19+0.5*x(6)))) (2*k(3)*p0*x(2)/(19+0.5*x(6))-2*k(4)*p0*x(4)/(19+0.5*x(6))*((1-2*x(1)-2*x(5)-x(6))*p0/(19+0.5*x(6)))) (k(4)*(2*(x(2)+x(3)+x(4))*p0/(19+0.5*x(6)))*((1-2*x(1)-2*x(5)-x(6))*p0/(19+0.5*x(6)))) (2*k(5)*((1-2*x(1)-2*x(5)-x(6))*p0/(19+0.5*x(6)))^2) ]; |
2樓2015-10-10 13:49:42
你的CEO
鐵蟲(chóng) (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 175.6
- 帖子: 71
- 在線(xiàn): 18小時(shí)
- 蟲(chóng)號(hào): 3071920
- 注冊(cè): 2014-03-19
- 性別: GG
- 專(zhuān)業(yè): 能源化工
|
??? Error using ==> feval Output argument "f" (and maybe others) not assigned during call to "C:\Users\cai\Desktop\Matlab模擬\PX程序\PX新動(dòng)力學(xué)\擬合程序\PX480602.m>ObjFunc4LNL". **********錯(cuò)誤提示***************** Error in ==> lsqnonlin at 203 initVals.F = feval(funfcn{3},xCurrent,varargin{:}); Error in ==> PX480602 at 30 [k,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@ObjFunc4LNL,k0,lb,ub,optimset('TolFun',1.0000e-12),x0,yexp); Caused by: Failure in initial user-supplied objective function evaluation. LSQNONLIN cannot continue. |
3樓2015-10-10 13:50:13
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專(zhuān)家經(jīng)驗(yàn): +1059 - 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線(xiàn): 3291.4小時(shí)
- 蟲(chóng)號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專(zhuān)業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
4樓2015-10-10 14:35:18
你的CEO
鐵蟲(chóng) (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 175.6
- 帖子: 71
- 在線(xiàn): 18小時(shí)
- 蟲(chóng)號(hào): 3071920
- 注冊(cè): 2014-03-19
- 性別: GG
- 專(zhuān)業(yè): 能源化工
5樓2015-10-18 16:59:01
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 你的CEO 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9 | Ymlll 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考研] 294求調(diào)劑材料與化工專(zhuān)碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)（081702）315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 材料專(zhuān)業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 材料專(zhuān)碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|