高阶统计量

说明：  9跨每跨35米桥梁响应的车桥耦合程序，使用Fortran编写(9 of 35 meters span bridge across every bridge coupled response procedures, the use of Fortran written)

一本介绍fortran语言的不错的书籍，希望你喜欢，版权归属原作者，仅供学习交流(A good book for learning Fortran language, if you read it carefully, you will get what you want.)

粒子群算法，也称粒子群优化算法或鸟群觅食算法（Particle Swarm Optimization），缩写为 PSO。它具有精度高，收敛快等优点，广泛用于系统辨识、函数优化等领域。文件pso.m是粒子群算法辨识程序，在此基础上，需要自行编写目标函数计算程序以实现不同系统的参数辨识。(Particle swarm optimization (PSO), also known as particle swarm optimization algorithm or bird swarm foraging algorithm (Particle Swarm Optimization), abbreviated as PSO. It has the advantages of high precision and fast convergence. It is widely used in system identification, function optimization and other fields. The document pso.m is a particle swarm algorithm identification program. On this basis, it is necessary to write an object function calculation program to realize the parameter identification of different systems.)

该损伤模型通过 商用有限元软件ABAQUS /Explicit 的用户子程序VUMAT 实现，并使用该模型对碳纤维增强环氧树脂 复合材料层合板在横向低速冲击作用下的损伤和变形行为进行预测分析． 数值仿真的结果与试验结 果进行了比较，吻合良好，验证了该模型的有效性．(Through coding user subroutine VUMAT of finite element software ABAQUS /explicit， and applied the model to an example，carbon fiber reinforced epoxy composite laminates under low velocity impact，the prediction for damage and deformation were good agreement with the experimental results．)

利用OPENMP实现矩阵相乘的C程序。Demonstrates a matrix multiply using OpenMP.()

T矩阵方法的FORTRAN环境下的程序，希望对大家有帮助！(T matrix method program FORTRAN environment, we want to help!)

程序是关于计算两相黎曼问题的一个求解器。可以直接设置。(The code is a solver about the two-phase Riemann problem. the user can use it directly.)

采用 Lax-Wendroff算法计算非粘Bergers方程(nonviscous Burgers equation, Lax-Wendroff method.)

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

使用匈牙利算法来解决指派问题，MATLAB代码(The Hungarian algorithm solves the assignment problem)