038] %%神经网络
[0039] clc
[0040] clearall
[0041] %%读入数据
[0042]t=xlsread('datal.xlsx'，'用于制备待测中药制剂的混合物料'，'Cl:Y1'）；％时 间
[0043]h=xlsread('datal.xlsx'，'用于制备待测中药制剂的混合物料'，' ）；％环 境湿度，其中：m=设置的环境湿度点个数+1
[0044]r=Xlsread('datal.xlsx'，'用于制备待测中药制剂的混合物料'，'C2:Yn'）；％吸 湿增重百分比，其中：n=测量所得吸湿增重百分比数值个数+1
[0045] %%构造输入和输出矩阵
[0046] [t_r,t_c]=size(t)
[0047] [h_r,h_c]=size(h)
[0048] fori=l:h_r
[0049] forj=l:t_c
[0050] X((i-l)*t_c+j, :) = [h(i)，t(j)];
[0051] Y((i_l)*t_c+j，：）=r(i，j);
[0052] end
[0053] end
[0054] X=X，；
[0055] Y=Y，；
[0056] %%BP神经网络
[0057] net=feedforwardnet(20)；
[0058] net=train(net;X，Y)%开始训练，其中pn，tn分别为输入输出样本
[0059] %%检验BP网络效果
[0060] Yl=sim(net,X)；
[0061] 三、将实际生产环境湿度值和最大吸湿增重值带入上述"环境湿度-吸湿时间-吸 湿增重"三维模型，所求算的吸湿时间即为安全生产时限。
[0062] 具体操作为：
[0063] 查阅药典附录制剂通则中所规定的待测剂型的最大含水量，将其设定为最大吸湿 增重值。将实际生产环境湿度值和最大吸湿增重值带入上述"环境湿度-吸湿时间-吸湿 增重"三维模型，所求算的吸湿时间即为安全生产时限。
[0064] 在上一段程序后面编写下列程序，即可求得吸湿时间（S卩：安全生产时限）：
[0065]%%xl中第一个输入数据是湿度值，第二个数据是吸湿增重百分比
[0066] xl=[0. 90; 100];
[0067] yl=sim(net,xl)
[0068] yl= ?
[0069] 实施例1 :预测环境湿度为50%条件下续断药材粉末制备散剂的安全生产时限
[0070] a)续断药材粉末等温吸湿数据的测定及其等温吸湿曲线的绘制
[0071] 在3个干燥器中分别配制相对湿度为3%、60%、95%的饱和盐溶液，将干燥器置于 25 °C的恒温条件下。
[0072] 取外径为50_、高为15_的干燥具塞称量瓶于实验前一天置于相应干燥器中平 衡24小时，精密称定其重量ml。将续断药材粉碎至120目，取粉末适量，平铺于上述称量瓶 中，控制续断粉末厚度为1_，精密称定重量m2。将称量瓶敞口，并与瓶盖同置于上述恒温 恒湿条件下，在1小时、2小时、3小时、4小时、6小时、8小时、12小时、24小时的时候盖好称 量瓶盖，精密称定重量m3,之后每隔24小时精密称定重量一次，直至吸湿增重百分率低于 0. 5%为止。
[0073] 按照下列公式计算增重百分率：
[0074]
[0075] 测量结果见表1所示。
[0076] 表1续断药材粉末等温吸湿数据测定结果
[0077]
[0078] 根据上述续断药材粉末等温吸湿数据的测定结果，绘制出续断药材粉末的3条等 温吸湿曲线，见图1。
[0079] b)建立续断药材粉末的"环境湿度-吸湿时间-吸湿增重"三维模型
[0080] 将上述3条等温吸湿曲线的数据存入Excel中，命名为datal.xlsx。进入Matlab 界面，编程如下程序，即可获得续断药材粉末的"环境湿度-吸湿时间-吸湿增重"的三维 模型Y1。
[0081] %%神经网络
[0082]clc
[0083]clearall
[0084] %%读入数据
[0085]t=xlsread('datal.xlsx'，' 续断药材粉末 '，'Cl:V1'）；％ 时间
[0086] h=xlsread('datal. xlsx'，' 续断药材粉末 '，' B2:B4'）；％ 环境湿度
[0087]r=xlsread(' datal. xlsx'，' 续断药材粉末 '，'C2:U4'）；％ 吸湿增重百分比
[0088] %%构造输入和输出矩阵
[0089] [t_r, t_c] =size (t)
[0090] [h_r, h_c] =size (h)
[0091] for i=l:h_r
[0092] for j=l:t_c
[0093] X((i-l)*t_c+j, :) = [h(i), t(j)];
[0094] Y((i-l)*t_c+j, :)=r (i, j);
[0095] end
[0096] end
[0097] X=X，；
[0098] Y=Y，；
[0099]%%BP神经网络
[0100] net=feedforwardnet (20);
[0101]net=train(net,X，Y)%开始训练，其中pn,tn分别为输入输出样本
[0102] %%检验BP网络效果
[0103] Yl=sim(net, X);
[0104]c)预测续断药材粉末在环境湿度为50%的条件下制备散剂时的安全生产时限
[0105] 查阅《中国药典》（2010版）一部附录制剂通则，确认散剂的最大含水量不得超过 9%，因此，可以认为，续断药材粉末在环境相对湿度为50%的条件下，制备的散剂具有质量 安全保障的最高含水量不应超过9%。所以，用续断药材粉末制备散剂时，必须在续断药材粉 末吸湿增重量达到含水量为9%之前完成全部的生产工序，这段时间即为续断药材粉末在 环境湿度为50%的条件下制备散剂时的安全生产时限。
[0106] 将设定的实际生产环境相对湿度值50%和续断药材粉末最大允许吸湿增重后的 含水量9%输入上述"环境湿度-吸湿时间-吸湿增重"三维模型中，运行下列程序，求解得 到的yl值即为安全生产时限。
[0107] 所述程序为：
[0108]xl=[50%;9%] ;%xl中第一个输入数据是环境湿度，第二个数据是含水量百分率
[0109] yl=sim(net, xl)
[0110] 程序输出结果为：
[0111] yl=67. 0
[0112] 即，50%湿度条件下，吸湿增重达到9%时所需的吸湿时间为67. 0小时。
[0113] 效果判断：
[0114] 通过本发明给出的"一种预测中药制剂安全生产时限的方法"对续断药材粉末安 全生产时限的预测，续断药材粉末在环境湿度为50%的条件下制备散剂时的安全生产时限 约为67. 0小时。
[0115]经测定，续断药材粉末在环境温度为25°C、相对湿度为50%时，达到含水量为9%时 需经时68小时，与预测结果极为接近，因此，可以认为本发明给出的技术方案具有实用性 和可操作性，对生产和科研具有指导意义。
[0116] 实施例2 :预测环境湿度为60%条件下续断垂盆草浸膏粉制备颗粒剂的安全生产 时限
[0117] a)垂盆草浸膏粉等温吸湿数据的测定及其等温吸湿曲线的绘制
[0118] 测定方法和使用的器具与实施例1相同，区别在于所用干燥器为6个，分别配制相 对湿度为3%、33%、43%、60%、75%、95%的饱和盐溶液，测定时间为1小时、2小时、3小时、4小 时、6小时、8小时、12小时、24小时，之后每隔24小时精密称定重量一次，直至达到平衡吸 湿量为止。
[0119] 测量结果见表2所示。
[0120] 表2垂盆草浸膏粉等温吸湿数