 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (989)
- >考研 (185)
- >導(dǎo)師招生 (132)
- >基金申請(qǐng) (19)
- >考博 (14)
- >蟲(chóng)友互識(shí) (10)
- >碩博家園 (8)
- >休閑灌水 (7)
- >博后之家 (5)
- >招聘信息布告欄 (4)
- >論文道賀祈福 (4)
- >復(fù)合材料 (2)
- >學(xué)術(shù)會(huì)議 (2)
- >教師之家 (2)
- >找工作 (2)
- >專業(yè)外語(yǔ) (1)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 366 | 回復(fù): 0 | ||
1135725495鐵桿木蟲(chóng) (著名寫(xiě)手)
|
[求助]
運(yùn)行matlab出現(xiàn)的問(wèn)題
|
|
Caused by: Failure in initial user-supplied objective function evaluation. FMINCON cannot continue. 誰(shuí)知道這是怎么回事啊？ function parafit % % r1 =(k(1)*c^2-0.5*k(2)*(9.404-c))/(1+k(3)*c+0.5*k(4)*(9.404-c)+2*k(5)*c^2/(k(4)*(9.404-c))); % % dCAdt = - r1; clear all clc % t/min CA / mol/L Kinetics=[0 7.9556 3 5.0759 6 4.3477 9 3.9301 12 3.6489 15 3.4439 18 3.2875 21 3.1645 31 2.8989 41 2.7540 51 2.6694 61 2.6182]; k0 = [0.1 0.1 0.01 0.01 0.01]; % 參數(shù)初值 lb = [0 0 0 0 0]; % 參數(shù)下限 ub = [100 100 100 100 100]; % 參數(shù)上限 x0 = [7.9556]; 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('\tk5 = %.11f\n',k(5)) 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('\tk5 = %.11f\n',k(5)) 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('\tk5 = %.11f\n',k(5)) 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('ca-pr','ca-real') 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.5*k(2)*(9.404-x))/(1+k(3)*c+0.5*k(4)*(9.404-x)+2*k(5)*x^2/(k(4)*(9.404-x))); |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有133人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有65人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
matlab實(shí)現(xiàn)牛頓迭代法求非線性方程組程序運(yùn)行出現(xiàn)問(wèn)題
已經(jīng)有8人回復(fù)
- COMSOL4.4不能成功啟用MATLAB，怎么辦？
已經(jīng)有5人回復(fù)
- matlab中斷循環(huán)
已經(jīng)有9人回復(fù)
- matlab提前設(shè)定語(yǔ)句運(yùn)行時(shí)間的問(wèn)題
已經(jīng)有4人回復(fù)
- 使用music_code運(yùn)行出現(xiàn)錯(cuò)誤
已經(jīng)有4人回復(fù)
- matlab 運(yùn)行ode45出錯(cuò)，不知道什么原因
已經(jīng)有4人回復(fù)
- linux 下 matlab 的調(diào)用問(wèn)題
已經(jīng)有3人回復(fù)
- matlab “Inner matrix dimensions must agree.”一直不知道哪里出錯(cuò)了 求助
已經(jīng)有6人回復(fù)
- matlab一運(yùn)行就出現(xiàn)的問(wèn)題
已經(jīng)有7人回復(fù)
- 請(qǐng)大俠指教matlab的這個(gè)錯(cuò)誤提示到底是個(gè)神馬意思？
已經(jīng)有12人回復(fù)
- 有個(gè)matlab程序，但是運(yùn)行不下來(lái)，求高手
已經(jīng)有23人回復(fù)
- 運(yùn)行matlab出現(xiàn)錯(cuò)誤
已經(jīng)有4人回復(fù)
- 關(guān)于matlab GUI打包生產(chǎn)exe文件，并使其在沒(méi)有安裝matlab的計(jì)算機(jī)上運(yùn)行的問(wèn)題。。
已經(jīng)有4人回復(fù)
- matlab的m文件關(guān)聯(lián)問(wèn)題
已經(jīng)有17人回復(fù)
- matlab解決重心法選址問(wèn)題，運(yùn)行不了啊，求助各位大神
已經(jīng)有10人回復(fù)
- matlab 運(yùn)行出錯(cuò)，求幫忙啦~~·
已經(jīng)有3人回復(fù)
- 問(wèn)下這個(gè)問(wèn)題怎么解決.....真心啊
已經(jīng)有4人回復(fù)
- matlab生成隨機(jī)行走網(wǎng)格 模擬材料生長(zhǎng)
已經(jīng)有26人回復(fù)
- 求助matlab-simulink運(yùn)行時(shí)報(bào)錯(cuò)問(wèn)題
已經(jīng)有3人回復(fù)
- 同樣的MATLAB程序，為何在別人的電腦上運(yùn)行不了
已經(jīng)有21人回復(fù)
- MATLAB安裝后出現(xiàn)一堆java的問(wèn)題
已經(jīng)有12人回復(fù)
- 【求助】誰(shuí)能幫我解決一下matlab安裝問(wèn)題【已解決】
已經(jīng)有8人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 1135725495 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
