MATLAB程序设计与应用刘卫国编课后实验答案81653Word文档下载推荐.docx
- 文档编号:21746205
- 上传时间:2023-02-01
- 格式:DOCX
- 页数:79
- 大小:556.48KB
MATLAB程序设计与应用刘卫国编课后实验答案81653Word文档下载推荐.docx
《MATLAB程序设计与应用刘卫国编课后实验答案81653Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《MATLAB程序设计与应用刘卫国编课后实验答案81653Word文档下载推荐.docx(79页珍藏版)》请在冰豆网上搜索。
1416i
0、7871+3、
7913+3、1416i
Columns5through8
0、7822+3、1416i
0、7602+3、
0、7254+3、
6784+3、1416i
Columns9through12
0、6196+3、1416i
0、5496+3、
0、4688+3、
3780+3、1416i
Columns13through16
0、2775+3、1416i
0、1680+3、
0、0497+3、
-0、
0771+3、1416i
Columns17through20
-0、2124+3、1416i+3、1416i
-0、3566+3、
-0、5104+3、1416i
-0、6752
Columns21through24
-0、8536+3、1416i+3、1416i
-1、0497+3、
-1、2701+3、1416i
-1、5271
Columns25through28
-1、8436+3、1416i
-2、2727+3、
-2、9837+3、1416i-37、0245
Columns29through32
-3、0017
-2、3085
-1、8971
-1、5978
Columns33through36
-1、3575
-1、1531
-0、9723
-0、8083
Columns37through40
-0、6567
-0、5151
-0、3819
-0、2561
Columns41through44
-0、1374-
■0、0255
0、0792
0、1766
Columns45through48
0>
2663
0、3478
4206
4841
Columns49through52
0、5379
0、5815
0、6145
6366
Columns53through56
0、6474
0、6470
6351
0、6119
Columns57through60
0、5777
0、5327
4774
4126
Column61
0、3388
z4=
00、2500
1、25001、0000
2、2500
2、已知:
12
34
4
A34
7
87,
B2
65
求下列表达式的值:
(1)A+6*B与A-B+l(其中I为单位矩阵)
(2)A*B与A、*B
(3)AA3与A、八3
⑷A/B及B\A
(5)[A,B]与[A([1,3],:
);
BA2]
A=[1234-4;
34787;
3657];
B=[13-1;
203;
3-27];
A+6、*B
A-B+eye(3)
A*B
A、*B
AA3
A、A3
A/B
B\A
[A,B]
[A([1,3],:
A+6、
*B
A-B+eye(3)
B\A[A,B]
ans=
18
52
-10
46
105
21
53
49
31
-3
32
8
84
67
68
44
62
309
-72
596
154
-5
241
102
261
9
-130
37226
233824
48604
247370
149188
600766
78688
454142
118820
1728
39304
-
64
343
658503
27
274625
16、
4000
-13、
6000
7、6000
35、
8000
-76、
2000
50、2000
67、
0000-134、
0000
68、0000
109、
4000-
131、
322、8000
-53、
85、
0000-171、0000
-61、
89、
8000-186、2000
-4
-1
87
-2
11
19
20
40
3、
设有矩阵A与B
16
6
10
17
A11
13
14
15,
B0
23
22
24
25
(1)
求它们的乘积
C。
(2)将矩阵C的右下角3X2子矩阵赋给D
(3)查瞧MATLAB工作空间的使用情况。
解:
、运算结果:
E=(reshape(1:
1:
25,5,5))'
;
F=[3016;
17-69;
023-4;
970;
41311];
C=E*F
H=C(3:
5,2:
3)
C=
93
150
77
258
335
237
423
520
397
588
705
557
753
890
717
H
4、完成下列操作:
(1)求[100,999]之间能被21整除的数的个数。
(2)建立一个字符串向量,删除其中的大写字母。
(1)结果:
m=100:
999;
n=find(mod(m,21)==0);
length(n)
43
(2)、建立一个字符串向量例如:
ch二'
ABC123d4e56Fg9:
则要求结果就是:
ch='
ABC123d4e56Fg9'
k=find(ch>
='
A'
&
ch<
Z'
ch(k)=[]ch=
123d4e56g9
实验二MATLAB矩阵分析与处理
1、设有分块矩阵AE'
R32
,其中
E、R、0、S分别为单位矩阵、随机矩阵、零矩
。
23
S22
阵与对角阵,试通过数值计算验证
A2
E
RRS
O
S2
M文件如下
输岀结果:
S=
10
02
A=
1、
00
、5383
、4427
、0000
、9961
、1067
0、0782
、9619
01
、
、0767
1、3280
、9923
0、3200
0、1564
2、8857
由ans,所以A
ERRS
Os2
2、产生5阶希尔伯特矩阵H与5阶帕斯卡矩阵P,且求其行列式的值Hh与Hp以及它们的条件数Th与Tp,判断哪个矩阵性能更好。
为什么?
M文件如下:
输出结果:
H=
1、
5000
3333
2500
1667
1429
0、2500
0、1667
0、1429
0、1250
0、2000
0、1250
0、1111
P=
11
12
13
15
14
35
15
70
Hh=
3、7493e-012
Hp=
Th=
4、7661e+005
Tp=
8、5175e+003
因为它们的条件数Th»
Tp,所以pascal矩阵性能更好。
3、建立一个5X5矩阵,求它的行列式值、迹、秩与范数。
解:
输出结果为
5070000
cl=
6、8500
c2=
5、4618
cinf=
4、已知
29
A20
求A的特征值及特征
向量
并分析其数学意义。
M文件如图:
7130
0、2803
0、2733
6084
-0、7867
0、8725
3487
0、5501
0、4050
-25、
3169
-10、5182
16、8351
数学意义:
V的3个列向量就是A的特征向量,D的主对角线上3个就是A的特征值,特别的,V
的3个列向量分别就是D的3个特征值的特征向量。
5、下面就是一个线性方程组
片
0.95
0.67
X2
X3
0.52
(1)求方程的解。
⑵将方程右边向量元素b3改为0、53再求解拼比较b3的变化与解的相对变化。
(3)计算系数矩阵A的条件数并分析结论。
输出结果:
X=
1、2000
0、6000
0、6000X2=
1、3533e+003
由结果,X与X2的值一样,这表示b的微小变化对方程解也影响较小,而A的条件数算得较小
所以数值稳定性较好,A就是较好的矩阵。
6、建立A矩阵,试比较sqrtm(A)与sqrt(A),分析它们的区别。
BEditor-胆]|乂|
EiltWitlest血tellItolED«
buED«
sktopItindo^MsId■JrX
scriptLn5Col3
运行结果有
bl
=
8891
-0、1102
2103
2917
2、1436
3698
3855
2、0760
7305
b2
4、
2、4495
2426
4721
2、2361
4641
2、8284
2、
2361
b=
6、0000
18、
20、
5、0000
12、
9、
8、0000
5、
分析结果知:
sqrtm(A)就是类似A的数值平方根(这可由b1*b1=A的结果瞧出),而sqrt(A)则就是对A中的每个元素开根号,两则区别就在于此。
实验三选择结构程序设计
、实验目的
1、掌握建立与执行M文件的方法。
2、掌握利用if语句实现选择结构的方法。
3、掌握利用switch语句实现多分支选择结构的方法。
4、掌握try语句的使用。
、实验内容
求分段函数的值。
2x
x6
x0且x3
2yx
5x6
0x5且x2及x3
x1
其他
用if语句实现,分别输出x=-5、0,-3、0,1、0,2、0,2、5,3、0,5、0时的y值。
解:
运算结果有
f(-5)
y=
14>
>
f(-3)y=
11>
f(i)y=
f⑵
f(2、5)
-0、2500
f(3)
f(5)
2、输入一个百分制成绩,要求输出成绩等级A、B、C、D、E。
其中90分~100分为A,80分~89分为B,79分~79分为C,60分~69分为D,60分以下为E。
要求:
(1)分别用if语句与switch语句实现。
(2)输入百分制成绩后要判断该成绩的合理性,对不合理的成绩应输出出错信息。
试算结果:
score=88
grade=
Bscore=123
错误:
输入的成绩不就是百分制成绩
3、硅谷公司员工的工资计算方法如下:
(1)工作时数超过120小时者,超过部分加发15%。
(2)工作时数低于60小时者,扣发700元。
(3)其余按每小时84元计发。
试编程按输入的工号与该号员工的工时数,计算应发工资。
M文件下
4、设计程序,完成两位数的加、减、乘、除四则运算,即产生两个两位随机整数,再输入一个运算符号,做相应的运算,并显示相应的结果。
M文件如下;
运算结果例:
a=
38
33
输入一个运算符:
人
c=
false
92
输入一个运算符:
+
5、建立5x6矩阵,要求输出矩阵第n行元素。
当n值超过矩阵的行数时,自动转为输出矩阵最后一行元素,并给出出错信息。
M文件如下:
运算结果如下
输入一个5行6列矩阵A=[123455;
234576;
222223;
1123973;
234567]输入一正整数n=4
1123973
234567]输入一正整数n=6
234567ans=
Errorusing==>
disp
Toomanyinputarguments、
实验四循环结构程序设计
1、掌握利用for语句实现循环结构的方法。
2、掌握利用while语句实现循环结构的方法。
3、熟悉利用向量运算来代替循环操作的方法。
1、根据一
III
~~2
n
,求n的近似值。
n分别取100、1000、10000
时,结果就是多少?
要求:
分别用循环结构与向量运算(使用sum函数)来实现。
运行结果如下
K>
%循环结构计算pi值
y=o;
n=input('
n='
fori=1:
n
y=y+1/i/i;
end
pi=sqrt(6*y)
n=100
pi=
3、1321
n=1000
3、1406
n=10000
3、1415
%向量方法计算Pi值
i=1、/(1:
n)、A2;
s=sum(i);
pi=sqrt(6*s)
pi=3、1415
111
2、根据y1,求:
35"
I2n1
⑴y<
3时的最大n值。
⑵与
(1)的n值对应的y值。
M—文件如下:
春-UntitledSr匚^叵|冈
y=0;
n=0;
whiley<
n=n+1;
y=y+1/(2*n-1);
y
ify>
n=n-1;
3、0033
n=
57n=
56
3、考虑以下迭代公式:
a
Xn1
bXn
其中a、b为正的学数。
(1)编写程序求迭代的结果,迭代的终止条件为|Xn+1-Xn|W10-5,迭代初值X0=1、0,迭代次数不超过500次。
(2)如果迭代过程收敛于r,那么r的准确值就是——一,当(a,b)的值取(1,1)、
(8,3)、(10,0、1)时,分别对迭代结果与准确值进行比较。
请输入正数a=1
请输入正数b=1
x=
0、6180
r=
-4、
7016
-1、
6180
s=
-0、0000
-2、
请输入正数a=8
请输入正数b=3
1、7016
0.0-6、4031
请输入正数a=10
请输入正数b=0
、1
3、1127
3、1127
-4、7016
-3、2127
-6、3254
f1
f2
f3
fn
fn12fn2
fn3n
求f1~f100中:
(1)最大值、最小值、各数之与。
(2)正数、零、负数的个数。
M—文件
以下就是运算结果
max(f)=437763282635min(f)=-899412113528sum(f)=-7427456019
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 程序设计 应用 卫国 课后 实验 答案 81653