matlab pdepe求解偏微分方程

時間 2021-08-11 23:09:35

1樓:喝酸奶的一花

%這個問題和matlab的示例問題差不多,只要對c,f,s修改一下,同時有一點,t不能從0開始,不然有bug,所以設定了一個很小的數開始,結果都是差不多的

function slove_pdepe

m = 0;

x = linspace(0,1,20);

t = linspace(8e-32,2,5);

sol = pdepe(m,@pdefun,@icfun,@bcfun,x,t);

u = sol(:,:,1);

surf(x,t,u)

title('numerical solution computed with 20 mesh points.')

xlabel('distance x')

ylabel('time t')

figure

plot(x,u(end,:))

title('solution at t = 2')

xlabel('distance x')

ylabel('u(x,2)')

function [c,f,s] = pdefun(x,t,u,dudx)

c = pi^2.*sin(x).*t;

f = t.*dudx;

s = u.*x;

function u0 = icfun(x)

u0 = sin(pi*x);

function [pl,ql,pr,qr] = bcfun(xl,ul,xr,ur,t)

pl = ul;

ql = 0;

pr = pi * exp(-t);

qr = 1;

2樓:匿名使用者

先標記下,有空在幫你編寫程式

3樓:大工竹子

function pdex1

m=0;

x=linspace(0,1,20);

t=linspace(0,2,5);

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

u=sol(:,:,1);

surf(x,t,u)

figure

plot(x,u(end,:))

function[c,f,s]=pdex1pde(x,t,u,dudx)

c=pi^2;

f=dudx;

s=u^2;

function u0=pdex1ic(x)u0=sin(pi*x);

function [pl,ql,pr,qr]=pdex1bc(x1,u1,xr,ur,t)

pl=u1;

ql=0;

pr=pi*exp(-t);

qr=1;

matlab怎麼求解偏微分方程組啊,先

4樓:展得囊繁

pdepe()函式的一般呼叫格式是:sol=pde...

1例題:解下列偏微分方程組,方程如下圖...

2pdefun的命令為:function[c,f,s]=pdef...

3初始條件:初始條件標準形式:[u1;u2]=...

4邊界條件:左邊界:[0;u2]+[1;0].*f=[0...

5求解偏微分方程組計算程式如下所示:>>...

急求用matlab求解一階偏微分方程組的程式! 100

5樓:匿名使用者

3樓只是拋物線方程組求解而已,樓主求雙曲的

6樓:匿名使用者

用parabolic函式、matlab本身的偏微分函式庫有很多的,這個最常用。

matlab求偏微分方程 20

7樓:linuxc語言

用 ode 系列求偏微分的數值解,,用dsolve求解析解這是ode的

建立m函式檔案

function dx=myfun(t,x)dx(1)=20/7*x(1)-x(2)*x(3);

dx(2)=-10*x(2)+x(1)*x(3);

dx(3)=-4*x(3)+x(1)*x(2);

dx=dx(:);

matlab視窗主程式

x0=[3,-4,2];

t0=0.001:0.001:20;

[t,x]=ode45('myfun',[0.001,20],x0); %ode45會自動調整步長

x(:,:,1)

plot(t,x)

legend('x','y','z')

dsovle不是很難,,自己好好看看,,,呵呵

為什麼多元函式的x,y偏導數連續就可微

為什麼偏導數連續是可微的充分不必要條件 1 偏導數連續是可微分充分條件,但不是必要條件。2 比如下面這個函式f x,y 函式的表示式為當x,y均為有理數時f x,y x 2 y 2 當x,y中有一個變數為無理數時f x,y 0。3 考慮這個函式在 0,0 處的微分,顯然 u f x,y f 0,0 ...

可導一定可微,可微一定可導嗎,為什麼偏導數存在不一定可微?

月似當時 可微一定可導,可導不一定可微,各變數在此點的偏導數存在為其必要條件,其充要條件還要加上在此函式所表示的廣義面中在此點領域內不含有 洞 存在,可含有有限個斷點。在一元函式中,可導與可微等價。一元函式中可導與可微等價,它們與可積無關。多元函式可微必可導,而反之不成立。即 在一元函式裡,可導是可...

高分求解關於二階連續偏導數問題的解法,有追加

關於那個古怪的符號怎麼打我們就不說了,呵呵,就用6表示嘛,大家知道意思就行了。其實那個 平方 很簡單的,估計您好像時間緊沒有看到教材在某一個角落中是說到這個問題了的,舉個例子,6z 6x表示函式z對其中的一個自變數x求導,這個時候把y看做常數就行了,想必方法你也很熟悉,那麼6 z 6y 就表示的是函...