[发明专利]一种睫状肌预训练方法

权利要求书

1.一种睫状肌预训练方法,其主要特征包括：

步骤1：采集近视患者视力分析样本，获取眼轴长度、角膜曲率、A计权声功率级噪声强度、下额托板高度等参数；按采集时间先后组成训练数据集，将矩阵中的后n-d个样本数据，，…，组成自回归模型的输出矩阵，自回归模型的输入矩阵Z构建方式如下所示：

其中，n为训练样本数、m为监测对象的测量变量数、R为实数集、表示n×m维的实数矩阵、d为延时测量数据的个数、N＝n-d、上标号T表示矩阵或向量的转置；

步骤2：用公式对输入矩阵Z与输出矩阵Y中的各个列向量分别作归一化处理从而消除量纲的影响，x表示矩阵Z或矩阵Y中任意一个列向量，与分别表示向量x的均值与标准差；

步骤3：训练RBF神经网络的参数，包括隐层神经元节点数k、各隐层神经元节点的径向基函数参数，，…，，具体的实施过程如下所示：

①设置隐层节点数k后，随机从输入矩阵Z中选取k个行向量分别做为各聚类簇的初始中心点向量；

②计算矩阵Z中各行向量与k个中心点之间的距离，并根据距离最小值将该行向量划分进相应的聚类簇中；

③计算各聚类簇中所有归属行向量的均值向量，该向量即为新的中心点位置；

④判断各中心点位置是否收敛，若否，则返回步骤②；若是，则将收敛后的中心点记录为，,...,，并执行步骤⑤；

⑤根据公式计算RBF参数，其中i＝1，2，…，k；

步骤4:计算输入矩阵Z中第j行的行向量经隐层第i个神经元节点转换后的输出，，其中，下标号i＝1，2，…，k，j＝1，2，…，N，相应的隐层输出向量为；

步骤5：重复步骤4直至得到Z中所有N个行向量对应的隐层输出矩阵；

步骤6：计算隐层输出S到输出层输出Y之间的回归系数矩阵B，；

步骤7：对RBF神经网络模型的输出误差E＝Y-ZB中的每列作归一化处理，得到均值为0、标准差为1的新数据矩阵；

步骤8：建立基于主成分分析算法的距离计算模型，保留模型参数集，其中P为投影变换矩阵、Λ为主成分的协方差矩阵、与分别为距离指标的控制上限D与下限Q，具体的实施过程如下：

算的协方差矩阵；

②求解C所有的特征值所对应的特征向量；

③设置保留的主成分个数为满足的最小值，并将对应的个特征向量组成载荷矩阵；

④将转变成对角矩阵，即为主成分的协方差矩阵；

⑤根据公式，分别计算距离指标的控制上限与下限，其中，表示自由度为k与N-d-k的F分布在置信度（取99%）下的取值、表示自由度为的咖方分布在置信度下的取值、加权系数、a与v分别表示Q距离指标的估计均值与估计方差；

步骤9：收集最新采样时刻的数据样本，并将其延时测量数据组成自回归模型的输入向量，其中下标号t表示当前最新采样时刻；

步骤10：对输入z与输出实施与步骤(2)中相同的归一化处理；

步骤11：将向量z输入进离线建模阶段训练得到的RBF神经网络，根据公式计算得到隐层神经元的输出向量s；

步骤12：根据公式计算RBF神经网络的输出估计误差e，并对误差e实施与步骤(7)中相同的归一化处理，得到数据向量；

步骤13：调用步骤8中保留的参数集实施在线预测，具体的实施过程包括：

①根据公式，计算距离指标D与Q的具体数值；

②根据D与Q的具体数值以及相应的距离指标的控制上下限与判断距离控制器的移动距离。