matlab课后习题答案Word下载.docx
- 文档编号:21117550
- 上传时间:2023-01-27
- 格式:DOCX
- 页数:13
- 大小:63.80KB
matlab课后习题答案Word下载.docx
《matlab课后习题答案Word下载.docx》由会员分享,可在线阅读,更多相关《matlab课后习题答案Word下载.docx(13页珍藏版)》请在冰豆网上搜索。
4.下列命令执行后,L1、L2、L3、L4的值分别是多少?
A=1:
9;
B=10-A;
...
L仁A==B;
L2=A<
=5;
L3=A>
3&
A<
7;
L4=find(A>
7);
L1的值为[0,0,0,0,1,0,0,0,0]
L2的值为[1,1,1,1,1,0,0,0,0]
L3的值为[0,0,0,1,1,1,0,0,0]
L4的值为[4,5,6]
5.已知
-
23
41
10
-45
-0.778
65
01
5
A=
32
6
-9.54
54
3.14一
完成下列操作:
(1)取出A的前3行构成矩阵B,前两列构成矩阵C,右下角32子矩阵构成矩阵D,B与C的乘积构成矩阵E。
B=A(1:
3,:
C=A(:
1:
2);
D=A(2:
4,3:
4);
6.当A=[34,NaN,Inf,-Inf,-pi,eps,0]时,分析下列函数的执行结果:
all(A),any(A),isnan(A),isinf(A),isfinite(A)。
all(A)的值为0
any(A)的值为1
isnan(A)的值为[0,1,0,0,0,0,0]isinf(A)的值为[0,0,1,1,0,0,0]isfinite(A)的值为[1,0,0,0,1,1,1]
7.用结构体矩阵来存储5名学生的基本情况数据,每名学生的数据包括学号、姓名、专业和6门课程的成绩。
student
(1).id='
0001'
;
student
(1).name='
Tom'
student
(1).major='
computer'
student
(1).grade=[89,78,67,90,86,85];
&
建立单元矩阵B并回答有关问题。
B{1,1}=1;
B{1,2}='
Brenden'
B{2,1}=reshape(1:
9,3,3);
B{2,2}={12,34,2;
54,21,3;
4,23,67};
(1)size(B)和ndims(B)的值分别是多少?
size(B)的值为2,2。
ndims(B)的值为2。
(2)B
(2)和B⑷的值分别是多少?
_14_[12][34][2]〕
B
(2)=258,B(4)=[54][21][3]
[369」J4][23][67]j
⑶B(3)=[]和B{3}=[]执行后,B的值分别是多少?
当执行B(3)=[]后,
B={1,[1,4,7;
2,5,8;
3,6,9],{12,34,2;
54,21,3;
4,23,67}}当执行B{3}=[]后,
B={1,[];
[1,4,7;
4,23,67}}
1.写出完成下列操作的命令。
(1)建立3阶单位矩阵Ao
A=eye(3);
⑵建立5X6随机矩阵A,其元素为[100,200]范围内的随机整数
round(100+(200-100)*rand(5,6));
⑶产生均值为1,方差为0.2的500个正态分布的随机数。
1+sqrt(0.2)*randn(5,100);
(4)产生和A同样大小的幺矩阵。
ones(size(A));
⑸将矩阵A对角线的元素加30。
A+eye(size(A))*30;
(6)从矩阵A提取主对角线元素,并以这些元素构成对角阵B。
B=diag(diag(A));
90°
或右旋90°
的功能。
例如,原矩阵为A,A左旋
B=rot90(A);
C=rot90(A,-1);
3.建立一个方阵A,求A的逆矩阵和A的行列式的值,并验证A与A-1是互逆的答:
A=rand(3)*10;
B=inv(A);
C=det(A);
先计算B*A,再计算A*B,由计算可知B*A=A*B,即A-A-1=A-1•A是互
逆。
4.求下面线性方程组的解。
4x12x2-x3=2
<
3冶_x2+2x3=10
12x1*3x?
=8
A=[4,2,-1;
3,-1,2;
12,3,0];
b=[2;
10;
8];
x=inv(A)*b
方程组的解为x=
27.3333
5.求下列矩阵的主对角线元素、上三角阵、下三角阵、秩、范数、条件数和迹
-1
2
31
(1)A=
1
-4
0.43
432]
(2)B=|
3
J-8.9
421一
11
15
9
(1)取主对角线元素:
diag(A);
上三角阵:
triu(A);
下三角阵:
tril(A);
秩:
rank(A);
范数:
norm(A,1);
或norm(A);
或norm(A,inf);
条件数:
cond(A,1);
或cond(A,2);
或cond(A,inf)迹:
trace(A);
⑵【请参考⑴】
6.求矩阵A的特征值和相应的特征向量。
110.5
A=110.25
]0.50.252一
[V,D]=eig(A);
习题四
1.从键盘输入一个4位整数,按如下规则加密后输出。
加密规则:
每位数字都加上7,然后用和除以10的余数取代该数字;
再把第一位与第三位交换,第二位与第四位交换。
a=input('
请输入4位整数:
'
A=[a/1000,a/100,a/10,a];
A=fix(rem(A,10));
A=rem(A+7,10);
b=A(3)*1000+A(4)*100+A
(1)*10+A
(2);
disp(['
加密后的值为:
num2str(b)]);
2.分别用if语句和switch语句实现以下计算,其中a、b、c的值从键盘输入'
,广
ax+bx+c,0.5兰xc1.5y=«
asincb+x,1.5^xv3.5
Inb+£
3.5^xc5.5
Ix
(1)用if语句实现计算:
请输入a的值:
b=input('
请输入b的值:
’);
c=input('
请输入c的值:
x=input('
请输入x的值:
ifx>
=0.5&
x<
1.5
y=a*xA2+b*x+c;
end
=1.5&
3.5
y=a*((sin(b)Fc)+x;
=3.5&
5.5
y=log(abs(b+c/x));
y='
num2str(y)]);
(2)用switch语句实现计算:
请输入a的值:
c=input('
请输入c的值:
switchfix(x/0.5)
case{1,2}
casenum2cell(3:
6)
y=a*((sin(b))Ac)+x;
casenum2cell(7:
10)
y=log(abs(b+c/x));
num2str(y)]);
3.产生20个两位随机整数,输出其中小于平均值的偶数。
A=fix(10+89*rand(1,20));
sum=0;
fori=1:
20
sum=sum+A(i);
B=A(find(A<
(sum/20)));
C=B(find(rem(B,2)==0));
disp(C);
4.
MATLAB
输入20个数,求其中最大数和最小数。
要求分别用循环结构和调用的max函数、min函数来实现。
(1)用循环结构实现:
v_max=0;
v_min=0;
x=input(['
请输入第'
num2str(i),'
数:
’]);
v_max
v_max=x;
end;
ifx<
v_min
v_min=x;
最大数为:
num2str(v_max)]);
最小数为:
num2str(v_min)]);
(2)用max函数、min函数实现:
A(i)=input(['
请输入第;
num2str(i),'
num2str(max(A))]);
num2str(min(A))]);
5.已知:
2222^11263,分别用循环结构和调用MATLAB的sum函数求s的值。
(1)用循环结构实现:
s=0;
fori=0:
63
s=s+2Ai;
s
(2)调用sum函数实现:
s=0:
63;
s=2.As;
sum(s)
6.当n分别取100、1000、10000时,求下列各式的值。
(1)1丄1丄III(T)n11川(Rn2)
234n
1」丄1川
(二)
3574
11111
416644n(3)
要求分别用循环结构和向量运算(使用sum或prod函数)来实现。
sum=0;
fork=1:
100sum=sum+(-1F(k+1)/k;
sum
使用sum函数:
x=[];
10000
x=[x,(-1)A(k+1)/k];
sum(x)
(2)用循环结构实现:
100sum=sum+(-1)A(k+1)/(2*k-1);
x=[];
100
x=[x,(-1)A(k+1)/(2*k-1)];
(3)用循环结构实现:
sum=O;
sum=sum+1/(4Ak);
使用sum函数实现:
x=[x,1/(4八©
];
(4)用循环结构实现:
t=1;
t=t*(((2*k)*(2*k))/((2*k-1)*(2*k+1)));
t
使用prod函数实现:
x=[x,((2*k)*(2*k))/((2*k-1)*(2*k+1))];
prod(x)
7.编写一个函数文件,求小于任意自然数n的斐波那契(Fibnacci)数列各项
斐波那契数列定义如下:
f^1,n=1
«
f2=1,n=2
乙=fjfz,n>
functionx=fibnacci(n)
n
ifi<
=2
x(i)=1;
else
x(i)=x(i-1)+x(i-2);
8.编写一个函数文件,用于求两个矩阵的乘积和点乘,然后在命令文件中调用该函数。
函数文件myfnc.m:
function[x,y]=myfnc(A,B)
try
x=A*B;
catch
y=A.*B;
命令文件myexe.m:
A=input('
请输入矩阵A:
B=input('
请输入矩阵B:
[x,y]=myfnc(A,B);
iflength(x)==0
display('
两矩阵的维数不匹配,无法进行乘积运算!
disp('
矩阵A和矩阵B的乘积为:
x
disp(矩阵A和矩阵B的点乘为:
y
n
9.先用函数的递归调用定义一个函数文件求
im,然后调用该函数文件求
i=1
k1kdkTk
functionsum=myfnc(n,m)
ifn<
=1
sum=1;
sum=myfnc(n-1,m)+nAm;
1005010I
在命令窗口中调用myfnc.m文件,计算k、k^-:
k4kmk£
k
sum=myfnc(100,1)+myfnc(50,2)+myfnc(10,-1)
10.写出下列程序的输出结果。
1s=0;
a=[12,13,14;
15,16,17;
18,19,20;
21,22,23];
fork=a
forj=1:
4
ifrem(k(j),2)~=0
s=s+k(j);
执行结果为
s=108
2命令文件exe.m执行后的结果为:
x=
41220
y=
246
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 课后 习题 答案