Matlab程序设计教程(第二版)课后参考
答案
八年级地理上册填图题岩土工程勘察试题省略号的作用及举例应急救援安全知识车间5s试题及答案
第二章
实验
1.(1) w=sqrt(2)*(1+0.34245*10^(-6))
w = 1.4142
(2) a=3.5;
b=5;
c=-9.8;
x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a)
x =0.9829
(3) a=3.32;
b=-7.9;
y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a)
y = -128.4271
(4) t=[2,1-3i;5,-0.65];
z=0.5*exp(2*t)*log(t+sqrt(1+t*t))
z =
1.0e+004 *
0.0048 + 0.0002i 0.0048 - 0.0034i
1.5899 2.0090 - 1.3580i
2.实验二 第1
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
3.实验二 第2题
4. H=hilb(5);
P=pascal(5);
Hh=det(H)
Hh = 3.7493e-012
Hp=det(P)
Hp = 1
Th=cond(H)
Th = 4.7661e+005
Tp=cond(P)
Tp = 8.5175e+003
条件数越趋近于1,矩阵的性能越好,所以帕斯卡矩阵性能更好。
5. A=[-29,6,18;20,5,12;-8,8,5]
A =
-29 6 18
20 5 12
-8 8 5
[V,D]=eig(A)
V =
0.7130 0.2803 0.2733
-0.6084 -0.7867 0.8725
0.3487 0.5501 0.4050
D =
-25.3169 0 0
0 -10.5182 0
0 0 16.8351
V为A的特征向量,D为A的特征值。它们之间满足A*V=V*D
思考练习
1. 在MATLAB中6+7i是一个复数常量,6+7*i则是一个表达式。i是虚数单位,而I是单位向量。
2. A.*B表示A和B单个元素之间对应相乘,A*B是按矩阵乘法的规则乘。
A./B表示A矩阵除以B矩阵的对应元素,B.\A等价于A./B。
A/B表示A *inv(B)而B\A表示inv(B)* A,对于矩阵运算一般两式不等。
3.(1) A(7)=[]
(2) abs(x)
(3) reshape(x,3,4)
(4) k=find(t==0); t(k)=0
(5) ones(size(A))
(6) diag(A)
B= diag(diag(A))
4. 3+randn(25,20)
5.实验二 第3题
6. A=[34,NaN,Inf,-Inf,-pi,eps,0];
all(A)
ans = 0
any(A)
ans = 1
isnan(A)
ans = 0 1 0 0 0 0 0
isinf(A)
ans = 0 0 1 1 0 0 0
isfinite(A)
ans = 1 0 0 0 1 1 1
7.结构矩阵建立P46
结构矩阵名.成员名=表达式
单元矩阵建立 P47
矩阵元素用大括号括起来。
8.用稀疏存储方式可以可以节省内存空间; 运算规则上,稀疏矩阵与普通矩阵一样,只是矩阵的存储方式不同。在运算过程中,稀疏存储矩阵可以直接参与运算。当参与运算的对象不全是稀疏存储矩阵时,所得结果一般是完全存储形式。
第三章
实验
1.提示: 先根据输入的3位整数分别求个、十、百位上的数字a、b、c,反向输出整数:100*a+10*b+c
2(1)
m=input('请输入一个百分制的成绩: ');
if m>=0&m<=100
switch fix(m/10)
case {9,10}
disp('这个成绩所属的成绩等级:A');
case {8}
disp('这个成绩所属的成绩等级:B');
case {7}
disp('这个成绩所属的成绩等级:C');
case {6}
disp('这个成绩所属的成绩等级:D');
otherwise
disp('这个成绩所属的成绩等级:E');
end
else
disp('错误:输入的成绩不合理');
end
3
a=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]
b=max(a)
c=min(a)
b=0
c=0
for i=1:20
if a(i)>b
b=a(i);
end
if a(i)
5
n=5
A(n,n)
else
A(n)
end
7
f1.m代码如下:
function m=f1(n)
m=n+10*log(n^2+5);
f2.m 代码如下:
function m=f2(n)
m=0;
for i=1:n
m=m+i*(i+1);
end
然后在命令文件p3_11.m中调用f1.m和f2.m
p3_11.m代码如下:
y1=f1(40)/(f1(30)+f1(20));
y2=f2(40)/(f2(30)+f2(20));
y1
y2
8
函数文件factor1.m代码如下:
function f=factor1(k)
if k==1
f=1;
else
f=factor1(k-1)+k;
end
函数文件factor2.m代码如下:
function f=factor2(k)
if k==1
f=1;
else
f=factor2(k-1)+^2;
end
函数文件factor3.m代码如下:
function f=factor3(k)
if k==1
f=1;
else
f=factor3(k-1)+1/k;
end
思考练习
1.见课本3.1
2
N=[1,2,3,4,5]
A=N.*2
B=N./4
C=1./N
D=1./(N./2)./(N./2)/4
3
B=0
for i=1:20
A(i)=fix((99-10+1)*rand(1)+10);
B=A(i)+B;
end
B=B/20
for i=1:20
if A(i)0.5
c=c+1;
end
end
b=mean(a)
b=std(a,0,2)
b=max(a)
b=min(a)
c/30000
5
a=fix((100-0+1)*rand(5,100));
[b,c]=max(a,[],2)
[b,c]=min(a,[],2)
mean(a,2)
std(a,0,2)
[b,c]=sort(a,2,'descend')
zcj=sum(a,1)
xsxh=1:100
PAGE
9