计算方法实验11
数值计算方法LU分解法实验
数值计算方法LU分解法实验LU分解法是一种常见的数值计算方法，用于解线性方程组或求解矩阵的逆。
该方法的核心思想是将一个矩阵分解为一个下三角矩阵和一个上三角矩阵的乘积。
在本篇文章中，我们将进行关于LU分解法的实验，并探讨其性能和应用。
首先，我们需要明确LU分解法的数学原理。
假设我们有一个n阶方阵A，LU分解法的目标是找到两个矩阵L和U，使得A=LU。
其中L是一个下三角矩阵，U是一个上三角矩阵。
然后，我们可以将原始的线性方程组Ax=b转化为两个新的方程组Ly=b和Ux=y。
通过求解这两个方程组，我们可以得到原始方程组的解。
接下来，我们将通过一个具体的例子来说明LU分解法的步骤和计算过程。
假设我们有以下方程组：2x+y+z=8-3x-y+2z=-11-2x+y+2z=-3首先，我们将系数矩阵A进行LU分解。
在这个例子中，我们可以得到下三角矩阵L和上三角矩阵U：L=100-1.510-1-11U=21101.52.5001然后，我们将方程组转化为Ly=b和Ux=y的形式。
解这两个方程组，可以得到y和x的值。
最终，我们可以得到方程组的解为x=2，y=3，z=-1通过以上的实例，我们可以看到LU分解法的步骤较为繁琐，但是它的结果是准确的。
那么，接下来我们将进行一系列实验，来评估LU分解法的性能和应用。
首先，我们将进行LU分解法的准确性测试。
我们将随机生成一组方程组，并使用LU分解法求解出它们的解。
然后，我们将使用该解验证原方程组，并计算出其误差。
我们重复这个过程多次，并计算平均误差。
通过这次实验，我们可以判断LU分解法的准确性。
其次，我们将评估LU分解法的计算效率。
我们将随机生成不同规模的方程组，并使用LU分解法求解它们。
然后，我们记录下求解所需的时间，并绘制出问题规模和求解时间的关系图。
通过这个实验，我们可以了解LU分解法在不同规模问题上的计算效率，从而评估其可行性和应用范围。
此外，我们还可以将LU分解法与其他数值计算方法进行比较。
实验11乙酸乙酯皂化反应速率常数的测定
1780 lg K 0.00754 T 4.54 T
乙酸乙酯皂化反应的活化能: Ea=27.3KJ/mol
实验结果与讨论
⑴结果：实测值为Ｅa= ⑵计算实验偏差： ⑶分析产生偏差的原因： ⑷有何建议与想法？
注意事项：
1.电导率仪要进行温度补偿及常数校正； 2.反应液在恒温时都要用橡胶塞子盖好； 3.混合过程既要快速，又要小心谨慎，不 要把溶液挤出反应管； 4.严格控制恒温的温度，因为反应过程温 度对反应速率常数影响很大； 5.严格配制溶液的浓度，保证氢氧化钠与 乙酸乙酯的浓度相等，否则反应速率常 数计算公式将发生变化。
积分得：
1 1 Kt c0 x c0
dx K (C0 x) 2 dt
-------------（1）
本实验采用电导法测量皂化反应中电 导L随时间t的变化。设L0、Lt、L∞分别代 表时间为0、t、∞时溶液的电导。在稀溶 液下，乙酸乙醋皂化反应，在不同反应时 刻，其电导率与浓度的满足以下关系：
实验原理
乙酸乙酯皂化反应是个典型的二级反 应。设反应物起始浓度均为C0，经时间t 后产物的浓度为x。
CH3COOC2 H5 NaOH CH3COONa C2 H5OH t=0 C0 C0 0 0 t=t C0-x C0 -x x x t=∞ 0 0 C0 C0
该反应的速率方程为：
0.924541.168103 t 1.95106 t 2
式中：密度ρ的单位为g/ml; 温度t的单位为℃。
用移液管量取20ml0.02mol/LNaOH溶液和 20ml蒸馏水放于100ml烧杯中，混均后倒 入大试管中（盖上橡皮塞）。 在反应管的小管中加入20ml0.02mol/L乙 酸乙醋溶液，大管中加入20ml0.02mol/L NaOH溶液（均盖上橡皮塞）。 把它们同时放入恒温槽中，恒温15分钟。
计算方法-插值方法实验
实验一插值方法一. 实验目的（1）熟悉数值插值方法的基本思想，解决某些实际插值问题，加深对数值插值方法的理解。
（2）熟悉Matlab 编程环境，利用Matlab 实现具体的插值算法，并进行可视化显示。
二. 实验要求用Matlab 软件实现Lagrange 插值、分段线性插值、三次Hermite 插值、Aitken 逐步插值算法，并用实例在计算机上计算和作图。
三. 实验内容1. 实验题目 （1）已知概率积分dxe y xx ⎰-=22π的数据表构造适合该数据表的一次、二次和三次Lagrange 插值公式，输出公式及其图形，并计算x =0.472时的积分值。
答：①一次插值公式：输入下面内容就可以得到一次插值结果 >> X=[0.47,0.48];Y=[0.4937452,0.5027498]; >> x=0.472;>> (x-X(2))/(X(1)-X(2))*Y(1)+(x-X(1))/(X(2)-X(1))*Y(2)ans =0.495546120000000>>②两次插值公式为：输入下面内容就可以得到两次插值结果>> X=[0.46,0.47,0.48];Y=[0.4846555,0.4937452,0.5027498]; >> x=0.472;>>(x-X(2))*(x-X(3))/((X(1)-X(2))*(X(1)-X(3)))*Y(1)+(x-X(1))*(x-X(3))/((X(2)-X(1))*(X(2)-X(3)))*Y(2)+(x-X(2))*(x-X(1))/((X(3)-X(2))*(X(3)-X(1)))*Y(3)i 0123x 0.46 047 0.48 0.49 y0.4846555 0.4937452 0.5027498 0.5116683ans =0.495552928000000>>③三次插值公式为：输入下面内容就可以得到三次插值结果>> X=[0.46,0.47,0.48,0.49];Y=[0.4846555,0.4937452,0.5027498,0.5116683];>> x=0.472;>>(x-X(2))*(x-X(3))*(x-X(4))/((X(1)-X(4))*(X(1)-X(2))*(X(1)-X(3)))*Y(1)+(x-X(4))*( x-X(1))*(x-X(3))/((X(2)-X(4))*(X(2)-X(1))*(X(2)-X(3)))*Y(2)+(x-X(4))*(x-X(2))*( x-X(1))/((X(3)-X(4))*(X(3)-X(2))*(X(3)-X(1)))*Y(3)+(x-X(3))*(x-X(2))*(x-X(1))/(( X(4)-X(1))*(X(4)-X(2))*(X(4)-X(3)))*Y(4)ans =0.495552960000000输入下面内容，绘出三点插值的图：>> X=[0.46,0.47,0.48,0.49];Y=[0.4846555,0.4937452,0.5027498,0.5116683];>> x=linspace(0.46,0.49);>>y=(x-X(2)).*(x-X(3)).*(x-X(4))/((X(1)-X(4))*(X(1)-X(2))*(X(1)-X(3)))*Y(1)+(x-X(4) ).*(x-X(1)).*(x-X(3))/((X(2)-X(4))*(X(2)-X(1))*(X(2)-X(3)))*Y(2)+(x-X(4)).*(x-X(2) ).*(x-X(1))/((X(3)-X(4))*(X(3)-X(2))*(X(3)-X(1)))*Y(3)+(x-X(3)).*(x-X(2)).*(x-X(1) )/((X(4)-X(1))*(X(4)-X(2))*(X(4)-X(3)))*Y(4);>>plot(x,y)（注意上面的“.*”不能用“*”替代）；（2）将区间[-5,5]分为10等份，求作211)(x x f +=的分段线性插值函数，输出函数表达式及其图形，并计算x =3.3152时的函数值。
实验十一蔗糖转化速率常数的测定
如能在一定的时间间隔t 内测得一系列数据，因为t 为定值，所以 ln(c - c ¢) 对 t 作图， 即可由直线的斜率求出 k 值。 这个方法的困难是必须使t 为一定值，而这通常不易直接求得，需要从 c ~ t 图上求 得，因而又多一步计算手续。
六、旋光仪的使用
1．旋光仪的构造及原理 光是电磁波，自然光在传播方向垂直的各个方向上都有振动，称为普通光。当光在 某个方向的振动占优势时，这种光称为偏振光。某些有机化合物和晶体，因为他们的分 子结构和晶体结构的特殊性，当有偏振光通过时，可使偏振光的振动平面向左或向右偏 转一定角度。这种使偏振光的振动面发生偏转的化合物称为旋光物质或光学活性物质， 这是因为分子内部有手性中心或晶体的各向异性造成的。旋光物质所具有的这种使偏振 面偏转的能力称为旋光性。一定温度下，偏转的角度和方向与所测物质的本性、溶液的 浓度、液层的厚度及偏振光的波长等因素有关。利用旋光性，可对旋光物质进行定性和
1
2
3
4
5
6
图 11-2 旋光仪的光学系统 1.光源 2.起偏棱镜 3.石英片 4.样品管 5. 检偏镜 6.目镜 Figure11-2 Oource 2. polarizer 3. quartz piece 4.sample compartment 5.analyzer 6.detector
四、实验内容
1．将恒温槽与旋光管保温套相接，将超级恒温水浴的温度调节为 25℃，开通超级 恒温水浴的循环水。用受皿天平称取约 10 g 蔗糖于碘量瓶内，并加蒸馏水 50 mL 使蔗糖 完全溶解。若溶液浑浊，则需过滤。用移液管移取蔗糖溶液 25.00 mL 放入一个清洁干燥 的碘量瓶中，在另一个清洁干燥的碘量瓶内移入 25.00 mL 0.5 mol/L 的 HCl 溶液。将这 两只碘量瓶一起放入恒温槽中恒温 5 min ~10 min。 2．将恒温槽中 25mL HCl 溶液迅速倾入 25mL 蔗糖溶液中，立即按下秒表，作为反 应起点 t0。再将溶液倒回盛 HCl 溶液的碘量瓶中，如此反复三次，摇匀。迅速取少许混 合溶液荡洗旋光管 3 ~ 4 次，然后将此混合溶液装入旋光管，盛满，盖好管盖（勿使管 内有气泡） ，擦净，立刻放进旋光仪内。测量各时间 t 时溶液的旋光度 t，测定时要迅速、 准确。剩余混合反应液连同碘量瓶一起放入 65℃ 左右的高温水浴内加速反应至水解完 全。 3．测量时间间隔大致如下：从 t0 开始每隔 5、5、10、10、20、20、30 min 各测一 次。 4．t 测定完毕后，取出高温水浴内的反应液冷却至实验温度，装入旋光管，测其旋 光度的值。 6．实验结束后，必须洗净旋光管和碘量瓶，并将碘量瓶放入烘箱烘干。
实验11 SGPT测定 2010
各管摇匀，37℃水浴箱预热10分钟 各管摇匀，37℃水浴箱预热10分钟 各管摇匀， ℃水浴箱加热20分钟 各管摇匀，37℃水浴箱加热 分钟 充分摇匀， 分钟内， 充分摇匀，在10分钟内，520nm处比色 分钟内 处比色 0.0 0 0.1 100 0.2 200 0.3 300 0.4 400 0.5 500
【实验方法】 实验方法】
（一）制作标准曲线
（A）方法一
管号 试剂 0.1M磷酸盐缓冲液（ml） 磷酸盐缓冲液（ ） 磷酸盐缓冲液 SGPT底物溶液（ml） 底物溶液（ ） 底物溶液 2µmol/ml丙酮酸标准液（ml） 丙酮酸标准（ ） 丙酮酸标准液 2, 4－二硝基苯肼溶液 － 0.4M NaOH (ml) 相当于丙酮酸含量（ 相当于丙酮酸含量（µmol） ） 相当于酶活力单位/100ml血清 血清 相当于酶活力单位 A520 A测定管—A空白管 0.10 0.50 0.00 0.5 5.0 0.10 0.45 0.05 0.5 5.0 0.10 0.40 0.10 0.5 5.0 0.10 0.35 0.15 0.5 5.0 0.10 0.30 0.20 0.5 5.0 0.10 0.25 0.25 0.5 5.0 0 1 2 3 4 5
【思考题】 思考题】
1. 什么叫转氨基作用? 在蛋白质代谢中有何重要作 什么叫转氨基作用 用? 2. 为什么制作标准曲线时,需要加入一定量SGPT底 为什么制作标准曲线时,需要加入一定量SGPT底 需要加入一定量 物溶液? 物溶液 3. 为什么测定酶活力时需要有对照 为什么测定酶活力时需要有对照? 4. 两种标准曲线测定方法在设计上有什么优缺点? 两种标准曲线测定方法在设计上有什么优缺点 哪一种更为合理?为什么 哪一种更为合理 为什么? 为什么
FTM11-残余粘着率试验计算方法
FTM11 残余粘着率试验计算方法
残余粘着率强度通常由百分比表示，是将测试样品的剥离力比上控制带的剥离力，再取平均值。
若两个比值不相近，可能的话，需要尽快安排重复实验。
若无法重新实验，需将两个比值同时记录在报告中。
残余粘着率％＝测试胶带的剥离力/控制带的剥离力X100
若测试结果超过100%，简单的将结果记录为≥100%覆盖性较差的涂硅和一些渗透性的基材也会对数值有一定的影响，那样的话，结果之间的比较应该限制在相同的基材之间（思泰仪器）。
实验11 粘度法测定高聚物的摩尔质量
实验十三粘度法测定高聚物的摩尔质量【目的要求】1．了解粘度法测定高聚物摩尔质量的基本原理和方法；2．掌握用乌氏(Ubbelohde)粘度计测定高聚物溶液粘度的原理和方法；3．测定右旋糖苷的摩尔质量。
【实验原理】高聚物是由单体分子经加聚或缩聚过程得到的。
在高聚物中，由于聚合度的不同，每个高聚物分子的大小并非都相同，致使高聚物的分子质量大小不一，参差不齐，且没有一个确定的值。
因此，高聚物的摩尔质量是一个统计平均值。
高聚物摩尔质量不仅反映了高聚物分子的大小，而且直接关系到它的物理性能，是一个重要的基本参数。
测定高聚物摩尔质量的方法很多，例如渗透压、光散射及超离心沉降平衡等方法。
但是不同方法所得平均摩尔质量也有所不同，比较起来，粘度法设备简单，操作方便，并有很好的实验精度，是常用的方法之一。
用此法求得的摩尔质量称为粘均摩尔质量。
粘度是液体流动时内摩擦力大小的反映。
高聚物溶液的特点是粘度特别大，原因在于其分子链长度远大于溶剂分子，加上溶剂化作用，使其在流动时受到较大的内摩擦力，粘性液体在流动过程中所受阻力的大小可用粘度系数（简称粘度）来表示。
纯溶剂粘度反映了溶剂分子间的内摩擦力，高聚物溶液的粘度则是高聚物分子间的内摩擦力、高聚物分子与溶剂分子间的内摩擦力及溶剂分子间内摩擦力三者之和。
在相同温度下，通常高聚物溶液的粘度大于纯溶剂粘度，即：﹥。
为了比较这两种粘度，引入增比粘度的概念，以表示（13.1）式中，称为相对粘度，定义为溶液粘度与纯溶剂粘度的比值，即（13.2）反映的也是粘度行为，而则表示已扣除了溶剂分子间的内摩擦效应。
高聚物的增比粘度往往随质量浓度的增加而增加。
为了便于比较，将单位浓度所显示的增比粘度称为比浓粘度，而称为比浓对数粘度。
当溶液无限稀释时，高聚物分子彼此相隔甚远，它们之间的相互作用可以忽略，此时有关系式图II-13-1 和图1 ；212（13.3）式中，称为特性粘度，它反映的是高分子与溶剂分子之间的内摩擦，其数值取决于溶剂的性质以及高聚物分子的大小和形态。
实验十一 配位化合物的组成和稳定常数的测定
配位化合物的组成和稳定常数的测定
随着科学技术的发展，配合物在科学研究和生产实践中显示出越来越重要的意义， 配合物不仅在化学领域里得到广泛的应用，并且对生命现象也具有重要的意义。人体内 各种酶（生物催化剂）的分子几乎都含有以配合状态存在的金属元素。配位化学与有 机、分析等化学领域以及生物化学、药物化学、化学工业都有着密切的关系，应 用非常广泛。 配位化学的研究有一整套的物理实验方法，如测定稳定常数可以采用 pH 法、电位 法、极谱法、分光光度法、溶剂萃取法、离子交换法等，这些方法对于配合物的组成、 结构和构象的研究是十分有效的。 我国化学家陈荣悌还报道了利用高效液相色谱法来测 定配位化合物稳定常数的新方法。 从上世纪 40 年代起，人们开始涉足于生物配合物的研究，陆续报道了各种生命金 属与各种 α-氨基酸、生物配体、羧肽酶和碳酸酐酶配合物的稳定常数。从这些研究中得 到了许多有关生物配体的新信息。到了 60 年代，配合物的研究方法又有了新的发展。 一方面出现了全新的物理方法， 如 M?ssbauer(穆斯堡尔)谱； 另一方面对原有的方法作了 改进，研究范围进一步扩大，如 Raman(拉曼)光谱采用了激光源，提高了光谱强度，使 有色配合物的研究成为可能；核磁共振谱因位移试剂的应用，提高了分辨率。这些方法 被应用于生物体内酶的结构和催化反应的研究以及化学模拟的研究。
3．测定吸光度 在 500nm 波长下分别测定上述溶液的吸光度，将所得数据记录于表 11-1。以吸光 度对磺基水杨酸的摩尔分数作图。 从所得的等摩尔系列图中找出最大吸收处的的配位体摩尔分数（XL）和金属离子摩 尔分数（1-XL） ，由公式（1）计算得配合物的组成；由公式（2）计算得配合物的解离 度；据公式（3）计算得配合物的表观稳定常数。
' K稳
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验十一一、编程并在计算机上调试修改运行1. P142例1、P145例2、P154例32. P155 2、P156 7、9二、阅读理解下列程序，并在计算机上运行.1. %利用e x的展开式计算函数e x的值(1) qexp.mfunction [E,k,F]=qexp(x,w,n)E=0;F=1;for k=1:nE=E+F;F=x*F/k;if abs(E-exp(x))<wbreakendend(2) qex.mfunction [E,k,F]=qex(x)E=0;F=1;k=1;while abs(F)>0E=E+F;F=x*F/k;k=k+1;end2.% 利用sin(x)的展开式计算其函数值(1) qsin.mfunction [s,k,T]=qsin(x)s=0;T=x;k=1;while abs(T)>0s=s+T;T=-x*x*T/(2*k*(2*k+1));k=k+1;end(2) qsn.mfunction [s,k,T]=qsn(x,w,n)s=0;T=x;for k=1:ns=s+T;T=-x*x*T/(2*k*(2*k+1));if abs(s-sin(x))<wbreakendend3. play.mwhile 1x=fix(100*rand);n=7;test=1;for k=1:7numb=int2str(n);disp(['You have a right to ',numb,' guesses'])disp(['A guess is a number between 0 and 100'])guess=input('Enter your guess:');if guess<xdisp('Low')elseif guess>xdisp('High')elsedisp('You won')test=0;break;endn=n-1;endif test==1disp('You lost')endm=input('继续？''Y''or''N''','s');if m=='Y'continueelsebreakendend一、1.P142例1function [y,k]=Jacobi(A,b,x,e,N)D=diag(diag(A));L=tril(A)-D;U=triu(A)-D;for k=1:Ny=-(D^(-1))*(L+U)*x+(D^(-1))*b;if norm(x'-y'，inf)>ex=y;elsebreakendendif k==Ndisp('µü´úʧ°Ü')end>> A=[10 -1 -2;-1 10 -2;-1 -1 5]; >> b=[7.2 8.3 4.2]';>> x=[0 0 0]';>> [y,k]=Jacobi(A,b,x,1e-3,100)y =1.09981.19981.2998k =8>> [y,k]=Jacobi(A,b,x,1e-4,100)y =1.10001.20001.3000k =10>> [y,k]=Jacobi(A,b,x,1e-6,100)y =1.10001.20001.3000k =15P145例2function [y,k]=GaussSeidel(A,b,x,e,N) D=diag(diag(A));L=tril(A)-D;U=triu(A)-D;for k=1:Ny=-((D+L)^(-1))*U*x+((D+L)^(-1))*b;if norm(x'-y'，inf)>ex=y;elsebreakendendif k==Ndisp('µü´úʧ°Ü')end>> A=[10 -1 -2;-1 10 -2;-1 -1 5];>> b=[7.2 8.3 4.2]';>> x=[0 0 0]';>> [y,k]=GaussSeidel(A,b,x,1e-3,100)y =1.09991.19991.3000k =5>> [y,k]=GaussSeidel(A,b,x,1e-4,100) y =1.10001.20001.3000k =7>> [y,k]=GaussSeidel(A,b,x,1e-6,100) y =1.10001.20001.3000k =9>> [y,k]=GaussSeidel(A,b,x,1e-7,100) y =1.10001.20001.3000k =10P154例3function [y,k]=SOR(A,b,x,w,e)%wΪËɳÚÒò×Ó£¬eΪ¾«¶ÈD=diag(diag(A));L=tril(A)-D;U=triu(A)-D;k=0;while 1y=((D+w*L)^(-1))*((1-w)*D-w*U)*x+w*((D+w*L)^(-1))*b;if norm(x'-y',inf)>ex=y;k=k+1;elsebreakendend>> A=[-4 1 1 1;1 -4 1 1;1 1 -4 1;1 1 1 -4];>> b=[1 1 1 1]';>> x=[0 0 0 0]';>> [y,k]=SOR(A,b,x,1,1e-5)y =-1.0000-1.0000-1.0000-1.0000k =20>> [y,k]=SOR(A,b,x,1.1,1e-5)y =-1.0000-1.0000-1.0000-1.0000k =16>> [y,k]=SOR(A,b,x,1.2,1e-5) y =-1.0000-1.0000-1.0000-1.0000k =11>> [y,k]=SOR(A,b,x,1.3,1e-5) y =-1.0000-1.0000-1.0000-1.0000k =11>> [y,k]=SOR(A,b,x,1.4,1e-5) y =-1.0000-1.0000-1.0000-1.0000k =14>> [y,k]=SOR(A,b,x,1.5,1e-5) y =-1.0000-1.0000-1.0000-1.0000k =17>> [y,k]=SOR(A,b,x,1.6,1e-5) y =-1.0000-1.0000-1.0000-1.0000k =23>> [y,k]=SOR(A,b,x,1.7,1e-5) y =-1.0000-1.0000-1.0000-1.0000k =34>> [y,k]=SOR(A,b,x,1.8,1e-5)y =-1.0000-1.0000-1.0000-1.0000k =542.P155 例2Jacobi>> A=[10 0 1 -5;1 8 -3 0;3 2 -8 1;1 -2 2 7]; >> b=[-7 11 23 17]';>> x=[0 0 0 0]';>> [y,k]=Jacobi(A,b,x,1e-3,100)y =1.00040.4998-1.99993.0002k =12>> [y,k]=Jacobi(A,b,x,1e-5,100)y =1.00000.5000-2.00003.0000k =19>> [y,k]=Jacobi(A,b,[1 0 1 0]',1e-5,100) y =1.00000.5000-2.00003.0000k =19>> [y,k]=Jacobi(A,b,[1 1 1 1]',1e-5,100) y =1.00000.5000-2.00003.0000k =18GaussSeild>> A=[10 0 1 -5;1 8 -3 0;3 2 -8 1;1 -2 2 7]; >> b=[-7 11 23 17]';>> x=[0 0 0 0]';>> [y,k]=GaussSeidel(A,b,x,1e-3,100)y =1.00020.4999-1.99992.9999k =8>> [y,k]=GaussSeidel(A,b,x,1e-5,100)y =1.00000.5000-2.00003.0000k =12>> [y,k]=GaussSeidel(A,b,[1 0 1 0]',1e-5,100) y =1.00000.5000-2.00003.0000k =12>> [y,k]=GaussSeidel(A,b,[1 1 1 1]',1e-5,100) y =1.00000.5000-2.00003.0000k =12P156 7(1)Jacobi>> A=[1 0 1;-1 1 0;1 2 -3];>> b=[5 -7 -17]';>> x=[0 0 0]';>> [y,k]=Jacobi(A,b,x,1e-3,100)迭代失败y =1.9957-5.00243.0011k =100>> [y,k]=Jacobi(A,b,x,1e-5,10000)y =2.0000-5.00003.0000k =209GaussSeild>> A=[1 0 1;-1 1 0;1 2 -3];>> b=[5 -7 -17]';>> x=[0 0 0]';>> [y,k]=GaussSeidel(A,b,x,1e-3,100) 迭代失败y =-1-8k =100>> [y,k]=GaussSeidel(A,b,x,1e-5,10000) 迭代失败y =-1-8k =10000(2)Jacobi>> A=[1 0.5 0.5;0.5 1 0.5;0.5 0.5 1]; >> b=[0 0.5 -2.5]';>> x=[0 0 0]';>> [y,k]=Jacobi(A,b,x,1e-3,100)迭代失败y =1.33332.3333-3.6667k =100>> [y,k]=Jacobi(A,b,x,1e-3,10000)迭代失败y =1.33332.3333-3.6667k =10000GaussSeild>> A=[1 0.5 0.5;0.5 1 0.5;0.5 0.5 1]; >> b=[0 0.5 -2.5]';>> x=[0 0 0]';>> [y,k]=GaussSeidel(A,b,x,1e-5,100) y =1.00002.0000-4.0000k =14>> [y,k]=GaussSeidel(A,b,x,1e-3,100) y =0.99972.0005-4.0001k =9P156 9>> A=[4 3 0;3 4 -1;0 -1 4];>> b=[16 20 -12]';>> x=[0 0 0]';>> [y,k]=SOR(A,b,x,1,0.5e-4)y =1.50013.3333-2.1667k =21>> [y,k]=SOR(A,b,x,1.5,0.5e-4) y =1.50003.3333-2.1667k =17>> [y,k]=SOR(A,b,x,1.9,0.5e-4) y =1.50003.3333-2.1666k =113二、1、（1）>> [E,k,F]=qexp(5,1e-5,1000)E =148.4132k =22F =2.1212e-006>> [E,k,F]=qexp(3,1e-4,1000)E =20.0855k =14F =5.4864e-005（2）>> [E,k,F]=qex(5)E =148.4132k =255F =>> [E,k,F]=qex(3)E =20.0855k =225F =2、（1）>> [s,k,T]=qsin(pi/3) s =0.8660k =91T =>> [s,k,T]=qsin(pi/2) s =1.0000k =98T =（2）>> [s,k,T]=qsn(pi/2,1e-4,100) s =1.0000k =5T =-3.5988e-006>> [s,k,T]=qsn(pi/3,1e-4,10) s =0.8660k =4T =4.1734e-0063、You have a right to 7 guessesA guess is a number between 0 and 100 Enter your guess:12LowYou have a right to 6 guessesA guess is a number between 0 and 100 Enter your guess:56HighYou have a right to 5 guessesA guess is a number between 0 and 100 Enter your guess:2LowYou have a right to 4 guessesA guess is a number between 0 and 100 Enter your guess:1LowYou have a right to 3 guessesA guess is a number between 0 and 100 Enter your guess:3LowYou have a right to 2 guessesA guess is a number between 0 and 100 Enter your guess:6LowYou have a right to 1 guessesA guess is a number between 0 and 100 Enter your guess:5LowYou lost继续？'Y'or'N'N21。