• 大小: 8KB
    文件类型: .zip
    金币: 1
    下载: 0 次
    发布日期: 2021-02-02
  • 语言: Matlab
  • 标签:   

资源简介

一、椭圆型偏微分方程 1. Helmholtz方程 及 特例(Possion方程, Laplace方程); Helmholtz.m possion 2. 满足牛顿边值条件的Helmholtz方程;Helmholtz_Newton.m 二、抛物线型偏微分方程 1. 显式前向欧拉法 EF_Euler.m 2. 隐式

资源截图

代码片段和文件信息

function [uxt] = ECD_Wave(AxfTit0i1t0bx0bxfMN)
%解方程a u_xx = u_tt for 0<=x<=xf 0<=t<=T
% 初始条件: u(x0) = it0(x) u_t(x0) = i1t0(x)
% 边界条件: u(0t)= bx0(t) u(xft) = bxf(t)
% M :沿x轴的等分段数
% N :沿y轴的等分段数
dx = xf/M; x = [0:M]‘*dx;
dt = T/N; t = [0:N]*dt;
for i = 1:M + 1
    u(i1) = it0(x(i)); 
end
for k = 1:N + 1
    u([1 M + 1]k) = [bx0(t(k)); bxf(t(k))];
end
r = A*(dt/dx)^ 2; r1 = r/2; r2 = 2*(1 - r);
u(2:M2) = r1*u(1:M - 11) + (1 - r)*u(2:M1) + r1*u(3:M + 11) + dt*i1t0(x(2:M)); %(11.3.4)
for k = 3:N + 1
    u(2:Mk) = r*u(1:M - 1k - 1) + r2*u(2:Mk-1) + r*u(3:M + 1k - 1)- u(2:Mk - 2); %(11.3.3)
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         667  2007-01-21 14:33  Numerical Methods for PDE\ECD_Wave.m
     文件         578  2007-01-20 14:51  Numerical Methods for PDE\EF_Euler.m
     文件         472  2007-01-24 13:47  Numerical Methods for PDE\fem_basis_ftn.m
     文件        1072  2007-01-24 14:50  Numerical Methods for PDE\fem_coef.m
     文件         767  2007-01-20 20:10  Numerical Methods for PDE\Grank_Nicholson.m
     文件        1310  2007-01-19 10:43  Numerical Methods for PDE\Helmholtz.m
     文件        1629  2007-01-19 19:48  Numerical Methods for PDE\Helmholtz_Newton.m
     文件         663  2007-01-20 19:41  Numerical Methods for PDE\IB_Euler.m
     文件        1304  2007-01-19 10:19  Numerical Methods for PDE\possion.m
     文件        1549  2007-01-24 14:50  Numerical Methods for PDE\Show_Basis.m
     文件        1696  2007-01-21 14:00  Numerical Methods for PDE\TDE.m
     文件        1435  2007-01-21 15:45  Numerical Methods for PDE\Wave2.m

评论

共有 条评论