計算方法上机题.docxVIP

计算,要求误差小于,给出实现算法。算法组织：利用绝对值的递增相加可以减小误差具体算法：%% 公式求和 %%S = 0;k = 100000;%%利用for循环对公式进行求和for n = 1: 100000S = S+1/(k*k);k = k-1;end%%精确到小数点后6位,以保证误差范围vpa(S,7)运行程序后，得到的结果是1.644924，精确到了小数点后6位。编写实现对N阶非奇矩阵A进行LU分解的程序。算法组织：根据定理可知若N阶矩阵的各阶顺序主子式行列式不为零，则存在唯一的单位下三角矩阵L和上三角矩阵U，满足A=LU，因此首先需要判断A中的元素akk（k=0,1,2…，n）不为零。根据如下的式子计算L和U中的元素当i=1时，得到u1j=a1j(j=1,2,3,…,n),而当j=1时，得li1= ai1/u11，由此可以计算出L的第一列元素和U的第一行元素。然后可以依据下式计算L和U中的其他元素。具体算法：[n,n]=size(A); %检测A的阶数L=eye(n); %初始化单位下三角阵LU=zeros(n,n);%初始化上三角阵Ufor i=1:n U(1,i)=A(1,i);%给上三角阵U的第一行赋值 L(i,1)=A(i,1)/U(1,1);%给单位下三角阵L的第一列赋值endfor i=2:n for j=i:n for k=1:i-1 M(k)=L(i,k