matlab课后习题答案

题二
1 理解矩阵MATLAB基数象？
答：量成仅行列矩阵单数（标量）成仅含元素矩阵量单数作矩阵特例处理
矩阵MATLAB基重数象

2 设AB两维矩阵问：
(1) A*BA*B值否相等？
答：相等
(2) ABB\A值否相等？
答：相等
(3) ABB\A值否相等？
答：相等
(4) ABB\A代表数学含义什？
答：AB等效B逆右A矩阵A*inv(B)B\A等效B矩阵逆左A矩阵inv(B)*A

3 写出完成列操作命令
(1) 矩阵A第2~5行中第1 3 5列元素赋矩阵B
答：BA(25125) BA(25[1 3 5])
(2) 删矩阵A第7号元素
答：A(7)[]
(3) 矩阵A元素值加30
答：AA+30
(4) 求矩阵A维数
答：size(A)
ndims(A)
(5) 量 t 0元素机器零代
答：t(find(t0))eps
(6) 含12元素量 x 转换成矩阵
答：reshape(x34)
(7) 求字符串ASCII码
答：abs(123’) double(123’)
(8) 求ASCII码应字符
答：char(49)

4 列命令执行L1L2L3L4值分少？
A19B10A
L1AB
L2A<5
L3A>3&A<7
L4find(A>3&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 已知

完成列操作：
(1) 取出A前3行构成矩阵B前两列构成矩阵C右角子矩阵构成矩阵DBC积构成矩阵E
答：BA(13)
CA(12)
DA(2434)
EB*C
(2) 分求E10&A<25)
答：E find(A>10&A<25)[1 5]

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[897867908685]

8 建立单元矩阵B回答关问题
B{11}1
B{12}'Brenden'
B{21}reshape(1933)
B{22}{123425421342367}
(1) size(B)ndims(B)值分少？
答：size(B) 值2 2
ndims(B) 值2
(2) B(2)B(4)值分少？
答：B(2)B(4)
(3) 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 2 5 8 3 6 9] {12 34 2 54 21 3 4 23 67}}

题三
1 写出完成列操作命令
(1) 建立3阶单位矩阵A
答：Aeye(3)
(2) 建立5×6机矩阵A元素[100200]范围机整数
答：round(100+(200100)*rand(56))
(3) 产生均值1方差02500正态分布机数
答：1+sqrt(02)*randn(5100)
(4) 产生A样幺矩阵
答：ones(size(A))
(5) 矩阵A角线元素加30
答：A+eye(size(A))*30
(6) 矩阵A提取角线元素元素构成角阵B
答：Bdiag(diag(A))

2 函数实现方阵左旋90o右旋90o功例原矩阵AA左旋B右旋C

答：
Brot90(A)
Crot90(A1)

3 建立方阵A求A逆矩阵A行列式值验证AA1互逆
答：
Arand(3)*10
Binv(A)
Cdet(A)
先计算B*A计算A*B计算知B*AA*BA·A1 A1·A互逆

4 求面线性方程组解

答：
A[4213121230]
b[2108]
xinv(A)*b
方程组解x
5 求列矩阵角线元素三角阵三角阵秩范数条件数迹
(1) (2)
答：
(1) 取角线元素：
diag(A)
三角阵：
triu(A)
三角阵：
tril(A)
秩：
rank(A)
范数：
norm(A1) norm(A) norm(Ainf)
条件数：
cond(A1) cond(A2) cond(Ainf)
迹：
trace(A)
(2)请参考(1)

6 求矩阵A特征值相应特征量

答：
[VD]eig(A)
题四
1 键盘输入4位整数规加密输出加密规：位数字加7然10余数取代该数字第位第三位交换第二位第四位交换
答：
ainput('请输入4位整数：')
A[a1000a100a10a]
Afix(rem(A10))
Arem(A+710)
bA(3)*1000+A(4)*100+A(1)*10+A(2)
disp(['加密值：'num2str(b)])

2 分if语句switch语句实现计算中abc值键盘输入

答：(1) if语句实现计算：
ainput('请输入a值：')
binput('请输入b值：')
cinput('请输入c值：')
xinput('请输入x值：')
if x>05 & x<15
ya*x^2+b*x+c
end
if x>15 & x<35
ya*((sin(b))^c)+x
end
if x>35 & x<55
ylog(abs(b+cx))
end
disp(['y'num2str(y)])

(2) switch语句实现计算：
ainput('请输入a值：')
binput('请输入b值：')
cinput('请输入c值：')
xinput('请输入x值：')
switch fix(x05)
case {12}
ya*x^2+b*x+c
case num2cell(36)
ya*((sin(b))^c)+x
case num2cell(710)
ylog(abs(b+cx))
end
disp(['y'num2str(y)])

3 产生20两位机整数输出中均值偶数
答：
Afix(10+89*rand(120))
sum0
for i120
sumsum+A(i)
end
BA(find(A<(sum20)))
CB(find(rem(B2)0))
disp(C)

4 输入20数求中数数求分循环结构调MATLABmax函数min函数实现
答：
(1) 循环结构实现：
v_max0
v_min0
for i120
xinput(['请输入第' num2str(i) '数：'])
if x> v_max
v_maxx
end
if x< v_min
v_minx
end
end
disp(['数：' num2str(v_max)])
disp(['数：' num2str(v_min)])

(2) max函数min函数实现：
for i15
A(i)input(['请输入第' num2str(i) '数：'])
end
disp(['数：' num2str(max(A))])
disp(['数：' num2str(min(A))])

5 已知：分循环结构调MATLABsum函数求s值
答：
(1) 循环结构实现：
s0
for i063
ss+2^i
end
s

(2) 调sum函数实现：
s063
s2^s
sum(s)

6 n分取100100010000时求列式值
(1)
(2)
(3)
(4)
求分循环结构量运算（sumprod函数）实现
答：
(1) 循环结构实现：
sum0
for k1100
sumsum+(1)^(k+1)k
end
sum
sum函数：
x[]
for k110000
x[x (1)^(k+1)k]
end
sum(x)

(2) 循环结构实现：
sum0
for k1100
sumsum+(1)^(k+1)(2*k1)
end
sum
sum函数：
x[]
for k1100
x[x (1)^(k+1)(2*k1)]
end
sum(x)

(3) 循环结构实现：
sum0
for k1100
sumsum+1(4^k)
end
sum
sum函数实现：
x[]
for k1100
x[x 1(4^k)]
end
sum(x)

(4) 循环结构实现：
t1
for k1100
tt*(((2*k)*(2*k))((2*k1)*(2*k+1)))
end
t
prod函数实现：
x[]
for k1100
x[x ((2*k)*(2*k))((2*k1)*(2*k+1))]
end
prod(x)

7 编写函数文件求意然数n斐波契（Fibnacci）数列项斐波契数列定义：

答：
function xfibnacci(n)
for i1n
if i<2
x(i)1
else
x(i)x(i1)+x(i2)
end
end

8 编写函数文件求两矩阵积点然命令文件中调该函数
答：
函数文件myfncm：
function [x y] myfnc(A B)
try
xA*B
catch
x[]
end
yA*B
命令文件myexem：
Ainput('请输入矩阵A：')
Binput('请输入矩阵B：')
[x y]myfnc(A B)
if length(x)0
display('两矩阵维数匹配法进行积运算')
else
disp('矩阵A矩阵B积：')
x
end
disp('矩阵A矩阵B点：')
y

9 先函数递调定义函数文件求然调该函数文件求
答：
函数文件myfncm：
function summyfnc(n m)
if n<1
sum1
else
sum myfnc (n1 m)+n^m
end
命令窗口中调myfncm文件计算：
summyfnc(100 1)+ myfnc(50 2)+myfnc(101)

10 写出列程序输出结果
① s0
a[121314151617181920212223]
for ka
for j14
if rem(k(j)2)~0
ss+k(j)
end
end
end
s
答：执行结果
s108
② 命令文件exem执行结果：
x
4 12 20
y
2 4 6


文档香网(httpswwwxiangdangnet)户传

《香当网》用户分享的内容，不代表《香当网》观点或立场，请自行判断内容的真实性和可靠性！
该内容是文档的文本内容，更好的格式请下载文档