欢迎来到知嘟嘟! 联系电话:13095918853 卖家免费入驻,海量在线求购! 卖家免费入驻,海量在线求购!
知嘟嘟
我要发布
联系电话:13095918853
知嘟嘟经纪人
收藏
专利号: 2014104988547
申请人: 杭州电子科技大学
专利类型:发明专利
专利状态:无效专利
专利领域: 测量;测试
更新日期:2023-07-20
缴费截止日期: 暂无
价格&联系人
年费信息
委托购买

摘要:

权利要求书:

1.一种基于非线性拟合方式的肽质谱峰特征参数提取方法,其特征在于:设质谱图中某离子的高斯峰由N个样点组成,N≥3;对样点按其丰度值从大到小排序后,其坐标构成集合A;

A={(m1,d1),(m2,d2),…(mN,dN)}其中,mi表示质荷比,di表示丰度值,i∈{1,2,3,…,N};准备通过样点拟合出的高斯曲线其函数形式设为:其中,函数f(x,P)代表理论丰度值,自变量x代表质荷比,p1、p2和p3为待求解的高斯曲线特征参数,分别表征缩放因子、质心、标准差,构成特征参数向量P=[p1 p2 p3];

具体步骤如下:

步骤(1)根据丰度值最大的3个样点数据,对高斯曲线特征参数赋初值;

其中,ln(·)表示取自然对数操作;

步骤(2)选择合适数值初始化迭代步长参数λ,该参数初始化数值的大小将影响迭代次数和收敛速度;

步骤(3)计算拟合结果误差Err,判定迭代过程是否结束;

设定判决门限ε1,如果Err≤ε1,则处理过程结束,当前向量P中的特征参数值即为求解的最终结果;反之,如果Err>ε1,则进入步骤(4);

步骤(4)根据当前特征参数向量P,构造矩阵J;

T

步骤(5)计算每次迭代过程中,特征参数向量P的更新向量H=[Δp1 Δp2 Δp3],Δp1,Δp2和Δp3分别为特征参数p1,p2和p3的待定更新值;构造误差矢量E;

T

E=[d1-f(m1,P),d2-f(m2,P),…dN-f(mN,P)]则:

T T -1 T

H=[J×J+λ×diag(J×J)] ×J×E其中,diag(·)表示矩阵对角元素提取和创建对角阵操作;

步骤(6)计算更新向量H的度量值ρ(H);

步骤(7)更新特征参数向量P和迭代步长参数λ;设定判决门限ε2,如果更新向量H的度量值ρ(H)>ε2,则当前特征参数向量P数值由P+H替代,即P←P+H,完成更新,同时当前迭代步长参数λ数值减小至λ/K,即λ←λ/K;反之,如果ρ(H)≤ε2,则当前特征参数向量P保持不变,同时迭代步长参数λ数值增加K倍,即λ←K×λ;K为比例因子,取值范围为5~20;在完成特征参数向量P和迭代步长参数λ更新后,返回至步骤(3),进行下一轮迭代。