目录
第一章 2
第1题 2
第3题 2
第4题: 3
第二章 6
第1题 6
第2题 7
第4题 8
第三章 11
第1题 11
第五章 13
第1题 13
第七章 16
第1题: 16
第八章 18
第1题 18
第一章
第1题
程序:
format long
i=1; y1=sqrt(5);y2=sqrt(y1);
while y1-y2>10^-8
y1=5.^(1./(2.^i));y2=5.^(1./(2.^(i+1)));
i=i+1;
end
y1,y2,i
运行结果:
y1 =
y2 =
i =28
第3题
程序:
x=-pi::pi;
y=sin(x); y2=x;y5=x;y10=x;
for i=1:2
y2=(-1)^i.*(x.^(2*i+1))./(factorial(2*i+1))+y2;
end
for i=1:5
y5=(-1)^i.*(x.^(2*i+1))./(factorial(2*i+1))+y5;
end
for i=1:10
y10=(-1)^i.*(x.^(2*i+1))./(factorial(2*i+1))+y10;
end
plot(x,y,'-b',x,y2,'-g',x,y5,':r',x,y10,':k');
legend('y','y2','y5','y10')
运行结果:
第4题:
程序:
①等分100份
x=-10::10;y=-10::10;
[x,y]=meshgrid(x,y);
z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);
z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);
mesh(x,y,z)
运行结果:
②等分200份
x=-10::10;y=-10::10;
[x,y]=meshgrid(x,y);
z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);
z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);
mesh(x,y,z)
运行结果:
③等分300份
x=-10::10;y=-10::10;
[x,y]=meshgrid(x,y);
z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);
z=exp(-abs(x))+cos(x+y)+1./(x.^2+y.^2+1);
mesh(x,y,z)
运行结果:
第二章
第1题
程序:
列主元素求解
A=[1 1 1 1 1;1 2 3 4 5;1 3 6 10 15;1 4 10 20 35;1 5 15 35 70];
[L,U,P]=lu(A)
L =
0 0 0 0
0 0 0
0 0
0
-
U =
0
0 0 - - -
0 0 0 - -
0 0 0 0 -
P =
1 0 0 0 0
0 0 0 0 1
0 0 1 0 0
0 0 0 1 0
0 1 0 0 0
Inv(A)求逆:
inv(A)
ans =
- -
- - -
- -
- - -
- -
左除命令A\E
E=[1 0 0 0 0;0 1 0 0 0;0 0 1 0 0;0 0 0 1 0;0 0 0 0 1];
A\E
ans =
5 -10 10 -5 1
-10 30 -35 19 -4
10 -35 46 -
数值分析matlab作业 来自淘豆网www.taodocs.com转载请标明出处.