CN1850122A - 一种勾兑中药材或天然产物使提取物成分含量稳定的方法 - Google Patents

Abstract

本发明涉及一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其步骤如下：1)确定中药材或天然产物标准提取物的指标成分及其含量；2)确定中药材或天然产物勾兑后提取物质量合格时，其指标成分含量的相对偏差范围；3)取样并提取待勾兑各批次中药材或天然产物，测定其指标成分含量；4)确定中药材或天然产物勾兑参数、质量控制参数、勾兑目标，形成约束最优化问题；5)求解约束最优化问题，得到最优勾兑比；6)按最优勾兑比混合中药材或天然产物，经提取得指标成分含量合格的提取物。本发明方法可应用于中药、天然产物、生物制品、酒、香料、饮料、保健品、调味品、食用油、燃油、燃气、饲料的加工生产过程，保证其指标成分含量稳定。

Description

一种勾兑中药材或天然产物使提取物成分含量稳定的方法

技术领域

本发明涉及一种勾兑方法，特别是关于一种勾兑中药材或天然产物使提取物成分含量稳定的方法。

背景技术

中药材及天然产物是食品、保健品、药品、香料和化工原料等的重要来源。中国多种多样的自然条件，孕育了丰富的中药材及天然产物资源。根据《中国中药及天然产物资源》统计，我国有药用植物11146种，药用动物1581种，药用矿物80种，其中常用中药材及天然产物600多种。中药材和天然产物主要有人工栽培和野生两种生长方式，常用中药材及天然产物中有栽培品种200多种，年产量30万吨，约占常用药材品种的30％，年产量占常用药材品种年收购总量的近50％；野生药材收购品种400多种，年收购量40万吨，品种约占常用药材的70％，年收购量占收购总量的50％～60％。全国经营的药材品种1200多种，中药材及天然产物年收购总值160亿元左右。全国已建立中药材及天然产物生产基地600多个，生产专业户达34万户，种植面积达38.67万公顷，产量达35万吨。

保证中药材或天然产物的成分含量稳定具有重要意义。中药材及天然产物所含成分非常复杂，一味中药材或天然产物中往往含有几十种到上百种不同的化合物。中药材及天然产物正是依靠其含有的多种化学成分协同作用而达到保健和治疗疾病等目的。为了保证用药安全、疗效稳定、质量可控，必须对中药材或天然产物的重要成分的含量进行控制，使之稳定在一个合理的水平上。如德国著名的EGb 761银杏叶提取物等对成分含量控制十分严格。但是，中药材及天然产物的活性成分大部分为药用植物的次生代谢产物，受种植培育的气候与地理特征、采收季节、存储条件等因素影响，其化学成分的含量分布具有不同程度的波动性。在中药材及天然产物提取活性成分等加工过程中，不同批次原料的活性成分含量往往不一致，既使是道地药材，这种差异仍然较大，由此导致中药材及天然产物最终产品的活性成分含量差异很大，从而导致药效不稳定。总之，如何确保中药材及天然产物提取物质量的稳定，使各种需要控制的成分含量在一个允许的范围内波动，是中药及天然产物现代科技面临的难题之一，这一问题对天然产物提取物、中药注射剂的研制生产而言尤为重要。

为了保证中药材及天然产物提取物内在质量稳定均一，人们提出了一种勾兑中药材或天然产物提取物的方法，例如公开号为CN1586509A的专利申请中，采取先提取中药材或天然产物得到提取物，再按适当比例勾兑提取物，使其内在质量稳定的方法。该方法存在以下不足：首先需要得到多批次提取物才能进行勾兑，这样会占据大量存储容积，积压资金，生产周期也长。如果某些批次的提取物始终不能用于勾兑，就容易造成提取物变质浪费，而且成品的生产计划受提取物状况影响很大。此外，这种勾兑方法优化目标不够灵活，没有考虑勾兑物质量合格时使成本最低或其它性能指标最优，定义勾兑物质量合格条件不够灵活，没有提供限制勾兑物用量比例的方法。

发明内容

针对上述问题，本发明的目的是提供一种不需要多批次提取，操作简单灵活，可以避免浪费，降低成本的勾兑中药材或天然产物使提取物成分含量稳定的方法。

为实现上述目的，本发明采取以下技术方案：一种勾兑中药材或天然产物使提取物成分含量稳定的方法，它包括以下步骤：1)确定中药材或天然产物的标准提取物的指标成分及其含量；2)确定中药材或天然产物勾兑后提取物质量合格时其指标成分实际含量的相对偏差范围，及中药材或天然产物勾兑后提取物的指标成分实际含量与标准提取物的的指标成分含量相似系数的域值；3)取样并提取待勾兑各批次中药材或天然产物，测定各批次提取物指标成分含量，得到指标成分含量矩阵；4)确定中药材或天然产物勾兑质量控制方法的原料勾兑参数、质量控制参数、勾兑目标，形成中药材或天然产物勾兑约束最优化问题；5)求解中药材或天然产物勾兑约束最优化问题，得到最优的中药材或天然产物勾兑比；6)按照最优中药材或天然产物勾兑比混合各批次中药材或天然产物，经提取得到指标成分含量合格、与标准提取物品质相近的提取物。

所述中药材或天然产物包括单味中药材或天然产物，以及多味中药材或天然产物的混合物。

所述中药材或天然产物的标准提取物为下述中的一种：1)符合国家药典标准、部颁标准的中药或天然产物提取物；2)由若干批次中药材或天然产物的提取物，按等比例混合而成；3)若干批次中药材或天然产物等比例混合物的提取物；4)根据文献报道、经验或计算所确定的各指标成分的含量而制成的中药材或天然产物标准提取物。

所述指标成分是指需要对其含量进行控制的成分，包括有效成分；所述标准提取物指标成分的含量为r个指标成分的含量，并用Cs＝(Cs₁，...，Cs_r)^T表示，其中r为提取物指标成分的个数。

所述确定中药材或天然产物勾兑后提取物质量合格时，其指标成分实际含量的相对偏差范围，是指Ca＝(Ca₁，...，Ca_r)^T的上下相对偏差范围为 $c^u={(c_1^u,...,c_r^u)}^T,$ $c^l={(c_1^l,...,c_r^l)}^T,$ 即满足 ${\mathrm{Cs}}_i(1-c_i^l)\≤{\mathrm{Ca}}_i\≤{\mathrm{Cs}}_i(1+c_i^u),0\≤c_i^l\≤\mathrm{1,0}\≤c_i^u\≤1,i=1,...,r;$ 所述中药材或天然产物勾兑后提取物的指标成分实际含量用Ca表示，Ca＝(Ca₁，...，Ca_r)^T，所述中药材或天然产物勾兑后提取物的指标成分实际含量Ca与所述标准提取物的指标成分含量Cs之间的相似系数用Sa表示；所述中药材或天然产物勾兑后提取物的指标成分实际含量Ca与所述标准提取物的指标成分含量Cs之间的相似系数Sa的域值为S₀，即满足Sa≥S₀。

所述中药材或天然产物勾兑后提取物与标准提取物的相似系数，由夹角余弦法或相关系数法计算而得。

所述指标成分含量矩阵是指待勾兑t个批次中药材或天然产物的r个指标成分含量，得到大小为r×t的含量矩阵C，其中t为中药材或天然产物的批次数。

所述原料勾兑参数包括：t个批次中药材或天然产物的价格p＝(p₁，...，p_t)^T；t个批次中药材或天然产物勾兑时勾兑比x＝(x₁，...，x_t)^T的上下限 $x^u={(x_1^u,...,x_t^u)}^T,$ $x^l={(x_1^l,...,x_t^h)}^T,$ 即 $x_j^l\≤x_j\≤x_j^u,$ 其中 $0\≤x_j^l\≤x_l^u\≤1,j=1,...,t.$

所述质量控制参数包括：中药材或天然产物勾兑后提取物的r个指标成分含量偏离标准提取物含量时惩罚的权重w＝(w₁..w_i..w_r)^T，w_i≥0；指标成分含量上下偏差控制因子 ${\mathrm{\&eta;}}_i^u={({\mathrm{\&eta;}}_1^u,...,{\mathrm{\&eta;}}_r^u)}^T,{\mathrm{\&eta;}}_i^l={({\mathrm{\&eta;}}_1^l,...,{\mathrm{\&eta;}}_r^l)}^T,$ 使 ${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u),$ 其中勾兑后提取物的指标成分理论含量Cm＝C·x，Cm_i是Cm的第i个分量， $0<{\mathrm{\&eta;}}_i^l\&le;1,$ $0<{\mathrm{\&eta;}}_i^u\&le;1,i=1,...,r;$ 相似系数偏差控制因子ξ，使Sm≥ξ·S₀，其中Sm是Cs和Cm之间的相似系数， $1\&le;\mathrm{\&xi;}\&le;\frac{1}{S_0}.$

所述勾兑目标包括平方偏差质量最优、绝对值偏差质量最优、成本最低、其中一种成分含量最高、其中一种成分含量最低、其中一批中药材或天然产物用量最大、其中一批中药材或天然产物用量最小；

所述平方偏差质量最优勾兑目标的函数为

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{w}_i{\left(\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right)}^2;$

所述绝对值偏差质量最优勾兑目标的函数为

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{w}_i\left|\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right|;$

所述成本最低勾兑目标的函数为

$\min \underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{p}_j{x}_j;$

所述其中一种成分含量最高勾兑目标的函数为maxCm_i；

所述其中一种成分含量最低勾兑目标的函数为minCm_i；

所述其中一批中药材或天然产物用量最大勾兑目标的函数为maxx_j；

所述其中一批中药材或天然产物用量最小勾兑目标的函数为minx_j。

所述形成中药材或天然产物勾兑约束最优化问题包括中药材或天然产物勾兑方法平方偏差质量最优的最优化问题：

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{w}_i{\left(\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right)}^2$

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法绝对值偏差质量最优的最优化问题：

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{w}_i\left|\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right|$

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=1,2,...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法成本最低的最优化问题：

$\min \underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{p}_j{x}_j$

$s.t.\mathrm{Cm}=C\&CenterDot;x$

$\underset{i=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=1,2,...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法其中一种成分含量最高的最优化问题：

maxCm_i

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法其中一种成分含量最低的最优化问题：

minCm_i

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法其中一批中药材或天然产物用量最大的最优化问题：

maxx_j

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法其中一批中药材或天然产物用量最小的最优化问题：

minx_j

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...r$

Sm≥ξ·S₀。

所述指标成分含量偏差控制因子 ${\mathrm{\&eta;}}_i^u={({\mathrm{\&eta;}}_1^u,...,{\mathrm{\&eta;}}_r^u)}^T,{\mathrm{\&eta;}}_i^l={({\mathrm{\&eta;}}_1^l,...,{\mathrm{\&eta;}}_r^l)}^T,$ 其取值的方法是 ${\mathrm{\&eta;}}_i^l=1-{\mathrm{RBEr}}_i/c_i^l,{\mathrm{\&eta;}}_i^u=1-\mathrm{RBE}{r}_i/c_i^u,$ 其中RBEr_i是|(Cm_i-Ca_i)|/Cm_i在各种可能情况下取到的最大值，所述RBEr_i通过勾兑实验测得或估计得到。

所述求解中药材或天然产物勾兑约束最优化问题，采取单纯形法、Karmarkar法、罚函数法、障碍函数法、二次规划法、逐步二次规划法、信赖域法、遗传算法、模拟退火算法、粒子群算法、蒙特卡罗模拟法中的一种或多种。

本发明方法的实施需要根据勾兑的实际要求确定一个约束最优化问题(包括目标函数及其约束条件两部分)，在计算机上运行求解该问题的数学模型从而得到最优的原料勾兑比。其核心是建立约束最优化问题。现将有关原理说明如下：

1、确定中药材或天然产物原料的最优勾兑比的最优化问题

本发明通过求解一个约束最优化问题来获取原料最优勾兑比。最优化问题由目标函数和约束条件两个部分组成。

(1)最优化问题的目标函数

本发明提出的目标函数包括平方偏差质量最优、绝对值偏差质量最优、成本最低、其中一种成分含量最高、其中一种成分含量最低、其中一批原料用量最大、其中一批原料用量最小等七种目标函数。

平方偏差质量最优的目标函数其基本形式如下：

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{w}_i{\left(\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right)}^2$ (式1)

在(式1)中，目标函数是最小化中药材或天然产物原料勾兑后提取物指标成分含量Cm与标准提取物指标成分含量Cs相对偏差的加权平方和，它反映原料勾兑后提取物的各指标成分的总体偏差水平。

绝对值偏差质量最优的目标函数其基本形式如下：

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{w}_i\left|\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right|$ (式2)

在(式2)中，目标函数是最小化中药材或天然产物原料勾兑后提取物指标成分含量Cm与标准提取物指标成分含量Cs相对偏差的绝对值的加权和，它反映原料勾兑后提取物的各指标成分的总体偏差水平。

(式1)、(式2)中的权重w_i≥0对各成分的含量偏差进行惩罚，其取值根据需要确定。当w中各w_i不等时，意味着给各个成分的含量偏差分配不同的权重。当取各权重w_i相等时，目标函数可简化为

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{\left(\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right)}^2$ (式3)

或

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}\left|\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right|$ (式4)

根据总体质量最优的勾兑目标函数(式1)、(式2)，配合权值w的取值，可以推导出个别指标成分最优的勾兑目标函数。其方法是：设置w中仅有一个分量w_i不为零，如w＝(0..w_i..0)^T，此时目标函数(式1)、(式2)与下面的目标函数等价：

$\min w_i\mathrm{\&Sigma;}\left|\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right|$ (式5)

该目标函数的含义是最小化第i个指标成分的含量偏差。

成本最低的目标函数最小化原料投料成本，其目标函数为：

$\min \underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{p}_j{x}_j$ (式6)

其中一种成分含量最高的目标函数最大化第i个成分的含量，其目标函数为：

maxCm_i                                                     (式7)

其中一种成分含量最低的目标函数最小化第i个成分的含量，其目标函数为：

minCm_i                                                     (式8)

其中一种原料用量最大的目标函数最大化第j批原料的用量，其目标函数为：

maxx_j                                                      (式9)

其中一批原料用量最小的目标函数最小化第j批原料的用量，其目标函数为：

minx_j                                                      (式10)

在实际应用中，除以上七种目标函数外，还可以根据实际需要对以上目标函数进行修改或定义其它目标函数。

(2)最优化问题的约束条件

对于平方偏差质量最优、绝对值偏差质量最优、成本最低、其中一种成分含量最高、其中一种成分含量最低、其中一批原料用量最大、其中一批原料用量最小等七种目标函数，其约束条件都具有如下基本形式：

Cm＝C·x                                                   (式11)

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$ (式12)

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$ (式13)

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$ (式14)

Sm≥ξ·S₀                                                 (式15)

约束(式11)是原料勾兑后提取物指标成分理论含量Cm的计算公式，(式11)为矩阵矩形，如果写成分量形式则为Cm_i＝C·x，i＝1，2，..，r。

约束(式12)要求各种中药材或天然产物原料勾兑比之和为1。

约束(式13)要求各中药材或天然产物原料勾兑比x_j均大于等于一个预先给定的最小值x_j ^l，但又小于一个预先给定的最大值x_j ^u，即 $x_j^l\&le;x_j\&le;x_j^u,$ 其中 $0\&le;x_j^l\&le;x_l^u\&le;1,j=1,...,t.$ 利用约束(式13)可灵活地指定参与勾兑的各批次中药材或天然产物原料的比例范围，从而有选择性地多使用或少使用某些特定批次的中药材或天然产物原料。

通过约束(式14) ${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)$ 及指标成分含量偏差控制因子η_i ^u、η_i ^l来保证勾兑后提取物实际含量Ca合格的方法，即满足 ${\mathrm{Cs}}_i(1-c_i^l)\&le;{\mathrm{Ca}}_i\&le;{\mathrm{Cs}}_i(1+c_i^u),$ 其中c_i ^l、c_i ^u分别是合格提取物相对于Cs_i变化范围的下、上限比例， $0\&le;c_i^l\&le;1,0\&le;c_i^u\&le;1.$ 为了保证勾兑后提取物实际含量Ca_i合格，必须在勾兑时将理论含量Cm_i的变化范围从 ${\mathrm{Cs}}_i(1-c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+c_i^u)$ 缩小为 ${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l{c}_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^l{c}_i^u),$ 其中指标成分含量偏差控制因子取值

${\mathrm{\&eta;}}_i^l=1-{\mathrm{RBEr}}_i/c_i^l$ (式16)

${\mathrm{\&eta;}}_i^u=1-{\mathrm{RBEr}}_i/c_i^u$ (式17)

其中RBEr_i是实际含量Ca_i与理论含量Cm_i相对偏差|(Cm_i-Ca_i)|/Cm_i在各种可能情况下取到的最大值，称之为原料勾兑基础偏差，RBEr_i通过勾兑实验测得或估计得到。

通过相似系数偏差控制因子ξ进一步提高最优化模型中勾兑后提取物与标准提取物的相似系数域值的方法，即Sm≥ξ·S₀。约束(式15)中Sm是中药材或天然产物原料勾兑后提取物Cm与标准提取物Cs之间的相似系数，其中相似系数可定义在指标成分含量上、色谱指纹图谱上、或多个来源数据融合后的数据上，S₀是预先设定的相似系数下限值，低于该阈值的中药材或天然产物原料勾兑后提取物被认为不合格。由于相似系数偏差控制因子满足 $1\&le;\mathrm{\&xi;}\&le;\frac{1}{S_0},$ 因此Sm≥ξ·S₀进一步提高了相似系数域值，从而为中药材或天然产物原料勾兑后的提取物与标准提取物在相似系数上合格留下一定的裕量。其中相似系数的计算方法有多种，常用的方法有夹角余弦法、相关系数法等。其中，夹角余弦法计算公式为：

S＝a^Tb/(|a|·|b|)                                           (式18)

相关系数法计算公式为：

S＝cov(a，b)/(|a|·|b|)                                     (式19)

其中，cov是协方差函数， $\left|a\right|=\sqrt{\underset{i=1}{\overset{n}{\mathrm{\&Sigma;}}}{\left(a_i\right)}^2}$ 表示向量a的模长。

2、最优化问题的求解

根据最优化理论，以平方偏差质量最优、绝对值偏差质量最优为目标函数的勾兑问题是约束非线性规划问题。以平方偏差质量最优为目标函数的勾兑问题在不含约束(式15)时是凸二次规划问题。以成本最低、其中一种成分含量最高、其中一种成分含量最低、其中一批原料用量最大、其中一批原料用量最小为目标函数的最优化问题，在含约束(式15)时是一约束非线性规划问题，在不含约束(式15)时是一线性规划问题。

求解线性规划问题的经典方法有单纯形法、Karmarkar法等。求解约束非线性规划问题的经典方法有罚函数法、障碍函数法、二次规划、逐步二次规划方法和信赖域方法等。这些方法一般基于迭代搜索进行优化，若在有限的迭代步内找不到可行解，则可认为在给定的约束条件下，用当前批次中药材或天然产物原料勾兑时，不存在可行解，应加入新的中药材或天然产物原料批次重新计算，或者适度放松约束条件，直到找到满足要求的中药材或天然产物原料勾兑配方为止。非线性最优化方法有时仅能得到局部最优解，此时应调整优化的初始点，多次重启优化过程，得到全局最优解。除经典方法外，求解约束非线性规划问题还有一些智能化方法，如遗传算法、模拟退火算法、群智能方法(如粒子群算法)、蒙特卡罗(Monte Carlo)模拟法等。运用智能方法可得到优化问题的全局最优解。

根据勾兑的最优化问题选定合适的求解方法，并通过计算机程序设计实现求解该最优化问题的模型软件，求解该问题，得到最优勾兑比。

为便于理解上述内容，现将有关符号及数据说明如下，见表1。

表1

本发明方法建立了一种通过勾兑中药材或天然产物原料来保证其提取物指标成分含量稳定的方法，而不是采用对中药材或天然产物的提取物进行勾兑，其核心有三，一是原料勾兑法的提出，二是建立中药材或天然产物原料勾兑的最优化问题，包括目标函数及其约束条件，三是提出通过指标成分含量偏差控制因子、相似系数偏差控制因子将原料勾兑后提取物指标成分的实际含量控制在合格范围内的方法。

本发明采取对中药材或天然产物原料进行勾兑再从中提取活性成分的方法，具有以下优点：(1)本发明可以对不同产地、不同等级的中药材或天然产物原料搭配使用，既可以充分利用现有的中药材及天然产物资源，又达到了中药材及天然产物提取物质量稳定均一的目的，既扩大了原料来源，又降低了原料成本，而且可以事先确定能够用于勾兑的中药材或天然产物原料批次，便于企业制订采购计划和生产计划，节约运行成本。(2)与提取物勾兑方法相比，本发明方法不需要多批次提取，也不需要将多批次的提取物进行存储，可大量节约存储容积，缩短生产周期，利用提取操作标准化，使操作简易方便，同时避免了利用提取物勾兑时某些批次的提取物因长期没有投入勾兑而引起变质，造成的浪费。(3)本发明方法有平方偏差质量最优、绝对值偏差质量最优、成本最低、其中一种成分含量最高、其中一种成分含量最低、其中一批原料用量最大、其中一批原料用量最小等七种目标函数，目标函数灵活、通用、合理，可以充分适应勾兑的各种不同情况。(4)本发明提供了勾兑时对原料用量(勾兑比x)进行限制的方法，为用户管理原料库存提供方便，现有技术没有提供这种机制。(5)本发明对勾兑后提取物指标成分含量Ca偏离标准提取物指标成分含量Cs的偏差的限制方法更为灵活，Ca各成分含量的允许偏差的上下限约束分开，可以不相同，并能够处理某些组分含量在一定范围内是越高越好、或越低越好的情况。因此与现有技术使用绝对值约束，要求上下偏差相同的方法相比更加灵活。(6)本发明提供了对勾兑后提取物指标成分的理论含量Cm与实际含量Ca的差异进行控制的机制，本发明充分考虑到实施过程中理论含量Cm与实际含量Ca存在偏差，特别引入了两种提供裕量的偏差控制因子，有力地保证原料勾兑后的提取物合格，而现有技术没有提供这种机制。本发明方法可以广泛应用于中药、天然产物、生物制品、酒、香料、饮料、保健品、调味品、食用油、燃油、燃气、饲料等行业的加工生产。

附图说明

图1是本发明的步骤方框示意图

具体实施方式

以下通过具体实施例并结合附图详细说明本发明方法的实施过程。

实施例1：采用本发明方法勾兑丹参原料并保证其提取物指标成分含量稳定

丹参是丹参注射用水针剂和粉针剂的主要原料，也是香丹注射剂等其它注射剂的重要原料。本实施例将丹酚酸B、丹参素、原儿茶醛列为需要进行含量控制的指标成分。在全国范围内采集、收购野生及栽培丹参药材共计58批，其中野生14批，种植44批，药材由中医研究院中药及天然产物所鉴定。

设定本实施例的目标是考察本发明原料勾兑方法能否保证原料勾兑后提取物指标成分含量Cm与标准提取物的含量Cs相等，并考察原料勾兑基础偏差，即|(Cm_i-Ca_i)|/Cm_i在各种可能情况下取到的最大值。

采用本发明的原料勾兑法，勾兑丹参原料并保证其提取物指标成分含量稳定的具体步骤如下(如图1所示)：

1.确定丹参标准提取物及其指标成分含量。对58批丹参药材，经常规工艺提取后，用高效液相色谱法(HPLC)分别测得其色谱指纹图谱，并分别用丹酚酸B、丹参素、原儿茶醛的标准品标定含量，得到58批原料3个指标成分的含量。选择以上58个丹参提取物的等比例均匀混合物为标准提取物，通过计算平均值得到丹酚酸B、丹参素、原儿茶醛3个指标成分的含量分别为3.63mg/mL、0.49mg/mL、0.074mg/mL，它们构成标准提取物的Cs向量，即Cs＝(3.63 0.49 0.074)^T。

2.设定丹参药材原料勾兑后提取物质量合格时3个指标成分含量Ca相对于标准提取物含量Cs的相对偏差的上下域值均为0.1，即c^l＝c^u＝(0.1 0.1 0.1)^T。设定丹参药材原料勾兑后的提取物与标准提取物在指标成分含量上的相似系数域值S₀＝0.95。

3.为减少丹参药材勾兑时的原料批数，在58批丹参药材中随机选定其中的7批，从中取样并提取，测定这7批原料提取物3个指标成分含量，其结果见表2。由表2可见不同批次的原料提取物其指标成分含量差异非常大。

表2：标准提取物、原料提取物的指标成分含量及其价格

	标准提取物	原料1	原料2	原料3	原料4	原料5	原料6	原料7
									成分1丹酚酸B含量Cs1(mg/mL)	3.63	3.20	4.41	2.21	5.58	2.15	3.97	2.28
成分2丹参素含量Cs2(mg/mL)	0.49	0.64	0.52	0.32	0.67	0.12	0.71	0.37
									成分3原儿茶醛含量Cs3(mg/mL)	0.074	0.066	0.041	0.028	0.089	0.004	0.127	0.062
原料价格p(元)	100	90	110	60	120	40	130	80

七批原料提取物3个指标成分的含量构成一个3×7的含量矩阵C：

$C=\left(\begin{array}{ccccccc}3.20& 4.41& 2.21& 5.58& 2.15& 3.97& 2.28\\ 0.64& 0.52& 0.32& 0.67& 0.12& 0.71& 0.37\\ 0.066& 0.041& 0.028& 0.089& 0.004& 0.127& 0.062\end{array}\right)$

4.确定丹参药材原料勾兑质量控制方法的原料勾兑参数、质量控制参数、勾兑目标，形成约束最优化问题。

根据事先设定的本实施例目标，设置参数如下。

4.1确定原料勾兑参数

设定7批丹参药材原料的单价(元)，其价格向量为

       p＝(90 110 60 120 40 130 80)^T

对7批丹参药材勾兑比x的上下限不进行限制，即x^l＝(0 0 0 0 0 0 0)^T，x^u＝(1 1 1 1 1 1 1)^T。

4.2确定质量控制参数

设定丹参药材原料勾兑后提取物的3个指标成分含量偏离标准提取物含量时，对3个偏差采取相同的惩罚，故应取权重w＝(1 1 1)^T。

由于勾兑目标是考察原料勾兑方法能否保证原料勾兑后提取物Cm与标准提取物Cs相等，故须设定指标成分含量偏差控制因子极小，η^l＝η^u＝(10^-610^-610^-6)^T，从而利用约束 ${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)$ 对Cm提供接近于零的偏差控制裕量。

在本实施例中，由于c^l、c^u取值较小，并且指标成分含量偏差控制因子极小，即约束条件使Cm与Cs接近相等，故不考虑相似系数域值条件，它可以自动满足。

4.3确定勾兑目标函数

本实施例考虑平方偏差质量最优、绝对值偏差质量最优、成本最低、其中一种成分含量最高、其中一种成分含量最低、其中一批原料用量最大、其中一批原料用量最小等七种目标函数，具体的目标函数见表3。表3中成分编号、原料编号参见表2。

表3：原料勾兑目标函数及其简称

4.4确定最优化问题

所述各最优化问题如下所示。平方偏差质量最优的最优化问题为：

$\min {\left(\frac{{\mathrm{Cm}}_1-3.63}{3.63}\right)}^2+{\left(\frac{{\mathrm{Cm}}_2-0.49}{0.49}\right)}^2+{\left(\frac{{\mathrm{Cm}}_3-0.074}{0.074}\right)}^2$

s.t.Cm＝C·x

x₁+x₂+x₃+x₄+x₅+x₆+x₇＝1

0≤x₁≤1

0≤x₂≤1

0≤x₃≤1

0≤x₄≤1

0≤x₅≤1

0≤x₆≤1

0≤x₇≤1

3.63(1-0.1×10^-6)≤Cm₁≤3.63(1+0.1×10^-6)

0.49(1-0.1×10^-6)≤Cm₂≤0.49(1+0.1×10^-6)

0.074(1-0.1×10^-6)≤Cm₃≤0.074(1+0.1×10^-6)

绝对值偏差质量最优的最优化问题为：

$\min \left|\frac{{\mathrm{Cm}}_1-3.63}{3.63}\right|+\left|\frac{{\mathrm{Cm}}_2-0.49}{0.49}\right|+\left|\frac{{\mathrm{Cm}}_3-0.074}{0.074}\right|$

s.t.Cm＝C·x

x₁+x₂+x₃+x₄+x₅+x₆+x₇＝1

0≤x₁≤1

0≤x₂≤1

0≤x₃≤1

0≤x₄≤1

0≤x₅≤1

0≤x₆≤1

0≤x₇≤1

3.63(1-0.1×10^-6)≤Cm₁≤3.63(1+0.1×10^-6)

0.49(1-0.1×10^-6)≤Cm₂≤0.49(1+0.1×10^-6)

0.074(1-0.1×10^-6)≤Cm₃≤0.074(1+0.1×10^-6)

成本最低的最优化问题为：

min 90x₁+110x₂+60x₃+120x₄+40x₅+130x₆+80x₇

s.t.Cm＝C·x

x₁+x₂+x₃+x₄+x₅+x₆+x₇＝1

0≤x₁≤1

0≤x₂≤1

0≤x₃≤1

0≤x₄≤1

0≤x₅≤1

0≤x₆≤1

0≤x₇≤1

3.63(1-0.1×10^-6)≤Cm₁≤3.63(1+0.1×10^-6)

0.49(1-0.1×10^-6)≤Cm₂≤0.49(1+0.1×10^-6)

0.074(1-0.1×10^-6)≤Cm₃≤0.074(1+0.1×10^-6)

成分1含量最高的最优化问题为：

max Cm₁

s.t.Cm＝C·x

x₁+x₂+x₃+x₄+x₅+x₆+x₇＝1

0≤x₁≤1

0≤x₂≤1

0≤x₃≤1

0≤x₄≤1

0≤x₅≤1

0≤x₆≤1

0≤x₇≤1

3.63(1-0.1×10^-6)≤Cm₁≤3.63(1+0.1×10^-6)

0.49(1-0.1×10^-6)≤Cm₂≤0.49(1+0.1×10^-6)

0.074(1-0.1×10^-6)≤Cm₃≤0.074(1+0.1×10^-6)

成分1含量最低的最优化问题为：

min Cm₁

s.t.Cm＝C·x

x₁+x₂+x₃+x₄+x₅+x₆+x₇＝1

0≤x₁≤1

0≤x₂≤1

0≤x₃≤1

0≤x₄≤1

0≤x₅≤1

0≤x₆≤1

0≤x₇≤1

3.63(1-0.1×10^-6)≤Cm₁≤3.63(1+0.1×10^-6)

0.49(1-0.1×10^-6)≤Cm₂≤0.49(1+0.1×10^-6)

0.074(1-0.1×10^-6)≤Cm₃≤0.074(1+0.1×10^-6)

原料4用量最大的最优化问题为：

maxx₄

s.t.Cm＝C·x

x₁+x₂+x₃+x₄+x₅+x₆+x₇＝1

0≤x₁≤1

0≤x₂≤1

0≤x₃≤1

0≤x₄≤1

0≤x₅≤1

0≤x₆≤1

0≤x₇≤1

3.63(1-0.1×10^-6)≤Cm₁≤3.63(1+0.1×10^-6)

0.49(1-0.1×10^-6)≤Cm₂≤0.49(1+0.1×10^-6)

0.074(1-0.1×10^-6)≤Cm₃≤0.074(1+0.1×10^-6)

原料4用量最小的最优化问题为：

minx₄

s.t.Cm＝C·x

x₁+x₂+x₃+x₄+x₅+x₆+x₇＝1

0≤x₁≤1

0≤x₂≤1

0≤x₃≤1

0≤x₄≤1

0≤x₅≤1

0≤x₆≤1

0≤x₇≤1

3.63(1-0.1×10^-6)≤Cm₁≤3.63(1+0.1×10^-6)

0.49(1-0.1×10^-6)≤Cm₂≤0.49(1+0.1×10^-6)

0.074(1-0.1×10^-6)≤Cm₃≤0.074(1+0.1×10^-6)

5.在微型个人计算机上，利用Matlab程序设计语言，编制模型软件，实现单纯形法、罚函数法、障碍函数法、二次规划法、遗传算法，使该软件能求解中药材或天然产物原料勾兑约束最优化问题。

6.通过模型软件求解最优化问题

运用模型软件求解步骤4.4中的各个最优化问题，得到最优勾兑比x＝(x₁，...，x₇)^T以及相应的勾兑后提取物指标成分理论含量Cm和按该勾兑比投料的原料成本Cost，结果见表4。

表4：原料勾兑优化方法最优勾兑比、成分理论含量及其成本

目标函数	x1	x2	x3	x4	x5	x6	x7	Cm1	Cm2	Cm3	成本
												minQ2	0.58	1.05	1.25	22.23	30.30	37.22	7.37	3.63	0.49	0.074	95.5
minQ1	0	4.65	0	16.78	34.65	43.91	0	3.63	0.49	0.074	96.2
												minCost	2.40	0	0	20.3	35.63	41.67	0	3.63	0.49	0.074	94.9
maxCm	0	0	0	29.28	20.12	24.26	26.34	3.63	0.49	0.074	95.8
												minCm1	0	0	0	29.28	20.12	24.26	26.34	3.63	0.49	0.074	95.8
maxx4	0	0	0	29.28	20.12	24.25	26.34	3.63	0.49	0.074	95.8
												minx4	0	4.65	0	16.78	34.65	43.91	0	3.63	0.49	0.074	96.2

表4中勾兑比x_j单位为％，Cm单位为mg/mL，成本单位为元。以下各表相同。

从表4可见，在上述七种优化目标和给定的约束条件下，都能使Cm和Cs相等，但勾兑比不完全相同，它们在相应的优化目标下是最优的。勾兑目标函数不同，对应的原料成本也不同。

7.按照最优原料勾兑比混合各批次中药材或天然产物原料，经提取得到指标成分含量合格、与标准提取物品质相近的提取物。

按上述七种最优化问题所得最优勾兑比混合7个批次丹参原药材，经提取，再用高效液相色谱法(HPLC)测定得到丹酚酸B、丹参素、原儿茶醛成分的实际含量Ca，计算Ca与Cs的相对偏差Er_i＝(Cs_i-Ca_i)/Cs_i，并计算Ca与Cs在三个指标成分上的相关系数Sa。结果见表5。

表5：原料勾兑后提取物的指标成分含量及其相对偏差

目标函数	Ca1	Ca2	Ca3	Er1	Er2	Er3	Sa
								minQ2	3.80	0.50	0.072	-0.047	-0.020	0.027	1.00
minQ1	3.55	0.47	0.076	0.022	0.041	-0.027	1.00
								minCost	3.43	0.52	0.071	0.055	-0.061	0.041	1.00
minCm	3.73	0.50	0.070	-0.028	-0.020	0.054	1.00
								maxCm	3.58	0.51	0.069	0.014	-0.041	0.068	1.00
minx4	3.76	0.47	0.077	-0.036	0.041	-0.041	1.00
								maxx4	3.65	0.46	0.075	-0.006	0.061	-0.014	1.00

由表5可见，尽管原料勾兑后提取物的理论含量Cm与Cs相等，但由于实际勾兑操作会引入多种误差，使得实际含量Ca与理论含量Cm有偏离，但偏差并不太大，完全满足本实施例步骤2中设定的丹参药材原料勾兑后提取物质量合格时3个指标成分含量Ca相对于标准提取物含量Cs的相对偏差的上下域值均为0.1，及设定丹参药材原料勾兑后的提取物与标准提取物在指标成分含量上的相似系数域值S₀＝0.95的要求。

8.确定指标成分含量偏差控制因子η^l、η^u

本实施例表明，在本发明提出的七种勾兑目标下均能使Cm与Cs相等；但对原料勾兑方式，勾兑后提取物的实际含量Ca_i偏离其理论含量Cm_i。下面考察Ca_i与Cm_i的相对偏差|(Cm_i-Ca_i)|/Cm_i可能的最大值，即该成分的原料勾兑基础偏差RBEr_i，由于Cm与Cs相等，据表5可得RBEr＝(0.055 0.061 0.068)^T，在本实施例中，为了提供足够的偏差裕量并简便起见，设RBEr＝(0.080.080.08)^T。由于原料勾兑基础偏差RBEr_i的存在，为了保证勾兑后的实际含量Ca_i合格，必须将Cm_i的含量变化范围 ${\mathrm{Cs}}_i(1-c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+c_i^u)$ 缩小为 ${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l{c}_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^l{c}_i^u),$ 其由 ${\mathrm{\&eta;}}_i^l=1-{\mathrm{RBEr}}_i/c_i^l,{\mathrm{\&eta;}}_i^u=1-{\mathrm{RBEr}}_i/c_i^u.$ 由于c^l＝c^u＝(0.1 0.1 0.1)^T，RBEr＝(0.08 0.08 0.08)^T，因此η^l＝η^u＝(0.2 0.2 0.2)^T。

令本实施例中的参数η^l＝η^u＝(0.2 0.2 0.2)^T，其余约束条件及参数不变，对表4中的7个优化问题重新优化，结果如表6。

表6：原料勾兑优化方法最优勾兑比、成分理论含量及其成本

目标函数	x1	x2	x3	x4	x5	x6	x7	Cm1	Cm2	Cm3	成本
												minQ2	0.33	0.78	1.52	23.42	28.55	35.24	10.15	3.63	0.49	0.074	95.5
minQ1	0	4.65	0	16.78	34.65	43.91	0	3.63	0.49	0.074	96.2
												minCost	2.69	0	0	17.95	37.41	41.95	0	3.56	0.48	0.073	93.5
maxCm1	0	0	0	43.11	0	0	56.89	3.70	0.50	0.074	97.2
												minCm1	0	0	0	35.68	0	5.92	58.4	3.56	0.50	0.075	97.2
maxx4	0	0	0	43.19	1.97	0	54.85	3.70	0.49	0.073	96.5
												minx4	0	22.97	0	0	28.23	48.79	0	3.56	0.50	0.073	100.0

由于引入了偏差控制因子η_i ^l、η_i ^u，按上述最优勾兑比混合原料并提取，提取物的指标成分的实际含量能满足要求。

由表6可见，偏差控制因子η^l、η^u放宽后，某些目标函数的值得到明显改善，如最低成本优化从94.9元降为93.5元。

实施例2：采用本发明方法勾兑丹参原料并保证其提取物指标成分含量稳定

为充分体现本发明的使用效果，设定本实施例的目标是说明本发明在约束条件中可以对原料用量作强制要求、允许成分含量上下限偏差不同的应用方式。

本实施例有关符号的含义与实施例1相同。本实施例步骤1与实施例1相同，即Cs＝(3.63 0.49 0.074)^T。本实施例步骤2设定c^l＝(0.08 0.10 0.16)^T，c^u＝(0.20 0.10 0.08)^T，此时三个成分含量的允许范围各不相同；不考虑相似度约束。本实施例步骤3与实施例1相同，含量矩阵C见实施例1。本实施例步骤4.1设定7批丹参药材原料的单价(元)也与实施例1相同，即p＝(90 110 60 120 40 130 80)^T；但本实施例设定原料3用量不少于10％，原料5用量不少于20％，原料6用量不超过30％，即x^l＝(0 0 0.1 0 0.2 0 0)^T，x^u＝(1 1 1 1 1 1 0.3 1)^T。本实施例步骤4.2设定w＝(1 1 1)^T；根据实施例1中所得的丹参原料勾兑基础偏差RBEr＝(0.08 0.08 0.08)^T，以及(式16)、(式17)，计算可得 ${\mathrm{\&eta;}}_i^l={\left(\mathrm{00.20.5}\right)}^T,{\mathrm{\&eta;}}_i^u={\left(\mathrm{0.60.20}\right)}^T.$ 本实施例中，其步骤4.3设定七种目标函数，具体目标函数见表3；其步骤4.4由表3中的7个目标函数及上述约束条件及参数，不考虑相似度约束，形成约束最优化问题；其步骤5与实施例1相同；其步骤6运用模型软件解上述最优化问题，得到最优勾兑比x以及相应的勾兑后提取物指标成分理论含量Cm和成本Cost。结果见表7。

表7：原料勾兑优化方法最优勾兑比、成分理论含量及其成本

目标函数	x1	x2	x3	x4	x5	x6	x7	Cm1	Cm2	Cm3	成本
												minQ2	0	0	10	26.56	20.38	30	13.06	3.63	0.5	0.073	95.5
minQ1	0	0	10	26.55	20	30	13.45	3.63	0.5	0.074	95.6
												minCost	3.45	0	10	25.96	29.6	30	0.99	3.63	0.48	0.068	91.9
maxCm1	0	0	10	44.94	25.97	19.09	0	4.04	0.5	0.068	95.1
												minCm1	0	0	23.18	26.82	20	30	0	3.63	0.49	0.069	93.1
maxx4	0	0	10	46.95	20	12.93	10.12	4.01	0.5	0.068	95.2
												minx4	4.99	14.08	10	16.06	20	30	4.86	3.63	0.5	0.068	96.1

由于引入了偏差控制因子η_i ^l、η_i ^u，按上述最优勾兑比混合原料并提取，提取物的指标成分的实际含量能满足要求。

实施例3：采用本发明方法勾兑丹参原料并保证其提取物指标成分含量稳定

为充分体现本发明的使用效果，设定本实施例的目标是说明本发明利用不同的权重w进行优化的应用方式。本实施例考察平方偏差质量最优、绝对值偏差质量最优这两种使用权重w的目标函数，考虑不同的w取值情况。

本实施例有关符号的含义与实施例1相同。本实施例步骤1与实施例1相同，即Cs＝(3.63 0.49 0.074)^T。本实施例步骤2与实施例2相同，即c^l＝(0.08 0.10 0.16)^T，c^u＝(0.20 0.10 0.08)^T；不考虑相似度约束。本实施例步骤3与实施例1相同，含量矩阵C见实施例1。本实施例步骤4.1与实施例2相同，即p＝(90 110 60 120 40 130 80)^T；x^l＝(0 0 0.1 0 0.2 0 0)^T，x^u＝(1 1 1 1 1 1 0.3 1)^T。

本实施例步骤4.2确定质量控制参数时，设定丹参原料勾兑后提取物的3个指标成分含量偏离标准提取物含量，对3个偏差采取不同的惩罚，即当w＝(1 0 0)^T时，意味着最小化第一个成分丹酚酸B的含量偏差，同理w＝(0 1 0)^T、w＝(0 0 1)^T分别最小化丹参素和原儿茶醛的含量偏差。在此三种w取值情况下，平方偏差质量最优minQ2、绝对值偏差质量最优minQ1等价，故不予区分。当w_i不等时，意味着给三个成分的含量偏差分配不同的权重。本实施例步骤4.2具体设定w的约束见表8。

本实施例步骤4.3确定目标函数有两种类别，一是平方偏差质量最优的目标函数：

$\min w_1{\left(\frac{{\mathrm{Cm}}_1-3.63}{3.63}\right)}^2+w_2{\left(\frac{{\mathrm{Cm}}_2-0.49}{0.49}\right)}^2+w_3{\left(\frac{{\mathrm{Cm}}_3-0.074}{0.074}\right)}^2$

二是绝对值偏差质量最优的目标函数：

${\min w}_1\left|\frac{{\mathrm{Cm}}_1-3.63}{3.63}\right|+w_2\left|\frac{{\mathrm{Cm}}_2-0.49}{0.49}\right|+w_3\left|\frac{{\mathrm{Cm}}_3-0.074}{0.074}\right|$

本实施例步骤4.3设定目标函数为5个，具体目标函数见表8，其步骤4.4由表8中的5个目标函数及上述约束条件及参数，不考虑相似度约束，形成约束最优化问题；其步骤5与实施例1相同；其步骤6运用模型软件解上述最优化问题，得到最优勾兑比x以及相应的勾兑后提取物指标成分理论含量Cm和成本Cost，结果见表8。

表8：原料勾兑优化方法各批次用量、成分理论含量及其成本

目标函数	w	x1	x2	x3	x4	x5	x6	x7	Cm1	Cm2	Cm3	成本
													minQ2	(1 0 0)T	3.99	3.07	12.46	25.19	21.97	26.86	6.45	3.63	0.49	0.069	93.5
minQ2	(0 1 0)T	3.21	2.30	11.73	28.18	23.09	25.57	5.93	3.68	0.49	0.068	93.5
													minQ2	(0 0 1)T	0	0	10	27.60	20.00	30	12.40	3.66	0.50	0.074	96.0
minQ2	(5 2 1)T	0	0	10	26.58	20.91	30	12.51	3.63	0.49	0.073	95.3
													minQ1	(5 2 1)T	0	0	10	26.66	22.79	30	10.56	3.63	0.49	0.072	94.5

由于引入了偏差控制因子η_i ^l、η_i ^u，按上述最优勾兑比混合原料并提取，提取物指标成分的实际含量能满足要求。

本发明中提供的目标函数、约束条件、勾兑参数可以灵活搭配，以适应工业应用中通过原料勾兑方法保证提取物指标成分含量稳定的要求。不同的目标函数、不同的约束条件、不同的参数将得到不同的最优勾兑比。

本发明的原料勾兑方法可以和本申请人中请的另一项提取物勾兑方法的专利申请联合使用，可进一步提高中药材及天然产物提取物的质量稳定性。

以上实施例用来解释说明本发明，而不是对本发明进行限制。在本发明权利要求的保护范围内，对本发明做出的任何修改和改变，都属于本发明的保护范围。

Claims (13)

1、一种勾兑中药材或天然产物使提取物成分含量稳定的方法，它包括以下步骤：1)确定中药材或天然产物的标准提取物的指标成分及其含量；2)确定中药材或天然产物勾兑后提取物质量合格时其指标成分实际含量的相对偏差范围，及中药材或天然产物勾兑后提取物的指标成分实际含量与标准提取物的指标成分含量相似系数的域值；3)取样并提取待勾兑各批次中药材或天然产物，测定各批次提取物指标成分含量，得到指标成分含量矩阵；4)确定中药材或天然产物勾兑质量控制方法的原料勾兑参数、质量控制参数、勾兑目标，形成中药材或天然产物勾兑约束最优化问题；5)求解中药材或天然产物勾兑约束最优化问题，得到最优的中药材或天然产物勾兑比；6)按照最优中药材或天然产物勾兑比混合各批次中药材或天然产物，经提取得到指标成分含量合格、与标准提取物品质相近的提取物。

2、如权利要求1所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于：所述中药材或天然产物包括单味中药材或天然产物，以及多味中药材或天然产物的混合物。

3、如权利要求1或2所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于，所述中药材或天然产物的标准提取物为下述中的一种：1)符合国家药典标准、部颁标准的中药或天然产物提取物；2)由若干批次中药材或天然产物的提取物，按等比例混合而成；3)若干批次中药材或天然产物等比例混合物的提取物；4)根据文献报道、经验或计算所确定的各指标成分的含量而制成的中药材或天然产物标准提取物。

4、如权利要求1所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于：所述指标成分是指需要对其含量进行控制的成分，包括有效成分；所述标准提取物指标成分的含量为r个指标成分的含量，并用Cs＝(Cs₁，...，Cs_r)^T表示，其中r为提取物指标成分的个数。

5、如权利要求1所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于：所述确定中药材或天然产物勾兑后提取物质量合格时，其指标成分实际含量的相对偏差范围，是指Ca＝(Ca₁，...，Ca_r)^T的上下相对偏差范围为 $c^u={(c_1^u,...,c_r^u)}^T,$ $c^l={(c_1^l,...,c_r^l)}^T,$ 即满足 ${\mathrm{Cs}}_i(1-c_i^l)\&le;{\mathrm{Ca}}_i\&le;{\mathrm{Cs}}_i(1+c_i^u),0\&le;c_i^l\&le;1,$ $0\&le;c_i^u\&le;1,i=1,...,r;$ 所述中药材或天然产物勾兑后提取物的指标成分实际含量用Ca表示，Ca＝(Ca₁，...，Ca_r)^T，所述中药材或天然产物勾兑后提取物的指标成分实际含量Ca与所述标准提取物的指标成分含量Cs之间的相似系数用Sa表示；所述中药材或天然产物勾兑后提取物的指标成分实际含量Ca与所述标准提取物的指标成分含量Cs之间的相似系数Sa的域值为S₀，即满足Sa≥S₀。

6、如权利要求1或5所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于：所述中药材或天然产物勾兑后提取物与标准提取物的相似系数，由夹角余弦法或相关系数法计算而得。

7、如权利要求1所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于：所述指标成分含量矩阵是指待勾兑t个批次中药材或天然产物的r个指标成分含量，得到大小为r×t的含量矩阵C，其中t为中药材或天然产物的批次数。

8、如权利要求1所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于，所述原料勾兑参数包括：t个批次中药材或天然产物的价格p＝(p₁，...，p_t)^T；t个批次中药材或天然产物勾兑时勾兑比x＝(x₁，...，x_t)^T的上下限 $x^u={(x_1^u,...,x_t^u)}^T,$ $x^l={(x_1^l,...,x_t^h)}^T,$ 即 $x_j^l\&le;x_j\&le;x_j^u,$ 其中 $0\&le;x_j^l\&le;x_l^u\&le;1,j=1,...,t.$

9、如权利要求1所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于，所述质量控制参数包括：中药材或天然产物勾兑后提取物的r个指标成分含量偏离标准提取物含量时惩罚的权重w＝(w₁..w_i..w_r)^T，w_i≥0；指标成分含量上下偏差控制因子 ${\mathrm{\&eta;}}_i^u={({\mathrm{\&eta;}}_1^u,...,{\mathrm{\&eta;}}_r^u)}^T,$ ${\mathrm{\&eta;}}_i^l={({\mathrm{\&eta;}}_1^l,...,{\mathrm{\&eta;}}_r^l)}^T,$ 使 ${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u),$ 其中勾兑后提取物的指标成分理论含量Cm＝C·x，Cm_i是Cm的第i个分量， $0<{\mathrm{\&eta;}}_i^l\&le;1,$ $0<{\mathrm{\&eta;}}_i^u\&le;1,$ i＝1，...，r；相似系数偏差控制因子ξ，使Sm≥ξ·S₀，其中Sm是Cs和Cm之间的相似系数， $1\&le;\mathrm{\&xi;}\&le;\frac{1}{S_0}.$

10、如权利要求1所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于：所述勾兑目标包括平方偏差质量最优、绝对值偏差质量最优、成本最低、其中一种成分含量最高、其中一种成分含量最低、其中一批中药材或天然产物用量最大、其中一批中药材或天然产物用量最小；

所述平方偏差质量最优勾兑目标的函数为

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{w}_i{\left(\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right)}^2;$

所述绝对值偏差质量最优勾兑目标的函数为

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{w}_i\left|\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right|;$

所述成本最低勾兑目标的函数为

$\min \underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{p}_j{x}_j;$

所述其中一种成分含量最高勾兑目标的函数为max Cm_i；

所述其中一种成分含量最低勾兑目标的函数为min Cm_i；

所述其中一批中药材或天然产物用量最大勾兑目标的函数为max x_j；

所述其中一批中药材或天然产物用量最小勾兑目标的函数为min x_j。

11、如权利要求1所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于，所述形成中药材或天然产物勾兑约束最优化问题包括中药材或天然产物勾兑方法平方偏差质量最优的最优化问题：

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{w}_i{\left(\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right)}^2$

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法绝对值偏差质量最优的最优化问题：

$\min \underset{i=1}{\overset{r}{\mathrm{\&Sigma;}}}{w}_i\left|\frac{{\mathrm{Cm}}_i-{\mathrm{Cs}}_i}{{\mathrm{Cs}}_i}\right|$

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法成本最低的最优化问题：

$\min \underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{p}_j{x}_j$

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法其中一种成分含量最高的最优化问题：

maxCm_i

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法其中一种成分含量最低的最优化问题：

minCm_i

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法其中一批中药材或天然产物用量最大的最优化问题：

maxx_j

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀；

中药材或天然产物勾兑方法其中一批中药材或天然产物用量最小的最优化问题：

minx_j

s.t.Cm＝C·x

$\underset{j=1}{\overset{t}{\mathrm{\&Sigma;}}}{x}_j=1$

$x_j^l\&le;x_j\&le;x_j^u,j=\mathrm{1,2},...,t$

${\mathrm{Cs}}_i(1-{\mathrm{\&eta;}}_i^l\&CenterDot;c_i^l)\&le;{\mathrm{Cm}}_i\&le;{\mathrm{Cs}}_i(1+{\mathrm{\&eta;}}_i^u\&CenterDot;c_i^u)i=\mathrm{1,2},...,r$

Sm≥ξ·S₀。

12、如权利要求1或9所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于：所述指标成分含量偏差控制因子 ${\mathrm{\&eta;}}_i^u={({\mathrm{\&eta;}}_1^u,...,{\mathrm{\&eta;}}_r^u)}^T,$ ${\mathrm{\&eta;}}_i^l={({\mathrm{\&eta;}}_1^l,...,{\mathrm{\&eta;}}_r^l)}^T,$ 其取值的方法是 ${\mathrm{\&eta;}}_i^l=1-\mathrm{RBE}{r}_i/c_i^l,$ ${\mathrm{\&eta;}}_i^u=1-\mathrm{RBE}{r}_i/c_i^u,$ 其中RBEr_i是|(Cm_i-Ca_i)|/Cm_i在各种可能情况下取到的最大值，所述RBEr_i通过勾兑实验测得或估计得到。

13、如权利要求1所述的一种勾兑中药材或天然产物使提取物成分含量稳定的方法，其特征在于：所述求解中药材或天然产物勾兑约束最优化问题，采取单纯形法、Karmarkar法、罚函数法、障碍函数法、二次规划法、逐步二次规划法、信赖域法、遗传算法、模拟退火算法、粒子群算法、蒙特卡罗模拟法中的一种或多种。

Images