数值分析-2120111..docVIP

数值分析实验 李红芳 2120111352 数值计算方法实验题 学 院： 机电 学院 班 级： 化学工程与技术11S3班 姓 名： 学 号： 第一章 第4题 a=-10;b=10;c=100;d=200;n=100;eps1=0.01; x=linspace(a,b,n);y=linspace(a,b,n); [X,Y]=meshgrid(x,y); for i=1:n %计算函数值z ，并作定义域裁剪 for j=1:n if (1-X(i,j))eps1|X(i,j)-Y(i,j)eps1 %if语句这样用 z(i,j)=NaN; %作定义域裁剪，定义域以外的函数值为NaN else z(i,j)=exp(-abs(X(i,j)))+cos(X(i,j)+Y(i,j))+1./(X(i,j)^2+Y(i,j)^2+1); end end end zz=-20*ones(1,n);plot3(x,x,zz),grid off,hold on %画定义域的边界线 mesh(X,Y,z) %绘图，读者可用meshz, surf, meshc在此替换之 xlabel(x),ylabel(y),zlabel(z), box on %把三维图形封闭在箱体里 第二章 第2题 function x=zhuigan(a,b,c,d,n) % % a为块矩阵A,b为块矩阵B,c为块矩阵C,d为等式右侧数据,n是电流数,x为所求个电路的电流量； format long; for i=2:n a(i)=a(i)./b(i-1); b(i)=b(i)-c(i-1).*a(i); d(i)=d(i)-a(i).*d(i-1); end d(n)=d(n)./b(n); n=n-1; for i=n:-1:1 d(i)=(d(i)-c(i).*d(i+1))./b(i); end x=d; 输入题中所给信息，调用函数 a=[0;-2;-2;-2;-2;-2;-2;-2]; b=[2;5;5;5;5;5;5;5]; c=[-2;-2;-2;-2;-2;-2;-2;0]; d=[220./27;0;0;0;0;0;0;0]; n=8; x=zhuigan(a,b,c,d,n) 可得 x = 8.147775166909105 4.073701092835030 2.036477565178471 1.017492820111148 0.507254485099400 0.250643392637350 0.119353996493976 0.047741598597591 即各段电路电流量：i1=8.148;i2=4.074;i3=2.036;i4=1.017;i5=0.507;i6=0.251;i7=0.119;i8=0.048; 第三章 第1题 （1）Jacobi迭代法 function jacob=jacob(a,b,n,x,e,N) %a矩阵A,b为矩阵B,x为初值,n为项数,e为误差,N为迭代次数,jacob为方程组解 k=1; x1=[0;0;0]; while kN k=k+1; %k为迭代次数 for i=1:n x1(i)=(b(i)-a(i,:)*x)./a(i,i)+x(i); end if norm(x1-x)e break; end x=x1; end if k=N fprintf(error!!) else jacob=x1; end k a=[10,1,2,3,4;1,9,-1,2,-3;2,-1,7,3,-5;3,2,3,12,-1;4,-3,-5,-1,15]; b=[12;-27;14;-17;12]; x=[0;0;0;0;0]; e=0.001; n=5; N=1