 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 688 | 回復(fù): 0 | ||
vas111新蟲 (初入文壇)
|
[求助]
為何decic函數(shù)提示Not enough input arguments.？
|
|
函數(shù)如下： function TwoDimensions_OCM_Main % 用對(duì)稱正交配置法求解固定床反應(yīng)器二維擬均相穩(wěn)態(tài)模型（二維穩(wěn)態(tài)PDE方程組） % (只在r方向無因次化) clear all clc format long global AA BB RR N Tw r DeltaH m n R v P T0 F0 R = 8.314; %摩爾氣體常數(shù)(單位：J/(mol*K)) %配置常數(shù)的選取 N = xlsread('Calculating_parameter.xlsx', 1, 'B6'); % 內(nèi)配置點(diǎn)個(gè)數(shù) if N == 3 AA = xlsread('Configuration_constants.xlsx', 2, 'B2:E5'); BB = xlsread('Configuration_constants.xlsx', 2, 'B11:E14'); RR = xlsread('Configuration_constants.xlsx', 2, 'B19:E19'); elseif N == 5 AA = xlsread('Configuration_constants.xlsx', 2, 'H2:M7'); BB = xlsread('Configuration_constants.xlsx', 2, 'H11:M16'); RR = xlsread('Configuration_constants.xlsx', 2, 'H19:M19'); elseif N == 7 AA = xlsread('Configuration_constants.xlsx', 2, 'P2:W9'); BB = xlsread('Configuration_constants.xlsx', 2, 'P11:W18'); RR = xlsread('Configuration_constants.xlsx', 2, 'P19:W19'); else fprintf('區(qū)域選取錯(cuò)誤！N應(yīng)為3、5或7的整數(shù)！'); return end %=========================================================================== % 輸入系數(shù) Er = xlsread('Calculating_parameter.xlsx', 1, 'B2'); % 徑向有效擴(kuò)散系數(shù)（單位： m^2/h） u = xlsread('Calculating_parameter.xlsx', 1, 'B3'); % 空塔氣速（單位： m/h） rhoB = xlsread('Calculating_parameter.xlsx', 1, 'B4'); % 催化劑堆積密度(單位：kg/m^3) lambda_er = xlsread('Calculating_parameter.xlsx', 1, 'B5');% 導(dǎo)熱系數(shù)(單位：W/(m*K)) n = xlsread('Calculating_parameter.xlsx', 1, 'B7');% 反應(yīng)組分?jǐn)?shù) hw = xlsread('Calculating_parameter.xlsx', 1, 'B8');% 壁面對(duì)流傳熱系數(shù)(單位：W/(m^2·K)) Tw = xlsread('Calculating_parameter.xlsx', 1, 'B9');% 壁面溫度（單位： K） T0 = xlsread('Calculating_parameter.xlsx', 1, 'B10');% 進(jìn)口溫度（單位： K） P = xlsread('Calculating_parameter.xlsx', 1, 'B11');% 進(jìn)口壓力（單位： Pa） r = xlsread('Calculating_parameter.xlsx', 1, 'B12');% 內(nèi)徑(單位：m) L = xlsread('Calculating_parameter.xlsx', 1, 'B13');% 管長(zhǎng)(單位：m) %=========================================================================== Cp = xlsread('Calculating_parameter.xlsx', 4, 'B2:Z2'); %各組分的比熱容(單位：kJ/(kmol*K)) DeltaH = xlsread('Calculating_parameter.xlsx', 4, 'B9:Z9'); %各個(gè)反應(yīng)的反應(yīng)熱(單位：kJ/mol) F0 = xlsread('Calculating_parameter.xlsx', 4, 'B3:Z3'); %各組分的進(jìn)口摩爾流速(單位：kmol/(m^2*h)) v = xlsread('Calculating_parameter.xlsx', 3, 'B2:Z21'); %反應(yīng)計(jì)量系數(shù)矩陣 m = rank(v); %獨(dú)立反應(yīng)數(shù) %=========================================================================== Xi(:,:,m) = 0; %對(duì)反應(yīng)進(jìn)度Xi初始化 [A,B,E] = AtoE(Er,u,rhoB,lambda_er,hw); %利用函數(shù)AtoE對(duì)假定的參數(shù)A,B,C,D,E賦值 y0(1:m*N) = 0; %狀態(tài)變量Xi的初值，設(shè)為0 y0(m*N+1: (m+1)*N) = T0; %狀態(tài)變量T的初值，設(shè)為0 y0 = transpose(y0); %對(duì)y0轉(zhuǎn)置 t0 = 0; fix_y0=ones((m+1)*N,1); dy0 = zeros((m+1)*N,1); fix_dy0=zeros((m+1)*N,1); [y0,dy0]=decic(@Equations,t0,y0,fix_y0,dy0,fix_dy0); [z,y] = ode15i(@Equations,[0 0.0023],y0,dy0,[],A,B,lambda_er,rhoB,E,Cp); %ode15i求解方程 %============================================================================ function dydx = Equations(t,y,dy,A,B,lambda_er,rhoB,E,Cp) global AA BB RR N Tw r DeltaH m n R v P F0 T= y(m*N+1: (m+1)*N); dTdl = dy(m*N+1: (m+1)*N); for j=1:m X(:,j) = y((j-1)*N+1:j*N);%j表示第j個(gè)獨(dú)立反應(yīng),其中X的前j列表示Xi的 dXdl(:,j) =dy((j-1)*N+1:j*N); end Rgene = ReactionRate(T(1:N),X)%普遍反應(yīng)速率,結(jié)果返回應(yīng)為矩陣，行代表對(duì)應(yīng)的反應(yīng)數(shù)m，列代表配置點(diǎn)個(gè)數(shù)N for i=1:n for j=1:N F(i,j) = F0(i)+sum(X(j,: )/(3.14*r^2).*v(:,i)'); dF(i,j) = sum(v(1:m,i).*dXdl(j,1:m)'/(3.14*r^2)); end end for i=1:m for j=1:N CA(j,i) = P/R/T(j)*F(i,j)/sum(F(:,j)); end end for i=1:m CAb(i)=-sum(AA(N+1,1:N).*CA(:,i)'/AA(N+1,N+1)); end for i=1:m % m表示獨(dú)立反應(yīng)數(shù)的個(gè)數(shù) n表示反應(yīng)組分?jǐn)?shù) for j=1:N % N表示正交配置點(diǎn)個(gè)數(shù) dC(j,i)=-P/R/(T(j)^2)*dTdl(j)*F(i,j)/sum(F(:,j))+(P/R/T(j)*(dF(i,j)*sum(F(:,j))-F(i,j)*(sum(dF(:,j))))/(sum(F(:,j)))^2)... -(A/r^2*(sum((BB(j,1:N)+AA(j,1:N)./RR(j)).*CA(:,i)')+(BB(j,N+1)+AA(j,N+1)./RR(j)).*CAb(i))-B*(-sum(v(1:m,i).*Rgene(1:m,j)))); end end Tb = (r*E*(T(N)-Tw)- sum(AA(N+1,1:N).*T(1:N)')/AA(N+1,N+1)); for j=1:N FiCpi(j) = sum(F(1:n,j)'.*Cp(1:n)); end for j=1:N dT(j) = dTdl(j)-((lambda_er/FiCpi(j))/r^2*(sum((BB(j,1:N)+AA(j,1:N)./RR(j)).*T(1:N)')+(BB(j,N+1)+AA(j,N+1)./RR(j))*Tb)+(rhoB/FiCpi(j))*sum(-DeltaH(1:m)'.*(-Rgene(1:m,j)))); end dC = reshape(dC,N*m,1);% 將dC變化為列向量 dydx = [dC;dT']; 然后錯(cuò)誤如下： Error using Equations (line 28) Not enough input arguments. Error in decic (line 67) res = feval(odefun,t0,y0,yp0,varargin{:}); Error in TwoDimensions_OCM_Main (line 62) [y0,dy0]=decic(@Equations,t0,y0,fix_y0,dy0,fix_dy0); |
回復(fù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(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分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(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ù)
找到一些相關(guān)的精華帖子，希望有用哦~
Undefined function 'mashgrid' for input arguments of type 'double'.
已經(jīng)有5人回復(fù)
- matlab Undefined function or method '～' for input arguments of type 'char'.
已經(jīng)有12人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲開始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 vas111 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +3 | www_q 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +5 | Doleres 2026-03-19 | 5/250 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|