2014计算方法参考答案及评分标准.docVIP

2013-2014学年第二学期计算方法（A卷）参考答案： 一、 1．，或 （4分）； （4分）。 2． ，得8分； ，得6分；其他有部分正确者得4分。 3． （4分）；收敛（2分）；，或，或原系数矩阵按行严格对角占优（2分）。 4． 优点：（1）收敛快/平方收敛（4分）；（2）适用面广（1分）；（3）可求重根和复根（1分）；缺点：（4）对初值要求高（1分）；（5）需要计算导函数的值（1分）。 5. （4分） ；（4分） 若计算错误，但能给出公式： 可得2分。 二、 6．命令或.m文件，例如： %定义函数 vpGauss.m如下： function b vpGauss A,b fprintf 列主元消元法\n ; [n,n] size A ; for k 1:n-1 [maxv,r] max abs A k:n,k ; if maxv 1e-12 fprintf 主元过小，停止计算\n ; break; return end q r+k-1; A [k q],: A [q k],: ; b [k q] b [q k] ; A k+1:n,k A k+1:n,k /A k,k ; A k+1:n,k+1:n A k+1:n,k+1:n - A k+1:n,k *A k,k+1:n ; b k+1:n b k+1:n -A k+1:n,k *b k ; end if abs A n,n 1e-12 %% 缺该判断，扣1分。 fprintf 主元过小，停止计算\n ; break; return end b n b n /A n,n ; for k n-1:-1:1 b k b k -A k,k+1:n *b k+1:n /A k,k ; end %%程序结束 %命令窗口运行： A [1，2，3; 4，5，6; 7，8，1]; b [0；2；-2]; b vpGauss A,b %或 vpGauss A,b （算法正确12分；执行正确3分。编程命令有误每处扣1分，有严重语法错误扣5分）。 7．命令或.m文件，例如： 1 程序一 %%计算第二类边界条件的三次c样条的结点上的二阶导数值M function M Mspline2_1 X,Y,y2_0,y2_1 %%根据X、Y计算n阶三对角系数矩阵的对角系数a、b、c和右端项f; n length X ; h 2:n X 2:n -X 1:n-1 ; n1 n-2; a zeros 1,n1 ;b a;c a;f a; b b+2; a 1:n1 h 2:n-1 ./ h 2:n-1 +h 3:n ; c 1:n1 1-a 1:n1 ; Yf 2:n Y 2:n -Y 1:n-1 ; f 1:n1 6* Yf 3:n ./h 3:n -Yf 2:n-1 ./h 2:n-1 ./ h 2:n-1 +h 3:n ; f 1 f 1 -a 1 *y2_0; f n1 f n1 -c n1 *y2_1; %% A zeros n1,n1 ; % 该段程序用三对角追赶法代替也可以： A 1,1 b 1 ; % M SDJ1 n1,a,b,c,f ; for k 2:n1 % 并给出追赶法的程序。 A k,k b k ; A k,k-1 a k ; A k-1,k c k-1 ; end %%然后调用上题的列主元法解线性方程组 M vpGauss A,f’ ; %% M [y2_0; M; y2_1]; %% %结束 2 程序二： %%计算三次c样条曲线在任意点x上的函数值y function y yMsx X,Y,M,x %%根据X、Y、M和x计算y; n length X ; h 2:n X 2:n -X 1:n-1 ; %%判断x落在哪个小区间[xj,xj+1];j n1 n1 1;n2 n; while n2~ n1+1 n5 fix n1+n2 /2 ; if x X n5 n1 n5; else n2 n5; end end %% %%计算y y M n1 * X n2 -x ^3/ 6*h n2 + M n2 * x-X n1 ^3/ 6*h n2 ; y y+ Y n1 -M n1 *h n2 *h n2 /6 * X n2 -x /h n2 ; y y+ Y n2 -M n2 *h n2 *h n2 /6 * x-X n1 /h n2 ; %% %结束 3 主执行程序： %命令窗口运行： 或写入一个函数进行执行 X [1 2 3 4 5 6 7 8 9 10]; Y [0.0000 0.6931 1.0986 1.3863 1.60