数值计算方法实验.docVIP

P54页例4—3 Jacobi迭代法 n=4; A=[10,-1,2,0;-1,11,-1,3;2,-1,10,-1;0,3,-1,8];b=[6,25,-11,15]; x0=[0,0,0,0]; x0=zeros(n,1);x=x0; epsilon=input(\n 精度=); N=input(\n 最大迭代次数N=); fprintf(\n %d:,0); for i=1:n fprintf(%f,x0(i)); end %一下是迭代过程 for k=1:N %这是第K步迭代，迭代前的向量在X0【】中，迭代后的再X【】中 normal=0; for i=1:n x(i)=b(i); for j=1:n if j~=i x(i)=x(i)-A(i,j)*x0(j); end end x(i)=x(i)/A(i,i); temp=abs(x(i)-x0(i));%求范数与迭代在同一个循环中 if tempnormal normal=temp;%这里用的是无穷范数 end end %第i步迭代结束 fprintf(\n %d:,k); for i=1:n x0(i)=x(i);%为下一步迭代准备初值 fprintf(%f,x(i));%输出迭代过程 end if normalepsilon return; end end fprintf(\n\n迭代%d次后仍未求得满足精度的解\n,N); Gauss-Seidel迭代法 n=4; A=[10,-1,2,0;-1,11,-1,3;2,-1,10,-1;0,3,-1,8];b=[6,25,-11,15]; x0=[0,0,0,0]; x0=zeros(n,1);x=x0; 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); 实验5 p83 第一题 n=3;x=[1,15,22]; y=[17-6+(58-9)/60.0,17-6+(36-23)/60.0,17-6+(26-31)/60.0]; xx=8; yy=0; for i=1:3 temp=y(i); for j=1:3 if j~=i temp=temp*(xx-x(j))/(x(i)-x(j)); end end yy=yy+temp; end fprintf(\n插值多项式在x=%f点的值为%f，,xx,yy); fprintf(而f(%f)=%f\n,xx,[17-6+(58-9)/60.0,17-6+(36-23)/60.0,17-6+(26-31)/60.0]); 第2题 f=fopen(FittingData.txt,r); n=fscanf(f,%d,1); for i=1:n tx=fscanf(f,%f,1); tz=fscanf(f,%f,1); x(i)=1/t; z(i)=log(y); end fclose(f); sum_x=0;sum_x2=0;sum_z=0;sum_xz=0; for i=1:n sum_x=sum_x+x(i); sum_x2=sum_x2+x(i)*x(i);