matlab

间断有限元方法（ Runge–Kutta discontinuous Galerkin ）求解N-S方程 算例为一维sod问题 流通量为BR1(Discontinuous Finite Element Method (Runge-Kutta discontinuous Galerkin) N-S equation solving a one-dimensional example problem sod Liquidity is BR1)

通过不精确的ALM算法，解析出张量矩阵中的低秩部分和稀疏部分。(via the Inexact ALM method will solve the programming.)

计算基尔霍夫近似下镜面与非镜面的反射系数，同时得到距离散射点2米的接收功率(Calculation Kirchhoff approximation of specular and non-specular reflection coefficient, while the reception power obtained from the scattering points of 2 meters)

松弛法解方程C程序(relaxation method of solving equations C program.)

说明：  Fractal分形，简短的代码，清晰可用(Fractal)

提供函数能进行1维2维的FFT，输入输出的序列长度应为2的整数次方倍(Provides functions can be one-dimensional two-dimensional FFT, sequence length input and output should be an integer power of 2 times)

burgers方程求解，采用fortran编写，有详细的说明，效果非常好(burgers equations, using fortran prepared a detailed description of the effect is very good)

二维粘弹性介质中波的传播的解析形式。应用Carcione等在1988年发表的《WAVE PROPAGATION SIMULATION IN A LINEAR VISCOELASTIC MEDIUM》一文中的方法，更加详细的介绍见英文描述。该代码对于做波场模拟的同学有很好的参考价值。( MATLABScript for obtaning the analytical solution for 2D wave propagation in a viscoelastic medium, based on the Appendix B of paper of Carcione et al. (1988)(1), corrected in (2). The source used is a Ricker wavelet in the y velocity component. A figure of the source temporal function and its inverse Fourier transform are provided to verify the good working of the method. Also works for elastic variables by substituting the values of vp and vs by real constant ones (vp_0 & vs_0). Attenuation is implemented by GMB-EK rheology presented in (4) and thoroughfully explained in (3). In order to make a proper comparison with numerical data it is recommended to either use the same rheology in both the numerical and the analytical. It is also possible to get an almost constant Q value by increasing a lot the number of mechanisms used in both the numerical and the analytical solutions, thus minimizing the effects of the Q fitting in the comparison.)

说明：  共轭梯度法求解线性方程组,主要采用FR共轭梯度法(solve liner solutions from CG)

本程序是根据jacobi过关法求实对称矩阵的全部特征值和特征向量(This procedure is based on realistic method jacobi clearance of all symmetric matrix eigenvalue and eigenvectors)