UPW_utux0_2.rar - function [ue,un]=UPW_utux0_2(v,dt,t)
一个简单的双曲型偏微分方程:
ut + ux = 0
初始条件为:
u(x,0) = exp[-10(4x-1)^2]
边界条件为:
u(-1,t)=0,u(1,t)=0.
本题要求:
使用迎风格式,选择 v=0.5, 计算并画出当dt=0.01和0.0025时,
方程在t=0.5,x在(-1,1)时的数值解和精确解
输入:
v--即a*dt/dx
dt--数值格式的时间步
t--要求解的时间
输出:
ue--在时间t时的1×N精确解矩阵
un--在时间t时的1×N数值解矩阵
输出图像:
精确解和数值解的图像
,2014-07-02 20:12:34,下载1次
UPW_utux0.rar - function [ue,un]=UPW_utux0(v,dt,t)
一个简单的双曲型偏微分方程:
ut + ux = 0
初始条件为:
u(x,0) = 1, x≤0
0, x>0.
边界条件为:
u(-1,t)=1,u(1,t)=0.
本题要求:
使用迎风格式,选择 v=0.5, 计算并画出当dt=0.01和0.0025时,
方程在t=0.5,x在(-1,1)时的数值解和精确解
输入:
v--即a*dt/dx
dt--数值格式的时间步
t--要求解的时间
输出:
ue--在时间t时的1×N精确解矩阵
un--在时间t时的1×N数值解矩阵
输出图像:
精确解和数值解的图像,2014-07-02 20:11:40,下载7次
LW_utux0_3.rar - function un=LW_utux0_3(dx,t)
Burgers equation:
ut + (1/2*u^2)x = 0
初始条件为:
u(x,0) = exp[-10(4x-1)^2]
边界条件为:
u(0,t)=0,u(1,t)=0
本题要求:
使用Lax-Windroff格式,选择 dx=0.01, 计算并画出当
t=0.15,和t=0.3时的数值解
输入:
dx--数值格式的x轴上的分割
r--r=dt/dx,本题预设r=0.5
t--要求解的时间
输出:
un--在时间t时的1×N数值解矩阵
输出图像:
数值解的图像,2014-07-02 20:09:57,下载7次
LW_utux0_2.rar - function [ue,un]=LW_utux0_2(v,dt,t)
一个简单的双曲型偏微分方程:
ut + ux = 0
初始条件为:
u(x,0) = exp[-10(4x-1)^2]
边界条件为:
u(-1,t)=0,u(1,t)=0.
本题要求:
使用Lax-Windroff格式,选择 v=0.5, 计算并画出当dt=0.01和0.0025时,
方程在t=0.5,x在(-1,1)时的数值解和精确解
输入:
v--即a*dt/dx
dt--数值格式的时间步
t--要求解的时间
输出:
ue--在时间t时的1×N精确解矩阵
un--在时间t时的1×N数值解矩阵
输出图像:
精确解和数值解的图像,2014-07-02 20:08:36,下载3次
LW_utux0.rar - function [ue,un]=LW_utux0(v,dt,t)
一个简单的双曲型偏微分方程:
ut + ux = 0
初始条件为:
u(x,0) = 1, x≤0
= 0, x>0.
边界条件为:
u(-1,t)=1,u(1,t)=0.
本题要求:
使用Lax-Windroff method,选择 v=0.5, 计算并画出当dt=0.01和0.0025时,
方程在t=0.5,x在(-1,1)时的数值解和精确解
输入:
v--即a*dt/dx
dt--数值格式的时间步
t--要求解的时间
输出:
ue--在时间t时的1×N精确解矩阵
un--在时间t时的1×N数值解矩阵
输出图像:
精确解和数值解的图像,2014-07-02 20:04:06,下载5次