Matlab 解 PDE 的問題

請問Matlab有辦法解有 x 和 y 方向的PDE嗎?! 我試過pdepe 不過他裡面是解 t 和 x 的 這樣有辦法改嗎?! 我的式子是 Vx(dCa/dx) = DA(d2Ca/dy2) kCa我為這個問題困擾很久了

有高手可以幫幫我嗎!! 感恩
我程式有寫出來囉~但發現你的k實在太大了

造成Ca瞬間達到10^178以上請問您的k是有有問題呢若沒問題

我把程式碼貼上

你再研究看看--------------------------------------------------------function dy = model(x

y

Cam)global n uy Ay By wy%% Parameter SettingVmax = 0.055; % (m/s)B = 1.5*(10^-3); % (m)k = 3808; D = 1.56*(10^-9); % (m2/s)H = 0.08317; %% Add Boundary Condition at y-DirectionCa = zeros(n 2

1);Ca(2:n 1) = y;Ca(n 2) = H*Cam;Ca(1) = -Ay(1

2:n 2)*Ca(2:n 2)/Ay(1

1);% Derivative of Ca% dCa/dydyCa = Ay(2:n 1

:)*Ca;% d2Ca/dy2 dy2Ca = By(2:n 1

:)*Ca;%% Equation Settingvx = 4*Vmax*uy(2:n 1).*(1.0 - uy(2:n 1));Deff = D/(2*B);dy = (Deff./vx).*dy2Ca k*Ca(2:n 1);--------------------------------------------------------clear allclcticglobal n uy Ay By wy%% Orthogonal Collocation Settingn = 12;n0 = 1;n1 = 1;al = 0;be = 0;[uy

Ay

By

Wy] = uAB(n

n0

n1

al

be);%% Boundary Condtion of x Direction% Initial at x = 0

Ca = 0.000308Ca0 = 0.000308;y0 = Ca0*ones(n

1);% ODE Solvertspan = linspace(0

1

10);Cam = [0.0138 0.0122 0.0107 0.0091 0.0075 0.006 0.004 0.0028 0.0013 0.000308];[x

y] = ode45(*model

[tspan(1):0.05:tspan(2)]

y0

[]

Cam(1));%%toc--------------------------------------------------------uAB的程式碼

可以在這裏找到http://tw.group.knowledge.yahoo.com/matlab-elisha0524/listitem/view?iid=1
我會用OC解

要嗎??
OC(orthogonal collocation)正交配置法

可惜不是內建Matlab內建也有解PDE的方法

例如BVP4C

或BVP5C...等或許你可以先把你的方程式貼上來大家一起研究看看喔
你的邊界是不是弄反了

應該是這樣x = L

Ca = 0y = B

Ca = HCamy=-B

dCa/dy = 0題目不難

應該很快就可以解出來了
你也把常數值告訴我吧還有Cam
我若要使用OC來解

方程式可能要改一下

你看看一下對不對喔設D = 2*B;方程式變更如下Vmax*(1-(2*y/D)^2)*(dCa/dx) = D*(d2Ca/dy2) kCaBoundary Conditionx = L

Ca = 0y = 0

dCa/dy = 0y = D

Ca = H*Cam
另外y = D

Ca = HCamCam = [0.0138 0.0122 0.0107 0.0091 0.0075 0.006 0.004 0.0028 0.0013 0.000308]其對應的x為?我僅知

x=L時 Cam = 0.0138

x=0 時 Cam = 0.000308

但其它的點也需知道喔~~
呵呵

不好意思

符號沒設好我設y = Dh

Ca = HCam好了
我沒說清楚我的Dh = 2*B;不是你方程式中的擴散係數

我重設了你的座標

所以你的邊界條件也必需一起更改
我重設的邊界條件有沒有設錯呢?x = L

Ca = 0y = 0

dCa/dy = 0y = 2*B

Ca = H*Cam

面試問題,沒問題先生,楓之谷的問題,星座 血型個性問題,今晚哪裡有問題,windows live communications platform發生問題 必須關閉 謹此致歉,感情問題,報稅問題,婆媳問題,數學問題問題,PDE,Matlab,Boundary Condition,Parameter,Derivative,程式碼,Orthogonal,Collocation,uAB

會員管理|運算|讀取|程式|字串|API|C|do-while|read|語法|迴圈|函數|堆疊|撰寫|for|進位|code|數值|題目|會員登入|C |線性|指令|矩陣|8051|陣列|連結|程式設計|