《数值分析报告上机题(matlab版)(东南大学)》由会员分享，可在线阅读，更多相关《数值分析报告上机题(matlab版)(东南大学)（23页珍藏版）》请在金锄头文库上搜索。
1、word 数值分析上机报告第一章一、题目准确值为。1) 编制按从大到小的顺序，计算SN的通用程序。2) 编制按从小到大的顺序，计算SN的通用程序。3) 按两种顺序分别计算,并指出有效位数。编制程序时用单精度4) 通过本次上机题，你明白了什么？二、通用程序clearN=input(Please Input an N (N1):);AccurateValue=single(0-1/(N+1)-1/N+3/2)/2);Sn1=single(0);for a=2:N; Sn1=Sn1+1/(a2-1);endSn2=single(0);for a=2:N; Sn2=Sn2+1/(N-a+2)2-1);
2、endfprintf(The value of Sn using different algorithms (N=%d)n,N);disp(_)fprintf(Accurate Calculation %fn,AccurateValue);fprintf(Caculate from large to small %fn,Sn1);fprintf(Caculate from small to large %fn,Sn2);disp(_)三、求解结果Please Input an N (N1):102The value of Sn using different algorithms (N=100
3、)_Caculate from large to small Caculate from small to large _Please Input an N (N1):104The value of Sn using different algorithms (N=10000)_Caculate from large to small Caculate from small to large _Please Input an N (N1):106The value of Sn using different algorithms (N=1000000)_Caculate from large
4、to small Caculate from small to large _四、结果分析有效位数 n 顺序 100 100001000000从大到小633从小到大566可以得出，算法对误差的传播又一定的影响，在计算时选一种好的算法可以使结果更为准确。从以上的结果可以看到从大到小的顺序导致大数吃小数的现象，容易产生较大的误差，求和运算从小数到大数算所得到的结果才比拟准确。第二章一、题目1给定初值与容许误差，编制牛顿法解方程f(x)=0的通用程序。2给定方程,易知其有三个根a) 由牛顿方法的局部收敛性可知存在当时，Newton迭代序列收敛于根x2*。试确定尽可能大的。b)试取假如干初始值，观察当
5、时Newton序列的收敛性以与收敛于哪一个根。3通过本上机题，你明白了什么？二、通用程序%寻找最大的delta值%clear%flag=1;k=1;x0=0;while flag=1 delta=k*10-6; x0=delta; k=k+1; m=0; flag1=1; while flag1=1 & m=103 x1=x0-fx(x0)/dfx(x0); if abs(x1-x0)=10-6 flag=0; endendfprintf(The maximun delta is %fn,delta); % 定义函数f(x)function Fx=fx(x) Fx=x3/3-x;% 定义导函数
6、df(x)function Fx=dfx(x) Fx=x2-1;% Newton法求方程的根%clear%ef=10-6; %给定容许误差10-6k=0;x0=input(Please input initial value Xo:);disp(k Xk);fprintf(0 %fn,x0); flag=1;while flag=1 & k=103 x1=x0-fx(x0)/dfx(x0); if abs(x1-x0)ef flag=0; end k=k+1; x0=x1;fprintf(%d %fn,k,x0); end 三、求解结果结果为：即得最大的为0.774597，Newton迭代序列
7、收敛于根=0的最大区间为-0.774597，0.774597。在区间上各输入假如干个数，计算结果如下：区间Please input initial value Xo:-30k XkPlease input initial value Xo:-10k XkPlease input initial value Xo:-10000k XkPlease input initial value Xo:-100k XkPlease input initial value Xo:-50k XkPlease input initial value Xo:-3k Xkk XkPlease input initi
8、al value Xo:-8k XkPlease input initial value Xo:-7k XkPlease input initial value Xo:-5k Xk结果显示，以上初值迭代序列均收敛于-1.732051，即根。在区间即区间-1，-0.774597上取-0.774598，-0.8，-0.85，-0.9，-0.99，计算结果如下：k Xkk Xkk Xkk Xkk Xk计算结果显示，迭代序列局部收敛于-1.732051，即根，局部收敛于1.730251，即根。在区间即区间-0.774597，0.774597上，由search.m的运行过程明确，在整个区间上均收敛于0，即根。k Xkk Xkk Xkk Xkk Xk在区间即区间0.774597，1上取0.774598，0.8，0.85，0.9，0.99，计算结果如下：计算结果显示，迭代序列局部收敛于-1.732051，即根，局部收敛于1.730251，即根。Please input initial value Xo:4k XkPlease input initial value Xo:3k Xkk Xk区间上取1
