《北理工数值分析大作业.docx》由会员分享，可在线阅读，更多相关《北理工数值分析大作业.docx（19页珍藏版）》请在金锄头文库上搜索。
1、北理工数值分析大作业数值剖析上机做业第 1 章1.1盘算积分，n=9。（请求盘算了局具备6位无效数字）步伐：n=1:19;I=zeros(1,19);I(19)=1/2*(exp(-1)/20)+(1/20);I(18)=1/2*(exp(-1)/19)+(1/19);for i=2:10I(19-i)=1/(20-i)*(1-I(20-i);endformat longdisp(I(1:19)了局截图及剖析：正在MATLAB中运转以上代码，患上到了局以下图所示：当盘算到数列的第10项时，所患上的了局即为n=9时的正确积分值。与6位无效数字可患上.1.2分手将区间-10.10分为100,200
2、,400等份，使用mesh或者surf下令绘出2元函数z=的3维图形。步伐： x = -10:0.1:10;y = -10:0.1:10;X,Y = meshgrid(x,y);Z = exp(-abs(X)+cos(X+Y)+1./(X.2+Y.2+1);subplot(2,2,1);mesh(X,Y,Z);title(步少0.1) x = -10:0.2:10;y = -10:0.2:10;X,Y = meshgrid(x,y);Z = exp(-abs(X)+cos(X+Y)+1./(X.2+Y.2+1);subplot(2,2,1);mesh(X,Y,Z);title(步少0.2)x
3、= -10:0.05:10;y = -10:0.05:10;X,Y = meshgrid(x,y);Z = exp(-abs(X)+cos(X+Y)+1./(X.2+Y.2+1);subplot(2,2,1);mesh(X,Y,Z);title(步少0.05)了局截图及剖析：由图可知，步少越小时，画患上的图形越粗确。 第 2 章 试用MATLAB 编程真现逃赶法供3对于角圆程组的算法，并思索梯形电路电阻成绩：电路中的电流128,i i i 谦足以下线性圆程组：121232343454565676787822/252025202520252025202520250i i V R i i i i
4、i i i i i i i i i i i i i i i i -=-+-=-+-=-+-=-+-=-+-=-+-=-+=设220,27V V R =，供各段电路的电流量。处置思绪：不雅察该圆程的系数矩阵可知，它是一个3对于角矩阵，故可使用逃赶法对于其举行供解。 步伐:for i=1:8a(i)=-2;b(i)=5;c(i)=-2;d(i)=0; enda(1)=0;b(1)=2;c(8)=0;d(1)=220/27; for i=2:8a(i)=a(i)/b(i-1); b(i)=b(i)-c(i-1)*a(i); d(i)=d(i)-a(i)*d(i-1); endd(8)=d(8)/b(
5、8); for i=7:-1:1d(i)=( d(i)-c(i)*d(i+1) )/b(i); end for i=1:8 x(i)=d(i); endx了局截图及剖析：正在MATLAB 中运转以上代码，患上到了局以下图所示：图中8个值挨次为128,i i i 的数值。 第 3 章 试分手用（1）Jacobi 迭代法；（2）Gauss-Seidel 解线性圆程组 1234510123412191232721735143231211743511512x x x x x ?-?-=?-?-?迭代初初背量与(0)(0,0,0,0,0)T x =. 3.1 Jacobi 迭代法 步伐： A=10 1
6、2 3 4;1 9 -12 -3; 2 -1 73 -5; 3 2 3 12 -1;4 -3 -5 -1 15; b=12;-27;14;-17;12; x0=0;0;0;0;0;D=diag(diag(A); I=eye(5); L=-tril(A,-1); B=I-DA; g=Db; y=B*x0+g; n=1;while norm(y-x0)=1.0e-6 x0=y;y=B*x0+g; n=n+1; endfprintf(%8.6fn,y); n患上到此了局时迭代次数为67次，到达粗度请求。3.2 Gauss-Seidel迭代法：步伐： A=10 1 2 3 4;1 9 -12 -3;2
7、 -1 73 -5;3 2 3 12 -1;4 -3 -5 -1 15;b=12;-27;14;-17;12;x0=0;0;0;0;0;D=diag(diag(A);U=-triu(A,1);L=-tril(A,-1);M=(D-L)U;g=(D-L)b;y=M*x0+g;n=1;while norm(y-x0)=1.0e-6x0=y;y=M*x0+g;n=n+1;endfprintf(%8.6fn,y);Gauss-Seidel迭代法只要要迭代38次便可谦足粗度请求。第 4 章 设A=?-162621666612，与先用幂法迭代3次，患上到A 的按模最年夜特性值的远似值，与为其整数全体，再用
8、反幂法盘算A的按模最年夜特性值的更粗确的远似值，请求偏差小于.步伐：A=12 6 -6;6 16 2; -6 2 16; x0=1;1;1;y=x0;b=max(abs(x0);k=1; while ( kx=A*y;b=max(abs(x);y=x./b; k=k+1;fprintf(eig1 equals %6.4fn,b); end bb0=fix(b); I=eye(3,3); x0=1;1;1;y=x0;l=0;bb=max(abs(x0);k=1; while ( abs(bb-l)=1.0e-10 ) l=bb;x=(A-bb0*I)y;bb=max(abs(x);y=x./bb
9、; eig=l+b; fprintf(eig2(%d) equals %12.10fn,k, eig); k=k+1; end真验截图及剖析：由图可知，由幂法3次迭代后患上到的特性值为19.4，而由反幂法患上到的特性值为20.3999999999.偏差小于第 5 章试编写MATLAB函数真现Newton插值，请求能输入插值多项式。对于函数f(x)=正在区间-5,5上真现10次多项式插值。请求：（1）输入插值多项式。（2）正在区间-5,5内匀称拔出99个节面，盘算那些节面上函数f（x）的远似值，并正在统一张图上绘出本函数以及插值多项式的图形。（3）不雅察龙格征象，盘算插值函数正在各节面处的偏差，
10、并绘出偏差图。5.1输入插值多项式步伐：x=-5:1:5;y=1./(1+4*(x.2);newpoly(x,y)function c,d=newpoly(x,y)n=length(x);d=zeros(n,n);d(:,1)=y;for j=2:nfor k=j:nd(k,j)=(d(k,j-1)-d(k-1,j-1)/(x(k)-x(k-j+1);endendc=d(n,n);for k=(n-1):-1:1c=conv(c,poly(x(k);m=length(c);c(m)=c(m)+d(k,k);endend了局及剖析：ans =Columns 1 through 2-0.00004
11、9595763049Columns 3 through 40.002740165908483 0.000000000000000Columns 5 through 6-0.051421507076720 0.000000000000000Columns 7 through 80.392014985282312 0.000000000000000Columns 9 through 10-1.143284048351025 0.000000000000001Column 111.00000000000000010次插值多项式由下到低系数为Columns 1至Column 11 5.2本函数取插值多
12、项式的图形步伐：x=-5:1:5;y=1./(1+4*(x.2);n=newpoly(x,y);x0=-5:0.1:5;y0=1./(1+4*(x0.2);vn=polyval(n,x0);plot(x0,vn,-r,x0,y0,-b);xlabel(x);ylabel(y);真验了局截图：yx本函数取插值多项式的图形如上图所示，蓝色为本函数的图形，白色为插值多项式的图形。5.3各节面的偏差及偏差图步伐：format long;x=-5:1:5;y=1./(1+4*(x.2);n=newpoly(x,y);x0=-5:0.1:5;y0=1./(1+4*(x0.2);vn=polyval(n,x
13、0);plot(x0,y0-vn,-r);xlabel(x);ylabel(y);真验了局截图：yx 偏差图如上图所示。第 6 章 炼钢厂出钢时所用的衰钢火的钢包，正在利用历程中因为钢液及炉渣对于包衬耐水质料的侵蚀，使其容积没有断减年夜。经实验，钢包的容积取响应的利用次数的数据列表以下：选用单直线xb a y11+=对于数据举行拟开，利用最小2乘法供出拟开函数，做出拟开直线图。 处置思绪：用Y 代替1/y ，用X 代替1/x ，本直线化为Y=a+bx ，单直线转化为一次线性圆程，利用最小2乘法供出该一次圆程的系数。步伐：x=2 3 5 6 7 9 10 11 12 14 16 17 19 20
14、;y=106.42 108.26 109.58 109.5 109.86 110 109.93 110.59 110.60 110.72 110.9 110.76 111.1 111.3;k1=0; k2=0; k3=0; k4=0;for i=1:14k1=k1+1/x(i);endfor i=1:14k2=k2+1/y(i);endfor i=1:14k3=k3+1/(x(i)2;endfor i=1:14k4=k4+1/(x(i)*y(i);endb=(k1*k2-14*k4)/(k12-14*k3)a=k2/14-k1*b/14plot(x,y,r*)hold onx=2:0.01:20;y=1./(a+b./x);plot(x,y)xlabel(x)ylabel(y)grid on真验了局截图取剖析：即最小2乘法供出拟开函数为：=0.008973+
