偏微分方程 python_有关一个偏微分方程组的求解？

题目写的看不懂。我也瞎答。

写了个MATLAB的小程序，用特征线法求解 偏微分方程组。

[XX,TT]=meshgrid(0:0.4:4,0:0.1:1);

N=size(XX,2);

T=size(XX,1);

u1=zeros(T,N);

u2=zeros(T,N);

dx=0.4;

dt=0.1;

ds=0.1;

u10=1:N;u20=(N:-1:1); % initial value

u10t=1:T;u20t=1:T; % value at a given node point

u1(1,:)=u10;

u2(1,:)=u20;

u1(:,1)=u10t;

u2(:,1)=u20t;

for tt=2:T % 特征线法求解线性偏微分方程组

for ii=2:N

u1(tt,ii)=(-400*u1(tt-1,ii-1)+400*u2(tt-1,ii-1))*ds;

u2(tt,ii)=(0.35*u1(tt-1,ii)-0.35*u2(tt-1,ii))*dt;

end

end

quiver(XX,TT,u1,u2)