Boundary+Processing 源代码

function BoundaryConditions(h,N) global rou T V; %Subsonic inflow boundary(point 1) conditions rou(h,1)=1; V(h,1)=2*V(h,2)-V…

function BoundaryConditions(h,N,PN) global rou T V; %Subsonic inflow boundary(point 1) conditions rou(h,1)=1; V(h,1)=2*V(h,2…

A=[5,-4,1,0;-4,6,-4,1;1,-4,6,-4;0,1,-4,5]; b=[2,-1,-1,2]'; n=length(b); w=10; D=diag(diag(A)); CL=-triu(A,1); CLZ=CL'; L=((…

function Ae=extend(A, flag) % extend.m extends data across the boundaries (to make other codes % uniform and simp…

function rbf_db(w, t, sigma, st_sz) % w - weights % t - centers % sigma - center spread % st_sz - step size to be used in…

function [xi,yi,zi] = gridder(x,y,z,xi,yi, bndyval, colbndyfac, rowbndyfac,... bndydist) % GRIDDER Data gridding. % % [xi,yi,zi…

function [snapshot]=afd_bc_inner(delx,delt,velocity,snap1,snap2,snapshot,boundary) % AFD_BC_INNER ... apply absorbin bcs to inner…

% CREWES finite difference seismic modelling % Acoustic finite difference modelling toolbox % % Sample scripts % HIGHV_WEDGE .…

function [snapshot]=afd_bc_outer(delx,delt,velocity,snap1,snap2,snapshot,boundary) % AFD_BC_OUTER ... apply absorbing boundary co…

function [snapshot,z,x]=afd_snap(delx,delt,velocity,snap1,snap2,laplacian,boundary) % AFD_SNAP ... take one finite difference tim…