西安交大计算方法考试

《计算方法B》上机

实验报告

学院：班级: 姓名：学号：

机械工程学院

2015年12月22日1

1.计算以下和式：S??1nn?016?211??4????8n?18n?48n?58n?6?，要求： ??(1)若保留11个有效数字，给出计算结果，并评价计算的算法；

(2)若要保留30个有效数字，则又将如何进行计算。

实现思想：

以上问题出现了近似数相减的问题，为了减小误差，可分别求得减数之和以及被减数之和，最后将两者相减。另外，减数与被减数求和均为同号计算，按照绝对值递增顺序相加可减小舍入误差。此题中对有效数字有要求，因而计算时首先需要根据有效数字位数计算得出迭代次数，以保证计算值的精度。

源程序：

m=input('输入有效数字个数m=');

s0=1;s1=0;s2=0;n=0; %判断迭代次数

while s0>=0.5*10^-(m-1)

s0=4/(16^n*(8*n+1))-2/(16^n*(8*n+4))-1/(16^n*(8*n+5))-1/(16^n*(8*n+6)); n=n+1; end

%分别求解各项并求和 for k=n-1:-1:0

a1=4

《计算方法B》上机

实验报告

学院：班级: 姓名：学号：

机械工程学院

2015年12月22日1

1.计算以下和式：S??1nn?016?211??4????8n?18n?48n?58n?6?，要求： ??(1)若保留11个有效数字，给出计算结果，并评价计算的算法；

(2)若要保留30个有效数字，则又将如何进行计算。

实现思想：

以上问题出现了近似数相减的问题，为了减小误差，可分别求得减数之和以及被减数之和，最后将两者相减。另外，减数与被减数求和均为同号计算，按照绝对值递增顺序相加可减小舍入误差。此题中对有效数字有要求，因而计算时首先需要根据有效数字位数计算得出迭代次数，以保证计算值的精度。

源程序：

m=input('输入有效数字个数m=');

s0=1;s1=0;s2=0;n=0; %判断迭代次数

while s0>=0.5*10^-(m-1)

s0=4/(16^n*(8*n+1))-2/(16^n*(8*n+4))-1/(16^n*(8*n+5))-1/(16^n*(8*n+6)); n=n+1; end

%分别求解各项并求和 for k=n-1:-1:0

a1=4

计算方法（C）

目 录

第1章 绪论

1.1 数值计算 1.2 数值方法的分析

1.2.1 计算机上数的运算 1.2.2 算法分析 第2章 线性代数方程组 2.1 Gauss消去法

2.1.1 消去法 2.1.2 主元消去法 2.2 矩阵分解

2.2.1 Gauss消去法的矩阵意义 2.2.2 矩阵的LU分解及其应用 2.2.3 其他类型矩阵的分解 2.2.4 解三对角矩阵的追赶法 2.3

线性方程组解的可靠性 2.3.1 向量和矩阵范数

2.3.2 残向量与误差的代数表征 2.4

解线性方程组解的迭代法 2.4.1 基本迭代法 2.4.2 迭代法的矩阵表示 2.4.3 收敛性

第3章 数据近似

3.1 多项式插值 3.1.1 插值多项式

3.1.2 Lagrange插值多项式 3.1.3 Newton插值多项式 3.1.4 带导数条件的插值多项式 3.1.5 插值公式的余项 3. 2 最小二乘近似

3.2.1 最小二乘问题的法方程 3.2.2 正交化算法 第4章 数值微积分

1

4.1 内插求积，Newton-Cotes公式 4.1.1 New

计算方法上机作业

1.对以下和式计算：S??1?4211?，要求： ????n?8n?48n?58n?6?n?016?8n?1?(1)若只需保留11个有效数字，该如何进行计算； (2)若要保留30个有效数字，则又将如何进行计算；

（1）解题思想和算法实现：

根据保留有效位数的要求，可以由公式

得出计算精度要求。只需要很

少内存，时间复杂度和d呈线性，不需要高浮点支持。先根据while语句求出符合精度要求的n值的大小，然后利用for语句对这n项进行求和，输出计算结果及n值大小即可。

（2）matlab源程序：

保留11位有效数字时； clear clc

format long n=0;

sum=1/(16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6)); while sum>=5*10^(-11); n=n+1;

sum=1/(16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6)); end

fori=0:n-1;

sum=sum+1/(16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)); end

vpa(sum,11) n

清洁验证残留限度的计算

根据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

计算机程序设计基础

（C++）

实 验 报 告

信息学院计算机教育中心

二零一四年三月

实验_0验报告

教学班级：___ 7班__学生学号：__20133936_学生姓名：____葛劲松 实验日期：_____2014.3.11______实验地点：______7307___(机房) 指导教师签名：___王威_______ 实验成绩：___________

一、实验目的

1．熟悉VS2010 集成开发环境；

2．学习新建和打开控制台应用程序的方法； 3．通过运行示例程序了解程序开发过程；

二、实验任务

1．新建源程序文件，编辑以下程序，编译，链接，得出运行结果，并保存程序： #include using namespace std; void main( )

{ cout<<'a'<

2．新建源程序文件，编辑以下程序，编译，链接，运行，输入36和64，得出运行结果，并保存程序：

#include using namespace std; intmain()

{ float num1,num2,sum;

cout<< \ cin>> num1;

cout<< \ cin>> num2;

sum=num1+num2;

cout<< \ return 0; }

三、实验结果

1.

哈密南-郑州±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

开始时间 2014年12月7日 星期日 20:57 完成于 耗时 分数 成绩

2014年12月7日 星期日 22:00 1 小时 3 分钟 39.00/40.00 97.50/满分100.00

窗体顶端 题目1 正确

获得1.00分中的1.00分

Flag question

题干

以下格式中，不属于视频文件格式的是（ ）。 选择一项：

a. WMV格式 b. MPG格式 c. MOV格式 d. MID格式

反馈

正确答案是：MID格式 题目2 正确

获得1.00分中的1.00分

Flag question

题干

以下对视频格式文件的描述中，不正确的是（ ）。 选择一项：

a. .ra、.rm 或 .rmvb 格式是Windows公司开发的一种静态图像文件格式

b. .ra、.rm 或 .rmvb 格式可以根据网络数据传输速率的不同而采用不同的压缩比率 c. .ra、.rm 或 .rmvb 格式可实现影像数据的实时传送和实时播放

d. .ra、.rm 或 .rmvb 格式是流式视频文件格式，用来在低速率的广域网上实时传输活动视频影像 反馈

正确答案是：.ra、.rm 或 .rmvb 格式是Windows公司开发的一种静态图像文

《数值计算方法》

邹昌文编

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为