太原理工大学计算方法上机题.docVIP

上机练习题 求非线性方程的根。 求方程在附近的是根，要求精度满足.(牛顿切线法) NewtonIteration x0=1.5 del=0.001 N=20 k x(k) 0 1.500000 1 0.784472 结果：0.739519 求方程在附近的是根，求出具有四位有效数字的根近似值..(简单迭代法) 程序 clear clc phi=inline((x^2+0.8)^(1/3)); %迭代函数 x0=input(x0=); del=input(del=); N=input(N=); n=1; fprintf(\n %2d %f ,0,x0); while nN x=phi(x0); if abs(x-x0)del fprintf(\n \n 近似解=%f \n,x); return end fprintf(\n %2d %f ,n,x); n=n+1; x0=x; end fprintf(\n \n %f d次迭代后未达到精度要求. \n,N); 运行结果 x0=1 del=0.00001 N=20 0 1.000000 1 1.216440 2 1.316116 3 1.363004 4 1.385180 5 1.395688 6 1.400671 7 1.403034 8 1.404155 9 1.404687 10 1.404939 11 1.405059 12 1.405116 13 1.405143 14 1.405155 近似解=1.405162 求解线性方程组（直接法或迭代法） 使用高斯-赛德尔迭代法求解 代码 clear clc n=input(n=);%矩阵的阶数 A=input(A=);%系数矩阵 b=input(b=); x=input(x=);%自变量 epsilon=input(\n 精度=); N=input(\n 最大迭代次数N=); fprintf(\n %d:,0); for i=1:n fprintf(%f,x(i)); end %以下是迭代过程 for k=1:N %这是第k步迭代，迭代前的向量在x0[]中，迭代后的向量在x[]中； normal=0; for i=1:n t=x(i); x(i)=b(i); for j=1:n if j~=i x(i)=x(i)-A(i,j)*x(j); end end x(i)=x(i)/A(i,i); temp=abs(x(i)-t);% 求范数于迭代在同一个循环中； if tempnormal normal=temp; %这里用的是无穷范数 end end %第i不迭代结束； fprintf(\n %d: ,k); for i=1:n fprintf(%f,x(i));%输出迭代过程 end if normalepsilon return end end fprintf(\n \n 迭代% d 次后仍未求得满足精度的解\n,N); 结果 n=4 A=[2,2,1,-3;-2,1,-1,-3;8,-1,3,8;10,4,4,3] b=[8,1,-1,8] x=[1,-1,2,-2] 精度=0.001 最大迭代次数N=10 0:1.000000-1.0000002.000000-2.000000 1: 1.000000-1.0000002.000000-2.000000 故原方程的解为 高斯-赛德尔 结果 n=6 A=[17.031,-0.615,-2.991,1.007,-1.006,0;-1,34.211,-1,-2.1,6.3,-1.7;0,0.5,13,-0.5,1,-1.5;4.501,3.11,-3.907,-61.705,12.17,8.99;0.101,-0.812,-0.017,-0.91,4.918,0.1;1,2,3,4.5,5,21.803] b=[0.23,-22.322,54,240.236,29.304,-117.818] x=