photonic_crystal_method_of_characteristic_matrix_1

说明：  MATLAB使用欧拉Euler法求解微分方程组 源程序代码(SOURCE PROGRAM CODE OF DIFFERENTIAL EQUATIONS SOURCE SOLUTION USING EULER METHOD IN MATLAB)

岩石损伤与断裂力学一直是有限元的弱势领域，必须依赖MATLAB一类的数学软件协助，本书籍提供了基础教程与相关案例。(Rock damage and fracture mechanics have been the weak field of finite element, must rely on MATLAB class of mathematical software to help, this book provides basic tutorials and related cases.)

matpower2机5节点程序进行潮流计算(matpower2 Machine 5 node program)

此程序用于计算矩形截面的弯矩曲率特性以及不同极限状态的识别(This program is used to calculate the moment curvature characteristic of the rectangular cross-section, and the identification of different limit state)

对信号的FFT（快速傅里叶变换）变换，针对LabView内部信号。(The signal FFT (Fast Fourier Transform) transformation, the internal signal for LabView.)

说明：  非常好的地形改正程序，只需要有点位的经纬度，重力值和高程，内附例子和说明。(Program to compute terrain effects on gravimetric quantities.various modes of terrain effects may be computed)

基于CPLEX半连续变量求解发电机组单时段经济调度问题，给定负荷需求，求每台发电机输出P=0或Pmin<=P<=Pmax,使发电成本最小。(Semi-continuous variable based on CPLEX to solve a single period of economic generator scheduling problem, given the load demand, seeking each generator output P = 0 or Pmin <= P <= Pmax, the minimum cost of power generation.)

test1.txt 第一列为时间t，0.001s为步长；第二列为系统的输入信号；第三列为系统的输出信号。 signal_analyzer.m首先将test1.txt 的数据读入workspace，然后进行FFT，分析出两信号的频率，输出对于输入的幅值衰减（db），输出对输入的相位差。(test1.txt the first time as t, 0.001s for the step the second as the system input signal the third as the system output signal. signal_analyzer.m first test1.txt data read workspace, and then proceed to FFT, analysis of the two signals in the frequency, the output for input amplitude attenuation (db), the output of the input phase difference.)

计算偏心圆导热计算以及同位网格上的SIMPLE算法 (Calculating thermal calculations and eccentric position on the grid with SIMPLE algorithm)

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