计算方法上机作业.docxVIP

目录 TOC \o 1-5 \h \z \o Current Document 计算方法A上机作业 1 上机练习目的 1 上机练习任务 1 计算方法A上机题目 1 程序设计要求 1 上机报告要求 1 \o Current Document QR分解法求解线性方程组 2 计算原理 2 程序框图 7 计算实习 8 Matlab 代码 8 \o Current Document 共轨梯度法求解线性方程组 10 计算原理 10 程序框图 11 计算实习 12 Matlab代码 12 \o Current Document 三次样条插值 14 计算原理 14 程序框图 16 计算实习 17 Matlab代码 17 \o Current Document 四阶龙格一库塔法求解常微分方程的初值问题 21 计算原理 21 程序框图 22 计算实习 23 Matlab 代码 23 计算方法A上机作业 上机练习目的 □ 复习和巩固数值计算方法的基本数学模型，全面掌握运用计算 机进行数值计算的具体过程及相关问题。 □利用计算机语言独立编写、调试数值计算方法程序，培养学牛 利用计算机和所学理论知识分析解决实际问题的能力。 上机练习任务 ?利用计算机语言编写并调试一系列数值方法计算通用程序，并 能正确计算给定题目，掌握调试技能。 ?掌握文件使用编程技能，如文件的各类操作，数据格式设计、 通用程序运行过程中文件输入输岀运行方式设计等。 ?写出上机练习报告。 计算方法A上机题目 QR分解方法求解线性方程组。（第二章） 共辘梯度法求解线性方程组。（第三章） 三次样条插值（第四章） 四阶龙格一库塔法求解常微分方程的初值问题 程序设计要求 程序要求是通用的，在程序设计吋要充分考虑哪些变量应该可 变的。 程序要求调试通过。 上机报告要求 报告内容包括： ? 每种方法的算法原理及程序框图。 ?程序使用说明。 ?算例计算结果。 QR分解法求解线性方程组 计算原理 当xeRft是任意给定的非零向量，veRfl是任意给定的单位向量，则存在初 等反射阵—几使得其如为常数，当取单位向量 时，由u确定的矩阵H必定满足Hx = G ,所以在计算过程屮取u的值为上述值。 设A是一个mx/7（mn）阶矩阵且它的列向量线性无关，则利用豪斯霍尔徳 变换可以把A逐步化为上梯形矩阵, a2n %丿 具体变换过程如下: 设q（i = l,2,???/）是m维单位坐标向量。 为把矩阵A的第一列q⑼化为（0,0,???,0）〃，取 x = ax⑼,v = q =（1,0,…,0）丁 w/?（或取卩=一弓），根据上式可得，取 其中 口0] —4弓「， 0\ -吗⑼ 一 九』_2吶_ _嗣 （ （0）\ 2““ -硕^不两-盲’） 令 H{=Im- 2u{u^ = Im - 2a；xa）^ ,用耳左乘 A（°）得, A⑴=耳屮）=（恥化厲才），V =（5弓4⑴4⑴,…，色⑴） 由此可知6 o? ? ?(113(123 a a(1 (2 a aVI/ 3aHnunXIV 1 /\ 2 a aa, =Hia『)= 由此可知 6 o ? ? ? (113(123 a a (1 (2 a a VI/ 3 a Hnun XIV 1 /\ 2 a a a, =Hia『)=(Zm - af 1wiw^)ajo) = aj0) — af1w7aj0)wi = a『)- 0ijWi = 4°)- 0i，(af)- 66) f a\j - 01 j(flu - 6)\ a2j - 0lja2l amj - ^baml / aij = a缪- 0血￥ - 6), j = 2,3,…,n, aij}=增 _0好朗,: = 2,3,??- ,m;J = 2,3,???“? 0=町1诃4°）=町半）一 o】6）Tq『） =町诃网）一赵））=时（孕碍-剧? 2,3,…， 2,3,…，n. i=2 第 2 步记 a$）= （/；）￡）,???,a舄）t = （aS?,a^l）T）T,显然* 0?取x = a^\v = ￡）（或取v = ））,其中f）为勺=（0,冷）t）t中的m-1维单位 坐标向址.这时，取 _ 玄$）- 72密） W2 血=11訝）一才）||2=両？ 其中 72 = -sgn（a￡））||a$）||2, w2 =叫）一 a2e^\ 而 11切||2 = |闘）- cr2e[1}h = \A（疋一巾遢））, 其中构造c T 2w2wJW2WJT他2奶OC2g =。2（6 - 其中 构造 c T 2w2wJ W2WJ T他2奶 OC2 g =。2（6 - a￡＞）? H2 = /m-1 - 2u2uJ = Zm-1 一 02 lw2w2 H2 = 1 0T 0 h2 用H2左乘4⑴得 4