An-Improved-Algorithm-for-Otsu

说明：  采用电磁仿真软件HFSS进行雷达散射截面(RCS)计算实例.(Monostatic Radar Cross Section)

利用c语言编写了迭代法解线性方程组的源程序，大大简化了计算量。(Use the C language source code of the iterative method for solving linear equations, which greatly simplify the calculation, suitable for solution of complex equations)

压缩包中包含五个CFD基础程序，其中MacCormack1DSolveRiemann.for是一维欧拉方程求解（黎曼问题），UpwindTVD_1D.for是一维可压粘性N-S方程求解，MAC-Chorin2D.for是利用MAC算法和Chorin压力迭代解法求解 二维不可压缩黏性平板间Poisuille流动问题，Lax-Webdroff2D.for是利用Lax-Webdroff差分格式求解二维平面激波反射问题，fvm_upwind_MAC_couette.for是以一阶迎风型离散格式和Chorin压力迭代求解couette二维流动问题。(CFD compression package contains five basic program, which is an Euler equation solving MacCormack1DSolveRiemann.for (Riemann problem), UpwindTVD_1D.for is one-dimensional compressible viscous NS equation solving, MAC-Chorin2D.for is the use of algorithms and MAC Chorin pressure iterative solution for solving the two-dimensional incompressible flow problems between viscous flat Poisuille, Lax-Webdroff2D.for using Lax-Webdroff difference scheme for solving the two-dimensional plane shock wave reflection problem, fvm_upwind_MAC_couette.for based on first order upwind scheme and Discrete Chorin pressure couette iterative solution of two-dimensional flows.)

contourlet变换工具箱，附带贝叶斯阈值去噪(Toolbox contourlet transform, with Bayesian threshold denoising)

提取nastran“.f06”结果中的节点频响函数数据，可简单修改提取模态频率和 其它数据(Extract nastran, " f06" results in the the node frequency response function data, can be simple modification to extract the modal frequencies and other data)

K均值算法实现二维点的聚类，可以手动选择分的类别数(K mean algorithm for clustering two dimensional points)

使用Matlab编程，对教材74页例3.3.1进行编程仿真。其中，正电压A的值、噪声方差值、每个码元周期内的采样点数N自行设定（可设置为可调的变量）。噪声值可使用Matlab中产生高斯随机数的函数来进行仿真。贝叶斯检测判决式中，先验概率P(H1)=P(H0)，判断错误的代价因子设为1，判断正确的代价因子设为0. 按照上述参数的设定进行仿真，实现对仿真数据的贝叶斯检测判决；循环生成仿真数据，并进行判决结果的统计，即记录判决正确的次数，估计正确判决的概率。(The use of Matlab programming, teaching materials 74 pages of 3.3.1 programming simulation. Among them, the value of positive voltage A, the variance of noise, the number of sampling points in each symbol cycle (N) can be set by itself (adjustable variable). Noise values can be simulated using the function of generating Gauss random numbers in Matlab. In the Bias detection decision, the prior probability P (H1) =P (H0), the factor of error judgment is set to 1, and the correct cost factor is set to 0. According to the simulation parameters and realize the judgment of Bayesian detection simulation data; generating cycle simulation data, statistics and decision, which records the correct number of judgment, the estimated probability of correct decision.)

首先根据高斯白噪声频率充满整个频带的特性，通过双树复小波包变换对高斯白噪声进行分解，利用频带能量泄漏的定量分析方法，验证了双树复小波包变换具有较低的频带能量泄漏特性；其次利用双树复小波包变换逐层分解信号，对每层分解所得分量求其FFT谱的峭度，得到基于双树复小波包变换的谱峭度图，根据图中峭度最大的原则，可以自动准确的选择信号分解最佳层数和最佳分量；最后将基于双树复小波包变换的谱峭度图的故障诊断方法应用于实际工程中，对齿轮故障振动信号进行分析，选择最佳分解层数和分量后利用希尔伯特包络解调，有效准确地提取了故障特征信息，验证了方法的可行性和有效性(The parameters of a filter were determined by experience, and that has a great influence on the results of signal processing. The discrete wavelet packet transform has a larger energy leakage of frequency band, which obviously affected the results of the envelope demodulation. It is necessary to have a method with a lower energy leakage of the frequency band before envelope demodulation. The dual tree complex wavelet packet transform (DT-CWPT) was a new signal processing method that had many good qualities. Because the energy leakage of the frequency band was smaller when the signal was decomposed by a dual tree complex wavelet packet transform, the dual tree complex wavelet packet transform was used to extract the fault feature information in the field of fault diagnosis. In this paper, first, according to the characteristics of Gaussian white noise, whose frequency was full of the whole frequency band, Gaussian white noise was decomposed by a dual-tree complex wavel)

control souce code with LQR

地震数据球面扩散振幅补偿程序，fortran源代码。(sphdiv for seismic explorer)