东南大学数值分析编程作业-20211119001719.docx

Word Word Word Word 数值分析编程作业 姓 名：卢一鸣 学 号：151106 Word nJttiS Word nJttiS Word PjttiW 指导老师：吴宏伟 原创力文档 原创力文档 下载高淸无水印 目录 P2。.第17题：舍入误差与有效数 3 P56.第20题：Newton迭代法 5 P126.第39题：列主元Gauss消去法13P127 .第40 P126.第39题：列主元Gauss消去法 13 P127 .第40题：逐次超松弛迭代法 16 P195.第37题：3次样条插值函数 21 P256.第23题：重积分的计算 23 P307.第23题：常微分方程初值问题数值解 25 P346.第 10 题：抛物方程 Crank-Nicolson 格式 30注，源程序放在对应的文件央里.如第一题放在“P20.第17題：舍入误差与有效数?文件 央里 P20.第17题：舍入误差与有效数 设京吉,其精确值为找土土 ⑴编制按从小到大的顺序s* =—!一十?4-+??…+ N 22-1 31 」一 ?计算S、,的通用程序： N2-l ? (2)编制按从小到大的顺序5“=白+浦=+. NF k W-1 -1 ????+二，计算L的通用程序: 2 -1 按两种顺序分别计算鸟砂，S时.并指出方效位数(编储程序时用单粘度)： 通过本上机題你明白了什么？ 答 ⑴⑵程序：程序放在Computation.m, goon.m文件里，运行程序时只需输入命令main 即叽 %此函数为main函散 缶此函数笫一个For循环N按从小到大顺序 disp ((blanks (0) z * 4^ S。神 SNtfi： 1 J); disp ((blanks (10),?N 按从小到大 J); disp ([blanks (10),，N 按从大到小?]); disp ([blanks (10) F f 真实值,]); for N-[100 10000 1000000] w=sprintf(%s%d,?N=,N); disp([blanks(15)Fw]); computation(N); Word Word Word Word 原创力文档 原创力文档 max.bookl18kom 下莪高清无 end goon (; function [ a ] = computation( N ) 吿本函数订算三种SN值：N按从小到大，N按从 大到小.真实值 半第二个for循环N按从大到小顺序 愚菊二个fodg孙N按从大到小顺序 SN=0; SN?single(SN); +N从小到大的顺序 for i-2:N SN-SN+1/; end a-sprintf CtS^Of^SN); disp{(blanks(10),a]); SN-0; SN-single(SN); 冬N从大到小的顺序 for i-N:-l:2 SN-SN+l/(iai-l); end a-sprintf(* %5.20f?r SN); disp([blanks(10)fa]); ，氏实值 b-0.5*(1.5-1/N-1/(N+1)); function [ a j ? goon() WUNTITLED3 Summary of Chis function goes here % Detailed explanation goes here while 1 disp (现在你可以尝试其他的Nffi,如想结 束程序?清输入N?1 -)； N=inpuW 诂输入 N： ?); if N=1 return; else w=sprintf ?N=,N); disp([blanks (15)rw]); computation(N); end end Word Word nJtftW Word Word nJtftW Word nJttiS Word nJttiS b-sprintf(?%5.20f?rSN); disp {[blanks(10)9b]); end （3） 运行结果， ? main 本函数三种SN值： N按M小到大 N按从大到小 直实值 N=100 0. 74004948139190674000 0.74004954099655151000 0.74004954099655151000 N= 10000 0.74985212087631226000 0.74989998340606689000 0.74989998340606689000 N= 1000000 0.74985212087631226000 0.74999898672103882000 0.74999898672103882000 现在你可以尝试其他的N鳳如想结束程序，请输入NSl 清输入N： 4 N=4 0,525