《数值计算方法》实验指导书(学生版)要点

合集下载

数值计算方法实验报告(含所有)

数值计算方法实验报告(含所有)

本科实验报告课程名称:计算机数值方法实验项目:计算机数值方法实验实验地点:专业班级:学号:学生姓名:xxx指导教师:xxx太原理工大学学生实验报告学院名称软件学院专业班级1217班学号201200xxxx 学生姓名xx 实验日期2014.05.21 成绩课程名称数值计算方法实验题目实验一方程求解一、实验目的和要求熟悉使用、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。

选择上述方法中的两种方法求方程:二分法f(x)=x3+4x2-10=0在[1,2]内的一个实根,且要求满足精度|x*-x n|<0.5×10-5二、主要设备笔记本 HP ProBook 6470b 一台编译软件:VC++6.0三、实验内容和原理函数f(x)在区间(x,y)上连续,先在区间(x,y)确定a与b,若f(a),f(b)异号,说明在区间(a,b)内存在零点,然后求f[(a+b)/2]。

假设F(a)<0,F(b)>0,a<b,①如果f[(a+b)/2]=0,该点即为零点;②如果f[(a+b)/2]<0,则区间((a+b)/2,b)内存在零点,(a+b)/2≥a;③如果f[(a+b)/2]>0,则区间(a,(a+b)/2)内存在零点,(a+b)/2≤b;返回①重新循环,不断接近零点。

通过每次把f(x)的零点所在区间收缩一半的方法,使区间内的两个端点逐步逼近函数零点,最终求得零点近似值。

四、操作方法与实验步骤1. 二分法:#include<stdio.h>#include<stdlib.h>#include<math.h>int main(){double a=1.0, b=2.0;double x,s;printf(" An\t\tBn\t\tF(Xn)\n");while(1){x=(a+b)/2;s=pow(x,3)+4*x*x-10;if (-0.000005 < s && s < 0.000005){break;}else if(s < 0){a=x;}else if(s > 0){b=x;}printf("%f\t%f\t%f\n",a,b,s);}printf("X的值为:%f\n",x);printf("误差:\t%f\n",s);return 0;}2. 割线法:#include"stdio.h"#include"math.h"int main(){float c,a=1.0,b=2.0;printf("每次得到的X的近似值:\n");while(1){c=b-(b*b*b+4*b*b-10)*(b-a)/(b*b*b+4*b*b-(a*a*a+4*a*a));if(fabs(b-c)<0.5*0.00001)break;b=c;printf("%f\n",b);}printf("X的值为:%f\n",c);}五、实验结果与分析二分法割线法分析:由程序知,使用二分法和割线法均能计算出方程的根,但利用割线法要比二分法计算的次数少,并且能够较早的达到精度要求。

数值计算方法学习指导

数值计算方法学习指导



x ∗ = ±10m × ( a1 × 10−1 + a2 × 10−2 + L + ak × 10− k + L + an × 10− n + L )
其中 m 是整数, a1 ≠ 0 , a1 , a2 , L , an 是 0 到 9 中的一个数字,若 x − x ≤

(1.1)
1 ×10m − n ,则 2
1 10 + 99
=
1 = 0.050125639L 10 + 9.94987
利用有效数字的多少来比较不同算法的优劣,说明了算法选取的重要性。 记x=
99, x ∗ = 9.94987 , e( x ∗ ) = x ∗ − x ,则 e( x ∗ ) ≤
1 × 10−5 , 2
由 e(10 − x ) ≈ − e( x ) 得
∗ ∗ ∗ ∗


运算得到的误差限分别为
∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ε ( x1 ± x2 ) ≈ ε ( x1 ) ± ε ( x2 ) ;ε ( x1 x2 ) ≈ x1 ε ( x2 ) + x2 ε ( x1 );
(1.2)
Hale Waihona Puke ∗ ∗ ∗ ∗ ∗ x1 ε ( x2 ) + x2 ε ( x1 ) ∗ x1 ε( ∗ ) ≈ ( x2 ≠ 0) 2 x2 x∗ 2
n
x 的相对误差限。
这是一元函数的误差传播问题,只需利用传播公式计算即可。 由 x = 10 ± 5% 知近似值为 x = 10 ,绝对误差限为 ε ( x ) = 5% 。
∗ ∗

1 −1 1 ∗ n (x ) = n

《数值计算方法实验》课程实验教学大纲

《数值计算方法实验》课程实验教学大纲

《数值计算方法实验》课程实验教学大纲课程名称:数值计算方法实验(Numerical Method Experiments)课程编码:10020400410 课程负责人:课程性质:独立设课课程属性:专业基础实验课学时学分:总学时18 总学分0.5 实验学时18 实验学分0.5开出时间:三年级上学期适用专业:信息与计算科学综合性、设计性实验项目数1个,总学时:4 其中课内学时:4课外学时:0主笔人:审核人:、批准人:一、课程简介《数值计算方法》在信息与计算科学领域有着非常重要的地位,为计算机编程提供算法;对培养学生的抽象思维能力,提高学生的编程能力有很重要的作用;是为我系信息与计算科学专业高年级学生开设的一门重要课程,它为计算机及其相关专业人员解决数值计算方面的问题提供方法,对提高学生的利用计算机解决实际问题的能力有很大帮助。

《数值计算方法实验》作为《数值计算方法》课程的必要实践环节,其主要目的是让学生在学习理论教学中关于典型数学问题的数值求解方法后,能够构造求解该类问题数值解的算法,并编程上机实现算法,在上机过程中加强对算法的理解,并应用算法去解决实际问题,另外通过编程练习提高学生的程序设计能力。

本实验课程中涉及MATLAB软件、插值、数据拟合、数值积分、线性方程(非线性方程)求解、矩阵特征值特征根的计算、微分方程求解等方面的内容。

二、实验目的及要求本课程一共安排了8个实验,要求学生能够依据课本提供的理论方法设计相应的算法,并利用Matlab数学软件平台编写程序求解特定问题的数值解,并在计算机上调试,进而验证算法,并可利用调试成功的程序解决实际问题。

1.验证性实验实验一和实验二主要使学生掌握Matlab的软件环境,并能应用Matlab编写数值计算方面的程序及绘图。

实验三,四,六,七,八是让学生针对理论课中学习过的不同问题编程求解他们的数值解并在计算机上验证。

2.设计性实验实验五为设计性实验,要求学生自行针对特定问题设计算法,根据算法编写程序,并引导学生对实验结果进行观察,分析和归纳,进而猜想出一般结果三、实验方式及要求1、验证性实验以传授知识为主,要求学生掌握基础知识、基本技能。

《数值计算方法》上机实验报告

《数值计算方法》上机实验报告

《数值计算方法》上机实验报告华北电力大学实验名称数值il•算方法》上机实验课程名称数值计算方法专业班级:电力实08学生姓名:李超然学号:200801001008 成绩: 指导教师:郝育黔老师实验日期:2010年04月华北电力大学实验报告数值计算方法上机实验报吿一.各算法的算法原理及计算机程序框图1、牛顿法求解非线性方程*对于非线性方程,若已知根的一个近似值,将在处展开成一阶xxfx ()0, fx ()xkk泰勒公式"f 0 / 2 八八,fxfxfxxxxx 0 0 0 0 0 kkkk2!忽略高次项,有,fxfxfxxx 0 ()()(),,, kkk右端是直线方程,用这个直线方程来近似非线性方程。

将非线性方程的**根代入,即fx ()0, X ,* fxfxxx 0 0 0 0, ,, kkkfx 0 fx 0 0,解出fX 0 *k XX,, k' fx 0 k水将右端取为,则是比更接近于的近似值,即xxxxk, Ik, Ikfx ()k 八XX, Ikk* fx()k这就是牛顿迭代公式。

,2,计算机程序框图:,见,,3,输入变量、输出变量说明:X输入变量:迭代初值,迭代精度,迭代最大次数,\0输出变量:当前迭代次数,当前迭代值xkl,4,具体算例及求解结果:2/16华北电力大学实验报吿开始读入l>k/fx()0?,0fx 0 Oxx,,01* fx ()0XX,,,?10kk, ,1,kN, ?xx, 10输出迭代输出X输出奇异标志1失败标志,3,输入变量、输出变量说明: 结束例:导出计算的牛顿迭代公式,并il •算。

(课本P39例2-16) 115cc (0), 求解结果:10. 75000010.72383710. 72380510. 7238052、列主元素消去法求解线性方程组,1,算法原理:高斯消去法是利用现行方程组初等变换中的一种变换,即用一个不为零的数乘 -个 方程后加只另一个方程,使方程组变成同解的上三角方程组,然后再自下而上 对上三角3/16华北电力大学实验报告方程组求解。

《数值计算方法》学习指南

《数值计算方法》学习指南

学习指南“数值计算方法”是以各类数学问题的数值解法作为研究对象,并结合现代计算机科学与技术为解决科学与工程中遇到的各类数学问题提供算法,其教学目的是:让学生掌握数值计算方法的基本概念和数值求解各类数学问题的基本方法,为学习者以后解决科学与工程中的实际问题打好基础。

在学习该课程应学习如下内容:一、误差分析1、要求准确掌握近似数绝对误差、相对误差和有效数字的基本概念,以及误差在近似数运算中的传播。

2、理解算法的数值稳定性问题。

二、多项式插值1、要求重点掌握构造插值多项式的两种方法:Lagrange插值法及Newton插值方法;2、掌握多项式插值的误差估计公式及它的导出;3、重点掌握构造低次Hermite插值多项式的基函数方法;4、了解关于样条函数的两种表示方法:截断幂函数表示法;B—样条函数表示法;5、掌握三次样条函数插值多项式构造的三弯矩方法;三、最佳逼近及其实现1、了解最佳一致逼近多项式的存在、唯一性及特征定理;2、掌握内积空间最佳逼近问题及最佳平方逼近多项式的计算问题;3、掌握切比雪夫多项式的性质及其应用问题;4、掌握离散数据拟合的最小二乘方法。

四、数值积分方法与数值微分1、重点掌握内插型求积公式的构造以及数值积分公式的代数精确度问题;2、了解外推法的思想和它在构造高精度数值求积公式中的应用;3、重点掌握Gauss型求积公式的性质和它的构造;4、理解数值微分公式构造的基本方法。

五、线性代数方程组的解法1、了解Gauss消去法的计算过程与Gauss消去法的实质——矩阵的三角分解;2、了解基于矩阵分解的Doolittle方法和Cholesky方法;3、重点掌握用追赶法求解三对角方程组;4、掌握简单迭代法收敛的充分必要条件和若干充分条件;5、重点掌握Jacobi迭代法与Gauss—Seidel迭代法的计算格式与收敛的充分条件;6、理解共轭向量系实对称矩阵的共轭梯度法;六、矩阵特征值问题的解法1、了解估计矩阵特征值的盖尔园方法与对称矩阵的极大、极小定理;2、重点掌握乘幂法的计算格式以及求矩阵特征值与特征向量的反乘幂法;3、掌握Householder矩阵的基本性质以及它在矩阵约化中的作用;4、理解求实对称矩阵特征值的二分法的基本思想;5、理解用QR方法求矩阵特征问题。

数值计算方法实验指导书

数值计算方法实验指导书

3) 特殊的三维图形函数 * [x,y,z]=sphere(n) % 画球,n 默认值 20 例: >> [a,b,c]=sphere(40); >> surf(a,b,c) >> axis('equal'); >> axis('square'); * [x,y,z]=cylinder(R,N) 格数 例: >> >> >> >> x=0:pi/20:pi*3; r=5+cos(x); [a,b,c]=cylinder(r,30); mesh(a,b,c) %R 母线向量,N 分
-1-
数值计算方法实验指导书
电子与信息工程系
-2-
数值分析实验< 数值分析实验<一> ---- Matlab 绪论 一\ 实验目的 1) 熟悉 Matlab 的运行环境及各种窗口 2) 掌握 Matlab 的矩阵变量类型,矩阵输入和矩阵的基本运算 3) 掌握命令及函数文件的作用及区别,并编写简单的 M 文件 4) 能熟练的向查寻目录中添加新目录,掌握常用的 Matlab 系统命令 二\ 实验内容 一> Matlab 启动与环境设置 1)启动 双击桌面图标 开始>程序>Matlab 安装目录>bin>matlab 2)环境设置 命令窗口(Command Window) 执行命令行,Matlab 主窗口; 窗口颜色及字体 File>Preferences.. 当前目录(Current Directory) File>Set Path 用于将新文件夹加入搜索路径,设置当前文件默 认目录; 3)Matlab 常用命令 上下箭头 调出最新用过的命令,重新执行 cd+目录名 改变当前目录 help 显示当前搜索路径中所有目录名称 help+函数(类)名 查找函数(类),给出函数用法及参数 lookfor+函数关键字 查询根据关键字搜索到的相关函数 exist+变量名 变量检验函数 what 目录中文件列表 who 内存变量列表 whos 内存变量详细信息 which 确定文件位置 clc 清屏 ! 调用 Dos 命令 4)联机演示系统 Help>Demos.. 输入命令:intro 二>Matlab 基本运算操作 1)数据类型 变量 区分大小写,长度不超过 31,字母开头 常量 i,j 虚单位,定义 sqrt(-1) pi 圆周率 eps 浮点运算的相对精度 exp(-52) NaN Not-a-Number,表示不定值

数值计算方法与算法第三版答案 数值计算方法学习指导书

数值计算方法与算法第三版答案 数值计算方法学习指导书

数值计算方法与算法第三版答案数值计算方法学习指导书数值计算方法学习指导书是怎么样的?以下是小编分享给大家的数值计算方法学习指导书简介的资料,希望可以帮到你!数值计算方法学习指导书内容简介《数字信号处理学习指导》是浙江省高等教育重点建设教材、应用型本科规划教材《数字信号处理》(唐向宏主编,浙江大学出版社出版,以下简称教材)的配套学习指导书,内容包括学习要求、例题分析、教材习题解答、自测练习以及计算机仿真实验等。

学习指导书紧扣教材内容,通过例题讲解,分析各章节的学习重点、难点以及需要理解、掌握和灵活运用的基本概念、基本原理和基本方法。

全书共有66例例题分析、121题题解、2套自测练习和6个MAT1AB计算机仿真实验。

数值计算方法学习指导书目录绪论第1章离散时间信号与系统1.1 学习要点1.2 例题1.3 教材习题解答第2章离散系统的变换域分析与系统结构2.1 学习要点2.2 例题2.3 教材习题解答第3章离散时间傅里叶变换3.1 学习要点3.2 例题3.3 教材习题解答第4章快速傅里叶变换4.1 学习要点4.2 例题4.3 教材习题解答第5章无限长单位冲激响应(IIR)数字滤波器的设计5.1 学习要点5.2 例题5.3 教材习题解答第6章有限长单位冲激响应(FIR)数字滤波器的设计6.1 学习要点6.2 例题6.3 教材习题解答第7章数字信号处理中的有限字长效应7.1 学习要点7.2 例题7.3 教材习题解答第8章自测题8.1 自测题(1)及参考答案8.2 自测题(2)及参考答案第9章基于MA TLAB的上机实验指导9.1 常见离散信号的MA TLAB产生和图形显示9.2 信号的卷积、离散时间系统的响应9.3 离散傅立叶变换9.4 离散系统的频率响应分析和零、极点分布9.5 IIR滤波器的设计9.6 FIR滤波器的设计数值计算方法学习指导书内容文摘第1章离散时间信号与系统1.1 学习要点本章主要介绍离散时间信号与离散时间系统的基本概念,着重阐述离散时间信号的表示、运算,离散时间系统的性质和表示方法以及连续时间信号的抽样等。

数值计算方法-实验指导书-8学时

数值计算方法-实验指导书-8学时

本科生实验报告实验课程数值计算方法学院名称机电工程学院专业名称机械工程学生姓名学生学号指导教师实验地点实验成绩二〇二三年十一月二〇二三年十二月目录实验1 MATLAB的基本运算 (1)实验2 MALTAB多项式运算 (14)实验3 图形的可视化及修饰处理 (16)实验4 SIMULINK仿真实验 (19)实验1 MATLAB的基本运算一、实验目的基本掌握MATLAB向量、矩阵、数组的生成及其基本运算(区分数组运算和矩阵运算)、常用的数学函数、以及字符串的操作。

二、实验内容1向量的生成和运算;2矩阵的创建、引用和运算;3多维数组的创建及运算;4字符串的操作。

三、实验步骤1.向量的生成和运算(1) 向量的生成①直接输入法:>> A=[1;2;3;4;5]②冒号生成发:>> A=1:2:10 ,B=1:10 ,C=10:-1:1函数生成法:linspace( )是线性等分函数,logspace( )是对数等分函数。

>> A=linspace(1,10)B=linspace(1,30,10)>> A=logspace(0,4,5)(2) 向量的运算维数相同的行向量只可以相加减,维数相同的列向量也可以相加减,标量可以与向量直接相乘除。

>> A=[1 2 3 4 5], B=3:7,>> AT=A', BT=B',>> E1=A+B, E2=A-B>> F=AT-BT,>> G1=3*A, G2=B/3,②向量的点积与叉积运算。

>> A=ones(1,10);B=1:10; BT=B';>> E1=A*BT2.矩阵的创建、引用和运算(1)矩阵的创建和引用m 元素构成的矩阵结构,行向量和列向量是矩阵的特殊形式。

矩阵是由n①直接输入法:>> A=[1 2 3;4 5 6]>> B=[ 1 4 72 5 83 6 9 ]>> A(1)>> A(4:end)>> B(:,1)>> B(:)>> B(5)抽取法>> clear>> A=[1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16] >> B=A(1:3,2:3)>> C=A([1 3],[2 4])>> A([1 3;2 4])③函数法:>> A=ones(3,4)>> B=zeros(3)>> C=eye(3,2)>> D=magic(3)(2) 矩阵的运算①矩阵的加减、数乘与乘法已知矩阵:>> A=[1 23 -1],>> B=[-1 01 2]>> A+B >> 2*A >> 2*A-3*B >> A*B②矩阵的逆矩阵>> format rat;A=[1 0 1;2 1 2;0 4 6] >> A1=inv(A)>> A*A1③矩阵的除法>> a=[1 2 1;3 1 4;2 2 1],b=[1 1 2],d=b'>> c1=b*inv(a),c2=b/a>> c3=inv(a)*d, c4=a\d3.多维数组的创建及运算数组运算用小圆点加在运算符的前面表示,以区分矩阵的运算。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

理学院《数值计算方法》实验指导书适合专业:信息与计算科学数学与应用数学统计学贵州大学二OO七年八月前言《数值计算方法》包括很多常用的近似计算的处理手段和算法,是信息与计算科学,数学与应用数学,统计学等专业的必修课程。

为了加强学生对该门课程的理解,使学生更好地掌握书中的数值计算方法、编制程序的能力,学习数值计算方法课程必须重视实验环节,即独立编写出程序,独立上机调试程序,必须保证有足够的上机实验时间。

在多年教学实践基础上编写了《数值计算方法》实验指导书,目的是通过上机实验,使学生能对教学内容加深理解,同时培养学生动手的能力。

本实验指导书,可与《数值计算方法》教材配套使用,但是又有独立性,它不具体依赖哪本教材,主要的计算方法在本指导书中都有,因此,凡学习数值计算方法课程的学生都可以参考本指导书进行上机实验。

上机结束后,按要求整理出实验报告。

实验报告的内容参阅《数值计算方法》实验指导书。

目录第一章函数基本逼近(一)——插值逼近实验一Lagrange插值法第二章函数基本逼近(二)——最佳逼近实验二数据拟合的最小二乘法第三章数值积分与数值微分实验三自适应复化求积法第四章线性代数方程组求解实验四Gauss列主元消去法实验五解三对角方程组的追赶法实验六Jacobi迭代法第五章非线性方程的数值解法实验七Newton迭代法第六章常微分方程数值解法实验八常微分方程初值问题的数值方法实验一 Lagrange 插值法实验学时:2 实验类型:验证 实验要求:必修一.实验目的010100,()()()().n nnnjn n ij i i jj iagrange x x x x yy y y x agrange x f x x x yL L x x ==≠-=≈-∑∏通过L 插值法的学习掌握如何根据已知函数表构造L 插值多项式用二.实验内容1.算法设计。

2.编写相应的程序上机调试。

3.已知下列函数表0.320.340.36sin 0.3145670.3334870.352274x x用上述程序验证用线性插值计算sin 0.3367的近似值为0.330365,用抛物插值计算sin 0.3367的近似值为0.330374。

4.已知下列函数表111213ln 2.3979 2.4849 2.5649x x用上述程序分别用线性插值与抛物插值计算ln11.75的近似值。

三.实验组织远行要求统一进行实验,一人一组。

四.实验条件PC 机,vc++6.0,Internet 网。

五.实验步骤1.根据算法事先写出相应程序。

2.启动PC 机,进入vc 集成环境,输入代码。

3.编译调试。

4.调试通过,计算出正确结果。

实验二 数据拟合的最小二乘法实验学时:2 实验类型:验证 实验要求:必修一.实验目的掌握数据拟合的思想,清楚数据拟合与插值法的区别及最小二乘原理在数据拟合中的重要作用及最小二乘解的求法。

二.实验内容1.算法设计。

2.编写相应的程序上机调试。

3.给定数据点0.00.20.40.60.8 1.0 1.20.9 1.9 2.8 3.3 4.0 5.7 6.5i ix y用上述程序验证用一次多项式拟合这组数据的结果为1()0.8428574.571429s x x =+。

4.给定数据点310135631013i ix y -----用上述程序分别用一次、二次和三次多项式拟合这组数据。

三.实验组织远行要求统一进行实验,一人一组。

四.实验条件PC 机,vc++6.0,Internet 网。

五.实验步骤1.根据算法事先写出相应程序。

2.启动PC 机,进入vc 集成环境,输入代码。

3.编译调试。

4.调试通过,计算出正确结果。

实验三 自适应复化求积法实验学时:2 实验类型:验证 实验要求:必修一.实验目的复合Simpson 公式是提高精度的行之有效的方法,但是在使用求积公式之前,必须先给出步长。

步长取得太大精度难以保证,步长太小则导致计算量的增加,而事先给出一个合适的步长往往是困难的,因此提出了在求积过程中根据精度要求自动确定积分步长的Simpson 公式的逐次分半加速法。

二.实验内容1.算法设计。

2.编写相应的程序上机调试。

3.用上述程序验证1sin x dx x ⎰的值为0.946083,61102ε-=⨯。

4、用上述程序计算2170,10x e dx ε--=⎰。

5、用上述程序计算60,10ε-=⎰。

三.实验组织远行要求统一进行实验,一人一组。

四.实验条件PC 机,vc++6.0,Internet 网。

五.实验步骤1.根据算法事先写出相应程序。

2.启动PC 机,进入vc 集成环境,输入代码。

3.编译调试。

4.调试通过,计算出正确结果。

实验四 Gauss 列主元消去法实验学时:2 实验类型:验证实验要求:必修一.实验目的用Gauss 列主元消去法解线性方程组111211*********2.n n n n nn n n a a a x b a a a x b a a a x b ⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭二.实验内容1.算法设计。

2.编写相应的程序上机调试。

3、用上述程序验证线性方程组12331831111-⎛⎫ ⎪-- ⎪ ⎪⎝⎭123x x x ⎛⎫ ⎪ ⎪ ⎪⎝⎭=15156⎛⎫⎪- ⎪ ⎪⎝⎭的解为1231.000000, 2.000000, 3.000000x x x ===。

4 、用上述程序解线性方程组123112*********x x x --⎛⎫⎛⎫⎛⎫⎪⎪ ⎪--= ⎪⎪ ⎪ ⎪⎪ ⎪--⎝⎭⎝⎭⎝⎭。

三.实验组织远行要求统一进行实验,一人一组。

四.实验条件PC 机,vc++6.0,Internet 网。

五.实验步骤1.根据算法事先写出相应程序。

2.启动PC 机,进入vc 集成环境,输入代码。

3.编译调试。

4.调试通过,计算出正确结果。

实验五 解三对角方程组的追赶法实验学时:2 实验类型:验证 实验要求:必修一.实验目的用追赶法解三对角线性方程组,Ax f =其中11112222211111,,0:,0(2,31)n n n n n nn i i i i i n n b c x f a b c x f A x f a b c x f a b A i n b c b a c a c b a ---⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪ ⎪ ⎪ ⎪⎪=== ⎪ ⎪ ⎪⎪ ⎪⎪⎝⎭⎝⎭ ⎪⎝⎭⎧>>⎪⎪≥+≠=-⎨⎪>>⎪⎩并且满足二.实验内容1.算法设计。

2.编写相应的程序上机调试。

3、用上述程序验证线性方程组3123123113⎛⎫ ⎪ ⎪ ⎪ ⎪⎝⎭1234x x x x ⎛⎫ ⎪ ⎪ ⎪ ⎪⎝⎭=1010⎛⎫⎪ ⎪ ⎪ ⎪⎝⎭的解为12340.552632,0.657895,0.868421,0.289474x x x x ==-==-。

4、用上述程序解线性方程组2112112112112-⎛⎫ ⎪-- ⎪⎪-- ⎪-- ⎪ ⎪-⎝⎭12345x x x x x ⎛⎫ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭10007⎛⎫⎪ ⎪⎪ ⎪⎪⎪⎝⎭。

三.实验组织远行要求统一进行实验,一人一组。

四.实验条件PC 机,vc++6.0,Internet 网。

五.实验步骤1.根据算法事先写出相应程序。

2.启动PC 机,进入vc 集成环境,输入代码。

3.编译调试。

4.调试通过,计算出正确结果。

实验六 Jacobi 迭代法实验学时:2 实验类型:验证 实验要求:必修一.实验目的.Jacobi Ax b =用迭代法解线性方程组二.实验内容1.算法设计。

2.编写相应的程序上机调试。

3、用上述程序验证线性方程组411141114⎛⎫ ⎪ ⎪ ⎪⎝⎭123x x x ⎛⎫ ⎪ ⎪ ⎪⎝⎭=7107⎛⎫⎪⎪ ⎪⎝⎭的解为4(0)1230.999980, 1.999980,0.999980,10,(0,0,0)T x x x x ε-=====。

4 、用上述程序解线性方程组152310127.211028.3,10,115 4.2x x x ε---⎛⎫⎛⎫⎛⎫⎪⎪ ⎪--== ⎪⎪ ⎪ ⎪⎪ ⎪--⎝⎭⎝⎭⎝⎭(0)(0,0,0)T x =。

三.实验组织远行要求统一进行实验,一人一组。

四.实验条件PC 机,vc++6.0,Internet 网。

五.实验步骤1.根据算法事先写出相应程序。

2.启动PC 机,进入vc 集成环境,输入代码。

3.编译调试。

4.调试通过,计算出正确结果。

实验七 Newton 迭代法实验学时:2实验类型:验证实验要求:必修一.实验目的利用N e w t o n 迭代公式1()(0,1,)()n n n n f x x x n f x +=-='求非线性方程()0f x =的近似根。

二.实验内容1.算法设计。

2.编写相应的程序上机调试。

3. 用上述程序验证10x xe -=在00.5x =附近的根为410.567143,102ε-=⨯,允许最大迭代次数为100。

3. 用上述程序求3310x x --=在02x =附近的根,610ε-=,允许最大迭代次数为200。

三.实验组织远行要求统一进行实验,一人一组。

四.实验条件PC 机,vc++6.0,Internet 网。

五.实验步骤1.根据算法事先写出相应程序。

2.启动PC 机,进入vc 集成环境,输入代码。

3.编译调试。

4.调试通过,计算出正确结果。

实验八 常微分方程初值问题的数值方法实验学时:4实验类型:综合实验要求:必修一.实验目的用数值微分法与数值积分法求一阶常微分方程初值问题0(,)[,]()y f x y x a b y a y '=∈⎧⎨=⎩在离散点01n a x x x b ≤<<<≤上的近似值12,,,n y y y 。

二.实验内容1.算法设计。

2.编写相应的程序上机调试,并对实验结果进行分析和比较。

3.取0.1h =,分别用E u l e r 法,改进E u l e r 法求解初值问题2[0,1](0)1x y y x y y ⎧'=-∈⎪⎨⎪=⎩,将其绝对误差限相比较,说明哪种方法精度高,准确解y =4.分别取0.025h =和0.005h =,用E u l e 法求解初值问题100[0,0.1](0)1y y x y '=-∈⎧⎨=⎩,输出准确解与数值解,并分析结果,准确解100x y e -=。

三.实验组织远行要求统一进行实验,一人一组。

四.实验条件PC 机,vc++6.0,Internet 网。

五.实验步骤1.根据算法事先写出相应程序。

2.启动PC 机,进入vc 集成环境,输入代码。

相关文档
最新文档