实验二色素的稳定性试验
一、概述
色素是对食品进行着色处理的添加剂，使用的目的是为了保持食品原料的原有颜色，弥补加工中天然色素在加工中的变化导致产品的变色现象；使产品的颜色与其风味保持一致；满足消费者的“心理”要求。
食用色素包括天然色素和人工合成色素两大类。无论哪类色素，在加工过程中，常常会由于各种因素而发生变色，衡量色素在加工中色调稳定性的指标有坚牢度和稳定性，前者是指在食品固体成分表面上的色素，对环境因素的抵抗作用，而后者是描述均匀分散在食品表面凝固的色素对环境的抵抗能力。
掌握影响色素稳定性的因素，对于正确使用色素具有重要意义。
二、实验目的
1.掌握食用色素坚牢度和稳定性的评价方法
2.掌握天然色素的提取方法
3.掌握影响色素坚牢度与稳定性的因素
三、实验原理
各种色素都有其固定的吸收波长和吸光度值，其坚牢度和稳定性可以用以下指标表示：耐热性、耐光性、耐酸性、耐碱性、耐氧化性、耐还原性、耐金属离子性和耐菌性。
将色素溶液加热、调节PH值、加入氧化剂和还原剂以及金属离子，通过测定上述各溶液的吸光度，并与之正常色素溶液的吸光度进行比较，就可以确定色素的坚牢度和稳定性。
四、实验原料与仪器
1.原料：新鲜花卉、柠檬黄， 1%高锰酸钾、1%Na
2S
2
O
3
、1%抗坏血酸、1%三氯化
铁、1%氯化铜、5%H
2O 2
2.仪器：烧杯、容量瓶（25ml、100ml、1000ml）、移液管（1ml）、试剂瓶、量筒、电炉、酸度计、分光光度计、水浴锅
五、实验方法
1.试剂的配制
（1）色素溶液的配制
称取柠檬黄1.25g，分别放入烧杯中，加入少量蒸馏水溶解后，移入1000ml 容量瓶，用蒸馏时定容至刻度。
取新鲜花卉100g，加水于60℃浸泡至花卉颜色变淡，将提取液移入1000ml 容量瓶定容。
（2）氧化剂、还原剂的配制
分别称取1g高锰酸钾、硫代硫酸钠、抗坏血酸放入烧杯，用少量蒸馏水溶解，移入100ml容量瓶中，以蒸馏水定容至刻度。
（3）金属离子溶液的配制
分别称取4.84g三氯化铁、2.116g氯化铜与烧杯中，用蒸馏水溶解，定容至100ml。
2.色素溶液特征吸光度的测定
（1）色素溶液吸光值的测定
取柠檬黄溶液0.2ml，加4.8 ml水，用分光光度计在 428nm波长下测定其吸光度。
（2）天然色素溶液吸光度值的测定
取天然色素溶液，用分光光度计在190－800nm波长下测定其吸光度，根据测定结果绘制吸光度值―吸光波长特征曲线，以最大吸光度的波长为该色素的特征吸光波长。
3.色素坚牢度与稳定性测定
（1）温度的影响
分别吸取天然色素、柠檬黄色素溶液0.2ml，加4.8 ml水，于70℃水浴锅内恒温水浴1h，冷却后于相应的特征波长下测定其吸光度值。
（2）氧化剂、还原剂的影响
分别吸取天然色素、柠檬黄色素溶液0.2ml，加4.8 ml水，加入1ml5%H
2O
2、
1ml1%Na
2S
2
O
3
溶液，在特征波长下测定其吸光度值。
（3）金属离子的影响
分别吸取天然色素、柠檬黄色素溶液0.2ml，加4.8 ml水，，分别加入0.5ml 的Fe3+（Cu2+）溶液，在特征波长下测定其吸光度值。
六、结果与分析
1、色素溶液的特征吸光值
柠檬黄λ428nm =
天然色素的特征波长λ＝ ；该波长下的吸光度值为А＝ 天然色素吸光值――吸光波长的曲线图
吸光度--吸光波长
0.2
0.4
0.6
0.8
1
1.2
12
3
吸光波长吸光度
2、温度、氧化剂、还原剂、金属离子对色素溶液的吸光度产生什么影响？
6、的影响
七、思考题：
1、影响食用色素坚牢度与稳定性的因素有那些？
2、比较天然色素与人工合成色素的稳定性
实验课程：数值计算方法专业：数学与应用数学班级：08070141 学号：37 姓名：汪鹏飞 中北大学理学院
实验1 赛德尔迭代法 【实验目的】 熟悉用塞德尔迭代法解线性方程组 【实验内容】 1.了解MATLAB 语言的用法 2.用塞德尔迭代法解下列线性方程组 1234123412341234 54 1012581034 x x x x x x x x x x x x x x x x ---=-??-+--=?? --+-=??---+=? 【实验所使用的仪器设备与软件平台】 计算机，MATLAB7.0 【实验方法与步骤】 1.先找出系数矩阵A ，将前面没有算过的x j 分别和矩阵的(,)A i j 相乘,然后将累加的和赋值给sum ，即(),j s u m s u m A i j x =+?.算 出()/(,) i i x b sum A i i =-,依次循环,算出所有的i x 。 2.若i x 前后两次之差的绝对值小于所给的误差限ε，则输出i x .否则重复以上过程，直到满足误差条件为止. 【实验结果】 (A 是系数矩阵，b 是右边向量，x 是迭代初值，ep 是误差限) function y=seidel(A,b,x,ep) n=length(b); er=1; k=0; while er>=ep
k=k+1; for i=[1:1:n] q=x(i); sum=0; for j=[1:1:n] if j~=i sum=sum+A(i,j)*x(j); end end x(i)=(b(i)-sum)/A(i,i); er=abs(q-x(i)); end end fprintf('迭代次数k=%d\n',k) disp(x') 【结果分析与讨论】 >> A=[5 -1 -1 -1;-1 10 -1 -1;-1 -1 5 -1;-1 -1 -1 10]; b=[-4 12 8 34]; seidel(A,b,[0 0 0 0],1e-3) 迭代次数k=6 0.99897849430002 1.99958456867649 2.99953139743435 3.99980944604109
实验一误差分析 实验1.1（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 其中ε（1.1）和（1.221,,,a a 的输出b ”和“poly ε。 （1（2 （3）写成展 关于α solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。 实验过程： 程序： a=poly(1:20); rr=roots(a); forn=2:21 n form=1:9 ess=10^(-6-m);
ve=zeros(1,21); ve(n)=ess; r=roots(a+ve); -6-m s=max(abs(r-rr)) end end 利用符号函数：（思考题一）a=poly(1:20); y=poly2sym(a); rr=solve(y) n
很容易的得出对一个多次的代数多项式的其中某一项进行很小的扰动，对其多项式的根会有一定的扰动的，所以对于这类病态问题可以借助于MATLAB来进行问题的分析。 学号：06450210 姓名：万轩 实验二插值法
附录ⅩⅠⅩ C原料药与药物制剂稳定性试验指导原则 稳定性试验的目的是考察原料药或药物制剂在温度、湿度、光线的影响下随时间变化的规律，为药品的生产、包装、贮存、运输条件提供科学依据，同时通过试验建立药品的有效期。 稳定性试验的基本要求是：（1）稳定性试验包括影响因素试验、加速试验与长期试验。影响因素试验用一批原料药或一批制剂进行。加速试验与长期试验要求用三批供试品进行。（2）原料药供试品应是一定规模生产的，供试品量相当于制剂稳定性试验所要求的批量，原料合成工艺路线、方法、步骤应与大生产一致。药物制剂供试品应是放大试验的产品，其处方与工艺应与大生产一致。药物制剂如片剂、胶囊剂，每批放大试验的规模，片剂至少应为10 000片，胶囊剂至少应为10 000粒。大体积包装的制剂如静脉输液等，每批放大规模的数量至少应为各项试验所需总量的10倍。特殊品种、特殊剂型所需数量，根据情况另定。（3）供试品的质量标准应与临床前研究及临床试验和规模生产所使用的供试品质量标准一致。（4）加速试验与长期试验所用供试品的包装应与上市产品一致。（5）研究药物稳定性，要采用专属性强、准确、精密、灵敏的药物分析方法与有关物质（含降解产物及其他变化所生成的产物）的检查方法，并对方法进行验证，以保证药物稳定性试验结果的可靠性。在稳定性试验中，应重视降解产物的检查。（6）由于放大试验比规模生产的数量要小，故申报者应承诺在获得批准后，从放大试验转入规模生产时，对最初通过生产验证的三批规模生产的产品仍需进行加速试验与长期稳定性试验。 本指导原则分两部分，第一部分为原料药，第二部分为药物制剂。 一、原料药 原料药要进行以下试验。 （一）影响因素试验 此项试验是在比加速试验更激烈的条件下进行。其目的是探讨药物的固有稳定性、了解影响其稳定性的因素及可能的降解途径与降解产物，为制剂生产工艺、包装、贮存条件和建立降解产物分析方法提供科学依据。供试品可以用一批原料药进行，将供试品置适宜的开口容器中（如称量瓶或培养皿），摊成≤5mm厚的薄层，疏松原料药摊成≤10mm厚薄层，进行以下试验。当试验结果发现降解产物有明显的变化，应考虑其潜在
L- 腈化物稳定性试验报告 一、概述 L-腈化物是L- 肉碱生产过程中的第一步中间体（第二步中间体： L-肉碱粗品；第三步中间体：L-肉碱潮品），由于L- 肉碱生产工艺为 间歇操作，即每生产一步中间体，生产完毕并出具合格检测报告后，存 入中间体仓库，以备下一步生产投料所需。根据本公司L- 肉碱产品的 整个生产周期，L- 腈化物入库后可能存放的最长时间为4 周（约28 天）。以此周期为时间依据制定了L- 腈化物稳定性试验方案，用于验 证L-腈化物在再试验期限内的各项质量指标数据的稳定性，并且能否符 合L- 腈化物的质量标准，此次稳定性试验的整个周期为28 天，具体 的稳定性试验方案以ICH 药物稳定性指导原则为基础制定，以确保L- 腈化化物稳定性试验的可操作性。 二、验证日期 2010 年1 月13 日- 2010 年2 月10 日 三、验证方案 1）样品储存和包装： 考虑到L- 腈化物今后的贮藏、使用过程，本次用于稳定性试验的样品 批次与最终规模生产所用的L- 腈化物的包装和放置条件相同。 2）样品批次选择：此次稳定性试验共抽取三批样品，且抽取样品的批次与 最终规模生产时的合成路线和生产工艺相同
3）抽样频率和日期：从2010.1.13 起，每隔7 天取样一次，共取五次，具体日期为：2010.1.13 、2010.1.20 、2010.1.27 、 2010.2.3 、2010.2.10 ，以确保试验次数足以满足L- 腈化物的稳 定性试验的需要。。 4）检测项目：根据L- 腈化物的质量标准的规定，此次稳定性试验的检测项目共五项，分别为外观、氯含量、熔点、比旋度、干燥失重。这 些指标在L- 腈化物的储存过程中可能会发生变化，且有可能影响 其质量和有效性。 5）试样来源和抽样：L- 腈化物由公司102 车间生产，经检测合格后储存于中间体仓库，本次稳定性试验的L- 腈化物均取自于该中间体仓 库，其抽样方法和抽样量均按照L- 腈化物抽样方案进行抽样。抽 样完毕后直接进行检测分析，并对检测结果进行登记，保存，作为稳 定性数据评估的依据。 四、稳定性试验数据变化趋势分析及评估 通过对三批L- 腈化物的稳定性试验，对其物理、化学方面稳定性资料进行评价，旨在建立未来相似情况下，大规模生产出的L- 腈化物是否适用 现有的再试验期（28天）。批号间的变化程度是否会影响未来生产的
本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 2016年 4 月 8 日
y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为：f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法： // 方程求根（割线法）.cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"
心得体会 使用不同的方法，可以不同程度的求得方程的解，通过二分法计算的程序实现更加了解二分法的特点，二分法过程简单，程序容易实现，但该方法收敛比较慢一般用于求根的初始近似值，不同的方法速度不同。面对一个复杂的问题，要学会简化处理步骤，分步骤一点一点的循序处理，只有这样，才能高效的解决一个复杂问题。
目的：考核药品质量稳定性，为确定药品的合理有效期提供科学依据。 范围：原料、中间产品、成品的稳定性试验 1．药物稳定性试验的目的是考察药物在规定的条件下随时间变化的规律，为药品的生产、包装、贮存、运输条件提供科学依据，同时通过试验建立药品的有效期。 2．稳定性试验的基本要求 2.1稳定性试验包括影响因素试验、加速试验和长期试验。其中影响因素试验适用于原料药，而加速试验和长期试验适用于原料药和制剂。 2.2 加速试验和长期试验样品所用容器、包装材料和包装方式应与上市产品一致。 2.3 应重视有关物质的检查。 3．影响因素试验 影响因素试验是在比加速试验更强烈的条件下进行，其目的是探讨药物的固有稳定性，了解影响其稳定性的因素及可能的降解途径与降解产物，为制剂生产工艺、包装、贮存条件与建立降解产物的分析方法提供科学依据。试验时，将样品置于适宜的容器中，摊成≤5mm厚的薄层，疏松原料药可摊成≤10mm厚的薄层，进行以下试验： 3.1高温试验 3.1.1将供试品置于适宜的容器中（如称量瓶或培养皿中）摊成≤5mm厚的薄层，疏松原料药可摊成≤10mm厚的薄层，在60℃下放置10天，于第5天和
第10天取样，按照稳定性重点考察项目进行检测，即：性状、熔点、含量、有关物质、吸湿性、以及根据药品性质选定的考察项目。 3.1.2 若检测结果与0天结果比较，有明显变化（如含量下降5%），则40℃下按同法进行试验。如60℃的结果无明显变化，不再进行40℃试验。 3.2 高湿度试验 在25℃将供试品开口置于相对湿度90±5%的恒湿密闭容器中（在密闭容器中放置硝酸钾饱和溶液），放置10天，于第5天和第10天取样，检验上述项目，重点考察供试品的吸湿性。如吸湿增重5%以上，则在相对湿度75±5%条件下（在密闭容器中放置氯化钠饱和溶液），同法进行试验。如吸湿增重5%以下，而且其他项目符合要求，则不再进行此项试验。 3.3 强光照射试验 将供试品开口置于装有日光灯的光照箱或适宜的光照箱中，照度为4500±500lx，放置10天，于第5天和第10天取样，检验上述项目，特别要注意外观的变化。光照箱要有恒定的照度，要注意防尘。 4．加速试验 4.1加速试验是在超常的条件下进行，其目的是通过加速药物的化学或物理变化，探讨药物的稳定性，为药品评审、包装、运输及贮存提供必要的资料。具体试验方法为： 取供试品三批，按市售包装，在温度40±2℃，相对湿度75±5%（在密闭容器中放置氯化钠饱和溶液）条件下，放置六个月。分别在第1个月、2个月、3个月、6个月末取样检测。按稳定性重点考察项目检测。如在6个月内供试品经检测不符合质量标准，则应在30±2℃，相对湿度60±5%（在密闭容器中放置亚硝酸钠饱和溶液）条件下，同法进行加速试验，时间仍为6个月。 4.2对温度特别敏感的药物，预计只能在冰箱中（4～8℃）保存的药物，其加速试验，可在25±2℃、相对湿度60±5%（在密闭容器中放置亚硝酸钠饱和溶液）条件下进行，时间为6个月。 4.3 对于包装在半透性容器的制剂，如塑料袋装溶液、塑料瓶装滴眼剂和滴
专业 序号 姓名 日期 实验1 算法的数值稳定性实验 【实验目的】 1．掌握用MATLAB 语言的编程训练，初步体验算法的软件实现； 2．通过对稳定算法和不稳定算法的结果分析、比较，深入理解算法的数值稳定性及其重要性。 【实验内容】 1．计算积分 ()dx a x x I n ?+=1 0) (n (n=0,1,2......，10) 其中a 为参数，分别对a=0.05及a=15按下列两种方案计算，列出其结果，并对其可靠性，说明原因。 2．方案一 用递推公式 n aI I n 1 1n + -=- (n=1,2,......,10) 递推初值可由积分直接得)1 ( 0a a In I += 3. 方案二 用递推公式 )1 (11-n n I a I n +-= (n=N,N-1,......,1) 根据估计式 ()()()11111+<<++n a I n a n 当1 n a +≥n 或 ()()n 1 111≤<++n I n a 当1 n n a 0+< ≤ 取递推初值为 ()()()() 11212])1(1111[21N +++=++++≈N a a a N a N a I 当1 a +≥ N N 或
()()]1111[21N N a I N +++= 当1 a 0+< ≤N N 计算中取N=13开始 【解】:手工分析怎样求解这题。 【计算机求解】:怎样设计程序？流程图？变量说明？能否将某算法设计成具有形式参数的函数 形式？ 【程序如下】: % myexp1_1.m --- 算法的数值稳定性实验 % 见 P11 实验课题(一) % function try_stable global n a N = 20; % 计算 N 个值 a =0.05;%或者a=15 % %-------------------------------------------- % % [方案I] 用递推公式 %I(k) = - a*I(k-1) + 1/k % I0 =log((a+1)/a); % 初值 I = zeros(N,1); % 创建 N x 1 矩阵(即列向量),元素全为零 I(1) =-a*I0+1; for k = 2:N I(k) =-a*I(k-1)+1/k; end % %--------------------------------------------
(此文档为word格式，下载后您可任意编辑修改！) 2012级6班###(学号)计算机数值方法 实验报告成绩册 姓名:宋元台 学号： 成绩：
数值计算方法与算法实验报告 学期： 2014 至 2015 第 1 学期 2014年 12月1日课程名称: 数值计算方法与算法专业:信息与计算科学班级 12级5班 实验编号： 1实验项目Neton插值多项式指导教师：孙峪怀 姓名：宋元台学号：实验成绩： 一、实验目的及要求 实验目的： 掌握Newton插值多项式的算法，理解Newton插值多项式构造过程中基函数的继承特点，掌握差商表的计算特点。 实验要求： 1. 给出Newton插值算法 2. 用C语言实现算法 二、实验内容 三、实验步骤(该部分不够填写.请填写附页)
1.算法分析： 下面用伪码描述Newton插值多项式的算法： Step1 输入插值节点数n，插值点序列{x(i),f(i)},i=1,2,……,n,要计算的插值点x. Step2 形成差商表 for i=0 to n for j=n to i f(j)=((f(j)-f(j-1)(x(j)-x(j-1-i)); Step3 置初始值temp=1,newton=f(0) Step4 for i=1 to n temp=(x-x(i-1))*temp*由temp(k)=(x-x(k-1))*temp(k-1)形成 (x-x(0).....(x-x(i-1)* Newton=newton+temp*f(i); Step5 输出f（x）的近似数值newton(x)=newton. 2.用C语言实现算法的程序代码 #include
范围：药物制剂。 责任：检验员、QA监控员、化验室主任、质保科科长、质量部负责人。 内容： 稳定性试验的目的是考察原料药或药物制剂在温度、湿度、光线的影响下随时间变化的规律，为药品的生产、包装、贮存、运输条件提供科学依据，同时通过试验建立药品的有效期。 稳定性试验的基本要求是：（1）稳定性试验包括影响因素试验、加速试验与长期试验。影响因素试验用1批原料药进行。加速试验与长期试验要求用3批供试品进行。（2）原料药供试品应是一定规模生产的。供试品量相当于制剂稳定性实验所要求的批量，原料药物合成工艺路线、方法、步骤应与大生产一致。药物制剂的供试品应是放大试验的产品其处方与生产工艺应与大生产一致。药物制剂如片剂、胶囊剂，每批放大试验的规模，片剂至少应为10 000片，胶囊剂至少应为10 000粒。大体积包装的制剂如静脉输液等，每批放大规模的数量至少应为各项试验所需总量的10倍。特殊剂型、特殊品种所需数量，根据具体情况另定。(3)供试品的质量标准应与临床前研究及临床试验和规模生产所使用的供试品质量标准一致。(4)加速试验与长期试验所用供试品的包装应与上市产品一致。(5)研究药物稳定性，要采用专属性强、准确、精密、灵敏的药物分析方法与有关物质(含降解产物及其他变化所生成的产物)的检查方法，并对方法进行验证，以保证药物稳定性试验结果的可靠性。在稳定性试验中，应重视降解产物的检查。（6）由于放大试验比规模生产的数量要小，故申报者应承诺在获得批准后，从放大试验转入规模生产时，对最初通过生产验证的3批规模生产的产品仍需进行加速试验与长期稳定性试验。 本指导原则分两部分，第一部分为原料药，第二部分为药物制剂。 1.原料药 原料药要进行以下试验。 1.1影响因素试验 此项试验是在比加速试验更激烈的条件下进行。其目的是探讨药物的固有稳定性、了解影响其稳定性的因素及可能的降解途径与降解产物，为制剂生产工艺、包装、贮存条件与建立降解产物的分析方法提供科学依据。供试品可以用一批原料药进行，将供试品置适宜的开口容器中(如称量瓶或培养皿)，摊成≤5mm厚的薄层，疏松原料药摊成≤10mm厚薄层，进行以下
附件3 特殊医学用途配方食品稳定性研究要求（试行） 一、基本原则 特殊医学用途配方食品稳定性研究是质量控制研究的重要组成部分，其目的是通过设计试验获得产品质量特性在各种环境因素影响下随时间 稳定性研究用样品应在满足《特殊医学用途配方食品良好生产规范》要求及商业化生产条件下生产，产品配方、生产工艺、质量要求应与注册申请材料一致，包装材料和产品包装规格应与拟上市产品一致。 影响因素试验、开启后使用的稳定性试验等采用一批样品进行；加速试验和长期试验分别采用三批样品进行。 （二）考察时间点和考察时间
稳定性研究目的是考察产品质量在确定的温度、湿度等条件下随时间变化的规律，因此研究中一般需要设置多个时间点考察产品的质量变化。考察时间点应基于对产品性质的认识、稳定性趋势评价的要求而设置。加速试验考察时间为产品保质期的四分之一，且不得少于3个月。长期试验总体考察时间应涵盖所预期的保质期，中间取样点的设置应当考虑产品的稳定性特点和产品形态特点。对某些环境因素敏感的产品，应适当增加考 3.检验方法：稳定性试验考察项目原则上应当采用《食品安全国家标准特殊医学用途配方食品通则》（GB 29922）、《食品安全国家标准特殊医学用途婴儿配方食品通则》（GB 25596）规定的检验方法。国家标准中规定了检验方法而未采用的，或者国家标准中未规定检验方法而由申请人自行提供检验方法的，应当提供检验方法来源和（或）方法学验证资料。检验方法应当具有专属性并符合准确度和精密度等相关要求。
四、试验方法 （一）加速试验 加速试验是在高于长期贮存温度和湿度条件下，考察产品的稳定性，为配方和工艺设计、偏离实际贮存条件产品是否依旧能保持质量稳定提供依据，并初步预测产品在规定的贮存条件下的长期稳定性。加速试验条件由申请人根据产品特性、包装材料等因素确定。 ％。如在6 温度 ％， 25℃±2℃ 长期试验是在拟定贮存条件下考察产品在运输、保存、使用过程中的稳定性，为确认贮存条件及保质期等提供依据。长期试验条件由申请人根据产品特性、包装材料等因素确定。 长期试验考察时间应与产品保质期一致，取样时间点为第一年每3个月末一次，第二年每6个月末一次，第3年每年一次。 如保质期为24个月的产品，则应对0、3、6、9、12、18、24月样品进行
实验一 误差分析 实验（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 )1.1() ()20()2)(1()(20 1∏=-=---=k k x x x x x p 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 )2.1(0 )(19=+x x p ε 其中ε是一个非常小的数。这相当于是对（）中19x 的系数作一个小的扰动。我们希望比较（）和（）根的差别，从而分析方程（）的解对扰动的敏感性。 实验内容：为了实现方便，我们先介绍两个Matlab 函数：“roots ”和“poly ”。 roots(a)u = 其中若变量a 存储n+1维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,,,+n a a a ，则输出u 的各分量是多项式方程 01121=+++++-n n n n a x a x a x a 的全部根；而函数 poly(v)b =
的输出b 是一个n+1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“poly ”是两个互逆的运算函数。 ;000000001.0=ess );21,1(zeros ve = ;)2(ess ve = ))20:1((ve poly roots + 上述简单的Matlab 程序便得到（）的全部根，程序中的“ess ”即是（）中的ε。 实验要求： （1）选择充分小的ess ，反复进行上述实验，记录结果的变化并分析它们。 如果扰动项的系数ε很小，我们自然感觉（）和（）的解应当相差很小。计算中你有什么出乎意料的发现表明有些解关于如此的扰动敏感性如何 （2）将方程（）中的扰动项改成18x ε或其它形式，实验中又有怎样的现象 出现 （3）（选作部分）请从理论上分析产生这一问题的根源。注意我们可以将 方程（）写成展开的形式， ) 3.1(0 ),(1920=+-= x x x p αα 同时将方程的解x 看成是系数α的函数，考察方程的某个解关于α的扰动是否敏感，与研究它关于α的导数的大小有何关系为什么你发现了什么现象，哪些根关于α的变化更敏感 思考题一：（上述实验的改进） 在上述实验中我们会发现用roots 函数求解多项式方程的精度不高，为此你可以考虑用符号函数solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。
1.1 术语 1.1.1 长期稳定性试验（long term stability testing） 在控制的环境条件下（接近实际条件）留样，并定期进行检验。其目的为了确 认适当的贮存条件及制定产品有效期或复验期提供依据。 1.1.2 加速稳定性试验（Accelerated stability testing） 通过使用超常的贮存条件来加速产品化学降解或物理变化。评估在运输期间短 期的贮存条件与标签上规定的条件不符时，可能会造成对产品质量的影响；是 对长期稳定性研究数据的补充。 1.1.3 影响因素试验 1.1.4 是在比加速试验更激烈的条件下进行。其目的是探讨药物的固有稳定性、了解 影响其稳定性的因素及可能的降解途径与降解产物，为制剂生产工艺、包装、 贮存条件和建立降解产物分析方法提供科学依据。 1.2 稳定性年度计划 1.2.1 QC经理根据生产品和QA的要求，种制订年度稳定性试验计划，包括：品名、 规格、批次、批号、试验方式、稳定性试验项目及方法、责任人等。 1.2.2 稳定性计划经QC经理审核，交质量部负责人批准方可生效。 1.3 原料药影响因素试验 1.3.1 原料药进行该试验。用一批代表性的原料药进行，将样品除去包装后，平放在 称量瓶或培养皿中摊成≤5 mm厚的薄层，疏松原料药摊成≤10 mm厚薄层， 进行以下实验，考察各项指标的变化情况。 1.3.2 高温试验 温度分别选择40℃、60℃两个温度水平。将样品在60℃温度下放置10天， 于第5天和第10天取样，按稳定性重点考察项目进行检测。若样品无明显变 化，则不进行40℃条件下的试验；若样品有明显变化(如含量下降5%，鉴别 不明显，外观色泽变化大等)，则须在40℃条件下同法进行试验。 1.3.3 高湿试验 选择相对湿度为75%±5%、90%±5%两个湿度水平（T=25℃）。 将样品置于相对湿度90%±5%条件下（装KNO3饱和溶液的干燥器中，用封口 胶密封）放置10天，于第5天和第10天取样，按稳定性重点考察项目要求 检测，同时准确称量试验前后样品的重量，若样品无明显变化，则不进行75%±5%
稳定性数据评价 1.介绍 1.1 指南的目的 该指南的目的是为了提供如何使用根据ICH指南Q1A（R）里详述的“新原料药和制剂稳定性试验”原则（以后提到即作为总指导原则）而产生的稳定性数据的介绍来建议再试验期或货架期。该指南描述了何时及如何使用有限外推法来建议关于原料药的再试验期或超出来自长期储存条件的数据的观测范围的原料药货架期。 1.2 背景 总指导原则提供的关于稳定性数据的评价和统计分析的指南是性质上简要和范围上有限制。尽管总指导原则指出回归分析是可接收的方法来分析关于再试验期或货架期评价的定量稳定性数据，并建议用0.25显著性水平操作合并批的统计测试，它很少包括细节。另外，总指导原则不包括当复合因素包含在全面或折合-设计调查的情况。当到该方针的第4步，总指导原则的评价部分将会重复，因此删去。 1.3 指南的范围 该指南，总指导原则的附件，目的是当基于定量和定性测试性质的稳定性数据评价而建议再试验期或货架期和贮存条件时提供预期值的清晰解释。该指南概括了基于单个或复合因素和全面或折合-设计调查得出的稳定性数据以确定再试验期或货架期的介绍。ICH Q6A 和Q6B提供了关于调整和证实认可标准的指南。 2. 指南 2.1 一般原则 正规稳定性调查的设计和实行应符合总指导原则列出的原则。稳
定性调查的目的是，在测试最少三批原料药或制剂基础上，确立适用于将来在相似环境下生产和包装批的再试验期或货架期和标签贮存说明。 在稳定性资料的说明和评价里应采用系统性方法，其中应包括，视情况而，从物理、化学、生物和微生物试验，包括从那些与剂型有关的特定性质（例如，固体口服剂型的溶解速率）的结果。如果合适，应注意回顾质量平衡的合适性。应该考虑能引起质量平衡明显不足的因素，例如，降解机理和稳定性-显示能力和分析方法内在可变性。单批的变化程度作用以后生产批次在其再试验期或货架期间仍保留在其认可标准内的信心。 该指南里关于统计法的介绍不意味着当统计计算被证明是多余时，用统计计算仍可取。但在一些情况下统计分析在再试验期或货架期的外推法里是有用的且在其它情况可能提倡将次用于核实再试验期或货架期。 稳定性数据测定的基本原则同于单个-与多个-因素调查和全面-与折合-设计调查。正规稳定性调查里的数据测定，并视情况而定，使用支持数据来确定可能作用原料药或制剂的质量和性能的关键质量性质。应各自评估每个性质和为了建议再试验期或货架期而由调查结果构成的全面评估。所提议的再试验期或货架期不应超过任何单个性质的预测。 附录A里提供的流程图和附录B里提供的关于如何分析和评价从多因素或折合设计得到的关于适当的定量试验性质的长期稳定性数据。用于数据分析的统计方法应该考虑稳定性调查为估计再试验期或货架期而提供有效统计结论。附录B也应该提供关于如何使用再试验
留样观察稳定性试验与持续稳定性考察的关系 新版GMP第225条明确了“留样”的概念：“企业按规定保存的、用于药品质量追溯或调查的物料、产品样品为留样。”并特别指出，“用于产品稳定性考察的样品不属于留样。”这说明企业在实际工作中容易将两者混淆。 对于“稳定性考察”，新版GMP专辟“持续稳定性考察”一节，详细规定了各项要求。在有关持续稳定性考察条款中，还有稳定性试验、长期稳定性试验等表述。留样与稳定性考察、稳定性试验等概念有必要进一步厘清，使其更具可操作性。 1998版GMP条款中，对于留样、稳定性考察等没有具体的要求，对其提出详细要求的是在新版GMP实施以后。也可以说，留样与稳定性考察的全面规范始于新版GMP。因此，企业在实施过程中，难免会出现一些认识与理解上的偏差。 留样功能只限于质量追溯与调查 留样的目的只限于物料与产品质量的追溯与调查。留样通常只是在以下情况时发挥作用：一是产品在流通、使用等环节，被监管部门依法抽样检测，发现质量问题后，与生产企业进行相关确认。此时，留样可以成为企业是否需要申请复检的依据之一。如果留样检测出现同样的质量缺陷，企业就应该转入相关的偏差调查与处理。 二是产品在流通、使用等环节，出现对病患者身体的疑似伤害甚至死亡，需要对生产企业及其涉及的产品进行追溯调查。 三是产品在流通、使用过程中出现质量争议时，留样可以作为相关证据，用于调查取证。 留样无论对企业还是监管者，或是对经营、使用环节，乃至使用该留样产品的病人，都具有非常重要的意义。任何环节出现的可能与产品质量相关的问题，都需要通过留样进行追溯与调查。在对与药品相关的违法犯罪行为实施惩处过程中，留样又具有重要的法律意义。在规定的保存期内，不允许灭失留样。否则，一旦发生法律问题，就会被认为是销毁证据。 鉴于留样的特定功能，在保存期限内，除了发生外观目测异常情况需要进行彻底调查外，其外包装完整性不应受到损坏。保存期内的留样，只进行一般的目检观察，不进行其它可能破坏原包装的任何检验检测。为了与需要进行产品内在质量检测的稳定性考察明确区分，留样与需要进行稳定性考察的样品应该分开存放，原则上不应保存在同一间库房内。留样观察记录与稳定性考察记录更应严格区分，不能混淆。 稳定性考察相关概念 在一些与药品质量管理相关的培训资料及专业书籍中，常出现与稳定性考察相关的概念，包括稳定性考察、持续稳定性考察、长期稳定性考察等。
9001 原料药物与制剂稳定性试验指导原则 稳定性试验的目的是考察原料药物或制剂在温度、湿度、光线的影响下随时间变化的规律，为药品的生产、包装、贮存、运输条件提供科学依据，同时通过试验建立药品的有效期。 稳定性试验的基本要求是：（1）稳定性试验包括影响因素试验、加速试验与长期试验。影响因素试验用1批原料药物或1批制剂进行；如果试验结果不明确，则应加试2个批次样品。生物制品应直接使用3个批次。加速试验与长期试验要求用3批供试品进行。（2）原料药物供试品应是一定规模生产的。供试品量相当于制剂稳定性试验所要求的批量，原料药物合成工艺路线、方法、步骤应与大生产一致。药物制剂供试品应是放大试验的产品，其处方与工艺应与大生产一致。药物制剂如片剂、胶囊剂，每批放大试验的规模，片剂至少应通常为100 000片，胶囊剂至少应为100 000粒。大体积包装的制剂如静脉输液等，每批放大规模的数量至少应为各项试验所需总量的10倍。特殊品种、特殊剂型所需数量，根据情况另定。（3）供试品的质量标准应与临床前研究及临床试验和规模生产所使用的供试品质量标准一致。（4）加速试验与长期试验所用供试品的包装应与上市产品一致。（5)研究药物稳定性，要采用专属性强、准确、精密、灵敏的药物分析方法与有关物质（含降解产物及其他变化所生成的产物）的检查方法，并对方法进行验证，以保证药物稳定性试验结果的可靠性。在稳定性试验中，应重视降解产物的检查。（6）由于放大试验比规模生产的数量要小，故申报者应承诺在获得批准后，从放大试验转入规模生产时，对最初通过生产验证的3 批规模生产的产品仍需进行加速试验与长期稳定性试验。（7）对包装在非渗透容器内的药物制剂可不考虑药物的湿敏感性或可能的溶剂损失，其稳定性研究可在任何湿度下进行。 制剂质量的“显著变化”通常定义为：①含量与初始值相差5%；或采用生物或免疫法测定时效价不符合规定；②任何降解产物超过标准限度要求；③外观、物理常数、功能试验（如颜色、相分离、再分散性、粘结、硬度、每揿剂量等）
L-腈化物稳定性试验报告 一、概述 L-腈化物是L-肉碱生产过程中的第一步中间体（第二步中间体：L-肉碱粗品；第三步中间体：L-肉碱潮品），由于L-肉碱生产工艺为间歇操作，即每生产一步中间体，生产完毕并出具合格检测报告后，存入中间体仓库，以备下一步生产投料所需。根据本公司L-肉碱产品的整个生产周期，L-腈化物入库后可能存放的最长时间为4周（约28天）。以此周期为时间依据制定T L-腈化物稳定性试验方案，用于验证L-腈化物在再试验期限内的各项质量指标数据的稳定性，并且能否符合L-腈化物的质量标准，此次稳定性试验的整个周期为28天，具体的稳定性试验方案以ICH药物稳定性指导原则为基础制定，以确保L-腈化化物稳定性试验的可操作性。 二、验证日期 2010 年1月13日----2010 年2月10日 三、验证方案 1）样品储存和包装： 考虑到L-腈化物今后的贮藏、使用过程，本次用于稳定性试验的样品批次与最终规模生产所用的L-腈化物的包装和放置条件相同。 2）样品批次选择：此次稳定性试验共抽取三批样品，且抽取样品的批次与 最终规模生产时的合成路线和生产工艺相同 3）抽样频率和日期：从2010.1.13起，每隔7天取样一次，共取五次，具体日期
为：2010.1.13、2010.1.20、2010.1.27、2010.2.3、2010.2.10，以确保试验 次数足以满足L-腈化物的稳定性试验的需要。。 4）检测项目：根据L-腈化物的质量标准的规定，此次稳定性试验的检测项目共五项，分别为外观、氯含量、熔点、比旋度、干燥失重。这些指标 在L-腈化物的储存过程中可能会发生变化，且有可能影响其质量和有效 性。 5）试样来源和抽样：L-腈化物由公司102车间生产，经检测合格后储存于中间体仓库，本次稳定性试验的L-腈化物均取自于该中间体仓库，其抽样方法和抽样量均按照L-腈化物抽样方案进行抽样。抽样完毕后直接进行检测分析，并对检测结果进行登记，保存，作为稳定性数据评估的依据。 四、稳定性试验数据变化趋势分析及评估 通过对三批L-腈化物的稳定性试验，对其物理、化学方面稳定性资料进行评价，旨在建立未来相似情况下，大规模生产出的L-腈化物是否适用现有的再试验期（28天）。批号间的变化程度是否会影响未来生产的L-腈化物在再试验期内是否仍符合其质量规范。本次试验数据以表格、图解的形式给出，从而对L-腈化物的稳定性数据进行有效的评估。
实验报告一 题目：非线性方程求解 摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。 前言：（目的和意义） 掌握二分法与Newton法的基本原理和应用。 数学原理： 对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。 对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b)<0，且f(x)在[a,b]内仅有一个实根x*，取区间中点c，若，则c恰为其根，否则根据f(a)f(c)<0是否成立判断根在区间[a,c]和[c,b]中的哪一个，从而得出新区间，仍称为[a,b]。重复运行计算，直至满足精度为止。这就是二分法的计算思想。
Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式 产生逼近解x*的迭代数列{x k}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为 其中r为要求的方程的根的重数，这就是改进的Newton法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。 程序设计： 本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x); y=-x*x-sin(x); 写成如上形式即可，下面给出主程序。 二分法源程序： clear %%%给定求解区间 b=1.5; a=0;
%%%误差 R=1; k=0;%迭代次数初值 while (R>5e-6) ; c=(a+b)/2; if f12(a)*f12(c)>0; a=c; else b=c; end R=b-a;%求出误差 k=k+1; end x=c%给出解 Newton法及改进的Newton法源程序：clear %%%% 输入函数 f=input('请输入需要求解函数>>','s') %%%求解f(x)的导数 df=diff(f);
***药物稳定性试验方案 一、试验目的 ***药物稳定性试验的目的是考察***药物在温度、湿度、光线的影响下随时间变化的规律，考察其在影响因素试验、加速试验的条件下各项指标是否符合***药物产品质量标准的要求。 二、试验要求 ***药物稳定性试验的基本要求有以下几个方面： 1、稳定性试验项目包括影响因素试验与加速试验。 2、***药物供试品可以从橡胶膏生产车间生产的同一批制剂中抽取，抽样数量符合稳定性试验要求。 3、加速试验与长期试验***药物所用供试品的容器和包装材料及包装方式应与上市产品一致。 4、***药物稳定性试验，要采用专属性强、准确、精密、灵敏的药物分析方法，以保证药物稳定性结果的可靠性。 三、试验方法 (一)影响因素试验 此项试验是在比加速试验更激烈的条件下进行。其目的是考察***药物的固有稳定性，了解其在高温、高湿及低温条件下各项质量指标的稳定性及变化情况。 1．高温试验 ***药物置药物稳定性检查仪中，60℃温度下放置10天，于第5天和第10天取样，按***药物成品质量标准进行全检。若***药物供试品有明显变化(如含量下降5％)，则在40℃条件下同法进行试验。若60℃无明显变化，不再进行40℃试验。 2．高湿度试验 ***药物置恒湿密闭容器中，在25℃分别于相对湿度90％±5％条件下放置10天，于第5天和第10天取样，按***药物成品质量标准进行全检，同时准确称量试验前后供试品的重量，以考察供试品的吸湿潮解性能。若吸湿增重5％以上，则在相对湿度75％±5％条件下，同法进行试验；若吸湿增重5％以下，且其他考察项目符合要求，则不再进行此项试验。恒湿条件可通过在密闭容器如干燥器下部
药物稳定性试验统计分析方法 在确定有效期的统计分析过程中，一般选择可以定量的指标进行处理，通常根据药物含 量变化计算，按照长期试验测定数值，以标示量％对时间进行直线回归，获得回归方程，求 出各时间点标示量的计算值（y'）,然后计算标示量（y'）95％单侧可信限的置信区间为y'±z，其中： zt N 2 1(X X) S 2() 2 NXiX （12-21） 式中，t N-2—概率0.05，自由度N-2的t单侧分布值（见表12-4），N为数组；X0—给定自变量；X—自变量X的平均值； Q S（12-22） N2 式中， 22 Q；L yy—y的离差平方和，LyyN L yy bLyy()/ xy ；L xy—xy的离差乘 积之和L xy xy(x)(y)/N；b—直线斜率。 将有关点连接可得出分布于回归线两侧的曲线。取质量标准中规定的含量低限（根据各品种实际规定限度确定）与置信区间下界线相交点对应的时间，即为药物的有效期。根据情况也可拟合为二次或三次方程或对数函数方程。 此种方式确定的药物有效期，在药物标签及说明书中均指明什么温度下保存，不得使用“室温”之类的名词。 例：某药物在温度25±2℃，相对温度60±10%的条件下进行长期实验，得各时间的标示量如表12-4。 表12-4供试品各时间的标示量 时间/月03691218 标示量/%99.397.697.398.496.094.0 以时间为自变量（x），标示量%（y）为因变量进行回归，得回归方程y=99.18－0.26x，r=0.8970，查T单侧分布表，当自由度为4，P=0.05得 tN-2=2.132 S Q N2 3.444 4 0.9279