版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時(shí)熱門版塊排行榜

返回列表

guanwei

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1075.1
散金: 195
紅花: 4
帖子: 247
在線: 53.7小時(shí)
蟲號(hào): 1638561
注冊: 2012-02-23
性別: GG
專業(yè): 化工系統(tǒng)工程

[求助] MATLAB擬合,正丁醇與乳酸酯化反應(yīng)PH模型的擬合程序出現(xiàn)問題已有1人參與

各位位，您好!
本人研究的方向是反應(yīng)精餾，最近在評(píng)估PH模型的擬合找到反應(yīng)速率常數(shù)，用最小二乘法進(jìn)行計(jì)算，看到文獻(xiàn)的編輯的程序，多次嘗試，有如下疑問，希望各位蟲友能幫忙解決

如下363.15K,正丁醇與乳酸酯化反應(yīng)PH模型的擬合程序，

function KineticsEst1_int11
clear all
clc
t=[0 20 40 100 270 450 750 1410 1590 1810];%時(shí)間
rA=[0.0004765 0.0004507  0.0004250  0.0003477 0.001549 0.0001117 0.0000396 0.0000018  0.0000038 0.0000062];%反應(yīng)速率
Pa=[0.042729739 0.041609737 0.040801805 0.038526838 0.033412812 0.030519156 0.028695168 0.02792077 0.027918742 0.0278887981]; %乳酸的摩爾分率
Pb=[0.123670897  0.122550895  0.121742963 0.119467996 0.11435397  0.111460314  0.109636326  0.108861928 0.1088599 0.108829956];%正T醇的摩爾分率
Pc=[0.833599364  0.834719366 0.835527299 0.837802265 0.842916291 0.845809947 0.847633933 0.848408333 0.848410361 0.848440305];%水的摩爾分率
Pd=[0.001120002 0.001927934 0.004202901 0.009316927 0.012210583 0.014034571 0.014808969 0.014810996  0.014840941  0.014850941];%乳酸正丁酯的摩爾分率
%線性擬合
P=2.337132745*Pa.*Pb-0.564655417*Pc.*Pd;y=rA';X=[ones(size(y))  P'];

b=X\y;k=b(2);
%非線性擬合
beta0=[k]
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunc,beta0,[],[],[],rA,Pa,Pb,Pc,Pd)

%擬合效果圖(實(shí)驗(yàn)與擬合的比較)
figure(1);plot(t,rA,'.')
r_poly=beta(1)*(2.337132745*Pa.*Pb-0.564655417*Pc.*Pd);
hold on;plot(t, r_poly,'g')
figure(2);plot(Pb,rA,'.'); hold on;plot(Pb, r_poly,'g')

%-----------------------------------------------------------------------------一
Function f=ObjFunc(beta,rA,Pa,Pb,Pc,Pd)
f=rA-beta(1)*(2.337132745*Pa.*Pb-0.564655417*Pc.*Pd);

出現(xiàn)以下錯(cuò)誤
??? Error using ==> feval
Undefined function or method 'ObjFunc' for input arguments of type 'double'.

Error in ==> lsqnonlin at 203
         initVals.F = feval(funfcn{3},xCurrent,varargin{:});

Error in ==> KineticsEst1_int11 at 33
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...

Caused by:
Failure in initial user-supplied objective function evaluation. LSQNONLIN
cannot continue.

回復(fù)此樓

» 猜你喜歡

請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有193人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有78人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)

1樓 2014-09-05 17:17:22

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

月只藍(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í)
蟲號(hào): 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

感謝參與，應(yīng)助指數(shù) +1

包括MATLAB在內(nèi)的4種常見軟件的程序代碼或者操作步驟參見：
http://www.gaoyang168.com/bbs/viewthread.php?tid=6572963&authorid=1122189

如果是代數(shù)方程，MATLAB程序編寫不方便，建議嘗試用1stopt。
如果是常微分方程，可在本版塊搜索相關(guān)MATLAB程序，或者求助專家dingd。

贊一下

回復(fù)此樓

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

2樓2014-09-05 18:18:29

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

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