Question

用matlab解偏微分方程组时程序怎么也调试不对，具体程序如下，望大神指点

用matlab解偏微分方程组时程序怎么也调试不对，具体程序如下，望大神指点function [ c,f,s ] = pdefun(x,t,u,Du)
k=0.07266;
E=5;
R=200;
%UNTITLED2 此处显示有关此函数的摘要
% 此处显示详细说明
c=[1;exp(2*k*x);14/17];
f=[-u(2);-u(2)^2/u(1)*exp(-2*k*x)*6/5+1/(45*7.266)*exp(-2*k*x)*(1/2*x)*Du(1);-u(2)*u(3)/u(1)];
s=[-2*k*u(2);-k/u(1)*(18/5*u(2)^2-155/126*u(3)^2)*exp(2*k*x)+3^(1/3)/(45*0.07266)*(-3*u(2)/u(1)^2+2*u(3)/E);-14/17*(34/7*k*u(2)*u(3)/u(1)+1/(45*0.07266)*(5/2*u(3)/u(1)^2+2*u(2)/E))];
end

function [u0]= pdeic(x)
%PDEIC 此处显示有关此函数的摘要
% 此处显示详细说明
k=0.07266;
u0=[3^(1/3)*exp(-2*k*x/3);exp(-2*k*x);0];
end
function [pa,qa,pb,qb] = pdebc(xa,ua,xb,ub,t)
%UNTITLED2 此处显示有关此函数的摘要
% 此处显示详细说明
R=200;
k=0.07266;
pa=[ua(1)-3^(1/3);ua(2)-1;ua(3)];
qa=[0;0;0];
pb=[ub(1)-3^(-1/3)*exp(-2*k*R/3);ub(2)-exp(-2*R*k);ub(3)-1+exp(-2*k*R)];
qb=[0,0,0];

end
clear;clc;

x=0:20:200;
t=0:0.5:60;
m=0;

sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);

u(1)=sol(:,:,1);

u(2)=sol(:,:,2);

figure;

surf(x,t,u(1))

title('u(1)(x,t)')

xlabel('Distance x')

ylabel('Time t')

figure;

surf(x,t,u(2))

title('u(2)(x,t)')

xlabel('Distance x')

ylabel('Time t')

Answer 1

用matlab求解偏微分方程组时，其程序存在几个小问题：

1、u变量不能以向量变量来表示。u(1)应改写为u1，u(2)应改写为u2

2、同理， surf(x,t,u(1))应改写为surf(x,t,u1)，surf(x,t,u(2))应改写为surf(x,t,u2)

3、更改后，可以得到其数值解

追问

你好，我尝试着做了一下，但是还是显示做不出

可以私信一下吗，如果可以，我先采纳了

或者，麻烦您把改过的代码发一下也行

追答

追问

pdefun中的u（1）需要替换吗

求全部代码

我照您的改了还是无法运行

还是报错

出错 pdepe (line 246)
[c,f,s] = feval(pde,xi(1),t(1),U,Ux,varargin{:});

出错 main (line 11)
sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);

追答

把全部命令和函数放在一个文件里，执行

你可以看看私信的内容

追问

好的