NEWMARKmethod

内点惩罚函数法C语言程序，适用于机械行业机械优化设计(SUMT)

说明：  一般的最小二乘通常是一次拿到全部的数据， 对所有数据进行统一优化计算得到模型系数。递推最小二乘是以一种递推的方式计算最小二乘，每次使用最新的测量值，来不断更新模型系数。(General least squares usually get all the data at one time, and perform unified optimization calculations on all the data to get the model coefficients. Recursive least squares is to calculate the least squares in a recursive way, each time the latest measured value is used to continuously update the model coefficients.)

地球科学（Fortran90）: 计算空间重力异常、扰动重力、垂线偏差或重力梯度的垂向梯度。 计算方式：数值积分、二维平面FFT、二维球面FFT或一维FFT。(Earth Science (Fortran90): calculate the gravity anomaly, gravity disturbance, the deflection of the vertical or gravity gradient vertical gradient. Calculation: numerical integration, FFT of the two-dimensional plane, two-dimensional spherical FFT or one-dimensional FFT.)

简单的平面有限元程序，采用三结点三角形单元；采用等带宽存贮技术；采用高斯消元法解线性方程组。(Simple plane finite element program, using three node triangular element use of such bandwidth, storage technology using Gauss elimination method to solve linear equations.)

二维粘弹性介质中波的传播的解析形式。应用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.)

潮流计算中拉丁超立方采样的算法，可利用其进行概率潮流计算(Flow calculation Latin hypercube sampling algorithm, can use its probabilistic load flow calculation)

很详细的一维FDTD介绍，以及MATLAB程序说明(Very detailed one-dimensional FDTD presentation and description of the MATLAB program)

适合于稀薄流动计算的气体动力学格式源码，该源码为xukun老是JCP论文中的算例一致(UGKS1D- shock structure calculation using Unified Gas-Kinetic Scheme)

三种帧同步算法的MATLAB代码,具有良好的算法效果，可以使用，请大家下载(Three kinds of frame synchronization algorithm MATLAB code, a good algorithm results, you can use, please download)

曲线的高斯拟合,去除趋势性,对于一般非平稳序列的分析.(Gaussian curve fitting and removal of the trend, and for most non-stationary sequence analysis.)