




全文预览已结束
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1矩阵与数值分析上机实习1. 考虑计算给定向量的范数:输入向量x = (x1; x2; ; xn)T ,输出 x 1, x 2, x 1。请编制一个通用程序,并用你编制的程序计算如下向量的范数:x = (1;111)T; ; y = (1; 2; ; n)T :23n对n = 10,100,1000甚至更大的n计算其范数,你会发现什么结果?你能否修改你的程序使得计算结果相对精确呢?2. 考虑y = f(x) = ln(1+x),其中定义f(0) = 1,此时f(x)是连续函数。用此公式计算x当x10 15; 10 15时的函数值,画出图像。另一方面,考虑下面算法:d = 1 + xif d = 1 then y = 1elsey = ln d=(d 1) end if用此算法计算x10 15; 10 15时的函数值,画出图像。比较一下发生了什么?3. 首先编制一个利用秦九韶算法计算一个多项式在给定点的函数值的通用程序,你的程序包括输入多项式的系数以及给定点,输出函数值。利用你编制的程序计算p(x) = (x 2)9 = x9 18x8 + 144x7 672x6 + 2016x5 4032x4 + 5376x3 4608x2 + 2304x 512在x = 2邻域附近的值。画出p(x)在x1:95; 20:5上的图像。4. 编制计算给定矩阵A的LU 分解和P LU 分解的通用程序,然后用你编制的程序完成下面两个计算任务:(1) 考虑10. . 01. .n n;A =.1 . . . . . .1111R.011111 自己取定x Rn,并计算b = Ax。然后用你编制的不选主元和列主元的Gauss消去法求解该方程组,记你计算出的解为x。对n从5到30估计计算解的精度。(2) 对n从5到30计算其逆矩阵。25. 编制计算对称正定阵的Cholesky分解的通用程序,并用你编制的程序计算Ax = b,其中A = (aij )Rn n,aij =1。b可以由你自己取定,对n从10到20验证程序的可靠性。i+j16. (1) 编制程序House(x),其作用是对输入的向量x,输出单位向量u使得(I2uuT )x =x 2e1。(2)编制Householder变换阵H = I 2uuT Rn n乘以A Rn m的程序HA,注意,你的程序并不显式的计算出H。(3)考虑矩阵12341323A =22e;102730275=2用你编制的程序计算H使得HA的第一列为 e1的形式,并将HA的结果显示。7. 用Jacobi和Gauss-Seidel迭代求解下面的方程组,输出迭代每一步的误差 xk x:5x1x2 3x3 = 2x1+ 2x2 + 4x3 = 13x1+ 4x2 + 15x3 = 108. 取不同的初值用Newton迭代以及弦截法求方程x3 + 2x2 + 10x 100 = 0的实根,列表或者画图说明收敛速度。9. 用二分法求方程ex cos x + 2 = 0在区间0; 4 上的所有根。10. 考虑函数f(x) = sin( x),x 0; 1。用等距节点作f(x)的Newton插值,画出插值多项式以及f(x)的图像,观察收敛性。11. 对函数f(x) = 1+1x2 ,x 5; 5,取不同的节点数n,用等距节点作lagrange插值,观察Runge现象。12. 令f(x) = e3x cos( x),考虑积分02 f(x)dx。区间分为50,100,200,500,1000等,分别用复合梯形以及复合Simpson积分公式计算积分值,将数值积分的结果与精确值比较,列表说明误差的收敛性。13. 分别用2点,3点以及5点的Gauss型积分公式计算如下定积分:1x2sin x2(1)dx;(2)0dx:1x1x214. 考虑微分方程初值问题:dx1dt=(ty y2);(t+1)2x(0) = 2:分别用Euler法,改进的Euler法,Runge-Kutta法求解该方程。分别取步长为0:1,0:01,0:001,计算到x(1),画图说明结果。3要求:1. 考试前提交作业(以Word形式提交结果,按照题目要求提交代码以及数值结果),主题写学号+姓名,发送至邮箱:张宏伟老师:zhuke 2015163.com孟兆良老师
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- DB14-T 1636-2025 甜樱桃蜜蜂授粉技术规程
- 拆分协议书:企业分家及业务拆分合同
- 车辆损伤赔偿与车辆维修进度协议
- 汽车后市场车辆挂靠业务合作协议
- 2025年初中物理八年级下册(沪科版)教学课件 第十章 第一节
- 2025年管理学选修课期末考试试题及答案
- 标杆企业高速公路工程施工标准化手册
- 智能制造厂房承包经营与人才培训合同
- 原粮仓储仓单质押贷款服务合同
- 场监管总局全国劳动合同示范文本库(外国子公司)
- TYNZYC 0104-2023 云南省中药材追溯体系 第4部分:追溯码编码规范
- 玉米套种红薯技术方案
- 楼梯踏步玻璃槽的施工方案
- 11测量放线复核记录报表(全自动计算-只需要输入坐标)
- 木偶奇遇记 阅读单
- 华为云桌面解决方案介绍
- 新时代中国特色社会主义理论与实践(2021版)课后思考题答案
- 骨痿临床路径及表单
- 上海市闵行区2022-2023学年高一下学期期末数学试题(无答案)
- 北京市智慧社区建设指导标准
- 2023年郴州市发改系统政策业务知识抢答赛题库
评论
0/150
提交评论