矩阵与数值分析上机实习题汇总.doc

矩阵与数值分析上机实习 设, 其精确值为. 编制按从大到小的顺序, 计算的通用程序 编制按从小到大的顺序, 计算的通用程序 按两种顺序分别计算并指出有效位数（编制程序时用单精度） 通过本上机题，你明白了什么 从小到大, 代码: %1---SN = % N = input(please input a number(N=2)) if(N 2) disp(wrong number) else S = 0; for j = 2:1:N S = S + 1/(j^2 -1); end disp(S:) disp(S) end 结果 please input a number(N=2)10^2 N = 100 S: 7.4005e-001 clear please input a number(N=2)10^4 N = 10000 S: 7.4990e-001 clear please input a number(N=2)10^6 N = 1000000 S: 7.5000e-001 从大到小 代码: %1---SN = % eps(single) N = input(please input a number(N=2)) if(N 2) disp(wrong number) else S = 0; for j = N:-1:2 S = S + 1/(j^2 -1); end disp(S:) disp(S) end 结果 please input a number(N=2)10^2 N = 100 S: 7.4005e-001 clear ans = 1.1921e-007 please input a number(N=2)10^4 N = 10000 S: 7.4990e-001 clear ans = 1.1921e-007 please input a number(N=2)10^6 N = 1000000 S: 7.5000e-001 (4)计算的顺序影响结果。 秦九韶算法。已知n次多项式，用秦九韶算法编写通用的程序计算函数在x0点的值，并计算在23点的值. (提示：编写程序时，输入系数向量和点x0，输出结果，多项式的次数可以通过向量的长度来判断 代码: A = input(请输入系数,由高次幂开始); n = input(请输入计算变量的值); len = length(A); val = zeros(len); val(1) = A(1); %%printf(len = %c, len) for i = 2: 1: len %disp(val(i-1)) %disp(n) val(i) = val(i-1)* n + val(i); end %printf(?????á1?￡o%f, val(len)) val(len) 结果: 请输入系数，由高次幂开始[7 3 -5 11] 请输入计算变量的值23 ans = 85169 分别用Gauss消元法和列主元消去法编程求解方程组Ax=b，其中 ， . 高斯消去法代码: A= [31 -13 0 0 0 -10 0 0 0 -13 35 -9 0 -11 0 0 0 0 0 -9 31 -10 0 0 0 0 0 0 0 -10 79 -30 0 0 0 -9 0 0 0 -30 57 -7 0 -5 0 0 0 0 0 -7 47 -30 0 0 0 0 0 0 0 -30 41 0 0 0 0 0 0 -5 0 0 27 -2 0 0 0 -9 0 0 0 -2 29 ]; B =[-15, 27, 23, 0, -20, 12, -7, 7, 10]; B = B.; %disp(A) %disp(B) C= [A,B]; n = size(A,1); ra = rank(A); rc= rank(C); x = zeros(1,n); if ra ~= rc disp(no solution) else for i = 1: 1: (n-1) for j = i + 1: 1: n temp = (-C(j,i))/ C(i,i); for k = i:1:(n+1)