数值分析.doc

数值计算方法实验题
学院: 机电学院
班级: 化学工程与技术11S3班
姓名:
学号:
第一章
第4题
>> a=-10;b=10;c=100;d=200;n=100;eps1=;
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 =








即各段电路电流量:i1=;i2=;i3=;i4=;i5=;i6=;i7=;i8=;
第三章
第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 k<N
k=k+1;