数值分析课程课程设计汇总

合集下载

高等数值分析课程设计

高等数值分析课程设计

高等数值分析课程设计一、题目背景高等数值分析是计算数学领域的一门重要课程,它主要研究数值计算中的算法、误差分析、收敛性和稳定性等基本问题,涵盖了线性代数、数值微积分、常微分方程数值解等数学分支学科。

本文将介绍一项高等数值分析课程的设计,以增强学生对课程的理解和能力。

二、设计目标2.1 教学目标本课程设计旨在帮助学生:•掌握常见的数值分析算法;•熟悉各种算法的误差分析和收敛性;•能够独立设计和实现数值计算程序;•培养学生解决实际问题的能力。

2.2 实现目标为了实现教学目标,本课程设计将遵循以下原则:•采用案例分析和实例演示的方式,将数学理论与实际应用相结合;•强调算法的实现方法和效率分析;•通过小组合作的方式完成实践任务,培养学生的团队合作能力;•开设课程论文撰写指导和实践报告撰写指导课程,提高学生的学术写作能力。

三、课程内容本课程的教学安排如下:3.1 理论讲授•数值线性代数•数值微积分•常微分方程数值解•偏微分方程数值解3.2 实践任务•实现线性方程组求解算法•实现求解非线性方程的算法•实现常微分方程数值解算法•实现偏微分方程数值解算法3.3 课程论文和实践报告撰写要求每个学生提交一篇课程论文和一份实践报告,内容包括理论和实践部分。

论文部分主要包括:•算法的理论分析和数学推导;•算法的实现方法和效率分析;•算法的收敛性和稳定性分析。

实践报告部分主要包括:•实践任务的设计和实现方法;•算法实现的过程与结果分析;•算法的应用和实用性分析。

四、教学评估本课程的教学评估主要包括以下几个方面:4.1 学生成绩评估学生成绩评估包括平时分、实验成绩、论文得分和考试成绩。

其中,实验成绩和论文得分占总成绩的比重大于考试成绩。

4.2 教学效果评估教学效果评估将从以下几个方面进行:•学生数学知识的掌握程度;•学生对数值计算的算法和方法的理解程度;•学生的编程能力和算法实现的水平;•学生实践能力和团队协作能力的培养。

大学数值分析课程设计

大学数值分析课程设计

大学数值分析课程设计一、课程目标知识目标:1. 理解数值分析的基本概念,掌握数值计算方法及其数学原理;2. 掌握线性代数、微积分等基本数学工具在数值分析中的应用;3. 学会分析数值算法的稳定性和误差,评估数值结果的正确性。

技能目标:1. 能够运用数值分析方法解决实际工程和科学研究问题;2. 掌握常用数值分析软件的使用,提高数据处理和问题求解的效率;3. 培养编程实现数值算法的能力,提高解决复杂问题的技能。

情感态度价值观目标:1. 培养学生对数值分析的浓厚兴趣,激发学习积极性;2. 培养学生的团队合作精神,提高沟通与协作能力;3. 增强学生的数学素养,使其认识到数学在科学研究和社会发展中的重要性。

课程性质分析:本课程为大学数值分析课程,旨在教授学生数值计算的基本理论和方法,培养学生解决实际问题的能力。

学生特点分析:学生具备一定的高等数学基础,具有较强的逻辑思维能力和抽象思维能力。

教学要求:1. 注重理论与实践相结合,提高学生的实际操作能力;2. 鼓励学生主动参与讨论,培养学生的创新意识和解决问题的能力;3. 结合实际案例,强化学生对数值分析在工程和科研中的应用认识。

二、教学内容1. 数值分析基本概念:包括误差分析、稳定性、收敛性等;教材章节:第一章 数值分析概述2. 数值线性代数:矩阵运算、线性方程组求解、特征值与特征向量计算等;教材章节:第二章 线性代数的数值方法3. 数值微积分:数值积分、数值微分、常微分方程数值解等;教材章节:第三章 微积分的数值方法4. 非线性方程与系统求解:迭代法、牛顿法、弦截法等;教材章节:第四章 非线性方程与系统的数值解法5. 优化问题的数值方法:线性规划、非线性规划、最小二乘法等;教材章节:第五章 优化问题的数值方法6. 数值模拟与数值实验:蒙特卡洛方法、有限元方法、差分方法等;教材章节:第六章 数值模拟与数值实验7. 数值软件应用:MATLAB、Python等数值计算软件在数值分析中的应用;教材章节:第七章 数值软件及其应用教学进度安排:第1-2周:数值分析基本概念第3-4周:数值线性代数第5-6周:数值微积分第7-8周:非线性方程与系统求解第9-10周:优化问题的数值方法第11-12周:数值模拟与数值实验第13-14周:数值软件应用及综合案例分析教学内容确保科学性和系统性,注重理论与实践相结合,提高学生的实际操作能力。

《数值分析》课程教案

《数值分析》课程教案

《数值分析》课程教案数值分析课程教案一、课程介绍本课程旨在介绍数值分析的基本概念、方法和技巧,以及其在科学计算和工程应用中的实际应用。

通过本课程的研究,学生将了解和掌握数值分析的基本原理和技术,以及解决实际问题的实用方法。

二、教学目标- 了解数值分析的基本概念和发展历程- 掌握数值计算的基本方法和技巧- 理解数值算法的稳定性和收敛性- 能够利用数值分析方法解决实际问题三、教学内容1. 数值计算的基本概念和方法- 数值计算的历史和发展- 数值计算的误差与精度- 数值计算的舍入误差与截断误差- 数值计算的有效数字和有效位数2. 插值与逼近- 插值多项式和插值方法- 最小二乘逼近和曲线拟合3. 数值微积分- 数值积分的基本原理和方法- 数值求解常微分方程的方法4. 线性方程组的数值解法- 直接解法和迭代解法- 线性方程组的稳定性和收敛性5. 非线性方程的数值解法- 迭代法和牛顿法- 非线性方程的稳定性和收敛性6. 数值特征值问题- 特征值和特征向量的基本概念- 幂迭代法和QR方法7. 数值积分与数值微分- 数值积分的基本原理和方法- 数值微分的基本原理和方法四、教学方法1. 理论讲授:通过课堂授课,讲解数值分析的基本概念、原理和方法。

2. 上机实践:通过实际的数值计算和编程实践,巩固和应用所学的数值分析知识。

3. 课堂讨论:组织学生进行课堂讨论,加深对数值分析问题的理解和思考能力。

五、考核方式1. 平时表现:包括课堂参与和作业完成情况。

2. 期中考试:对学生对于数值分析概念、原理和方法的理解程度进行考查。

3. 期末项目:要求学生通过上机实验和编程实践,解决一个实际问题,并进行分析和报告。

六、参考教材1. 《数值分析》(第三版),贾岩. 高等教育出版社,2020年。

2. 《数值计算方法》,李刚. 清华大学出版社,2018年。

以上是《数值分析》课程教案的概要内容。

通过本课程的研究,学生将能够掌握数值分析的基本原理和技术,并应用于实际问题的解决中。

数值分析课程设计

数值分析课程设计

数值分析课程设计一、题目描述在本次数值分析课程设计中,我们需要实现下列内容:给定一个函数f(x),任取一个初值x0,使用牛顿法求出f(x)=0的一个根。

二、算法实现在数值计算中,牛顿法(Newton’s method) 是一种迭代算法,可以快速地求解方程的数值解,对于一般的实数函数,牛顿法可以用来求方程f(x)=0的根。

设x n是f(x)的根的一个近似值,y=f(x n)是对应的函数值,则用f(x)的一阶泰勒展开式$$ f(x) \\approx f(x_n)+f'(x_n)(x-x_n) $$且令上式等于零,得到牛顿迭代公式:$$ x_{n+1}=x_n-\\frac{f(x_n)}{f'(x_n)} $$若x0是f(x)的一个根的初始近似值,则$$ x_{n+1}=x_n-\\frac{f(x_n)}{f'(x_n)}, \\ n=0,1,2,\\cdots $$是迭代序列,如果 $\\lim\\limits_{n\\rightarrow \\infty}x_n=\\alpha$,且 $f(\\alpha)=0$,则 $\\alpha$ 是方程的一个根。

三、实验步骤1.确定初始值x0,计算f(x0)和f′(x0)。

2.按照牛顿法迭代公式计算x n+1。

3.如果满足指定的条件,则停止迭代,并输出x n+1。

4.否则,返回第二步迭代计算x n+2,直至满足指定的条件。

四、实验代码def newton_method(f, df, x0, eps=1e-8, max_iter=1000):'''利用牛顿法求解非线性方程f(x)=0的根。

:param f: 函数:param df: 导函数:param x0: 初值:param eps: 容差:param max_iter: 最大迭代次数:return:近似解'''n =1while True:x1 = x0 - f(x0) / df(x0)if abs(x1 - x0) < eps or n > max_iter:return x1x0 = x1n +=1五、实验结果我们使用上述实现的牛顿法来解决如下问题:$$ f(x) = x^2-3, \\ x_0=2 $$则f′(x)=2x。

数值分析方法课程设计

数值分析方法课程设计

数值分析方法课程设计背景介绍数值分析是一门研究求解各种数学问题的有效数值计算方法的学科,其应用广泛,如科学计算、工程设计和金融计算等领域。

在数值分析中,许多方法依赖于计算机的计算能力。

此外,数值分析还需要对数学理论和计算机科学两方面的知识有较深的理解。

本课程设计旨在通过实践,帮助学生深入了解数值分析方法及其应用,并提高学生的计算机编程能力。

课程设计目标•熟练掌握数值分析中的基本算法和方法,如插值法、数值积分等•能够将所学算法应用于实际问题,并编写可靠的程序解决问题•加深对计算机编程的理解和掌握,增强编程实践和创新能力•提高对数值分析和计算机科学交叉领域的理解课程内容第一部分:基本算法和方法1.数值微积分基本概念和原理2.插值法及其在实际中的应用3.数值积分的基本方法和理论基础4.常微分方程常用数值解法第二部分:实践应用与编程实现1.利用插值法和数值积分求解实际问题2.实现数值微积分和常微分方程的求解程序3.利用现有的数值分析软件解决实际问题,如 MATLAB 和 Python 等课程设计方案1.向学生介绍数值分析基本算法和方法,并讲解其理论基础和实际应用。

2.向学生提供一些实际问题,引导学生根据所学算法和方法进行求解。

3.给予学生一定的编程实践机会,让他们能够将所学算法实现为程序,并运用到具体的问题中。

4.通过课程作业、仿真实验等形式对学生进行考核和评价,确保学生能够有效掌握所学知识和能力。

评价标准1.学生掌握数值分析基本算法和方法的程度2.学生在实际问题中应用所学算法的能力3.学生编程实践和创新能力的水平4.学生对数值分析和计算机科学交叉领域的理解总结本课程设计旨在培养学生的数值分析和计算机编程实践,通过课程作业和编程实践等形式将理论知识与实际问题相结合,提高学生的实践应用能力。

同时,本课程设计也为学生未来的研究和工作提供了一定的基础。

数值分析课程设计(最终版)

数值分析课程设计(最终版)

数值分析课程设计(最终版)本⽂主要通过Matlab 软件,对数值分析中的LU 分解法、最⼩⼆乘法、复化Simpon 积分、Runge-Kutta ⽅法进⾏编程,并利⽤这些⽅法在MATLAB 中对⼀些问题进⾏求解,并得出结论。

实验⼀线性⽅程组数值解法中,本⽂选取LU 分解法,并选取数据于《数值分析》教材第5章第153页例5进⾏实验。

所谓LU 分解法就是将⾼斯消去法改写为紧凑形式,可以直接从矩阵A 的元素得到计算L 、U 元素的递推公式,⽽不需要任何步骤。

⽤此⽅法得到L 、U 矩阵,从⽽计算Y 、X 。

实验⼆插值法和数据拟合中,本⽂选取最⼩⼆乘拟合⽅法进⾏实验,数据来源于我们课堂学习该章节时的课件中的多项式拟合例⼦进⾏实验。

最⼩⼆乘拟合是⼀种数学上的近似和优化,利⽤已知的数据得出⼀条直线或者曲线,使之在坐标系上与已知数据之间的距离的平⽅和最⼩。

利⽤excel 的⾃带函数可以较为⽅便的拟合线性的数据分析。

实验三数值积分中,本⽂选取复化Simpon 积分⽅法进⾏实验,通过将复化Simpson 公式编译成MATLAB 语⾔求积分∫e ;x dx 10完成实验过程的同时,也对复化Simpon 积分章节的知识进⾏了巩固。

实验四常微分⽅程数值解,本⽂选取Runge-Kutta ⽅法进⾏实验,通过实验了解Runge-Kutta 法的收敛性与稳定性同时学会了学会⽤Matlab 编程实现Runge-Kutta 法解常微分⽅程,并在实验的过程中意识到尽管我们熟知的四种⽅法,事实上,在求解微分⽅程初值问题,四阶法是单步长中最优秀的⽅法,通常都是⽤该⽅法求解的实际问题,计算效果⽐较理想的。

实验五数值⽅法实际应⽤,本⽂采⽤最⼩⼆乘法拟合我国2001年到2015年的⼈⼝增长模型,并预测2020年我国⼈⼝数量。

关键词:Matlab ;LU 分解法;最⼩⼆乘法;复化Simpon 积分;Runge-Kutta⼀.LU分解法 (1)1.1实验⽬的 (1)1.2基本原理 (1)1.3实验内容 (2)1.4数据来源 (3)1.5实验结论 (3)⼆.Lagrange插值 (4)2.1实验⽬的 (4)2.2基本原理 (5)2.3实验内容 (5)2.4数据来源 (6)2.5实验结论 (6)三.复化simpon积分 (7)3.1实验⽬的 (7)3.2基本原理 (7)3.3实验内容 (7)3.4数据来源 (8)3.5实验结论 (8)四.Runge-Kutta⽅法 (9)4.1实验⽬的 (9)4.2基本原理 (9)4.3实验内容 (10)4.4数据来源 (11)4.5实验结论 (11)五.数值⽅法实际应⽤ (11)5.1实验⽬的 (11)5.2基本原理 (12)5.3实验内容 (12)5.4数据来源 (13)5.5实验结论 (13)总结 (16)参考⽂献 (17)⼀.LU 分解法1.1实验⽬的[1] 了解LU 分解法的基本原理和⽅法;[2] 通过实例掌握⽤MATLAB 求线性⽅程组数值解的⽅法; [3] 编程实现LU 分解1.2基本原理对于矩阵A ,若存在⼀个单位下三⾓矩阵L 和⼀个上三⾓U ,使得A =LU (1.1)。

[调研报告]《数值分析》课程设计报告范文

[调研报告]《数值分析》课程设计报告范文

课程设计报告课程设计题目: 非线性方程求解2011年 11月 27 日题目:用二分法,简单迭代法、牛顿迭代法以及弦截法求非线性方程2sin 02-=x x误差不超过10-4,输出迭代次数,初始值和根的近似值。

一、摘要在matlab 环境下运用熟悉的计算机编程语言结合二分法、简单迭代法、牛顿迭代法以及弦截法求解非线性方程,在运行完程序后,对运行结果做出了各方面的分析和比较。

最终得出二分法迭代次数最多,需14次,而简单迭代法、牛顿迭代法以及弦截法的迭代次数都较少,只需4—5次。

由于方程有多个解,所以当赋的初始值不同或给定的区间不同时,根的近似值也会有所不同。

二、设计目的用熟悉的计算机语言编程,上机完成用二分法、简单迭代法、牛顿迭代法以及弦截法求解非线性方程,掌握各种方法的理论依据及求解思路,了解各种迭代方法的异同。

三、理论基础二分法:二分法就是将方程根所在的区间平分为两个小区间,再判断根属于哪个小区间;把有根的小区间再平分为二,再判断根所在的更小的区间,对分;重复这一过程,最后求出所要的近似值。

简单迭代法:简单迭代法是将方程()0f x =化为一个等价的方程:()(1)x x ϕ=从而构成序列:1()0,1,2...(2)k k x x k ϕ+==即给定一个初值0x ,由(2)可算得10()x x ϕ=,再将1x 带入(2)的右端,又可得21()x x ϕ=,…。

我们{k x }为迭代序列,而称(1)式中的()x ϕ为迭代函数,(2)为迭代格式。

如果()x ϕ连续,迭代序列{k x }收敛于*x ,则*x 就是方程(1)的解。

事实上*1()()()lim lim lim k k k k k k x x x x ϕϕϕ+→∞→∞→∞===,又*1,k k x x +→∞→,亦即:**()x x ϕ=或*()0f x =所以,如果迭代序列收敛,总能收敛于原方程的解。

实际计算中,无穷过程不可能实现,只迭代到一定程度,取1k x +作为原方程的近似根。

数值分析导论第二版课程设计

数值分析导论第二版课程设计

数值分析导论第二版课程设计概述数值分析是一门应用数学学科,主要研究计算数学方法和算法,用于求解科学和工程问题中的数学问题。

本课程设计主要参考《数值分析导论第二版》一书,通过对各类数值算法的编程实现,深入理解数值分析中的基本概念和原理。

实验内容本课程设计分为三个部分,每个部分包含若干个小题目,学生需要根据教师布置的题目,自行编写相应的数值算法代码,并进行调试和验证。

第一部分问题描述利用梯形公式、Simpson公式和复合梯形公式计算以下函数的定积分:$$ \\int_0^1 x^2 e^{-x}dx $$其中,梯形公式和Simpson公式的误差要求必须小于10−4,复合梯形公式的误差要求必须小于10−6。

实验要求1.分别编写梯形公式、Simpson公式和复合梯形公式的算法;2.对每个算法进行误差分析,并计算误差;3.输出每个算法的结果以及误差。

第二部分问题描述实现三次样条插值算法,对以下函数进行插值计算:$$ f(x) = \\frac{1}{1+x^2}, ~ x\\in[-5,5] $$实验要求1.编写三次样条插值算法的代码;2.对插值结果进行可视化,画出原始函数及插值函数的曲线图;3.计算插值函数在x=0处的值,并与真实值进行比较。

第三部分问题描述实现共轭梯度法,对以下线性方程组进行求解:$$ \\begin{cases} 3x_1-x_2+x_3=6\\\\ x_1+4x_2-x_3=5\\\\2x_1-x_2+3x_3=2 \\end{cases} $$实验要求1.编写共轭梯度法的代码;2.对求解结果进行可视化,绘制解的分量x1,x2,x3随迭代次数的变化曲线;3.计算解的二范数,与真实解进行比较。

实验环境本课程设计建议使用Python语言进行编程实现,推荐使用Jupyter Notebook进行代码编写和操作演示。

相关库或工具包的版本要求如下:1.Python 3.0及以上版本;2.NumPy 1.0及以上版本;3.Matplotlib 2.0及以上版本;4.Jupyter Notebook5.7及以上版本。

(完整word版)数值分析课程设计实验二

(完整word版)数值分析课程设计实验二

实验二2.1一、题目:用高斯消元法的消元过程作矩阵分解。

设20231812315A ⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦消元过程可将矩阵A 化为上三角矩阵U ,试求出消元过程所用的乘数21m 、31m 、31m 并以如下格式构造下三角矩阵L 和上三角矩阵U(1)(1)212223(2)313233120231,1L m U a a m m a ⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦验证:矩阵A 可以分解为L 和U 的乘积,即A =LU 。

二、算法分析:设矩阵111213212223313233a a a A a a a a a a ⎛⎫ ⎪= ⎪ ⎪⎝⎭,通过消元法可以将其化成上三角矩阵U ,具体算法如下: 第1步消元:111111(1)22112(1)331130,0;;2,3;i i i i i i i i a m a a a a m a i a a m a +=≠⎧⎪=+=⎨⎪=+⎩ 得到111213(1)(1)12223(1)(1)323300a a a A a a a a ⎛⎫ ⎪= ⎪ ⎪⎝⎭第2步消元:(1)(1)(1)32322222(2)(1)(1)333332230,0;;a m a a a a m a ⎧+=≠⎪⎨=+⎪⎩ 得到的矩阵为111213(1)(1)22223(2)33000a a a A a a a ⎛⎫ ⎪= ⎪ ⎪⎝⎭三、程序及运行结果b1.mA=[20 2 3;1 8 1;2 -3 15];for i=1:2M(i)=A(i+1,1)/A(1,1);endfor j=2:3A1(j,2)=A(j,2)-M(j-1)*A(1,2);A1(j,3)=A(j,3)-M(j-1)*A(1,3);endM(3)=A1(3,2)/A1(2,2);A1(3,2)=0;A1(3,3)=A1(3,3)-M(3)*A1(2,3);M,A1运行结果为:M =0.0500 0.1000 -0.4051A1 =0 0 00 7.9000 0.85000 0 15.0443所以:10020230.051007.90.850.10.405110015.0443L U ⎛⎫⎛⎫ ⎪ ⎪== ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭验证:L=[1 0 0;0.05 1 0;0.1 -0.4051 1];U=[20 2 3;0 7.9 0.85;0 0 15.0443];A1=L*UA1 =20.0000 2.0000 3.00001.0000 8.0000 1.00002.0000 -3.0003 15.0000四、精度分析因为根据LU 的递推公式可知,L ,U 分别为下三角和上三角矩阵,其中L 不在对角线上的元素值为111()k ik ik is sk s kk l a l u u -==-∑,在计算每个系数时会产生相应的计算误差。

数值分析第七版课程设计

数值分析第七版课程设计

数值分析第七版课程设计一、实验背景数值分析是计算数学的重要分支,是研究利用计算机求解数值问题的方法和理论的一门学科。

本课程设计旨在通过实验,加深对数值分析相关算法的理解,提高数学建模和计算机编程的能力。

二、实验内容本次课程设计包括以下两个实验:实验一:插值与逼近1.将函数$f(x)=\\dfrac{1}{x}$在区间[1,2]上进行等距节点插值,节点数分别为5、10、15和20,误差使用最大误差和平均误差来比较。

2.使用Newton插值法和Lagrange插值法对于函数$f(x)=\\sin x$进行插值,比较两种方法的误差。

3.对于函数f(x),给定节点x0,x1,x2,x3,计算出f(x)在x=1.5处的三次Hermite插值。

4.对于函数$f(x)=\\dfrac{1}{1+x^2}$,使用最小二乘法对其进行多项式逼近,比较多项式次数为1、2、3和4时的逼近结果。

实验二:数值微积分1.使用五点中心公式,计算f″(x)的近似值,并比较二、四、六、八次公式的精度。

2.使用梯形公式和Simpson公式分别求解函数$f(x)=\\cos(x^2)$在区间[0,1]上的定积分,比较两种方法的精度。

3.使用数值微积分方法计算曲线y=x3+2x+1在区间[0,1]上的弧长,步长分别为0.2、0.1、0.05和0.025,并比较不同步长对计算结果的影响。

三、实验要求1.使用MATLAB或Python等编程语言完成实验,并提交完整的程序代码以及实验报告。

2.实验报告应包括实验的目的、原理、过程、结果及其分析等内容。

3.程序代码应具有较好的结构性、可读性和可复用性,其中涉及到的算法应有详细的注释。

四、实验评分1.实验报告50分,其中内容占30分,格式和排版占20分。

2.程序代码50分,其中正确性占30分,可读性和可复用性占20分。

3.本次课程设计总成绩为实验报告分数和程序代码分数的加权平均分。

数值分析课程设计目录+正文

数值分析课程设计目录+正文

一、目录一、目录 (1)二、设计内容 (2)1.Gauss列主元素消去法解线性方程组的算法设计 (2)1.1 模块设计 (2)1.2主要程序代码 (2)1.3 运行结果 (5)2.Romberg算法 (6)2.1模块设计 (6)2.2 主要程序代码 (6)2.3 运行结果 (7)三、小结 (8)四、参考文献 (8)二、设计内容1.Gauss列主元素消去法解线性方程组的算法设计1.1 模块设计求任意阶(n<10)线性方程组的能力,主要包括一下模块:(1)在源程序中定义求解的方程阶数;(2)输入方程的系数矩阵和常数矩阵;(3)获取增广矩阵(4)Gauss列主元素消去:找主元、消元(分解)、回代、输出结果。

1.2主要程序代码#include <iostream.h>#include <iomanip.h>#include <math.h>#define n 4class Guss{private:float a[n][n];float b[n];float Array[n][n+1];public:Guss() //构造矩阵{cout<<"请输入"<<n<<"*"<<n<<"维系数矩阵:"<<endl;for(int i = 0;i<n;i++) // 输入系数矩阵{for(int j = 0;j<n;j++)cin>>a[i][j];}cout<<"请输入"<<1<<"*"<<n<<"维常数矩阵:"<<endl;for(i = 0;i<n;i++) // 输入常数矩阵{cin>>b[i];}/* for(int i1 = 0;i1<n;i1++){for(int j1 =0;j1<n+1;j1++){cout<<a[i1][j1]<<setw(15);}cout<<b[i1];cout<<endl;}*/}void zgjz() //获得增广矩阵{// float c[n][n+1]; //定义增广矩阵// float *p;// p = &a[0][0];for(int i = 0;i<n;i++) //有系数矩阵和常数矩阵获得增广矩阵205 {for(int j = 0;j<n+1;j++){if(j == n)Array[i][j] = b[i];elseArray[i][j] = a[i][j];}}// return &c[0][0];// return *p; //返回增广矩阵cout<<"增广矩阵为:"<<endl;for(int i1 = 0;i1<n;i1++){for(int j1 =0;j1<n+1;j1++){cout<<Array[i1][j1]<<setw(15);}cout<<endl;}}void Gjfc() //Guss列主元素消元法{float max,u,z; //定义三个float变量,max标记最大值,u用来交换,z用来计算最后的结果int flag; //flag用来标记下标for(int i = 0;i<n-1;i++){max = Array[i][i];flag = i;for(int k = i+1;k<n;k++){if(fabs(max) < fabs(Array[k][i])) //找主元{max = Array[k][i];flag = k;}}if(flag != i){for(int j = i;j<n+1;j++){u = Array[i][j];Array[i][j] = Array[flag][j];Array[flag][j] = u;}}for(int x = i+1;x<n;x++) //消元{u = Array[x][i]/max;for(int y = i;y<n+1;y++){Array[x][y] = Array[x][y]-u*Array[i][y];}}}for(int x=n-1;x>=0;x--) //回代{z=0;for(int y=x+1;y<n;y++)z=z+Array[x][y]*Array[y][n];Array[x][n]=(Array[x][n]-z)/(Array[x][x]);//计算结果}cout<<"Guss消元所得方程组的解为:"<<endl;for(i = 0;i<n;i++) //输出结果cout<<Array[i][n]<<setw(15);cout<<endl;}};//类的结束int main(){Guss G;// float Array[n][n+1];G.zgjz();G.Gjfc();return 0;}1.3 运行结果2.Romberg算法2.1模块设计利用递推的复合梯形公式,并结合外推加速公式的Romberg算法。

数值分析课程设计反幂法

数值分析课程设计反幂法

数值分析课程设计反幂法一、教学目标本节课的教学目标是让学生掌握反幂法的原理及其在数值分析中的应用。

知识目标要求学生了解反幂法的定义、性质及其在求解非线性方程中的应用;技能目标要求学生能够运用反幂法求解实际问题,并能够对结果进行分析和评价;情感态度价值观目标则是培养学生的探究精神、合作意识以及对于数学问题的兴趣。

二、教学内容本节课的教学内容主要包括反幂法的原理及其在数值分析中的应用。

首先,介绍反幂法的定义和性质,通过具体的例子让学生理解反幂法的含义和作用。

然后,讲解反幂法在求解非线性方程中的应用,并通过实际问题让学生掌握反幂法的具体操作步骤。

最后,对反幂法的优缺点进行总结,并引导学生思考如何选择合适的数值方法。

三、教学方法为了达到本节课的教学目标,将采用多种教学方法进行教学。

首先,运用讲授法为学生讲解反幂法的原理和性质,让学生掌握基本知识。

其次,通过讨论法让学生探讨反幂法在实际问题中的应用,培养学生的动手能力和解决问题的能力。

此外,还可以采用案例分析法和实验法,让学生在实际操作中感受反幂法的应用,并能够对其进行评价。

四、教学资源为了支持本节课的教学内容和教学方法的实施,将准备以下教学资源。

首先,教材和相关参考书,为学生提供理论知识的学习;其次,多媒体资料,包括PPT、视频等,用于直观展示反幂法的原理和应用;最后,实验设备,如计算机、计算器等,让学生能够在实际操作中掌握反幂法。

通过这些教学资源的使用,丰富学生的学习体验,提高学生的学习效果。

五、教学评估本节课的教学评估将采取多元化方式进行,以确保评估的客观性和公正性,并全面反映学生的学习成果。

评估方式包括平时表现、作业、小测验和期末考试。

平时表现主要评估学生的课堂参与度和提问回答情况;作业则是对学生掌握反幂法知识的检验,要求学生独立完成相关练习题;小测验则是在课程中期进行,用以检查学生对反幂法的理解和应用能力;期末考试则是对学生整个学期学习成果的全面考核,包括理论知识的理解和实际问题的解决。

合工大数值分析课程设计

合工大数值分析课程设计

合工大数值分析课程设计一、课程目标知识目标:1. 掌握数值分析的基本概念、原理及方法,如插值、数值微积分、常微分方程数值解等;2. 理解数值算法的稳定性、收敛性等性能指标,并能够分析给定数值问题的适用算法;3. 了解数值分析在工程、物理及计算机科学等领域的应用,并能运用所学知识解决实际问题。

技能目标:1. 能够运用数值分析方法解决实际工程问题,具备数值计算编程能力;2. 能够运用所学软件(如MATLAB等)进行数值实验,分析实验结果,优化算法;3. 能够对给定数值问题进行误差分析,提出改进措施,提高计算精度。

情感态度价值观目标:1. 培养学生严谨的科学态度,认识到数值分析在工程技术领域的重要性;2. 激发学生对数值分析的兴趣,培养其主动探索、创新的精神;3. 增强学生的团队协作意识,提高沟通与交流能力。

本课程针对合肥工业大学数值分析课程设计,结合大三年级学生特点,注重理论与实践相结合,培养学生的数值计算能力和实际应用能力。

课程目标旨在使学生在掌握基本理论知识的基础上,能够解决实际问题,提高学生的综合素质,为未来的学术研究或工程实践打下坚实基础。

通过对课程目标的分解,教师可以更好地进行教学设计和评估,确保学生达到预期学习成果。

二、教学内容本章节教学内容主要包括以下几部分:1. 数值分析基本概念:介绍数值分析的定义、研究内容及其在工程中的应用。

- 教材章节:第1章 数值分析引论2. 插值法:讲解拉格朗日插值、牛顿插值、样条插值等方法。

- 教材章节:第2章 插值法3. 数值微积分:介绍数值积分和数值微分的基本原理及方法。

- 教材章节:第3章 数值微积分4. 常微分方程数值解:讲解初值问题和边值问题的数值解法。

- 教材章节:第4章 常微分方程数值解5. 线性方程组的迭代法:介绍雅可比迭代、高斯-赛德尔迭代等方法。

- 教材章节:第5章 线性方程组迭代法6. 数值算法性能分析:分析算法的稳定性、收敛性等性能指标。

数学分析课程设计

数学分析课程设计

实验一1.1水手、猴子和椰子问题:五个水手带了一只猴子来到南太平洋的一个荒岛上,发现那里有一大堆椰子。

由于旅途的颠簸,大家都很疲惫,很快就入睡了。

第一个水手醒来后,把椰子平分成五堆,将多余的一只给了猴子,他私藏了一堆后便又去睡了。

第二、第三、第四、第五个水手也陆续起来,和第一个水手一样,把椰子分成五堆,恰多一只猴子,私藏一堆,再去入睡,天亮以后,大家把余下的椰子重新等分成五堆,每人分一堆,正好余一只再给猴子,试问原先共有几只椰子?试分析椰子数目的变化规律,利用逆向递推的方法求解这一问题。

解:问题分析:假设第n个人分椰子前,有An个椰子,即第n-1醒来面对的椰子数为An-1,第n个水手醒来后,把椰子平分成五堆,将多余的一只给了猴子,他私藏了一堆后便又去睡了则易得:*51A-14Ann+=(1)设最后一堆为a个,则第五个人分完后余下A(6)=5*a+1;由(1)式可得出A(1),A(2)...A(5);由于椰子数目为整数,有a=1开始,建立循环,可以由matlab求得最小的a=1023 ,五次数量变化如下:Matlab窗口输入如下:>> a=1;A=zeros(1,6);while a>0,A(6)=5*a+1;if A(6)~=fix(A(6)),a=a+1;continue;endi=6;while i>=2,A(i-1)=A(i)*5/4+1;if A(i-1)~=fix(A(i-1)),break;endi=i-1;endif (A(1)~=0)&&(A(1)==fix(A(1))),break;enda=a+1;end>> aa =1023>> A(1)ans15621>>plot(A,’k’)1.2当0,1,2,,100n=时,选择稳定的算法计算积分1d10nxn xeI xe--=+⎰解:用蒙特卡罗方法求解该问题,在(0,1)上随机选取一定数量的点,以落在积分区域内的频率近似所求积分,matlab 程序如下: mtj.mfunction a=mtj(n)f=@(x)exp(-n*x)/(exp(-x)+10); N=10000; r=0;for i=1:N,A=rand(1,2); if A(2)<=f(A(1)) r=r+1; end end a=r/N;1.3 绘制静态和动态的Koch 分形曲线问题描述:从一条直线段开始,将线段中间的三分之一部分用一个等边三角形的另两条边代替,形成具有5个结点的新的图形;在新的图形中,又将图中每一直线段中间的三分之一部分都用一个等边三角形的另两条边代替,再次形成新的图形,这时,图形中共有17个结点。

数值分析课程设计95分

数值分析课程设计95分

数值分析课程设计95分一、教学目标本节课的教学目标是让学生掌握数值分析的基本概念和方法,培养学生运用数值分析解决实际问题的能力。

具体目标如下:1.知识目标:(1)了解数值分析的基本概念和研究对象;(2)掌握数值逼近、数值积分和数值解方程等基本方法;(3)了解数值分析在工程和科学计算中的应用。

2.技能目标:(1)能够运用数值分析方法解决实际问题;(2)能够正确选择合适的数值方法并分析其优缺点;(3)能够编写简单的数值计算程序。

3.情感态度价值观目标:(1)培养学生对数值分析的兴趣和好奇心;(2)培养学生勇于探索、严谨治学的科学态度;(3)培养学生团队协作、交流分享的合作精神。

二、教学内容本节课的教学内容主要包括数值分析的基本概念、数值逼近、数值积分和数值解方程等。

具体安排如下:1.数值分析的基本概念:介绍数值分析的研究对象、特点和方法。

2.数值逼近:包括插值法、函数逼近法和数值微积分。

3.数值积分:介绍数值积分的基本方法和误差分析。

4.数值解方程:包括线性方程组的求解、非线性方程的求解和最优化问题。

三、教学方法为了提高学生的学习兴趣和主动性,本节课将采用多种教学方法相结合的方式。

具体方法如下:1.讲授法:用于讲解数值分析的基本概念、原理和方法。

2.案例分析法:通过分析实际案例,让学生了解数值分析在工程和科学计算中的应用。

3.实验法:让学生动手编写数值计算程序,培养实际操作能力。

4.讨论法:学生分组讨论,促进学生间的交流与合作。

四、教学资源为了支持教学内容和教学方法的实施,本节课将采用以下教学资源:1.教材:《数值分析教程》;2.参考书:《数值分析与应用》;3.多媒体资料:数值分析相关视频、动画等;4.实验设备:计算机、编程环境等。

通过以上教学资源的使用,丰富学生的学习体验,提高学生的学习效果。

五、教学评估本节课的评估方式包括平时表现、作业和考试等,以全面反映学生的学习成果。

具体评估方式如下:1.平时表现:通过观察学生在课堂上的参与程度、提问回答等情况,评估学生的学习态度和理解程度。

数值分析原理课程设计

数值分析原理课程设计

数值分析原理课程设计背景数值分析是一门通过数值计算方法来解决数学问题的学科,它涉及到算法的设计和实现,同时结合计算机技术,对各种数学方法进行实现和应用。

在现代科学和工程领域中,数值分析已经成为一种必备的工具,旨在提高数学问题的求解效率和准确性。

本文将介绍数值分析原理的课程设计,并描述其基本原则和实践。

目的本次课程设计的目的是让学生能够掌握数值分析的基本原理和常用算法,实现并验证数学模型。

此外,通过具体的问题案例,帮助学生深刻理解数值分析的重要性和应用价值。

内容数值算法课程设计的第一部分是数值算法。

这一部分的学习内容包括数值微积分、求解线性方程组、差分方程、插值与拟合、非线性方程求解、最小二乘问题等。

其中数值微积分是数值分析的基础,包括求导和积分等,掌握了这个基础之后,学生才能更好地理解后续的数值算法。

数值计算模型的建立第二个部分是关于数值计算模型的建立。

本部分的主要内容是介绍如何构建数学模型和求解数学问题的方法。

具体来说,需要涉及的内容包括数据结构、矩阵的表示、求解常微分方程和偏微分方程等。

在学习这个部分时,学生将会学到如何将实际问题抽象成数学模型,以及如何用数值方法去解决这些问题。

应用实践最后一部分是应用实践。

在这个部分,学生需要根据老师提供的问题案例,自己构建数学模型并使用数值算法来解决问题。

例如,可以探讨一些经济学、物理学或工程学领域的问题,并在此基础上开展课程设计实践。

最终,学生需要提交一个课程设计报告,详细说明数值算法的应用和解决方案,并进一步总结和探讨所学知识。

教学方法本课程设计将采用线上学习和线下实践相结合的教学方法。

学生需要在课堂上自学,并与老师进行必要的讲解和互动。

在本课程设计中,老师将会给出一系列的问题案例,并提供相关的实践指导,鼓励学生们进行自主探究和思考。

结论数值分析原理的课程设计具有一定的难度和挑战,但是对于学生的数学素养和强化实践能力都非常有帮助。

通过学习本课程设计,学生不仅能够掌握数值分析的核心知识和基本技能,同时也能够深入了解数值方法的应用场景和实际应用。

数值分析课程设计课程设计

数值分析课程设计课程设计

数值分析课程设计一、课程设计目的和意义数值分析课程设计是通过选择数值分析中的一些基础算法,设计并编写计算机程序来解决实际的算法问题。

课程设计有助于学生更好地理解和掌握数值分析的基础理论知识,同时对于提高学生的编程能力以及培养学生解决实际问题的能力都具有很大的意义。

二、课程设计内容1.矩阵求解矩阵运算是数值分析中的一项基础知识,但学生在初学阶段往往会遇到矩阵运算方面的问题。

因此,在本课程的矩阵求解部分,学生将会学会如何利用数值分析算法对矩阵进行求解。

2.牛顿迭代法牛顿迭代法是数值分析中常用的一种迭代算法,它可以用来求解函数的根。

在课程设计的牛顿迭代法部分,学生将会深入学习该算法的理论知识,并通过编程实践来深化对其的理解。

3.插值和拟合对于实际问题中的数据,我们需要通过插值和拟合等方法来求取相关的函数。

因此,在课程设计的插值和拟合部分,学生将会学习到常用的插值和拟合算法,并通过实现相应的程序来加深对该算法的理解。

4.数值微积分数值微积分是数值分析中的一项基础知识,它是计算机科学中的一个重要组成部分。

在课程设计的数值微积分部分,学生将会在学习理论知识的基础上,通过编写相应的程序来巩固和加深对该算法的理解。

三、课程设计流程1.熟悉课程设计要求在开始课程设计之前,学生应该熟悉课程设计的要求和流程,明确自己需要完成的任务,并制定相应的计划。

2.确定课程设计题目根据课程设计的要求和个人兴趣,学生可以选择一些自己感兴趣的题目,并请教老师和同学进行相关意见的讨论和确认。

3.学习相关理论知识学生在开始进行课程设计之前,需要对所选择的算法进行深入的学习,并完成必要的理论知识的掌握。

4.开始进行编程在掌握相关的理论知识之后,学生开始进行计算机程序的编写,并不断尝试改进和优化。

5.进行结果验证在完成计算机程序的编写之后,学生需要对其进行一定程度的结果验证,并分析测试结果。

6.撰写课程设计报告在完成验证工作之后,学生需要根据要求撰写课程设计报告,并逐步改进报告的质量和便于理解程度。

《数值分析》课程设计

《数值分析》课程设计

问题的提出3.3 用SOR 方法解下列方程组(去松弛因子w=1.2),要求14||||10k k X X +-∞-<。

12142145x x x x +=⎧⎨-=⎩ 3.4 设 411011A ⎛⎫= ⎪⎝⎭,计算()cond A ∞。

3.5 用选列主元Gauss 消元法求解方程组12312312334721320x x x x x x x x x -+=⎧⎪-+-=-⎨⎪-+=⎩3.6 用追赶法解三对角方程组12345210001121000012100001210000120x x x x x -⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪-- ⎪ ⎪ ⎪⎪ ⎪ ⎪=-- ⎪ ⎪ ⎪-- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭ 3.7 用三角分解法解方程组123248541816862207x x x -⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪--= ⎪⎪ ⎪ ⎪⎪ ⎪--⎝⎭⎝⎭⎝⎭3.8 用选主元消元法计算下列行列式126324951。

一、问题分析1. 超松弛法是迭代方法的一种加速方法,其计算公式简单,但需要选择合适的松弛因子,以保证迭代过程有较快的收敛速度。

2. A 的条件数计算首先要获得A 的逆,而求A 的逆可以转化为求n 个方程组。

3. 完全主元消元法在计算过程中花费了大量的时间用于寻找主元。

同时,各变量的位置在消元过程中也可能会发生变化。

而列选主元法则可消除这个弊病。

4. 追赶法主要是解三对角方程组。

所谓追指消元过程,赶指回代过程。

5. Gauss 消元法是通过逐步消元过程,将方程组的系数矩阵A 转变为一个上三角矩阵。

三角分解法,就是把系数矩阵分解为两个三角阵。

6.将某一向量坐标同乘以某非零实数,加到另一向量上,行列式的值不变。

用选主元法将行列式矩阵变为三角阵,对角线上的数值相乘即为行列式的值。

二、编程解决3.3Sor法c语言编程:#include<stdio.h>#include<math.h>#include<stdlib.h>#define omega 1.2 //取值不合适结果可能发散void main(){double a[5][5];double b[5],x[5],f,t,y[5]={0,0,0,0,0};int i,j,n,cnt=0;printf("阶数:");scanf("%d",&n);printf("请输入%d阶的A矩阵\n",n);for(i=0;i<n;i++)for(j=0;j<n;j++)scanf("%lf",&a[i][j]);printf("请输入B矩阵\n");for(i=0;i<n;i++)scanf("%lf",&b[i]);printf("count\t");for(i=0;i<n;i++)printf("x[%d]\t\t",i);printf("收敛程度\n");do{for(i=0;i<n;i++)x[i]=y[i];for(i=0;i<n;i++){t=0;for(j=0;j<n;j++)t=t+a[i][j]*(j<i?y[j]:x[j]);y[i]=x[i]+omega*(b[i]-t)/a[i][i];printf("%d",cnt++);for(i=0;i<n;i++)printf("\t%lf",x[i]);f=0;for(i=0;i<n;i++)f+=fabs(y[i]-x[i]);printf("\t%g\n",f);}while(f>1e-4 && cnt<100);}所得结果:3.4 求逆、算条件数编程:#include <stdio.h>#include <math.h>#include <stdlib.h>#define N 5 //可修改,以改变可解决的最大维数。

12级数值分析课程设计

12级数值分析课程设计

12级数值分析课程设计一、课程目标知识目标:1. 理解数值分析的基本概念,掌握数值计算的基本方法,如插值、数值微积分、常微分方程数值解等;2. 掌握误差分析的基本理论,了解数值稳定性和收敛性的概念;3. 掌握线性代数、微积分等数学基础知识在数值分析中的应用。

技能目标:1. 能够运用数值分析方法解决实际问题,如求解非线性方程、线性方程组、优化问题等;2. 能够运用编程语言(如MATLAB、Python等)实现数值算法,并进行调试和优化;3. 能够运用误差分析理论评估数值算法的准确性和稳定性。

情感态度价值观目标:1. 培养学生对数值分析的兴趣和热情,激发学生主动探索数值方法解决实际问题的欲望;2. 培养学生的团队协作意识,学会与他人合作共同解决问题;3. 培养学生的创新精神和批判性思维,敢于对现有数值方法提出质疑,勇于尝试改进和创新。

课程性质:本课程为专业基础课,旨在培养学生运用数值方法解决实际问题的能力。

学生特点:12级学生已具备一定的数学基础和编程能力,具有较强的逻辑思维和分析能力。

教学要求:结合课程性质和学生特点,注重理论与实践相结合,强化上机实践和案例分析,提高学生的实际操作能力和问题解决能力。

通过本课程的学习,使学生能够熟练运用数值分析方法,为后续专业课程学习和未来从事相关工作打下坚实基础。

二、教学内容1. 数值分析基本概念:包括数值计算的误差、稳定性、收敛性等;教材章节:第一章 数值分析引论内容:误差分析、稳定性与收敛性、数值方法的分类与特点。

2. 数值微积分:包括数值积分和数值微分;教材章节:第二章 数值微积分内容:梯形公式、辛普森公式、高斯公式、数值微分方法。

3. 插值与逼近:包括拉格朗日插值、牛顿插值、样条插值等;教材章节:第三章 插值与逼近内容:线性插值、多项式插值、样条插值、最佳逼近问题。

4. 解非线性方程:包括二分法、牛顿法、弦截法等;教材章节:第四章 解非线性方程内容:迭代法、牛顿法、弦截法、非线性方程组的求解。

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

课程 设 计我再也回不到大二了,大学是那么短暂设计题目 数值分析 学生姓名 李飞吾 学 号 x x x x x x x x 专业班级 信息计x x x x x 班 指导教师设 计 题 目共15题如下成绩数值分析课程设计1.1 水手、猴子和椰子问题:五个水手带了一只猴子来到南太平洋的一个荒岛上,发现那里有一大堆椰子。

由于旅途的颠簸,大家都很疲惫,很快就入睡了。

第一个水手醒来后,把椰子平分成五堆,将多余的一只给了猴子,他私藏了一堆后便又去睡了。

第二、第三、第四、第五个水手也陆续起来,和第一个水手一样,把椰子分成五堆,恰多一只猴子,私藏一堆,再去入睡,天亮以后,大家把余下的椰子重新等分成五堆,每人分一堆,正好余一只再给猴子,试问原先共有几只椰子?(15621) 试分析椰子数目的变化规律,利用逆向递推的方法求解这一问题 解:算法分析:解该问题主要使用递推算法,关于椰子数目的变化规律可以设起初的椰子数为0p ,第一至五次猴子在夜里藏椰子后,椰子的数目分别为01234,,,,p p p p p 再设最后每个人分得x 个椰子,由题:14(1)5k k p p +=- (k=0,1,2,3,4)51(1)5x p =-所以551p x =+,11k k p p +=+利用逆向递推方法求解151,4k k p p +=+ (k=0,1,2,3,4)MATLAB 代码: n=input('n= '); n= 15621 for x=1:n p=5*x+1; for k=1:5 p=5*p/4+1; endif p==fix(p), break end enddisp([x,p])1.2 设,15nn x I dx x=+⎰ (1)从0I 尽可能精确的近似值出发,利用递推公式:115(1,2,20)n n I I n n-=-+=计算机从1I 到20I 的近似值;(2)从30I 较粗糙的估计值出发,用递推公式:111(30,29,,3,2)55n n I I n n-=-+=计算从1I 到20I 的近似值;解:首先第一步,估计0I 和30I 的值:syms x n;int (x^0/(5+x),0,1) ans=log(2)+log(3)-log(5) eval(ans) ans= 0.1823则取0I 为0.18 syms x n;int(x^30/(5+x),0,1) ans =931322574615478515625*log(2)+931322574615478515625*log(3)-931322574615478515625*log(5)-79095966183067699902965545527073/465817912560 eval(ans) ans = 0MATLAB 中小数点后保留四位,由上面计算知道积分的值不为了零。

所以30I 的取值为0.00001-0.0001 MATLAB 代码: i=input('i='); i=0.18;>> if i>=0.1&&i<=0.2 for n=1:1:20 i=-5*i+1/n endelseif i>0&&i<=0.0001 for n=30:-1:2i=(-1/5)*i+1/(5*n) end end i =1.1336e+005 i =-5.6679e+005 i =2.8339e+006 i =-1.4170e+007 i =7.0848e+007 i =-3.5424e+008 i =1.7712e+009 i =-8.8560e+009 i =4.4280e+010 i =-2.2140e+011同理输入积分初始值i=0时可以得 i=0.0884结果分析:第二种方法所得的结果相对来说比较精确一些,也比较可靠因为第一种方法每一迭代都将最初的误差放大了五倍,使得最终的误差越来越大;而第一种方法经过每一次迭代都将误差缩小为初始误差的五分之一,使得最终的误差越来越小,因此相对来说比较可靠,性能较好。

1.3 绘制Koch 分形曲线问题描述:从一条直线段开始,将线段中间的三分之一部分用一个等边三角形的另两条边代替,形成具有5个结点的新的图形(图1-4);在新的图形中,又将图中每一直线段中间的三分之一部分都用一个等边三角形的另两条边代替,再次形成新的图形(图1-5),这时,图形中共有17个结点。

这种迭代继续进行下去可以形成Koch 分形曲线。

问题分析:考虑由直线段(2个点)产生第一个图形(5个点)的过程,设1P 和5P 分别为原始直线段的两个端点。

现在需要在直线段的中间依次插入三个点234,,P P P 产生第一次迭代的图形(图1-4)。

显然,2P 位于1P 点右端直线段的三分之一处,4P 点绕2P 旋转60度(逆时针方向)而得到的,故可以处理为向量24P P 经正交变换而得到向量23P P ,形成算法如下: (1)2151()/3P P P P =+-; (2)41512()/3P P P P =+-;(3)T3242()P P P P A =+-⨯;在算法的第三步中,A 为正交矩阵。

这一算法将根据初始数据(1P 和5P 点的坐标),产生图1-4中5个结点的cos sin 33sin cos 33A ππππ⎡⎤-⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦坐标。

这5个结点的坐标数组,组成一个5×2矩阵。

这一矩阵的第一行为为1P 的坐标,第二行为1P 的坐标,第二行为2P 的坐标……第五行为5P 的坐标。

矩阵的第一列元素分别为5个结点的x 坐标 ,第二列元素分别为5个结点的y 坐标。

问题思考与实验:(1)考虑在Koch 分形曲线的形成过程中结点数目的变化规律。

设第k 次迭代产生结点数为k n ,第1k +迭代产生结点数为1k n +,试写出k n 和1k n +之间的递推关系式;(2)参考问题分析中的算法,考虑图1-4到图1-5的过程,即由第一次迭代的5个结点的结点坐标数组,产生第二次迭代的17个结点的结点坐标数组的算法;(3)考虑由第k 次迭代的k n 个结点的结点坐标数组,产生第1k +次迭代的1k n +个结点的结点坐标数组的算法;(4)设计算法用计算机绘制出如下的Koch 分形曲线(图1-6) 解:(1) 143k k n n +=- (2)(3)算法及(4)代码分析:p=[0 0;10 0]; %P 为初始两个点的坐标,第一列为x 坐标,第二列为y 坐标 n=2; %n 为结点数A=[cos(pi/3) -sin(pi/3);sin(pi/3) cos(pi/3)]; %旋转矩阵 for k=1:4d=diff(p)/3; %diff 计算相邻两个点的坐标之差,得到相邻两点确定的向量%则d 就计算出每个向量长度的三分之一,与题中将线段三等分对应 m=4*n-3; %迭代公式q=p(1:n-1,:); %以原点为起点,前n-1个点的坐标为终点形成向量p(5:4:m,:)=p(2:n,:); %迭代后处于4k+1位置上的点的坐标为迭代前的相应坐标 p(2:4:m,:)=q+d; %用向量方法计算迭代后处于4k+2位置上的点的坐标 p(3:4:m,:)=q+d+d*A'; %用向量方法计算迭代后处于4k+3位置上的点的坐标 p(4:4:m,:)=q+2*d; %用向量方法计算迭代后处于4k 位置上的点的坐标 n=m; %迭代后新的结点数目 endplot(p(:,1),p(:,2)) %绘出每相邻两个点的连线 axis([0 10 0 3.5])2.1 用高斯消元法的消元过程作矩阵分解。

设20231812315A ⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦消元过程可将矩阵A 化为上三角矩阵U ,试求出消元过程所用的乘数21m 、31m 、31m 并以如下格式构造下三角矩阵L 和上三角矩阵U(1)(1)212223(2)313233120231,1L m U a a m m a ⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦验证:矩阵A 可以分解为L 和U 的乘积,即A=LU 。

矩阵LU 分解MATLAB 代码: function hl=zhjLU(A)[n,n]=size(A);RA=rank(A); if RA~=ndisp('因为A 的n 阶行列式hl 等于零,所以A 不能进行LU 分解.A 的秩RA 如下:'); RA,hl=det(A); return endif RA==nfor p=1:nh(p)=det(A(1:p,1:p)); endhl=h(1:n); for i=1:nif h(1,i)==0disp('因为A 的各阶主子式等不等于零,所以A 能进行LU 分解.A 的秩RA 和各阶顺序主子式值hl 依次如下:'); RA,hl return end endif h(1,i)~=0disp('因为A 的各阶主子式都不等于零,所以A 能进行LU 分解.A 的秩RA 和各阶顺序主子式值hl 如下:');for j=1:nU(1,j)=A(1,j);endfor k=2:n for i=2:n for j=2:nL(1,1)=1;L(i,i)=1; if i>jL(1,1)=1;L(2,1)=A(2,1)/U(1,1);L(i,1)=A(i,1)/U(1,1);L(i,k)=(A(i,k)-L(i,1:k-1)*U(1:k-1,k))/U(k,k); elseU(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j); end end end endRA,hl,U,L · end end以上上代码保存为M 文件,并在命令窗口输入 A=[20 2 3;1 8 1; 2 -3 15]; b=[0 0 0]'; h=zhjLU(A)2.2 用矩阵分解方法求上题中A 的逆矩阵。

记1231000,1,0001b b b ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦ 分别求解方程组 123,,AX b AX b AX b ===由于三个方程组系数矩阵相同,可以将分解后的矩阵重复使用。

对第一个方程组,由于A=LU ,所以先求解下三角方程组1b =LY ,再求解上三角方程组=UX Y ,则可得逆矩阵的第一列列向量;类似可解第二、第三方程组,得逆矩阵的第二列列向量的第三列列向量。

由三个列向量拼装可得逆矩阵1-A 。

解:MATLAB 代码如下:b1=[1;0;0]; b2=[0;1;0]; b3=[0;1;1]; A=[20,2,3;1,8,1;2,-3,15];L=[1,0,0;0.05,1,0;0.1,-0.4051,1]; U=[20 2 3;0 7.9 0.85;0 0 15.0443]; Y1=L\b1X1=U\Y1 Y2=L\b2 X2=U\Y2 Y3=L\b3 X3=U\Y3 Y1 =1.0000 -0.0500 -0.1203 X1 =0.0517 -0.0055 -0.0080 Y2 =0 1.0000 0.4051 X2 = -0.0164 0.1237 0.0269 Y3 = 0 1.0000 1.4051 X3 =-0.0257 0.11650.09342.3 验证希尔伯特矩阵的病态性:对于三阶矩阵11/21/31/21/31/41/31/41/5H ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦取右端向量T[11/613/1247/60]b =,验证:(1)向量T T123[][111]X x x x ==是方程组b =HX 的准确解;(2)取右端向量b 的三位有效数字得T[1.83 1.080.783]b =,求方程组的准确解*X ,并与X 的数据[1,1,1]T作比较 。

相关文档
最新文档