数值计算方法实践作业.docx

用Gauss-Seidel迭代法求解方程组。解:function [x,n]=gauseidel(A,b,x0,eps,M)if nargin==3eps= 1.0e-6;M = 200;elseif nargin == 4M = 200;elseif nargin<3errorreturn;endD=diag(diag(A));L=-tril(A,-1);U=-triu(A,1);G=(D-L)\U;f=(D-L)\b;x=G*x0+f;n=1;while norm(x-x0)>=epsx0=x;x=G*x0+f;n=n+1;if(n>=M)disp('Warning: μü′ú′?êyì??à￡??é?ü2?ê?á2￡?');return;endend>> A=[1.446 0.7948 0.8801;0.6946 1.9568 0.1730;0.6213 0.5526 1.9797];>> b=[1 0 1]';>> x0=zeros(3,1);>> [x,n]=gauseidel(A,b,x0)x = 0.5891-0.24340.3882n =11用Jacobi迭代法求解方程组。解:function [x,n]=jacobi(A,b,x0,eps,varargin)if nargin==3 eps= 1.0e-6; M = 200;elseif nargin<3 errorreturnelseif nargin ==5 M = varargin{1};endD=diag(diag(A)); L=-tril(A,-1); U=-triu(A,1); B=D\(L+U);f=D\b;x=B*x0+f;n=1; while norm(x-x0)>=eps x0=x; x=B*x0+f; n=n+1;if(n>=M) disp('Warning: μü′ú′?êyì??à￡??é?ü2?ê?á2￡?');return;endend>> A=[0.9989-0.0005 -0.0002;-0.0046 0.9946 0.0077;-0.0002 0.0092 0.9941];>> b=[1 0 1]';>> x0=zeros(3,1);>> [x,n]=jacobil(A,b,x0)x = 1.0114-0.00311.0062n =4求出下表所列各数据的拉格朗日插值多项式，并计算x=1.6时y的值。X11.21.82.54y0.84150.93200.97380.5985-0.7568解：function f = Language(x,y,x0)syms t;if(length(x) == length(y)) n = length(x); else disp('xoíyμ???êy2??àμè￡?');return;end%?ì′íf = 0.0;for(i = 1:n) l = y(i); for(j = 1:i-1) l = l*(t-x(j))/(x(i)-x(j)); end;for(j = i+1:n) l = l*(t-x(j))/(x(i)-x(j)); %????à-??àêè??ùoˉêyend; f = f + l; %????à-??àêè?2??μoˉêy simplify(f); %?ˉ?òif(i==n)if(nargin == 3) f = subs(f,'t',x0); %????2??μμ?μ?oˉêy?μelse f = collect(f); %??2??μ?à??ê??1?a f = vpa(f,6); %??2??μ?à??ê?μ??μêy?ˉ3é6?????èμ?D?êyendendend>> x=[1 1.2 1.8 2.5 4];>> y=[0.8415 0.9320 0.9738 0.5985 -0.7568];>> f=Language(x,y)f =0.0328112*t^4 - 0.204917*t^3 - 0.0145485*t^2 + 1.05427*t - 0.0261189>> f=Language(x,y,1.6)f = 0.9992求出下表所列各数据的均差形式的牛顿插值多项式，并计算x=1.