CN111125625A - 一种基于嵌入式***的光谱基线计算方法 - Google Patents
一种基于嵌入式***的光谱基线计算方法 Download PDFInfo
- Publication number
- CN111125625A CN111125625A CN201911206716.6A CN201911206716A CN111125625A CN 111125625 A CN111125625 A CN 111125625A CN 201911206716 A CN201911206716 A CN 201911206716A CN 111125625 A CN111125625 A CN 111125625A
- Authority
- CN
- China
- Prior art keywords
- vector
- matrix
- value
- spectrum
- baseline
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于嵌入式***的快速光谱基线计算方法,包括如下步骤:(1)权重向量和调节系数初始化;(2)更改背景拟合算法;(3)引入新矩阵存储2n+1对角矩阵;(4)对矩阵进行压缩分解;(5)求解Lx=b;(6)求解Uz=x;(7)计算背景基线相关信息;(8)计算新的权重向量;(9)权重向量判断;(10)背景基线提取失效处理。本发明通过更改背景拟合算法及采用压缩空间下的矩阵分解,节省了数据存储空间并提高了计算速度,使一种无需人工干预的快速光谱基线计算方法能够在嵌入式***中得到应用。本发明也可适用于惩罚最小二乘法及在其基础上衍生的其他算法在嵌入式***中的应用。
Description
技术领域
本发明涉及一种光谱基线计算方法,尤其涉及一种数据存储空间有限情况下的快速光谱基线计算方法。
背景技术
随着光谱检测技术越来越广泛的应用于分析化学、物理、材料等领域。各种背景、噪声等引起的光谱干扰,使光谱信号信噪比和测量精度下降,已严重影响到物质的定量分析及设备标定。为降低测量光谱的背景噪声干扰,常需对光谱信号进行背景基线修正,以得到更好的定量分析结果。为实现背景基线的自动计算减少人工干预,基于惩罚最小二乘法衍生出的不对称最小二乘法和自适应迭代重加权惩罚最小二乘法等背景基线修正方法,已广泛应用于各种光谱分析。
惩罚最小二乘法是在标准最小二乘基础上加上平滑度惩罚项,公式如下:
S(z)=(y-z)T(y-z)+λzTDTDz
其中,y为长度为N的等间隔光谱采样信号;z为经过平滑滤波后的光谱拟合信号;y和z均为N×1的列向量;式中第一项(y-z)T(y-z)为标准最小二乘,通过误差平方和来表示z相对于y的保真度;式中第二项zTDTDz通过差分平方和来表示z的平滑度;λ作为调节系数,用来保证拟合信号z在保真度和平滑度之间达到平衡;λ值越大拟合信号z越平滑。D为微分矩阵,阶数可根据需要设置;常为一阶或二阶微分矩阵。
z=(I+λDTD)-1y
惩罚最小二乘法常用于信号的平滑滤波。为实现对带有多个峰值点的光谱信号进行背景基线提取,在惩罚最小二乘法基础上,通过引入权重向量w,来实现背景基线提取。引入权重向量后,S(z)公式如下:
S(z)=(y-z)TW(y-z)+λzTDTDz
z=(W+λDTD)-1Wy
其中,w为N×1的权重向量;W为基于w的N×N对角矩阵。可通过改变权重向量值,使背景基线不随信号中的峰值点变化。
不对称最小二乘法和自适应迭代重加权惩罚最小二乘法提供了两种不同的权重向量设置法。
在不对称最小二乘法中,权重向量w设置如下:
其中,p常设为0.001~0.1。该方法在纯背景区域权重值一样,无法根据拟合值和原始数据的差异灵活设定。当峰值点较小时,造成背景扣除误差过大。
在自适应迭代重加权惩罚最小二乘法中,权重向量w设置如下:
其中,t为当前的迭代次数,d为y-z中负元素的集合。当迭代次数达到设置值或满足|d|<0.001×|y|时,迭代停止,得到背景基线。
在嵌入式***进行光谱基线计算时,若直接对矩阵求逆,其时间复杂度为 O(N·N!);若采用高斯消去法,其时间复杂度为O(N3/3);其存储空间复杂度为 O(N2)。由于嵌入式***的数据存储空间和工作频率常受硬件限制,随着光谱测量分析设备测量种类增多,扫描覆盖的波长范围变宽,采集的光谱数据增加,其时间和空间复杂度成倍增加;如光谱数据为2500,其计算所需存储空间至少为6.25兆,而常用嵌入式***内部数据空间不超过1兆,使上述光谱基线计算方法很难在嵌入式***中实现。
发明内容
本发明的技术解决问题是:针对现有技术的使用局限,本发明提供了一种基于嵌入式***的快速光谱基线计算方法,基于算法本身特点,采用压缩空间下的矩阵分解法,解决光谱基线自动计算方法在嵌入式***的应用,解决在数据存储空间和计算时间受限情况下,光谱基线的实时计算。
本发明的技术解决方案是:一种基于嵌入式***的快速光谱基线计算方法,该方法包括如下步骤:
(1)、构建光谱基线拟合信号z表达式为z=(Q+T)-1Qy,并建立光谱基线拟合方程z=A-1b,其中,A=Q+T为N×N的上下带宽为n的2n+1对角矩阵, b=Qy为N×1的列向量,为N×N的对角矩阵,W为基于wt的N×N对角矩阵,T=DTD为N×N的上下带宽为n的2n+1对角矩阵,D为n阶微分矩阵;λ为调节系数,wt为权重向量;
(2)、对光谱信号进行等间隔采样,得到长度为N的光谱数字采样信号y;
(3)、初始化调节系数λ,令迭代次数t为0,权重向量wt初始化为全1的 N×1列向量;
(4)、求解光谱基线拟合方程z=A-1b,得到光谱基线拟合信号z,若出现光谱基线拟合方程无法求解的情况时,则回到步骤(2),重新获取光谱数字采样信号y;否则,进入步骤(5);
(5)、假设d=y-z,取向量d中小于0的元素组成向量d-,对向量d-求平均m和标准差σ;
(6)、采用如下公式得到更新后的权重向量wt+1,公式如下:
(7)、若权重向量wt+1满足如下权重向量判别条件,则完成光谱基线的提取;否则,将迭代次数t加1更新迭代次数,同时更新后的权重向量wt+1替代之前的权重向量wt,跳至步骤(4)重新执行;
权重向量判别条件如下:
θ为权重比例系数;
(8)、若步骤(4)~(7)迭代次数t达到预设的门限,仍无法提取出有效的光谱基线,降低调节系数λ值,令迭代次数t为0,重新执行步骤(4)~(7),直至提取出有效的光谱基线。
所述步骤(4)求解光谱基线拟合方程的具体实现为:
(4.1)、将矩阵A的2n+1对角数据按如下规则转换为矩阵C的元素,矩阵C 和A中的元素对应关系为ci-j+n,j=aij,i=0~N-1,j=0~N-1:
(4.2)、将矩阵A做压缩空间下的三角分解A=LU,L为下半带宽为n的单位下三角矩阵,U为上半带宽为n的上三角矩阵;
(4.3)、将b=Qy代入方程Lx=b,求解该方程得到向量x;
(4.4)、将步骤(4.3)得到的向量x代入方程Uz=x,求解方程得到拟合信号向量z。
分解后的矩阵L和U存储在矩阵C中:
其中,li,j为L的元素,其和C中元素对应关系为ci-j+n,j=lij;ui,j为U的元素,其和C中元素对应关系为ci-j+n,j=uij。
所述步骤(4.2)的具体实现为:
(4.2.1)、令k=1;
(4.2.2)、采用如下公式计算矩阵U的第k-1行数据,并将矩阵U中的数据放入矩阵C(2n+1)×N中的前n+1行,公式如下:
每计算完成一个c′k-jj+n,jj-1值,将c′k-jj+n,jj-1值赋值给ck-jj+n,jj-1;
(4.2.3)、采用如下公式计算矩阵L的第k-1列数据,并将矩阵L中的数据放入矩阵C(2n+1)×N中的最后n行,公式如下:
每计算完成一个c′ii-k+n,k-1值,将c′ii-k+n,k-1值赋值给cii-k+n,k-1;当cn,k-1为0时,说明光谱基线拟合方程无法求解,直接结束矩阵分解,不再执行后续步骤(4.2.4);
(4.2.4)、k值加1,若k>N时,则完成矩阵分解;否则,跳至步骤(4.2.2) 继续执行。
所述n取值为1~N-1。
所述步骤(4.3)中所述向量x保存在向量b的存储空间内,求解向量x的具体实现为:
(4.3.1)、通过公式x0=b0计算向量x的元素第一行元素x0;
(4.3.2)、令i′=2;
(4.3.3)、通过如下公式计算向量x的第i′-1个元素xi′-1:
(4.3.4)、将向量x的第i′-1个元素xi′-1存储至向量b的第i′-1个元素;
(4.3.5)、将i′值加1,更新;
(4.3.6)、若i′>N,则完成向量x的计算;否则,跳至步骤(4.3.3)继续执行。
所述步骤(4.4)中所述拟合信号向量z保存在向量b的存储空间内,求解拟合信号向量z的具体实现为:
(4.4.1)、计算拟合信号向量z的第N-1个采样值zN-1,具体计算公式为: zN-1=bN-1/cn,N-1;
(4.4.2)、将拟合信号向量z的第N-1个采样值zN-1存储至向量b的第N-1个元素的存储空间中,即:bN-1=zN-1;
(4.4.3)、令i″=N-1;
(4.4.4)、采用下列公式计算拟合信号向量第i″-1个元素zi″-1的计算公式为:
(4.4.5)、将拟合信号向量第i″-1个元素zi″-1存储至向量b的第i″-1个元素对应的存储空间bi″-1;
(4.4.6)、i″值减1;
(4.4.7)、若i″<1,则完成向量z的计算;否则,跳至步骤(4.4.4)继续执行。
所述的步骤(1)中调节系数λ的取值范围为104~105。
所述的步骤(7)中权重比例系数θ的取值范围为10-5~10-6。
本发明与现有技术相比的有益效果是:
(1)、本发明通过更改背景拟合算法,将2n+1对角矩阵乘法改为对角矩阵乘,在一次迭代过程中减少了n(2N-n-1)次乘法计算,节省运算时间。
(2)、本发明通过引入矩阵C将N×N的2n+1对角矩阵按照对应关系压缩放入(2n+1)×N的C矩阵,节省了(N-2n-1)×N的数据存储空间。通过采用压缩分解法,在拟合信号计算过程中仅占用了(2n+2)×N的数据存储空间,即矩阵C和向量b的数据空间,大大节省了拟合信号计算过程中的数据存储空间占用。同时,该算法仅对对角内的数据进行计算处理,使时间复杂度呈线性,使计算速度得到极大提升。通过节省数据存储空间和提高计算速度使一种快速光谱基线计算方法能够在嵌入式***中得到应用。
(3)、本发明中的拟合信号计算方法,也可适用于惩罚最小二乘法及在其基础上衍生的其他算法在嵌入式***中的应用。
附图说明
图1为本发明实施例的基于嵌入式***的快速光谱基线计算方法流程图;
图2为本发明实施例矩阵分解的流程图;
图3为本发明实施例的光谱信号、背景基线和去背景后的光谱信号示意图。
具体实施方式
以下结合附图和具体实施例对本发明进行详细说明。
本发明提供了一种基于嵌入式***的快速光谱基线计算方法,该方法如图1所示,包括如下步骤:
(1)、根据公式z=(W+λDTD)-1Wy,构建光谱基线拟合信号z表达式为 z=(Q+T)-1Qy,以减少矩阵的运算量,并建立光谱基线拟合方程z=A-1b,其中, A=Q+T为N×N的上下带宽为n的2n+1对角矩阵,b=Qy为N×1的列向量,为N×N的对角矩阵,W为基于wt的N×N对角矩阵,T=DTD为N×N 的上下带宽为n的2n+1对角矩阵,D为n阶微分矩阵;λ为调节系数,取值范围为104~105。wt为权重向量;所述n取值范围为1~N-1,一般取1或者2,N取250~2500。
(2)、对光谱信号进行等间隔采样,得到长度为N的光谱数字采样信号y;
(3)、初始化调节系数λ,令迭代次数t为0,权重向量wt初始化为全1的 N×1列向量;
(4)、求解光谱基线拟合方程z=A-1b,得到光谱基线拟合信号z,若出现光谱基线拟合方程无法求解的情况时,则回到步骤(2),重新获取光谱数字采样信号y;否则,进入步骤(5);
所述步骤(4)求解光谱基线拟合方程的具体实现为:
(4.1)、将矩阵A的2n+1对角数据按如下规则转换为矩阵C的元素,矩阵C 和A中的元素对应关系为ci-j+n,j=aij,i=0~N-1,j=0~N-1:
即
这样,矩阵C(2n+1)×N用于存储矩阵A的2n+1对角数据,对角以外数据不作保存,以节省数据的存储空间。
(4.2)、将矩阵A做压缩空间下的三角分解A=LU,L为下半带宽为n的单位下三角矩阵,U为上半带宽为n的上三角矩阵;
为节省数据存储空间,本发明直接将分解后的矩阵L和U存储在矩阵C中:
其中,li,j为L的元素,其和C中元素对应关系为ci-j+n,j=lij;ui,j为U的元素,其和C中元素对应关系为ci-j+n,j=uij;
如图2所示,具体的LU分解过程为:
(4.2.1)、令k=1;
(4.2.2)、采用如下公式计算矩阵U的第k-1行数据,并将矩阵U中的数据放入矩阵C(2n+1)×N中的前n+1行,公式如下:
每计算完成一个c′k-jj+n,jj-1值,将c′k-jj+n,jj-1值赋值给ck-jj+n,jj-1;
(4.2.3)、采用如下公式计算矩阵L的第k-1列数据,并将矩阵L中的数据放入矩阵C(2n+1)×N中的最后n行,公式如下:
每计算完成一个c′ii-k+n,k-1值,将c′ii-k+n,k-1值赋值给cii-k+n,k-1;当cn,k-1为0时,说明光谱基线拟合方程无法求解,直接回到步骤(2),重新获取光谱数字采样信号y;
(4.2.4)、k值加1,若k>N时,则完成矩阵分解;否则,跳至步骤(4.2.2) 继续执行。
(4.3)、将b=Qy代入方程Lx=b,求解该方程得到向量x;所述向量x保存在向量b的存储空间内,求解向量x的具体实现为:
(4.3.1)、通过公式x0=b0计算向量x的元素第一行元素x0;
(4.3.2)、令i′=2;
(4.3.3)、通过如下公式计算向量x的第i′-1个元素xi′-1:
(4.3.4)、将向量x的第i′-1个元素xi′-1存储至向量b的第i′-1个元素;
(4.3.5)、将i′值加1,更新;
(4.3.6)、若i′>N,则完成向量x的计算;否则,跳至步骤(4.3.3)继续执行。
(4.4)、将步骤(4.3)得到的向量x代入方程Uz=x,求解方程得到拟合信号向量z。
所述拟合信号向量z保存在向量b的存储空间内,求解拟合信号向量z的具体实现为:
(4.4.1)、计算拟合信号向量z的第N-1个采样值zN-1,具体计算公式为: zN-1=bN-1/cn,N-1;
(4.4.2)、将拟合信号向量z的第N-1个采样值zN-1存储至向量b的第N-1个元素的存储空间中,即:bN-1=zN-1;
(4.4.3)、令i″=N-1;
(4.4.4)、采用下列公式计算拟合信号向量第i″-1个元素zi″-1的计算公式为:
(4.4.5)、将拟合信号向量第i″-1个元素zi″-1存储至向量b的第i″-1个元素对应的存储空间bi″-1;
(4.4.6)、i″值减1;
(4.4.7)、若i″<1,则完成向量z的计算;否则,跳至步骤(4.4.4)继续执行。
(5)、假设d=y-z,取向量d中小于0的元素组成向量d-,对向量d-求平均m和标准差σ;
(6)、采用如下公式得到更新后的权重向量wt+1,公式如下:
(7)、若权重向量wt+1满足如下权重向量判别条件,则完成光谱基线的提取;否则,将迭代次数t加1更新迭代次数,同时更新后的权重向量wt+1替代之前的权重向量wt,跳至步骤(4)重新执行;
权重向量判别条件如下:
θ为权重比例系数,取值范围为10-5~10-6。
(8)、若步骤(4)~(7)迭代次数t达到预设的门限,仍无法提取出有效的光谱基线,降低调节系数λ值,令迭代次数t为0,重新执行步骤(4)~(7),直至提取出有效的光谱基线。
实施例:
一种基于嵌入式***的快速光谱基线计算方法,通过采用压缩空间矩阵分解法,在(2n+2)×N的数据存储空间内实现光谱基线计算;通过迭代计算,不断更新权重向量wt,最终得到光谱基线拟合信号z。
如图3所示,假设y为长度N=250的等间隔采样光谱信号,采用n=2的二阶微分矩阵,进行基线拟合方程求解。如图1所示,其具体步骤如下:
(1)、根据公式z=(W+λDTD)-1Wy,构建光谱基线拟合信号z表达式为 z=(Q+T)-1Qy,以减少矩阵的运算量,并建立光谱基线拟合方程z=A-1b,其中, A=Q+T为250×250的上下带宽为2的五对角矩阵,b=Qy为250×1的列向量,为250×250的对角矩阵,W为基于wt的250×250对角矩阵,T=DTD为250×250的上下带宽为2的五对角矩阵,D采用二阶微分矩阵,即所述n取值为 2;λ为调节系数,取值范围为104~105;wt为权重向量。
(2)、如图3所示,对光谱信号进行等间隔采样,得到长度N=250的光谱数字采样信号y;
(3)、调节系数λ初始化为1×105,令迭代次数t为0,权重向量wt初始化为全1的250×1列向量;
(4)、求解光谱基线拟合方程z=A-1b,得到光谱基线拟合信号z,若出现光谱基线拟合方程无法求解的情况时,则回到步骤(2),重新获取光谱数字采样信号y;否则,进入步骤(5);
所述步骤(4)求解光谱基线拟合方程的具体实现为:
(4.1)、将矩阵A的五对角数据按如下规则转换为矩阵C的元素,矩阵C和 A中的元素对应关系为ci-j+2,j=aij,i=0~249,j=0~249:
即
这样,矩阵C5×250用于存储矩阵A的五对角数据,对角以外数据不作保存,以节省数据的存储空间。
(4.2)、将矩阵A做压缩空间下的三角分解A=LU,L为下半带宽为2的单位下三角矩阵,U为上半带宽为2的上三角矩阵;
为节省数据存储空间,本发明直接将分解后的矩阵L和U存储在矩阵C中:
其中,li,j为L的元素,其和C中元素对应关系为ci-j+2,j=lij;ui,j为U的元素,其和C中元素对应关系为ci-j+2,j=uij;
如图2所示,具体的LU分解过程为:
(4.2.1)、令k=1;
(4.2.2)、采用如下公式计算矩阵U的第k-1行数据,并将矩阵U中的数据放入矩阵C5×250中的前三行,公式如下:
每计算完成一个c′k-jj+n,jj-1值,将c′k-jj+n,jj-1值赋值给ck-jj+n,jj-1;
(4.2.3)、采用如下公式计算矩阵L的第k-1列数据,并将矩阵L中的数据放入矩阵C5×250中的最后两行,公式如下:
每计算完成一个c′ii-k+n,k-1值,将c′ii-k+n,k-1值赋值给cii-k+n,k-1;当cn,k-1为0时,说明光谱基线拟合方程无法求解,直接回到步骤(2),重新获取光谱数字采样信号y;
(4.2.4)、k值加1,若k>N时,则完成矩阵分解;否则,跳至步骤(4.2.2) 继续执行。
(4.3)、将b=Qy代入方程Lx=b,求解该方程得到向量x;所述向量x保存在向量b的存储空间内,求解向量x的具体实现为:
(4.3.1)、通过公式x0=b0计算向量x的元素第一行元素x0;
(4.3.2)、令i′=2;
(4.3.3)、通过如下公式计算向量x的第i′-1个元素xi′-1:
(4.3.4)、将向量x的第i′-1个元素xi′-1存储至向量b的第i′-1个元素;
(4.3.5)、将i′值加1,更新;
(4.3.6)、若i′>N,则完成向量x的计算;否则,跳至步骤(4.3.3)继续执行。
(4.4)、将步骤(4.3)得到的向量x代入方程Uz=x,求解方程得到拟合信号向量z。
所述拟合信号向量z保存在向量b的存储空间内,求解拟合信号向量z的具体实现为:
(4.4.1)、计算拟合信号向量z的第N-1个采样值zN-1,具体计算公式为: zN-1=bN-1/cn,N-1;
(4.4.2)、将拟合信号向量z的第N-1个采样值zN-1存储至向量b的第N-1个元素的存储空间中,即:bN-1=zN-1;
(4.4.3)、令i″=N-1;
(4.4.4)、采用下列公式计算拟合信号向量第i″-1个元素zi″-1的计算公式为:
(4.4.5)、将拟合信号向量第i″-1个元素zi″-1存储至向量b的第i″-1个元素对应的存储空间bi″-1;
(4.4.6)、i″值减1;
(4.4.7)、若i″<1,则完成向量z的计算;否则,跳至步骤(4.4.4)继续执行。
(5)、假设d=y-z,取向量d中小于0的元素组成向量d-,对向量d-求平均m和标准差σ;
(6)、采用如下公式得到更新后的权重向量wt+1,公式如下:
(7)、若权重向量wt+1满足如下权重向量判别条件,则完成光谱基线的提取;否则,将迭代次数t加1更新迭代次数,同时更新后的权重向量wt+1替代之前的权重向量wt,跳至步骤(4)重新执行;
权重向量判别条件如下:
权重比例系数θ为10-5。
(8)、若步骤(4)~(7)迭代次数t达到50次后,仍无法提取出有效的光谱基线,降低调节系数λ值,令迭代次数t为0,重新执行步骤(4)~(7),直至提取出有效的光谱基线。
本实施例中,步骤(4)~(7)迭代23次后,提取出光谱基线拟合信号z。如图3所示,z为背景基线,d为去背景后的光谱信号。在主频为100MHz、嵌入式芯片DSP TMS320C6701上进行背景提取时,迭代一次所需时间为13.58ms;计算出本实施例的背景基线,所需时间为0.32s。
本发明说明书中未详细描述的内容为本领域技术人员公知技术。
Claims (9)
1.一种基于嵌入式***的快速光谱基线计算方法,其特征在于包括如下步骤:
(1)、构建光谱基线拟合信号z表达式为z=(Q+T)-1Qy,并建立光谱基线拟合方程z=A- 1b,其中,A=Q+T为N×N的上下带宽为n的2n+1对角矩阵,b=Qy为N×1的列向量,为N×N的对角矩阵,W为基于wt的N×N对角矩阵,T=DTD为N×N的上下带宽为n的2n+1对角矩阵,D为n阶微分矩阵;λ为调节系数,wt为权重向量;
(2)、对光谱信号进行等间隔采样,得到长度为N的光谱数字采样信号y;
(3)、初始化调节系数λ,令迭代次数t为0,权重向量wt初始化为全1的N×1列向量;
(4)、求解光谱基线拟合方程z=A-1b,得到光谱基线拟合信号z,若出现光谱基线拟合方程无法求解的情况时,则回到步骤(2),重新获取光谱数字采样信号y;否则,进入步骤(5);
(5)、假设d=y-z,取向量d中小于0的元素组成向量d-,对向量d-求平均m和标准差σ;
(6)、采用如下公式得到更新后的权重向量wt+1,公式如下:
(7)、若权重向量wt+1满足如下权重向量判别条件,则完成光谱基线的提取;否则,将迭代次数t加1更新迭代次数,同时更新后的权重向量wt+1替代之前的权重向量wt,跳至步骤(4)重新执行;
权重向量判别条件如下:
θ为权重比例系数;
(8)、若步骤(4)~(7)迭代次数t达到预设的门限,仍无法提取出有效的光谱基线,降低调节系数λ值,令迭代次数t为0,重新执行步骤(4)~(7),直至提取出有效的光谱基线。
2.根据权利要求1所述的一种基于嵌入式***的快速光谱基线计算方法,其特征在于所述步骤(4)求解光谱基线拟合方程的具体实现为:
(4.1)、将矩阵A的2n+1对角数据按如下规则转换为矩阵C的元素,矩阵C和A中的元素对应关系为ci-j+n,j=aij,i=0~N-1,j=0~N-1:
(4.2)、将矩阵A做压缩空间下的三角分解A=LU,L为下半带宽为n的单位下三角矩阵,U为上半带宽为n的上三角矩阵;
(4.3)、将b=Qy代入方程Lx=b,求解该方程得到向量x;
(4.4)、将步骤(4.3)得到的向量x代入方程Uz=x,求解方程得到拟合信号向量z。
4.根据权利要求2所述的一种基于嵌入式***的快速光谱基线计算方法,其特征在于所述步骤(4.2)的具体实现为:
(4.2.1)、令k=1;
(4.2.2)、采用如下公式计算矩阵U的第k-1行数据,并将矩阵U中的数据放入矩阵C(2n+1)×N中的前n+1行,公式如下:
每计算完成一个c'k-jj+n,jj-1值,将c'k-jj+n,jj-1值赋值给ck-jj+n,jj-1;
(4.2.3)、采用如下公式计算矩阵L的第k-1列数据,并将矩阵L中的数据放入矩阵C(2n+1)×N中的最后n行,公式如下:
每计算完成一个c'ii-k+n,k-1值,将c'ii-k+n,k-1值赋值给cii-k+n,k-1;当cn,k-1为0时,说明光谱基线拟合方程无法求解,直接结束矩阵分解,不再执行后续步骤(4.2.4);
(4.2.4)、k值加1,若k>N时,则完成矩阵分解;否则,跳至步骤(4.2.2)继续执行。
5.根据权利要求1所述的一种基于嵌入式***的快速光谱基线计算方法,其特征在于所述n取值为1~N-1。
7.根据权利要求2所述一种基于嵌入式***的快速光谱基线计算方法,其特征在于所述步骤(4.4)中所述拟合信号向量z保存在向量b的存储空间内,求解拟合信号向量z的具体实现为:
(4.4.1)、计算拟合信号向量z的第N-1个采样值zN-1,具体计算公式为:zN-1=bN-1/cn,N-1;
(4.4.2)、将拟合信号向量z的第N-1个采样值zN-1存储至向量b的第N-1个元素的存储空间中,即:bN-1=zN-1;
(4.4.3)、令i”=N-1;
(4.4.4)、采用下列公式计算拟合信号向量第i”-1个元素zi”-1的计算公式为:
(4.4.5)、将拟合信号向量第i”-1个元素zi”-1存储至向量b的第i”-1个元素对应的存储空间bi”-1;
(4.4.6)、i”值减1;
(4.4.7)、若i”<1,则完成向量z的计算;否则,跳至步骤(4.4.4)继续执行。
8.根据权利要求1所述的一种基于嵌入式***的快速光谱基线计算方法,其特征在于:所述的步骤(1)中调节系数λ的取值范围为104~105。
9.根据权利要求1所述的一种基于嵌入式***的快速光谱基线计算方法,其特征在于:所述的步骤(7)中权重比例系数θ的取值范围为10-5~10-6。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911206716.6A CN111125625B (zh) | 2019-11-29 | 2019-11-29 | 一种基于嵌入式***的光谱基线计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911206716.6A CN111125625B (zh) | 2019-11-29 | 2019-11-29 | 一种基于嵌入式***的光谱基线计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111125625A true CN111125625A (zh) | 2020-05-08 |
CN111125625B CN111125625B (zh) | 2023-06-06 |
Family
ID=70496218
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911206716.6A Active CN111125625B (zh) | 2019-11-29 | 2019-11-29 | 一种基于嵌入式***的光谱基线计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111125625B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113673409A (zh) * | 2021-08-16 | 2021-11-19 | 点靓纳谱(上海)生物医药科技有限公司 | 一种用于光谱分析的自动收敛光谱校正方法及装置 |
CN117633423A (zh) * | 2024-01-26 | 2024-03-01 | 苏州简测科技有限公司 | 一种自适应光谱去基线算法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090001262A1 (en) * | 2003-10-22 | 2009-01-01 | Erik Visser | System and Method for Spectral Analysis |
US20090285463A1 (en) * | 2008-04-18 | 2009-11-19 | Ricardo Otazo | Superresolution parallel magnetic resonance imaging |
US20100100336A1 (en) * | 2008-10-21 | 2010-04-22 | Wright David A | Methods of automated spectral peak detection and quantification without user input |
US20100108895A1 (en) * | 2008-11-05 | 2010-05-06 | Siemens Aktiengesellschaft | Background signal suppression in pet spectrums |
CN101763445A (zh) * | 2008-12-23 | 2010-06-30 | 北京理工大学 | 一种高光谱图像降维芯片 |
CN102542547A (zh) * | 2011-12-29 | 2012-07-04 | 北京航空航天大学 | 一种基于光谱约束的高光谱图像融合方法 |
CN103018194A (zh) * | 2012-12-06 | 2013-04-03 | 江苏省质量安全工程研究院 | 基于背景估计的非对称最小二乘基线校正方法 |
JP2013101608A (ja) * | 2011-10-17 | 2013-05-23 | Mizuho Information & Research Institute Inc | 非線形構造解析計算装置、非線形構造解析計算方法及び非線形構造解析計算プログラム |
CN103217409A (zh) * | 2013-03-22 | 2013-07-24 | 重庆绿色智能技术研究院 | 一种拉曼光谱预处理方法 |
CN106529563A (zh) * | 2016-09-19 | 2017-03-22 | 西安电子科技大学 | 基于双图稀疏非负矩阵分解的高光谱波段选择方法 |
CN106596506A (zh) * | 2016-12-16 | 2017-04-26 | 温州大学 | 一种基于压缩存储和列选主元高斯消去法的airPLS实现方法 |
CN108844939A (zh) * | 2018-03-14 | 2018-11-20 | 西安电子科技大学 | 基于非对称加权最小二乘的拉曼光谱检测基线校正方法 |
CN109709062A (zh) * | 2018-12-29 | 2019-05-03 | 广东工业大学 | 一种物质识别方法、装置和计算机可读存储介质 |
-
2019
- 2019-11-29 CN CN201911206716.6A patent/CN111125625B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090001262A1 (en) * | 2003-10-22 | 2009-01-01 | Erik Visser | System and Method for Spectral Analysis |
US20090285463A1 (en) * | 2008-04-18 | 2009-11-19 | Ricardo Otazo | Superresolution parallel magnetic resonance imaging |
US20100100336A1 (en) * | 2008-10-21 | 2010-04-22 | Wright David A | Methods of automated spectral peak detection and quantification without user input |
US20100108895A1 (en) * | 2008-11-05 | 2010-05-06 | Siemens Aktiengesellschaft | Background signal suppression in pet spectrums |
CN101763445A (zh) * | 2008-12-23 | 2010-06-30 | 北京理工大学 | 一种高光谱图像降维芯片 |
JP2013101608A (ja) * | 2011-10-17 | 2013-05-23 | Mizuho Information & Research Institute Inc | 非線形構造解析計算装置、非線形構造解析計算方法及び非線形構造解析計算プログラム |
CN102542547A (zh) * | 2011-12-29 | 2012-07-04 | 北京航空航天大学 | 一种基于光谱约束的高光谱图像融合方法 |
CN103018194A (zh) * | 2012-12-06 | 2013-04-03 | 江苏省质量安全工程研究院 | 基于背景估计的非对称最小二乘基线校正方法 |
CN103217409A (zh) * | 2013-03-22 | 2013-07-24 | 重庆绿色智能技术研究院 | 一种拉曼光谱预处理方法 |
CN106529563A (zh) * | 2016-09-19 | 2017-03-22 | 西安电子科技大学 | 基于双图稀疏非负矩阵分解的高光谱波段选择方法 |
CN106596506A (zh) * | 2016-12-16 | 2017-04-26 | 温州大学 | 一种基于压缩存储和列选主元高斯消去法的airPLS实现方法 |
CN108844939A (zh) * | 2018-03-14 | 2018-11-20 | 西安电子科技大学 | 基于非对称加权最小二乘的拉曼光谱检测基线校正方法 |
CN109709062A (zh) * | 2018-12-29 | 2019-05-03 | 广东工业大学 | 一种物质识别方法、装置和计算机可读存储介质 |
Non-Patent Citations (1)
Title |
---|
姜安;彭江涛;谢启伟;魏纪平;彭思龙;: "一种改进的非对称最小二乘基线校正算法" * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113673409A (zh) * | 2021-08-16 | 2021-11-19 | 点靓纳谱(上海)生物医药科技有限公司 | 一种用于光谱分析的自动收敛光谱校正方法及装置 |
CN117633423A (zh) * | 2024-01-26 | 2024-03-01 | 苏州简测科技有限公司 | 一种自适应光谱去基线算法 |
CN117633423B (zh) * | 2024-01-26 | 2024-04-05 | 苏州简测科技有限公司 | 一种自适应光谱去基线方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111125625B (zh) | 2023-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109543763B (zh) | 一种基于卷积神经网络的拉曼光谱分析方法 | |
CN110045321B (zh) | 基于稀疏和低秩恢复的稳健doa估计方法 | |
CA2537328C (en) | Method of processing and storing mass spectrometry data | |
CN111125625A (zh) | 一种基于嵌入式***的光谱基线计算方法 | |
CN106649026B (zh) | 适用于运维自动化***的监测数据压缩方法 | |
CN109946253B (zh) | 一种光谱去噪方法 | |
CN110865357A (zh) | 一种基于参数优化vmd的激光雷达回波信号降噪方法 | |
CN110503156B (zh) | 一种基于最小相关系数的多变量校正特征波长选择方法 | |
CN114611582B (zh) | 一种基于近红外光谱技术分析物质浓度的方法及*** | |
CN112766484A (zh) | 浮点神经网络模型量化***和方法 | |
CN117236084A (zh) | 一种木工机械加工生产智能管理方法及*** | |
US7656978B2 (en) | Fast Q-filter | |
CN108469609B (zh) | 一种用于雷达目标跟踪的检测信息滤波方法 | |
US7865322B2 (en) | Relative noise | |
CN115902391A (zh) | 基于形态学滤波和特征值分解的动态门限信号检测方法 | |
CN114627963B (zh) | 蛋白数据填充方法、***、计算机设备及可读存储介质 | |
CN115561282A (zh) | 一种地下水重金属的检测方法及*** | |
CN112801936B (zh) | 一种x射线荧光光谱自适应本底扣除方法 | |
CN115541021A (zh) | 拉曼光谱特征峰的定位方法、电子设备及存储介质 | |
CN114965425A (zh) | 基于迭代自适应加权算法的绝缘油拉曼光谱基线扣除方法和*** | |
CN109062861B (zh) | 一种基于滑动递推限幅滤波的数据处理方法 | |
CN114067169A (zh) | 一种基于卷积神经网络的拉曼光谱分析方法 | |
CN115136140A (zh) | 采集阶段峰宽确定和基线偏移估计的技术 | |
CN115102643B (zh) | 噪声基底跟踪方法、装置、设备及存储介质 | |
CN111949920A (zh) | 一种提取海冰观测值分解的空间多尺度信息的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |