 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 285 | 回復(fù): 0 | ||
鄭美琴琴金蟲(chóng) (著名寫(xiě)手)
|
[求助]
matlab程序出錯(cuò)
|
大神幫幫忙，感激不盡�。�！ M文件 function yt=ped_a(t,y) global nr nz dr dz drs dzs... r z Ds DL c q... cs v e R kf K... c0 p n ncall %1D to 2D矩陣 for i=1:nz for j=1:nr ij=(i-1)*nr+j; c(i)=y(i); q(i,j)=y(ij+nz); end end %r,z網(wǎng)格化 for i=1:nz for j=1:nr %2/r*qr if(j==1) %r=0 qr(i,j)=2.0*(q(i,j+1)-q(i,j))/drs; elseif(j==nr) %r=R qr(i,j)=2/R*kf/(p*Ds)*(c(i)-cs); else qr(i,j)=2/r(j)*(q(i,j+1)-q(i,j-1))/(2*dr); end %qrr if(j==1) qrr(i,j)=2.0*(q(i,j+1)-q(i,j))/drs; elseif(j==nr) qf(i,j)=q(i,j-1)+2.0*dr*kf/(p*Ds)*(c(i)-cs); qrr(i,j)=(qf(i,j)-2.0*q(i,j)+q(i,j-1))/drs; else qrr(i,j)=(q(i,j+1)-2.0*q(i,j)+q(i,j-1))/drs; end %cz if(i==1) %z=0 cz(i)=(c(i)-c0)/dz; elseif(i==nz) %z=L cz(i)=0; else cz(i)=(c(i)-c(i-1))/dz; end %czz czz(i)=(c(i+1)-2.0*c(i)+c(i-1))/dzs; %PDEs qt(i,j)=Ds*(qrr(i,j)+qr(i,j)); ct(i)=DL*czz(i)-v*cz(i)-3*(1-e)/(e*R)*kf*(c(i)-cs(i,j)); q(i,j)=K*cs(i,j)^(1/n); end end %2D to 1D for i=1:nz for j=1:nr ij=(i-1)*nr+j; yt(i)=ct(i); yt(ij+nz)=qt(i,j); end end %轉(zhuǎn)置和計(jì)數(shù) yt=yt'; ncall=ncall+1; Command window： %全部區(qū)域 >> global nr nz dr dz drs dzs... r z Ds DL c q... cs v e R kf K... p n ncall >> %模型參數(shù) >> c0=250; >> L=0.15; >> R=0.001; >> p=436.8; >> v=5.8*10^(-5); >> e=0.363; >> DL=1.66*10^(-7); >> Ds=2.4281*10^(-12); >> K=0.14; >> n=1; >> %z方向網(wǎng)格 >> nz=20; >> dz=L/nz; >> for i=1:nz z(i)=i*dz; end >> dzs=dz^2; >> %r方向網(wǎng)格 >> nr=7; >> dr=R/(nr-1); >> for j=1:nr r(j)=(j-1)*dr; end >> drs=dr^2; %ODE集成 >> tf=200; >> tout=[0.0:50:tf]; >> nout=5; >> ncall=0; >> %初始條件 >> for i=1:nz for j=1:nr if(i==1) c(i)=c0 else c(i)=0 end q(i,j)=0 y0(i)=c(i); y0((i-1)*nr+j+nz)=q(i,j) end end %ODE集成 >> reltol=1.0e-04; abstol=1.0e-04; options=odeset('RelTol',reltol,'AbsTol',abstol); [t,y]=ode15s(@ped_a,tout,y0,options); 模擬結(jié)果： ??? In an assignment A(I) = B, the number of elements in B and I must be the same. Error in ==> ped_a at 36 cz(i)=(c(i)-c0)/dz; Error in ==> funfun\private\odearguments at 110 f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ==> ode15s at 227 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ... |
回復(fù)此樓» 猜你喜歡
070300一志愿211，312分求調(diào)劑院校
已經(jīng)有4人回復(fù)
-
332求調(diào)劑
已經(jīng)有10人回復(fù)
-
279求調(diào)劑
已經(jīng)有6人回復(fù)
-
310求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料專(zhuān)碩 085600求調(diào)劑
已經(jīng)有3人回復(fù)
-
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可
已經(jīng)有12人回復(fù)
-
生物技術(shù)與工程
已經(jīng)有8人回復(fù)
-
085602 化學(xué)工程專(zhuān)碩 340分求調(diào)劑
已經(jīng)有4人回復(fù)
-
292求調(diào)劑
已經(jīng)有8人回復(fù)
-
297求調(diào)劑
已經(jīng)有13人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
請(qǐng)問(wèn)這段小matlab代碼怎么修改？
已經(jīng)有10人回復(fù)
- matlab程序錯(cuò)誤，提示如下【弄了半夜還是不行】
已經(jīng)有10人回復(fù)
- matlab 運(yùn)行ode45出錯(cuò)，不知道什么原因
已經(jīng)有4人回復(fù)
- matlab求函數(shù)積分出錯(cuò)
已經(jīng)有12人回復(fù)
- matlab一運(yùn)行就出現(xiàn)的問(wèn)題
已經(jīng)有7人回復(fù)
- 有個(gè)matlab程序，但是運(yùn)行不下來(lái)，求高手
已經(jīng)有23人回復(fù)
- matlab程序中調(diào)用m文件的問(wèn)題
已經(jīng)有13人回復(fù)
- 請(qǐng)教MATLAB中的LMI求解問(wèn)題，急�。�！
已經(jīng)有7人回復(fù)
- matlab-出錯(cuò)提示Undefined method for input arguments type 'double'.
已經(jīng)有5人回復(fù)
- 求助多輸出支持向量機(jī)的matlab程序
已經(jīng)有6人回復(fù)
- 大家愛(ài)幫我看看這段MATLAB程序?yàn)樯冻鲥e(cuò)了
已經(jīng)有5人回復(fù)
- 模式搜索MATLAB程序
已經(jīng)有9人回復(fù)
- matlab程序問(wèn)題，總是出錯(cuò)，麻煩大家看看問(wèn)題在那里？
已經(jīng)有6人回復(fù)
- 求MATLAB程序包SDELab
已經(jīng)有8人回復(fù)
- matlab賦值問(wèn)題
已經(jīng)有5人回復(fù)
- 同樣的MATLAB程序，為何在別人的電腦上運(yùn)行不了
已經(jīng)有21人回復(fù)
- 看不懂matlab程序該怎么辦？
已經(jīng)有7人回復(fù)
- 求助matlab編程
已經(jīng)有13人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 鄭美琴琴 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 282求調(diào)劑 +4 | wcq131415 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
[考研] 289求調(diào)劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
|
[考研] 343求調(diào)劑 +6 | 愛(ài)羈絆 2026-03-29 | 6/300 |
|
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲(chóng)085600 2026-03-27 | 8/400 |
|
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專(zhuān)業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 材料專(zhuān)碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚(yú)567 2026-03-23 | 4/200 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|