功能: 创建2D振动-多峰模型

格式: PeakVibrationCreate(n, p)

n : 一个元素个数为5的矩阵变量,里面存储多峰模型对应的[n1,n2,n3,n4,n5,n6]系数。此参数可以通过{PeakVibrationGetN<矩阵运算/PeakVibrationGetN>}获取
p : 多峰模型系数。此参数可以通过{PeakVibrationGetParam<矩阵运算/PeakVibrationGetParam>}获取。

说明:
1、本函数主要初始化如下数学模型,一般此参数最先由2D振动-多峰模型工具箱计算得来。
$$
\begin{aligned}
x(t)&=\underbrace{\alpha+\sum_{i=1}^{n_1}a_it^i}_{趋势项}
\\&+\underbrace{\sum_{i=1}^{n_2}\exp(-\beta_it)(b_i+g_it)+\sum_{i=1}^{n_3}(d_i\exp(-\gamma_it)+e_i)\cos(w_it+r_i)}_{振动项}
\\&+\underbrace{\sum_{i=1}^{n_4}\dfrac{2h_i}{\pi}\dfrac{p_i}{4(t-c_i)^2+ p_i^2} + \sum_{i=1}^{n_5}\dfrac{q_i}{\sigma_i\sqrt{2\pi}}\exp(-\dfrac{(t-u_i)^2}{2\sigma_i^2})+\sum_{i=1}^{n_6}\mathbf{weibull}(t,k_i,\lambda_i,s_i)}_{多峰项}
\\
\\
\mathbf{weibull}(t,k,\lambda,s)&=\begin{cases}
\dfrac{k}{\lambda}(\dfrac{t-s}{\lambda})^{k-1}\exp(-(\dfrac{t-s}{\lambda})^k) &,t>s
\\
0&,t\leq s
\end{cases}
\end{aligned}
$$
例子:
n = [ 1 1 2 1 1 0 ];
p =
[ 295845834.351268
 -168499.433596370
  0.08015859458975
 -462016.749564439
 -9768.24709729634
  263.611207456596
 -10.5633972790198
  0.86124003876391
  1.15360561518295
 -61.1969436802180
 -0.66865654417505
  0.32975728694554
  1.00526588654014
  2.36972517149763
  2.61681327858656
 -300418369509.681
  477.776139642666
  5.19585763105866
  23986519890.6085
  89.8494991514260
  16.4920376966502 ];
obj = PeakVibrationCreate(n, p);//创建模型
{PeakVibrationShow<矩阵运算/PeakVibrationShow>}(obj);//显示模型

y = {PeakVibrationCalcExpress<矩阵运算/PeakVibrationCalcExpress>}(obj, 1.26)//计算模型在1.26处的值
y =
[ 0.85776749253273 ]