1、专业 序号 姓名 日期 实验 1 算法的数值稳定性实验【实验目的】1掌握用 MATLAB 语言的编程训练，初步体验算法的软件实现；2通过对稳定算法和不稳定算法的结果分析、比较，深入理解算法的数值稳定性及其重要性。【实验内容】1 计算积分 (n=0,1,2，10) dxaIn10)(其中 a 为参数，分别对 a=0.05 及 a=15 按下列两种方案计算，列出其结果，并对其可靠性，说明原因。2 方案一 用递推公式 (n=1,2,10) naI1n递推初值可由积分直接得 )(0I3. 方案二 用递推公式 (n=N,N-1,1) )1(1-nnIaI根据估计式 当I 1na或当n11Ia 1na0取
2、递推初值为当12)1(12N NaNaaI 1aN或当121NaIN 1a0N计算中取 N=13 开始【程序如下】:% myexp1_1.m - 算法的数值稳定性实验% 见 P11 实验课题(一) %function yyjjglobal n aN = 20; % 计算 N 个值a =0.05;%或者 a=15% %-% % 方案 I 用递推公式 %I(k) = - a*I(k-1) + 1/k% I0 =log(a+1)/a); % 初值I = zeros(N,1); % 创建 N x 1 矩阵(即列向量),元素全为零I(1) =-a*I0+1;for k = 2:NI(k) =-a*I(k
3、-1)+1/k;end% %-% % 方案 II 用递推公式%I(k-1) = ( - I(k) + 1/k ) / a% II = zeros(N,1);if a = N/(N+1)II(N)=(2*a+1)/(2*a*(a+1)*(N+1);elseII(N) =(1/(a+1)/(N+1)+1/N)/2;endfor k = N:-1:2II(k-1) =(-II(k)+1/k)/a;end% %-% % 调用 matlab 高精度数值积分命令 quadl 计算以便比较III = zeros(N,1);for k = 1:Nn = k;III(k) = quadl(f,0,1);end%
4、 %-% % 显示计算结果clc fprintf(n 方案 I 结果 方案 II 结果 精确值)for k = 1:N,fprintf(nI(%2.0f) %17.7f %17.7f %17.7f,k,I(k),II(k),III(k)end% %-function y = f(x) % 定义函数global n a % 参量 n 为全局变量y =x.n./(a+x); % 注意:这里一定要 点 运算return% %-【运行结果如下】：当 a=0.05方案 I 结果 方案 II 结果 精确值I( 1) 0.8477739 -919648916620722180000.0000000 0.84
5、77739I( 2) 0.4576113 45982445831036109000.0000000 0.4576113I( 3) 0.3104528 -2299122291551805700.0000000 0.3104528I( 4) 0.2344774 114956114577590290.0000000 0.2344776I( 5) 0.1882761 -5747805728879515.0000000 0.1882761I( 6) 0.1572529 287390286443975.9400000 0.1572529I( 7) 0.1349945 -14369514322198.654
6、0000 0.1349945I( 8) 0.1182503 718475716110.0577400 0.1182503I( 9) 0.1051986 -35923785805.3917770 0.1051986I(10) 0.0947401 1796189290.3695889 0.0947401I(11) 0.0861721 -89809464.4275704 0.0861724I(12) 0.0790247 4490473.3047119 0.0790247I(13) 0.0729718 -224523.5883125 0.0729718I(14) 0.0677800 11226.250
7、8442 0.0677800I(15) 0.0632777 -561.2458755 0.0632777I(16) 0.0593361 28.1247938 0.0593361I(17) 0.0558567 -1.3474162 0.0558567I(18) 0.0527627 0.1229264 0.0527627I(19) 0.0499934 0.0464853 0.0499934I(20) 0.0475003 0.0476757 0.0475003当 a=15方案 I 结果 方案 II 结果 精确值I( 1) 0.0319222 0.0319222 0.0319222I( 2) 0.02
8、11673 0.0211673 0.0211673I( 3) 0.0158245 0.0158245 0.0158245I( 4) 0.0126326 0.0126326 0.0126326I( 5) 0.0105112 0.0105112 0.0105112I( 6) 0.0089993 0.0089993 0.0089993I( 7) 0.0078674 0.0078674 0.0078674I( 8) 0.0069883 0.0069883 0.0069883I( 9) 0.0062862 0.0062859 0.0062859I(10) 0.0057064 0.0057117 0.00
9、57117I(11) 0.0053136 0.0052336 0.0052337I(12) 0.0036289 0.0048293 0.0048296I(13) 0.0224896 0.0044830 0.0044838I(14) -0.2659159 0.0041831 0.0041831I(15) 4.0554050 0.0039207 0.0039207I(16) -60.7685756 0.0036893 0.0036893I(17) 911.5874579 0.0034837 0.0034837I(18) -13673.7563129 0.0033002 0.0032998I(19)
10、 205106.3973248 0.0031283 0.0031344I(20) -3076595.9098724 0.0030754 0.0029847 【结果分析】：1、综上所述，当 a=0.05 的时候，方案二算法的结果从 I(20)开始计算，刚开始的时候与精确解相差不大，但是随着计算的进行，误差变得越来越大，最终与原来的精确解相差十分巨大，而方案一算法的数值结果始终与精确解相差不大，是稳定的算法。2、当 a=15 的时候，反而是方案二的算法的数值结果与精确解更为接近，方案一算法的结果随着算法运算的进行，与精确解相差变大了。3、以上的实验说明了我们在进行数值分析时一定要选择合适的算法，而不能盲目的选择单一的算法，因为随着数值的变化，可能稳定的算法也会出现很大的误差。所以我们要根据实际问题来确定合适的算法，才能尽可能的减小误差。
