波偏微分方程的求解（2）

t = 0 :0.1:1;
x = 0:0.1:1;
L = 1
a = 1
for xi = x
u= sin(pi * xi / L) * cos(pi * a * t / L) + 0.5 *sin(3 * pi * xi / L) * cos(3 * pi * a * t / L);
plot (t, u);
hold on;
end

t = 0 :0.1:2;
x = 0:0.1:1;
L = 1
a = 1
for xi = x
u= sin(pi * xi / L) * cos(pi * a * t / L) + 0.5 *sin(3 * pi * xi / L) * cos(3 * pi * a * t / L);
plot (t, u);
hold on;
end

t = 0 :0.1:3;
x = 0:0.1:1;
L = 1
a = 1
for xi = x
u= sin(pi * xi / L) * cos(pi * a * t / L) + 0.5 *sin(3 * pi * xi / L) * cos(3 * pi * a * t / L);
plot (t, u);
hold on;
end

t = 0 :0.1:1;
x = 0:0.1:1;
L = 2
a = 1
for xi = x
u= sin(pi * xi / L) * cos(pi * a * t / L) + 0.5 *sin(3 * pi * xi / L) * cos(3 * pi * a * t / L);
plot (t, u);
hold on;
end

t = 0 :0.1:2;
x = 0:0.1:1;
L = 2
a = 1
for xi = x
u= sin(pi * xi / L) * cos(pi * a * t / L) + 0.5 *sin(3 * pi * xi / L) * cos(3 * pi * a * t / L);
plot (t, u);
hold on;
end

t = 0 :0.1:4;
x = 0:0.1:1;
L = 2
a = 1
for xi = x
u= sin(pi * xi / L) * cos(pi * a * t / L) + 0.5 *sin(3 * pi * xi / L) * cos(3 * pi * a * t / L);
plot (t, u);
hold on;
end

$ThesolutionisperiodicintimewithaperiodofT=\frac{2L}{\alpha }$

clear all,clc;
a=0:.01:1;
[x,t]=meshgrid(a);
u=0.5.*(sin(pi * (x+t)) + 0.5 * sin(3 * pi * (x+t))+sin(pi * (x-t)) + 0.5 * sin(3 * pi * (x-t)));
figure(1)
mesh(x,t,u);
axis([0 1 0 1 -1 1]);
title('u=0.5.*((x+t).*exp(-(x+t).^2)+(x-t).*exp(-(x-t).^2)); mesh')
colormap cool
colorbar

clear all,clc;
a=0:.01:2;
[x,t]=meshgrid(a);
u=0.5.*(sin(pi * (x+t)) + 0.5 * sin(3 * pi * (x+t))+sin(pi * (x-t)) + 0.5 * sin(3 * pi * (x-t)));
figure(1)
mesh(x,t,u);
axis([0 2 0 2 -2 2]);
title('u=0.5.*((x+t).*exp(-(x+t).^2)+(x-t).*exp(-(x-t).^2)); mesh')
colormap cool
colorbar

波偏微分方程的求解（1）
Octave online

欢迎关注二幺子的知识输出通道：