核酸报告计算方法

西安交通大学

计算方法B上机报告

班级：XXXXXXXX 姓名：XXX

学号：XXXXXXXX

2015/12/13 Sunday

目录

目录

1 题目一 ............................................................................................................................... 1 1. 数值计算........................................................................................................................ 1 2. 实现思想........................................................................................................................ 1 3. 源程序........................................................

《计算方法》实验报告

学号 实验项目名称 一、实验名称 计算方法实验 姓名 班级 实验一 插值与拟合 二、实验目的： (1)明确插值多项式和分段插值多项式各自的优缺点； (2)编程实现拉格朗日插值算法，分析实验结果体会高次插值产生的龙格现象； (3)运用牛顿插值方法解决数学问题。 三、实验内容及要求 1(1) 对于f(x)?,?5?x?5 21?x要求选取11个等距插值节点，分别采用拉格朗日插值和分段线性插值，计算x为0.5, 4.5处的函数值并将结果与精确值进行比较。 输入：区间长度，n(即n+1个节点)，预测点 输出：预测点的近似函数值，精确值，及误差 （2）已知1?1，4?2，9?3，用牛顿插值公式求5的近似值。 输入：数据点集，预测点。 输出：预测点的近似函数值 四、实验原理及算法描述 算法基本原理： (1)拉格朗日插值法

(2) 牛顿插值法 算法流程 五、程序代码及实验结果 （1） 输出： A．拉格朗日插值法 B.分段线性插值 X y(精确) y(拉格朗日) y(分段线性) 误差(拉) 误差(分) 0.500000 0.800000

西安交通大学

计算方法B上机报告

班级：XXXXXXXX 姓名：XXX

学号：XXXXXXXX

2015/12/13 Sunday

目录

目录

1 题目一 ............................................................................................................................... 1 1. 数值计算........................................................................................................................ 1 2. 实现思想........................................................................................................................ 1 3. 源程序........................................................

《计算方法》实验报告

学号 实验项目名称 一、实验名称 计算方法实验 姓名 班级 实验一 插值与拟合 二、实验目的： (1)明确插值多项式和分段插值多项式各自的优缺点； (2)编程实现拉格朗日插值算法，分析实验结果体会高次插值产生的龙格现象； (3)运用牛顿插值方法解决数学问题。 三、实验内容及要求 1(1) 对于f(x)?,?5?x?5 21?x要求选取11个等距插值节点，分别采用拉格朗日插值和分段线性插值，计算x为0.5, 4.5处的函数值并将结果与精确值进行比较。 输入：区间长度，n(即n+1个节点)，预测点 输出：预测点的近似函数值，精确值，及误差 （2）已知1?1，4?2，9?3，用牛顿插值公式求5的近似值。 输入：数据点集，预测点。 输出：预测点的近似函数值 四、实验原理及算法描述 算法基本原理： (1)拉格朗日插值法

(2) 牛顿插值法 算法流程 五、程序代码及实验结果 （1） 输出： A．拉格朗日插值法 B.分段线性插值 X y(精确) y(拉格朗日) y(分段线性) 误差(拉) 误差(分) 0.500000 0.800000

江西科技师范大学 计算方法 实验报告

江 西 科 技 师 范 学 院

实 验 报 告

课 程

系 别

班 级

学 号

姓 名

江西科技师范大学 计算方法 实验报告

目录

实验一 误差的传播与估计…………………………………… 实验二 拉格朗日插值多项式………………………………… 实验三 变步长复合梯形求积公式…………………………… 实验四 解非线性方程二分法………………………………… 实验五 一元非线性方程的迭代解法………………………… 实验六 列主元高斯消去法……………………………………

每次实验课必须带上此本子，以便教师检查预习情况和记录实验原始数据。实验时必须遵守实验规则。用正确的理论指导实践袁必须人人亲自动手实验，但反对盲目乱动，更不能无故损坏仪器设备。这是一份重要的不可多得的自我学习资料袁它将记录着你在大学生涯中的学习和学习成果。请你保留下来，若干年后再翻阅仍将感到十分新鲜，记忆犹新。它将推动你在人生奋斗的道路上永往直前！

江西科技师范大学 计算方法 实验报告

实验一 误差的传播与估计

一、 实验课程名称 数学实验

二、 试验项目名称 误差的传播与估计 三、 实验目的和要求

理解误差在算术运算中的传播方式及如何在算术运算中控制误差的传播。

四、

清洁验证残留限度的计算

根据GMP实施指南和相关要求，我们控制原料药（乙酰螺旋霉素）残留限度的计算依据如下：

计算方法：10ppm法、日剂量的千分之一、下批批量的0.1%（基于低毒性原料的杂质限度标准）

1、10ppm法：乙酰螺旋霉素批量为260kg，因残留物浓度最高为10*10-6，即10mg/kg，则残留物总量最大为：260*10*10-6=2600mg。则设备内表面残留物允许的限度为：

2600g?1000?100cm2?10%（保险系数）?70%（取样回收率） 残留限量A? 289.7m?10000＝20.31㎎/100㎝2

残留限度定为：20.31㎎/100㎝2/25ml=0.8124mg/ml

2、日剂量的千分之一：由于原料药生产清洁后用于生产药用辅料（醋酸钠），其为无活性物质，因此暂无法用此公式计算。

3、下批批量的0.1%（基于低毒性原料的杂质限度标准）

原料药（乙酰螺旋霉素）的最小批产量为260㎏，下批批量的0.1%，则乙酰螺旋霉素最大残留物为260g。

擦拭测试：擦拭面积以10㎝×10㎝的区域计 残留限量A?260g?1000?100cm2?10%（保险系数）?70%（取样回收率） 289.7m?10

哈密南-郑州±800kV特高压直流输电线路工程 晋1标段施工项目部 1.放线牵张力计算

（1）模拟放线弧垂，选取控制档、放线模板K值。 （2）计算控制档水平张力： Tn? 式中：

w2 2KTn——控制档水平张力，t ；

w2——导线单位重量，t ； K——模板K值。 （3）计算张力机出口张力：

ε(εn0?1) T0?n[Tn?w2?h0]

0n0(ε?1)ε1 式中：

T0——张力机出口张力，t ；

n——放线段内滑车数；

n0——张力场与控制档间滑车数；

ε——滑车摩擦系数；

?h0——控制档与张力场累计高差，m，控高为“+”。

（4）计算初始牵引力：

ε(εn?1) p0?k0[NTε?w1?h] 0n(ε?1)n

哈密南-郑州±800kV特高压直流输电线路工程 晋1标段施工项目部 1.放线牵张力计算

（1）模拟放线弧垂，选取控制档、放线模板K值。 （2）计算控制档水平张力： Tn? 式中：

w2 2KTn——控制档水平张力，t ；

w2——导线单位重量，t ； K——模板K值。 （3）计算张力机出口张力：

ε(εn0?1) T0?n[Tn?w2?h0]

0n0(ε?1)ε1 式中：

T0——张力机出口张力，t ；

n——放线段内滑车数；

n0——张力场与控制档间滑车数；

ε——滑车摩擦系数；

?h0——控制档与张力场累计高差，m，控高为“+”。

（4）计算初始牵引力：

ε(εn?1) p0?k0[NTε?w1?h] 0n(ε?1)n

《数值计算方法》

邹昌文编

2009年10月

上机实验指导书

“数值计算方法”上机实验指导书

实验一 误差分析

实验1.1（病态问题）

实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。

数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。

问题提出：考虑一个高次的代数多项式

p(x) (x 1)(x 2) (x 20) (x k)

k 120

(1.1)

显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动

p(x) x19 0

(1.2)

其中 是一个非常小的数。这相当于是对（1.1）中x19的系数作一个小的扰动。我们希望比较（1.1）和（1.2）根的差别，从而分析方程（1.1）的解对扰动的敏感性。

实验内容：为了实现方便，我们先介绍两个MATLAB函数：“roots”和“poly”。

u roots(a)

其中若变量a存储n+1维的向量，则该函数的输出u为

重大事故后果分析方法：泄漏

事故后果分析是安全评价的一个重要组成部分，其目的在于定量地描述一个可能发生的重大事故对工厂、厂内职工、厂外居民，甚至对环境造成危害的严重程度。分析结果为企业或企业主管部门提供关于重大事故后果的信息，为企业决策者和设计者提供关于决策采取何种防护措施的信息，如防火系统、报警系统或减压系统等的信息，以达到减轻事故影响的目的。火灾、爆炸、中毒是常见的重大事故，可能造成严重的人员伤亡和巨大的财产损失，影响社会安定。世界银行国际信贷公司(IFC)编写的《工业污染事故评价技术手册》中提出的易燃、易爆、有毒物质的泄漏、扩散、火灾、爆炸、中毒等重大工业事故的事故模型和计算事故后果严重度的公式，主要用于工业污染事故的评价。该方法涉及内容，也可用于火灾、爆炸、毒物泄漏中毒等重大事故的事故危险、危害程度的评价。

由于设备损坏或操作失误引起泄漏从而大量释放易燃、易爆、有毒有害物质，可能会导致火灾、爆炸、中毒等重大事故发生。 1 泄漏情况

1．1 泄漏的主要设备

根据各种设备泄漏情况分析，可将工厂(特别是化工厂)中易发生泄漏的设备分类，通常归纳为：管道、挠性连接器、过滤器、阀门、压力容器或反应器、泵、压缩机、储罐、