pdepe帮助!在PDE(边界条件解算器)

56 Ansichten(30天)的

Altere Kommentare anzeigen

沂州杜我2019年1月16日

0
Verknupfen

这位Frage祖茂堂Direkter链接

https://de.mathworks.com/matlabcentral/answers/440084-pdepe-help-boundry-condition-in-pde-solver

Kommentiert: 乔什·迈耶我2019年11月19日

Akzeptierte Antwort: 乔什·迈耶

pdex1.m

你好,我需要一些帮助在MATLAB解决微分方程。

我需要解决热传导方程:Cρ* dT (x, t) / dT = df (dT (x, t) / dx) / dx。

边界条件:T (x, 0) = T0;dT (L、t) / dx = 0;T (0, T) = Tg (T)[上边界条件,(在这里,Tg是instrument-recorded温度))

然而,非周期边界条件Tg (t)(温度随时间变化和与时间)没有固定的变换关系。

可以使用pdepe如何解决吗?

有人能帮助我吗? ?

这是我的代码:

1 Kommentar
Keine anzeigenKeine ausblenden

沂州杜我2019年1月16日

这是我的代码;

函数pdex1

m = 0;

x = [0 0.005 0.01 0.05 0.1 0.2 0.5 0.7 0.9 0.95 0.99 0.995 1];

t = [0 1 0.005 0.01 0.05 0.1 0.5 1.5 - 2);

C = 4.0;

k = 419;

密度= 8960;

TT = [25 30 35 40 45 50 55 60 65];

索尔= pdepe (m, @pdex1pde, @pdex1ic @pdex1bc x, t);

u1 =索尔(:,:1);

图

冲浪(x, t, u1)

标题(u1 (x, t))

包含(' x '距离)

ylabel (“t”)

% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

函数[c、f、s] = pdex01pde (x, t, u, DuDx)

c =密度* c;

f = k * DuDx;

y = u (1);

F = 0 * y;

s = F;

% % - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

函数情况= pdex1ic (x);

情况= 25;

% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

函数(pl ql,公关,qr] = pdex1bc (xl, ul, xr, ur, t)

pl = ul-TT(我);

ql = 0;

公关= 0;

qr = 1;

Melden您西奇,嗯祖茂堂kommentieren。

Melden您西奇一个,嗯这Frage祖茂堂beantworten。

Akzeptierte Antwort

乔什·迈耶我2019年2月15日

2
Verknupfen

这位Antwort祖茂堂Direkter链接

https://de.mathworks.com/matlabcentral/answers/440084-pdepe-help-boundry-condition-in-pde-solver answer_361115

你的边界条件包括线函数

pl = ul-TT(我);

TT(我) 不是内定义的函数,这就是导致错误。你需要通过工作 TT 作为一个额外的参数。但是,尽管如此,该指数我是未定义的。所以我建议做某事:

                             函数[pl, ql,公关,qr] = pdex1bc (xl, ul, xr, ur, t)
                            
                             Tg = 25;%或Tg = h (t)如果你能表达t的函数
                            
                             pl = ul-Tg;
                            
                             ql = 0;
                            
                             公关= 0;
                            
                             qr = 1;