CN108918758A - 化学反应光谱的cielab色空间色差△e与物质量计算方法 - Google Patents

化学反应光谱的cielab色空间色差△e与物质量计算方法 Download PDF

Info

Publication number
CN108918758A
CN108918758A CN201810840688.2A CN201810840688A CN108918758A CN 108918758 A CN108918758 A CN 108918758A CN 201810840688 A CN201810840688 A CN 201810840688A CN 108918758 A CN108918758 A CN 108918758A
Authority
CN
China
Prior art keywords
value
uproar
parameter
chemical reaction
spread out
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.)
Pending
Application number
CN201810840688.2A
Other languages
English (en)
Inventor
王飞
葛娜
何浩
王立杉
邹明强
崔宗岩
王涛
王宇曦
张昂
张妍妍
宁建华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201810840688.2A priority Critical patent/CN108918758A/zh
Publication of CN108918758A publication Critical patent/CN108918758A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N31/00—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods
    • G01N31/16—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods using titration
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/12—Circuits of general importance; Signal processing
    • G01N2201/129—Using chemometrical methods

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Molecular Biology (AREA)
  • Investigating Or Analysing Materials By The Use Of Chemical Reactions (AREA)

Abstract

本发明公开了化学反应光谱的CIELAB色空间色差△E与物质量计算方法，对化学反应容器内的化学反应溶液在可见光波长范围内依次测量加入的不同反应物体积V值对应的吸光度值，并计算为CIELAB色空间的色差△E同步测量加入的不同反应物体积V值对应的pH值、氢离子浓度c[H+]值，将色差△E值、pH值、c[H+]值经过数据降噪后计算出对应的△E噪值、pH噪值、c[H+]噪值，将该△E噪值参数与加入的不同反应物体积V值计算出△E衍值，并与对应的pH噪值、c[H+]噪值建立平面直角坐标曲线，实现了化学反应同步测量光光谱谱曲线测量的计算方法，可实现化学反应测量的自动化、批量化的检测规模检测。

Description

化学反应光谱的CIELAB色空间色差△E与物质量计算方法
技术领域
本发明涉及测量技术领域，尤其涉及化学反应光谱的CIELAB色空间色差△E与物质量计算方法。
背景技术
化学反应离不开溶液中离子变化，不同化学反应的成分对可见光光谱的吸收与反射不同，这种变化与溶液中化学成分的化学结构有关。这种专利技术包含了可见光范围内的连续光谱、溶液离子、试剂加入试剂的体积的精确测量与计算、试剂加入试剂体积对应的pH值、c[H+]值，经过数学方程去除噪音的干扰，建立直观的平面直角坐标的化学反应光谱变化曲线。用实时光谱曲线参数的变化，为了解、观察化学反应的进程带来了一种新颖且直观的检测手段。
目前，化学界常用的方法是通过单次测量结果估算化学反应进程，测量精度差、劳动强度大、方法简陋且繁琐，不能用反应曲线计量反应突变点。国内外尚无用于实时色度值参数的建立平面直角坐标化学反应光谱变化曲线对化学反应溶液进程检测、分析的方法，该方法将对光谱产生影响的化学反应用直观、即时的曲线在平面坐标系中表示出来，采用衍生参数△E衍作为纵坐标与pH值、c[H+]、加入试剂的体积V为横坐标，建立色空间△E衍生参数与pH值或c[H+]或者加入试剂的体积V的化学反应光谱反应曲线，该曲线可以实现实时的测量和反应化学反应变化，进行化学形态的表征。
发明内容
本发明的目的是为了解决现有技术中存在的缺点，而提出的化学反应光谱的CIELAB色空间色差△E与物质量计算方法。
为了实现上述目的，本发明采用了如下技术方案：
化学反应光谱的CIELAB色空间△E*值与物质量计算方法，计算步骤包括：对化学反应溶液在可见光波长范围内测量加入的不同反应物体积V值对应的一组波长的吸光度值，计算色差△E(CIELAB色空间的色差△E* ab，按照公式或者计算,CIE DE 2000总色差△E00按照公式△E* ab和△* 00统称为本发明的色差△E)，将该组△E值数据经过数学方程降噪后，计算出该点的△E噪参数，将△E噪与加入的不同反应物体积V值计算出△E噪的衍生参数△E衍，同步测量加入的不同反应物体积V值对应的pH值，将pH值计算为c[H+]，将pH值和c[H+]数据经过数学方程降噪后为pH噪值和c[H+]噪，再将衍生参数△E衍按加入的不同反应物体积V值顺序与对应的pH值噪、c[H+]噪建立平面直角坐标系，该平面直角坐标系中的曲线即为化学反应CIELAB色空间△E衍参数与加入反应物体积V同步测量的体积V值、pH值、离子浓度c[H+]的物质特征量的坐标曲线；
S01：根据加入的不同反应物体积V值的第x点测量的吸光度值，计算该测量点的参数的△Ex值，其中x＝1,2,3，…；
S02：根据加入的不同反应物体积V值的第y点测量的吸光度值，计算该测量点的的△Ey值，其中y＝2,3，4，…；
S03：根据第x点测量的pHx值，其中x＝1,2,3，…，计算该测量点的c[H+](mo/)x值，其中x＝1,2,3，…；
S04：根据第y点测量的pHy值，其中y＝2,3，4，…，计算该测量点的c[H+]y值，其中y＝2,3，4，…；
S05：依次将获得的△Ex值、△Ey值、pHx值、pHy值、c[H+]x值、c[H+]y值通过数学方程降噪计算为降噪后的△Ex噪值、△Ey噪值、pHx噪值、pHy噪、c[H+]x噪值、c[H+]y噪值。
优选的，根据下式计算△E参数中衍生参数△E衍的△En/Vm，△En/Vm＝│△E噪x n/Vm│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；
优选的，根据下式计算△E参数中衍生参数△E衍的△En/△Vm，△En/△Vm＝│△E噪x n/(Vx-Vy)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
优选的，根据下式计算△E参数中衍生参数△E衍的△En/pHm，△En/pHm＝│△E噪x n/pHm│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；
优选的，根据下式计算△E参数中衍生参数△E衍的△En/△pHm，△En/△pHm＝│△E噪x n/(pHx-pHy)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
优选的，根据下式计算△E参数中衍生参数△E衍的△En/c[H+] m，△En/c[H+] m＝│△E噪x n/c[H+]x噪 m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；
优选的，根据下式计算△E参数中衍生参数△E衍的△En/△c[H+] m，△△En/△c[H+] m＝│△E噪x n/(c[H+]x噪-c[H+]x噪值)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，
优选的，根据下式计算△E参数中衍生参数△E衍的△△En/Vm，△△En/Vm＝│(△E噪x-△E噪y)n/Vm│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
优选的，根据下式计算△E参数中衍生参数△E衍的△△En/△Vm，△△En/△Vm＝│(△E噪x-△E噪y)n/(Vx-Vy)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
优选的，根据下式计算△E参数中衍生参数△E衍的△△En/pHm，△△En/pHm＝│(△E噪x-△E噪y)n/pHm│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
优选的，根据下式计算△E参数中衍生参数△E衍的△△En/△pHm，△△En/△pHm＝│(△E噪x-△E噪y)n/(pHx-pHy)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
优选的，根据下式计算△E参数中衍生参数△E衍的△△En/c[H+] m，△△En/c[H+] m＝│(△E噪x-△E噪y)n/c[H+]x噪 m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
优选的，根据下式计算△E参数中衍生参数△E衍的△△En/△c[H+] m，△△En/△c[H+] m＝│(△E噪x-△E噪y)n/(c[H+]x噪-c[H+]x噪值)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，。
本发明的有益效果是：本发明实现了化学反应光谱曲线测量的计算方法，实现了化学反应光谱曲线的建立，其测量精度高，计算步骤简单，可实现化学反应测量的自动化、批量化的检测规模检测。
附图说明
图1为△E*3/V3-V的化学反应光谱曲线；
图2为△E*3/V3-pH的化学反应光谱曲线；
图3为△E*3/V3-c[H+]的化学反应光谱曲线；
图4为△E*3/△V3-V的化学反应光谱曲线；
图5为△E*3/△V3值-pH值的化学反应光谱曲线；
图6为△E*3/△V3值--c[H+]值的化学反应光谱曲线；
图7为△E*3/pH3-V的化学反应光谱曲线；
图8为△E*3/△pH3-pH的化学反应光谱曲线；
图9为△E*3/pH3--c[H+]的化学反应光谱曲线；
图10为△E*3/△pH3-V的化学反应光谱曲线；
图11为△E*3/△pH3值-pH值的化学反应光谱曲线；
图12为△E*3/△pH3值-c[H+]值的化学反应光谱曲线；
图13为△E*3/c[H+] 3-V的化学反应光谱曲线；
图14为△E*3/c[H+] 3-pH的化学反应光谱曲线；
图15为△E*3/c[H+] 3-c[H+]的化学反应光谱曲线；
图16为△E*3/△c[H+] 3-V的化学反应光谱曲线；
图17为△E*3/△c[H+] 3-pH的化学反应光谱曲线；
图18为△E*3/△c[H+] 3-c[H+]的化学反应光谱曲线；
图19为△△E*3/V3-V的化学反应光谱曲线；
图20为△△E*3/V3-pH的化学反应光谱曲线；
图21为△△E*3/V3-c[H+]的化学反应光谱曲线；
图22为△△E*3/△V3-V的化学反应光谱曲线；
图23为△△E*3/△V3值-pH值的化学反应光谱曲线；
图24为△△E*3/△V3值--c[H+]值的化学反应光谱曲线；
图25为△△E*3/pH3-V的化学反应光谱曲线；
图26为△△E*3/△pH3-pH的化学反应光谱曲线；
图27为△△E*3/pH3--c[H+]的化学反应光谱曲线；
图28为△△E*3/△pH3-V的化学反应光谱曲线；
图29为△△E*3/△pH3值-pH值的化学反应光谱曲线；
图30为△△E*3/△pH3值-c[H+]值的化学反应光谱曲线；
图31为△△E*3/c[H+] 3-V的化学反应光谱曲线；
图32为△△E*3/c[H+] 3-pH的化学反应光谱曲线；
图33为△△E*3/c[H+] 3-c[H+]的化学反应光谱曲线；
图34为△△E*3/△c[H+] 3-V的化学反应光谱曲线；
图35为△△E*3/△c[H+] 3-pH的化学反应光谱曲线；
图36为△△E*3/△c[H+] 3-c[H+]的化学反应光谱曲线。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
参照图1-图36，化学反应光谱的CIELAB色空间△E值与物质量计算方法，其计算步骤包括：
S01：对化学反应容器内的化学反应溶液在可见光波长范围内测量加入的不同反应物体积V值顺序n、n+1(n＝∞)、相同波长间隔△λ的一组波长的吸光度值Abs；
S02：将该组波长的吸光度值Abs计算为色空间的△E值；
S03：将该组波长的吸光度值计算为色空间的△E值；
S04：将该△E值组数据均经过数学方程降噪后，计算出该点的△E噪参数；
S05：依次测量加入的不同反应物体积V值n、n+1(n＝∞)测量点的加入试剂的体积V值；
S06：依次测量加入的不同反应物体积V值n、n+1(n＝∞)测量点的pH值；
S07：将该pH值组数据均经过数学方程降噪后，计算出该点的pH噪参数；
S08：将pH噪参数换算为c[H+]值的摩尔浓度moL/L；
S09：依测量加入的不同反应物体积V值的顺序将△E噪n、△E噪n+1(n＝∞)的△E噪参数与加入试剂的体积V值计算化学光谱衍生参数△E衍；
S10：将按测量加入的不同反应物体积V值的顺序n、n+1、n+x(x＝∞)的计算的△E噪的衍生参数△E衍与对应的pH值，建立平面平面直角坐标系，该平面直角坐标系中的曲线即为横坐标为pH值的化学反应光谱曲线；
S11：将按测量加入的不同反应物体积V值的顺序n、n+1、n+x(x＝∞)的计算的△E噪的衍生参数与对应的体积V值，建立平面直角坐标系，该平面直角坐标系中的曲线即为横坐标为加入试剂体积V的化学反应光谱曲线。
S12：将按测量加入的不同反应物体积V值的顺序n、n+1、n+x(x＝∞)的计算的pH噪的衍生参数与对应的体积V值，建立平面直角坐标系，该平面直角坐标系中的曲线即为横坐标为加入试剂体积V的化学反应光谱曲线。
S13：将按测量加入的不同反应物体积V值的顺序n、n+1、n+x(x＝∞)的计算的c[H+]噪的衍生参数与对应的体积V值，建立平面直角坐标系，该平面直角坐标系中的曲线即为横坐标为加入试剂体积V的化学反应光谱曲线。
本实施例中，测量时，用水作为空白参比溶液进行吸光度测量装置的空白校正，而后将加入化学反应容器中置于吸光度测量装置中，加入化学反应溶液，而后以波长间隔△λ为5nm、测量波长380nm～760nm的一组波长。
利用计算机进行计算，所述计算过程利用程序语言编为计算机程序，并作为应该程序安装于计算机中。所述计算过程包括:
SO1：用水作为空白参比溶液进行吸光度测量装置的空白校正，设置波长间隔△λ为5nm、测量范围为波长380nm～760nm；
SO2：称取110g氢氧化钠，溶于100mL无二氧化碳的水中，摇匀，注入聚乙烯容器中，密闭放置至溶液清亮。按表1的规定量，用塑料管量取上层清液5.4mL于试剂瓶中，用无二氧化碳的水稀释至1000mL，摇匀；
SO3：称取0.7523g于105℃～110℃电烘箱中干燥至恒量的工作基准试剂邻苯二甲酸氢钾，加50ml无二氧化碳的水溶解至化学反应容器中，加入浓度为10g/L的酚酞指示液2滴；
SO4：将化学反应容器送入测量仪器，测量加入的氢氧化钠溶液不同体积的对应的吸光度值Abs；
SO5：将该组测量的加入的氢氧化钠溶液不同体积V值顺序的波长的吸光度值Abs计算为对应顺序的CIELAB色空间的△E值；
SO6：将该组测量的不同体积V值顺序的波长的吸光度值Abs计算为对应顺序的色空间的△E值；
S07：将该不同体积V值顺序的△E值组数据经过数学方程降噪后，计算出对应顺序的△E噪参数；
S08：将对应顺序的△E噪参数计算为对应顺序的化学光谱衍生参数△E衍；
S09：将不同体积V值顺序对应的该组pH值数据均经过数学方程降噪后，计算对应的pH噪参数；
S10：将pH噪参数换算为对应的c[H+]值的摩尔浓度moL/L；
S11：将不同体积V值n、n+1、n+x(x＝∞)顺序计算的pH噪的衍生参数与对应的体积V值，建立平面直角坐标系，该平面直角坐标系中的曲线即为横坐标为加入试剂体积V值、纵坐标为化学光谱衍生参数△E衍的化学反应光谱CIELAB色空间曲线。
S12：将不同体积V值n、n+1、n+x(x＝∞)顺序计算的pH噪的衍生参数与对应的pH噪参数值，建立平面直角坐标系，该平面直角坐标系中的曲线即为横坐标为加入试剂体积V对应的pH噪参数值、纵坐标为化学光谱衍生参数△E衍的化学反应光谱CIELAB色空间曲线。
S13：将不同体积V值n、n+1、n+x(x＝∞)顺序计算的pH噪的衍生参数与对应的c[H+]值，建立平面直角坐标，该平面直角坐标中的曲线即为横坐标为加入试剂体积V对应的c[H+]值、纵坐标为化学光谱衍生参数△E衍的化学反应光谱CIELAB色空间曲线。
S14：在平面直角坐标系上色空间曲线上的凸出峰，即为氢氧化钠标准溶液标定中的滴定终点，该实验滴定的最大凸出峰对应的体积10.016ml即为滴定终点体积V；
S15：与实验样品同时做空白试验，该空白试验终点的最大凸出峰对应的体积0.100mL即为空白试验滴定终点体积V0；
S16：将空白试验的滴定终点体积V0(0.100mL)、实验样品滴定终点体积V(10.016mL)、邻苯二甲酸氢钾质量m(0.7525g)、邻苯二甲酸氢钾的摩尔质量(单位为克每摩尔(g/mol)[M(KHC8H4O4)＝204.22])带入计算公式cNaOH(mol/L)计算氢氧化钠的浓度cNaOH：cNaOH(mol/L)＝(m×1000)/((V-V0)×M),式中：m——邻苯二甲酸氢钾质量，单位为克(g)；V——氢氧化钠溶液体积，单位为毫升(mL)；V0——空白试验消耗氢氧化钠溶液体积，单位为毫升(mL)；M——邻苯二甲酸氢钾的摩尔质量，单位为克每摩尔(g/mol)[M(KHC8H4O4)＝204.22]，cNaOH＝0.3715mol/L。
加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的物质特征量的体积V值、pH值、c[H+]值、△V值、pH△值、△c[H+]值建立化学反应CIELAB色空间的△E衍参数与反应物体积V或反应物体积V对应的pH值、c[H+]值的平面坐标系曲线；
S01：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值，建立△En/Vm-V平面坐标系曲线，纵坐标参数为△En/Vm＝│△E噪x n/Vx m│(其中n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S02：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的增量△V，建立△En/△Vm-V平面坐标系曲线,纵坐标参数为△En/△Vm＝│△E噪x n/(Vx-Vy)m│(其中n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S03：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的pH值，建立△En/pHm-V平面坐标系曲线,纵坐标参数为△En/pHm＝│△E噪x n/pHx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S04：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的增量△pH值，建立△En/△pHm-V平面坐标系曲线,纵坐标参数为△En/△pHm＝│△E噪x n/(pHx-pHy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S05：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的c[H+]值，建立△En/c[H+] m-V平面坐标系曲线,纵坐标参数为△En/c[H+] m＝│△E噪x n/c[H+]x m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S06：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的增量△c[H+]值，建立△En/△c[H+] m-V平面坐标系曲线,纵坐标参数为△En/△c[H+] m＝│△E噪x n/(c[H+]x-c[H+]y)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S07：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值，建立△En/Vm-pH平面坐标系曲线,纵坐标参数为△En/Vm＝│△E噪x n/Vx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；)，横坐标为加入的不同反应物体积V值对应的pH；
S08：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的增量△V，建立△En/△Vm-pH平面坐标系曲线,纵坐标参数为△En/△Vm＝│△E噪x n/(Vx-Vy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的pH；
S09：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的pH值，建立△En/pHm-pH平面坐标系曲线,纵坐标参数为△En/pHm＝│△E噪x n/pHx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；)，横坐标为加入的不同反应物体积V值对应的pH；
S10：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的pH值的增量△pH值，建立△En/△pHm-pH平面坐标系曲线,纵坐标参数为△En/△pHm＝│△E噪x n/(pHx-pHy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的pH；
S11：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的c[H+]值，建立△En/c[H+] m-pH平面坐标系曲线,纵坐标参数为△En/c[H+] m＝│△E噪x n/c[H+]x m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；)，横坐标为加入的不同反应物体积V值对应的pH；
S12：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的c[H+]值的增量△c[H+]值，建立△En/△c[H+] m-pH平面坐标系曲线,纵坐标参数为△En/△c[H+] m＝│△E噪x n/(c[H+]x-c[H+]y)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的pH；
S13：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值，建立△En/Vm-c[H+]平面坐标系曲线,纵坐标参数为△En/Vm＝│△E噪x n/Vx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S14：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的增量△V，建立△En/△Vm-c[H+]平面坐标系曲线,纵坐标参数为△En/△Vm＝│△E噪x n/(Vx-Vy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S15：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的pH值，建立△En/pHm-c[H+]平面坐标系曲线,纵坐标参数为△En/pHm＝│△E噪x n/pHx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S16：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的增量△pH值，建立△En/△pHm-c[H+]平面坐标系曲线,纵坐标参数为△En/△pHm＝│△E噪x n/(pHx-pHy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S17：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的c[H+]值，建立△En/c[H+] m-c[H+]平面坐标系曲线,纵坐标参数为△En/Vm＝│E* 噪x n/c[H+]x m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S18：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△E衍与加入的不同反应物体积V值对应的增量△c[H+]值，建立△En/△c[H+] m-c[H+]平面坐标系曲线,纵坐标参数为△En/△Vm＝│△E噪x n/(c[H+]x-c[H+]y)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S19：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值，建立△△En/Vm-V平面坐标系曲线,纵坐标参数为△△En/Vm＝│(△E噪x-△E噪y)n/Vx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S20：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的增量△V，建立△△En/△Vm-V平面坐标系曲线,纵坐标参数为△△En/△Vm＝│(△E噪x-△E噪y)n/(Vx-Vy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S21：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的pH值，建立△△En/pHm-V平面坐标系曲线,纵坐标参数为△△En/pHm＝│(△E噪x-△E噪y)n/pHx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S22：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的增量△pH值，建立△△En/△pHm-V平面坐标系曲线,纵坐标参数为△△En/△pHm＝│(△E噪x-△E噪y)n/(pHx-pHy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S23：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的c[H+]值，建立△△En/c[H+] m-V平面坐标系曲线,纵坐标参数为△△En/Vm＝│(△E噪x-△E噪y)n/c[H+]x m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S24：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的增量△c[H+]值，建立△△En/△c[H+] m-V平面坐标系曲线,纵坐标参数为△△En/△c[H+] m＝│(△E噪x-△E噪y)n/(c[H+]x-c[H+]y)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值或其代表的质量数值；
S25：根据加入的不同反应物体积V值的顺序对应的E*参数中衍生参数△E衍与加入的不同反应物体积V值，建立△△En/Vm-pH平面坐标系曲线,纵坐标参数为△△En/Vm＝│(△E噪x-△E噪y)n/Vx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的pH；
S26：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的增量△V，建立△△En/△Vm-pH平面坐标系曲线,纵坐标参数为△△En/△Vm＝│(△E噪x-△E噪y)n/(Vx-Vy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的pH；
S27：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的pH值，建立△△En/pHm-pH平面坐标系曲线,纵坐标参数为△△En/pHm＝│(△E噪x-△E噪y)n/pHx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的pH；
S28：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的pH值的增量△pH值，建立△△En/△pHm-pH平面坐标系曲线,纵坐标参数为△△En/△pHm＝│(△E噪x-△E噪y)n/(pHx-pHy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的pH；
S29：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的c[H+]值，建立△△En/c[H+] m-pH平面坐标系曲线,纵坐标参数为△△En/Vm＝│(△E噪x-△E噪y)n/c[H+]x m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的pH；
S30：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的c[H+]值的增量△c[H+]值，建立△△En/△c[H+] m-pH平面坐标系曲线,纵坐标参数为△△En/△Vm＝│(△E噪x-△E噪y)n/(c[H+]x-c[H+]y)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的pH；
S31：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值，建立△△En/Vm-c[H+]平面坐标系曲线,纵坐标参数为△△En/Vm＝│(△E噪x-△E噪y)n/Vx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S32：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的增量△V，建立△△En/△Vm-c[H+]平面坐标系曲线,纵坐标参数为△△En/△Vm＝│(△E噪x-△E噪y)n/(Vx-Vy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S33：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的pH值，建立△△En/pHm-c[H+]平面坐标系曲线,纵坐标参数为△△En/pHm＝│(△E噪x-△E噪y)n/pHx m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S34：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的增量△pH值，建立△△En/△pHm-c[H+]平面坐标系曲线,纵坐标参数为△△En/△pHm＝│(△E噪x-△E噪y)n/(pHx-pHy)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S35：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的c[H+]值，建立△△En/c[H+] m-c[H+]平面坐标系曲线,纵坐标参数为△△En/Vm＝│(△E噪x-△E噪y)n/c[H+]x m│(其中，n＝1,2,3，…；m＝0,1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的c[H+]；
S36：根据加入的不同反应物体积V值的顺序对应的△E参数中衍生参数△△E衍与加入的不同反应物体积V值对应的增量△c[H+]值，建立△△En/△c[H+] m-c[H+]平面坐标系曲线,纵坐标参数为△△En/△Vm＝│(△E噪x-△E噪y)n/(c[H+]x-c[H+]y)m│(其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，)，横坐标为加入的不同反应物体积V值对应的c[H+]；
△E参数中衍生参数△E衍的平面坐标系曲线上最大凸出峰值对应的横坐标体积参数代表的物质特征量的值为化学反应突变点；
S01：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V或反应物体积V代表的物质质量参数的△En/Vm-V平面坐标系曲线,△En/Vm-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S02：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V或反应物体积增量△V代表的物质质量参数的△En/△Vm-V平面坐标系曲线,△En/△Vm-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S03：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH参数的△En/pHm-V平面坐标系曲线,△En/pHm-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S04：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH增量的△pH参数的△En/△pHm-V平面坐标系曲线,△En/△pHm-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S05：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]参数的的△En/c[H+] m-V平面坐标系曲线,△En/c[H+] m-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S06：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]增量的△c[H+]参数的的△En/c[H+] m-△V平面坐标系曲线,△En/△c[H+] m-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S07：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V或反应物体积V代表的物质质量参数的△En/Vm-pH平面坐标系曲线,△En/Vm-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S08：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V或反应物体积增量△V代表的物质质量参数的△En/△Vm-pH平面坐标系曲线,△En/△Vm-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S09：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH参数的△En/pHm-pH平面坐标系曲线,△En/pHm-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S10：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH增量的△pH参数的△En/△pHm-pH平面坐标系曲线,△En/△pHm-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S11：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]参数的的△En/c[H+] m-pH平面坐标系曲线,△En/c[H+] m-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S12：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]增量的△c[H+]参数的的△En/△c[H+] m-pH平面坐标系曲线,△En/△c[H+] m-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S13：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V或反应物体积V代表的物质质量参数的△En/Vm-c[H+]平面坐标系曲线,△En/Vm-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S14：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积增量△V代表的物质质量参数的△En/△Vm-c[H+]平面坐标系曲线,△En/△Vm-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S15：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH参数的△En/pHm-c[H+]平面坐标系曲线,△En/pHm-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S16：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH增量的△pH参数的△En/△pHm-c[H+]平面坐标系曲线,△En/△pHm-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S17：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]参数的的△En/c[H+] m-c[H+]平面坐标系曲线,△En/c[H+] m-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S18：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]参数的增量△c[H+]的△En/△c[H+] m-c[H+]平面坐标系曲线,△En/△c[H+] m-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S19：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V或反应物体积V代表的物质质量参数的△△En/Vm-V平面坐标系曲线,△△En/Vm-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S20：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V或反应物体积增量△V代表的物质质量参数的△△En/△Vm-V平面坐标系曲线,△△En/△Vm-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S21：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH参数的△△En/pHm-V平面坐标系曲线,△△En/pHm-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S22：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH增量的△pH参数的△△En/△pHm-V平面坐标系曲线,△△En/△pHm-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S23：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]参数的的△△En/c[H+] m-V平面坐标系曲线,△△En/c[H+] m-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S24：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]增量的△c[H+]参数的的△△En/c[H+] m-△V平面坐标系曲线,△△En/△c[H+] m-V平面坐标系曲线上最大凸出峰值对应的横坐标体积V代表的物质特征量的值为化学反应突变点；
S25：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V或反应物体积V代表的物质质量参数的△△En/Vm-pH平面坐标系曲线,△△En/Vm-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S26：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V或反应物体积增量△V代表的物质质量参数的△△En/△Vm-pH平面坐标系曲线,△△En/△Vm-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S27：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH参数的△△En/pHm-pH平面坐标系曲线,△△En/pHm-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S28：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH增量的△pH参数的△△En/△pHm-pH平面坐标系曲线,△△En/△pHm-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S29：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]参数的的△△En/c[H+] m-pH平面坐标系曲线,△△En/c[H+] m-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S30：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]增量的△c[H+]参数的的△△En/△c[H+] m-pH平面坐标系曲线,△△En/△c[H+] m-pH平面坐标系曲线上最大凸出峰值对应的横坐标上pH值代表的物质特征量的值为化学反应突变点；
S31：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V或反应物体积V代表的物质质量参数的△△En/Vm-c[H+]平面坐标系曲线,△△En/Vm-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S32：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积增量△V代表的物质质量参数的△△En/△Vm-c[H+]平面坐标系曲线,△△En/△Vm-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S33：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH参数的△△En/pHm-c[H+]平面坐标系曲线,△△En/pHm-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S34：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的pH增量的△pH参数的△△En/△pHm-c[H+]平面坐标系曲线,△△En/△pHm-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S35：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]参数的的△△En/c[H+] m-c[H+]平面坐标系曲线,△△En/c[H+] m-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点；
S36：根据建立化学反应光谱CIELAB色空间△E衍参数与反应物体积V对应的c[H+]参数的增量△c[H+]的△△En/△c[H+] m-c[H+]平面坐标系曲线,△△En/△c[H+] m-c[H+]平面坐标系曲线上最大凸出峰值对应的横坐标上c[H+]代表的物质特征量的值为化学反应突变点。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在发明的保护范围之内。

Claims (14)

1.化学反应光谱的CIELAB色空间△E值与物质量计算方法，其特征在于，计算步骤包括：对化学反应溶液在可见光波长范围内测量加入的不同反应物体积V值对应的一组波长的吸光度值，计算色差△E(CIELAB色空间的色差△E* ab，按照公式或者计算,CIE DE 2000总色差△E00按照公式计算，△E* ab和△* 00统称为本发明的CIELAB色空间的色差△E)，将该组△E值的数据经过数学方程降噪后，计算出该点的△E噪参数，将△E噪与加入的不同反应物体积V值计算出△E噪的衍生参数△E衍，同步测量加入的不同反应物体积V值对应的pH值，将pH值计算为氢离子浓度c[H+]，将pH值和c[H+]数据经过数学方程降噪后为pH噪值和c[H+]噪，再将衍生参数△E衍按加入的不同反应物体积V值顺序与对应的pH值噪、c[H+]噪建立平面直角坐标系，该平面直角坐标系中的曲线即为化学反应CIELAB色空间△E衍参数与加入反应物体积V同步测量的体积V值、pH值、离子浓度c[H+]的物质特征量的坐标曲线；
S01：根据加入的不同反应物体积V值的第x点测量的吸光度值，计算该测量点的参数的△Ex值，其中x＝1,2,3，…；
S02：根据加入的不同反应物体积V值的第y点测量的吸光度值，计算该测量点的的△Ey值，其中y＝2,3，4，…；
S03：根据第x点测量的pHx值，其中x＝1,2,3，…，计算该测量点的c[H+](mo/)x值，其中x＝1,2,3，…；
S04：根据第y点测量的pHy值，其中y＝2,3，4，…，计算该测量点的c[H+]y值，其中y＝2,3，4，…；
S05：依次将获得的△Ex值、△Ey值、pHx值、pHy值、c[H+]x值、c[H+]y值通过数学方程降噪计算为降噪后的△Ex噪值、△Ey噪值、pHx噪值、pHy噪、c[H+]x噪值、c[H+]y噪值。
2.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△En/Vm，△En/Vm＝│△E噪x n/Vm│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；
3.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△En/△Vm，△En/△Vm＝│△E噪x n/(Vx-Vy)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
4.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△En/pHm，△En/pHm＝│△E噪x n/pHm│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；
5.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△En/△pHm，△En/△pHm＝│△E噪x n/(pHx-pHy)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
6.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△En/c[H+] m，△En/c[H+] m＝│△E噪x n/c[H+]x噪 m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；
7.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△En/△c[H+] m，△△En/△c[H+] m＝│△E噪x n/(c[H+]x噪-c[H+]x噪值)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，
8.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△△En/Vm，△△En/Vm＝│(△E噪x-△E噪y)n/Vm│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
9.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△△En/Vm，△△En/Vm＝│(△E噪x-△E噪y)n/Vm│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
10.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△△En/△Vm，△△En/△Vm＝│(△E噪x-△E噪y)n/(Vx-Vy)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
11.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△△En/pHm，△△En/pHm＝│(△E噪x-△E噪y)n/pHm│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
12.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△△En/△Vm，△△En/△pHm＝│(△E噪x-△E噪y)n/(pHx-pHy)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
13.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△△En/c[H+] m，△△En/c[H+] m＝│(△E噪x-△E噪y)n/c[H+]x噪 m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，；
14.根据权利要求1所述的化学反应光谱的CIELAB色空间色差△E与物质量计算方法，其特征在于，根据下式计算△E参数中衍生参数△E衍的△△En/△c[H+] m，△△En/△c[H+] m＝│(△E噪x-△E噪y)n/(c[H+]x噪-c[H+]x噪值)m│，其中，n＝1,2,3，…；m＝1,2,3，…；x＝1,2,3，…；y＝2,3，4，…，。
CN201810840688.2A 2018-07-27 2018-07-27 化学反应光谱的cielab色空间色差△e与物质量计算方法 Pending CN108918758A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810840688.2A CN108918758A (zh) 2018-07-27 2018-07-27 化学反应光谱的cielab色空间色差△e与物质量计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810840688.2A CN108918758A (zh) 2018-07-27 2018-07-27 化学反应光谱的cielab色空间色差△e与物质量计算方法

Publications (1)

Publication Number Publication Date
CN108918758A true CN108918758A (zh) 2018-11-30

Family

ID=64417164

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810840688.2A Pending CN108918758A (zh) 2018-07-27 2018-07-27 化学反应光谱的cielab色空间色差△e与物质量计算方法

Country Status (1)

Country Link
CN (1) CN108918758A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021114109A1 (zh) * 2019-12-10 2021-06-17 深圳市能源环保有限公司 一种基于cielab色彩空间快速测算污泥热值的方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0658882A (ja) * 1992-08-05 1994-03-04 Nakano Vinegar Co Ltd 自動滴定装置
CN102253164A (zh) * 2011-04-26 2011-11-23 东北电力大学 基于溶液图像技术的水质碱度在线测量装置及其测量方法
CN105717108A (zh) * 2016-02-07 2016-06-29 王飞 化学分析液体颜色CIE1976L*a*b*色空间测定方法
CN105738357A (zh) * 2016-02-07 2016-07-06 王飞 化学分析用氢氧化钠溶液配置的CIE1976L*a*b*色空间方法
CN107421948A (zh) * 2017-07-25 2017-12-01 秦皇岛水熊科技有限公司 基于CIELab色空间的检测溶液颜色变化的方法
CN108051434A (zh) * 2017-12-13 2018-05-18 深圳市在田翊方科技有限公司 一种基于颜色识别的待测液浓度定量检测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0658882A (ja) * 1992-08-05 1994-03-04 Nakano Vinegar Co Ltd 自動滴定装置
CN102253164A (zh) * 2011-04-26 2011-11-23 东北电力大学 基于溶液图像技术的水质碱度在线测量装置及其测量方法
CN105717108A (zh) * 2016-02-07 2016-06-29 王飞 化学分析液体颜色CIE1976L*a*b*色空间测定方法
CN105738357A (zh) * 2016-02-07 2016-07-06 王飞 化学分析用氢氧化钠溶液配置的CIE1976L*a*b*色空间方法
CN107421948A (zh) * 2017-07-25 2017-12-01 秦皇岛水熊科技有限公司 基于CIELab色空间的检测溶液颜色变化的方法
CN108051434A (zh) * 2017-12-13 2018-05-18 深圳市在田翊方科技有限公司 一种基于颜色识别的待测液浓度定量检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李运奎等: "基于CIELAB 色空间的红葡萄酒颜色直观表征", 《农业机械学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021114109A1 (zh) * 2019-12-10 2021-06-17 深圳市能源环保有限公司 一种基于cielab色彩空间快速测算污泥热值的方法

Similar Documents

Publication Publication Date Title
Zhu et al. A colorimetric method for fluoride determination in aqueous samples based on the hydroxyl deprotection reaction of a cyanine dye
López‐de‐Alba et al. Simultaneous determination of synthetic dyes tartrazine, allura red and sunset yellow by differential pulse polarography and partial least squares. A multivariate calibration method
Shamsipur et al. Novel fluorimetric bulk optode membrane based on a dansylamidopropyl pendant arm derivative of 1-aza-4, 10-dithia-7-oxacyclododecane ([12] aneNS2O) for selective subnanomolar detection of Hg (II) ions
CN105738357A (zh) 化学分析用氢氧化钠溶液配置的CIE1976L*a*b*色空间方法
CN106290334B (zh) 钴锰原料中钴和锰含量的化学测定方法
CN100456023C (zh) 快速检测食品中含吊白块方法
Goddu et al. Photometric titrations
CN103994996B (zh) 钒铁中磷的分析方法
CN103760160A (zh) 一种过氧化氢的比色检测方法
CN104792775A (zh) 一种经济快捷型比色测定水体中cod的测试包及测定方法
CN105954341A (zh) 一种石灰石中氧化钙、氧化镁的连续分析测定方法
CN105717108A (zh) 化学分析液体颜色CIE1976L*a*b*色空间测定方法
Barek et al. Polarographic and voltammetric determination of selected triazine-based azo dyes with different reactive groups
CN110018157A (zh) 一种纺织品甲醛快速测试方法
CN108918760A (zh) 化学反应光谱的CIELAB色空间的a*值与物质量关系计算方法
CN108918758A (zh) 化学反应光谱的cielab色空间色差△e与物质量计算方法
CN110320315A (zh) 三元正极材料的残碱含量检测方法及其应用
Safavi et al. Kinetic spectrophotometric determination of traces of sulfide
CN108918757A (zh) 化学反应光谱cielab色空间的l*值与物质量关系计算方法
CN108918759A (zh) 化学反应光谱的CIELAB色空间的b*值与物质量关系计算方法
CN108918756A (zh) 化学反应光谱的cielab色空间c*值与物质量计算方法
CN106706603A (zh) 一种检测生铁中元素含量的方法
CN104076080A (zh) Sbs改性沥青中改性剂含量的测定方法
CN110208260A (zh) 一种面粉中溴酸钾快速检测试剂、试剂盒以及检测方法
Safavi et al. Kinetic spectrophotometric determination of traces of sulphite

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20181130

RJ01 Rejection of invention patent application after publication