MATALB与信号处理
实 验 报 告
实验项目名称: MATALB运算基础
学 院:____交通运输工程学院___
专 业:_________道路工程_______
学 号:___201610010415_____
班 级:____道路4班____________
报 告 人:____田然_____________
指 导 老 师:___ 胡双红________
实 验 时 间:_________2017.9.16_____
实验目的:
1. 熟悉启动和退出MATLAB的方法。
2. 掌握MATLAB命令窗口的组成。
3. 掌握建立矩阵的方法。
4. 掌握MATLAB各种表达式的书写规则以及常用函数的使用。
实验内容与要求:
实验1:
求下列表达式的值,然后显示MATLAB工作空间的使用情况(whos函数)并保存全部变量(save函数)。
(1)
提示:自然对数e用exp(1)函数产生
(2) ,其中
提示:ln运算用log()函数计算 根号运算用sqrt()函数计算
(3) ,
提示:利用冒号表达式生成a向量,求各点的函数值时用点乘运算
(4) ,其中
提示:用逻辑表达式求分段函数值。
实验2:
已知 ,
求下列表达式的值:
(1) A+6*B 和 A-B+I(其中I为单位矩阵)
(2) A*B 和A.*B
(3) A^3和A.^3
(4) A/B 及 B/A
(5) [A,B]和[A([1,3],:); B^2]
提示:次方运算用^
注意点运算和矩阵运算的区别
实验3:
设有矩阵A和B
,
(1) 求它们的乘积C。
(2) 将矩阵C的右下角32子矩阵赋给D(P25)。
(3) 查看MATLAB工作空间的使用情况。
实验4:
完成下列操作:
(1) 求[100,999]之间能被21整除的数的个数。
提示:先利用冒号表达式(P23),
再利用find(P33)和length(P24)函数。
(2) 建立一个字符串向量,删除其中的大写字母。
提示:利用find函数和空矩阵(P33)。
实验过程或源代码:
实验1:
代码: 1) z1=(2*sin(85*pi/180))/(1+exp(2))
z1 =
0.2375
2)x=[2,1+2i;-0.45,5]
2) x =
2.0000 + 0.0000i 1.0000 + 2.0000i
-0.4500 + 0.0000i 5.0000 + 0.0000i
>> z2=log(x+sqrt(1+x^2))/(2)
z2 =
0.7114 - 0.0253i 0.8968 + 0.3658i
0.2139 + 0.9343i 1.1541 - 0.0044i
3) a=-3:0.1:3
a =
Columns 1 through 13
-3.0000 -2.9000 -2.8000 -2.7000 -2.6000 -2.5000 -2.4000 -2.3000 -2.2000 -2.1000 -2.0000 -1.9000 -1.8000
z31=(exp(0.3*a)-exp(-0.3*a))/(2)
z31 =
Columns 1 through 13
-1.0265 -0.9840 -0.9423 -0.9015 -0.8615 -0.8223 -0.7838 -0.7461 -0.7090 -0.6725 -0.6367 -0.6014 -0.5666
z32=sin(a+0.3)
z32 =
Columns 1 through 13
-0.4274 -0.5155 -0.5985 -0.6755 -0.7457 -0.8085 -0.8632 -0.9093 -0.9463 -0.9738 -0.9917 -0.9996 -0.9975
z33=log((0.3+a)/2)
z33 =
Columns 1 through 7
0.3001 + 3.1416i 0.2624 + 3.1416i 0.2231 + 3.1416i 0.1823 + 3.1416i 0.1398 + 3.1416i 0.0953 + 3.1416i 0.0488 + 3.1416i
z3=z32.*z32+z33
z3 =
Columns 1 through 7
0.4828 + 3.1416i 0.5281 + 3.1416i 0.5813 + 3.1416i 0.6386 + 3.1416i 0.6958 + 3.1416i 0.7490 + 3.1416i 0.7939 + 3.1416i
Columns 8 through 14
0.8268 + 3.1416i 0.8442 + 3.1416i 0.8430 + 3.1416i 0.8209 + 3.1416i 0.7760 + 3.1416i 0.7073 + 3.1416i 0.6144 + 3.1416i
Columns 15 through 21
0.4977 + 3.1416i 0.3579 + 3.1416i 0.1964 + 3.1416i 0.0149 + 3.1416i -0.1849 + 3.1416i -0.4017 + 3.1416i -0.6348 + 3.1416i
Columns 22 through 28
-0.8852 + 3.1416i -1.1564 + 3.1416i -1.4578 + 3.1416i -1.8098 + 3.1416i -2.2631 + 3.1416i -2.9858 + 3.1416i -37.0245 + 0.0000i
Columns 29 through 35
-2.9858 + 0.0000i -2.2631 + 0.0000i -1.8098 + 0.0000i -1.4578 + 0.0000i -1.1564 + 0.0000i -0.8852 + 0.0000i -0.6348 + 0.0000i
Columns 36 through 42
-0.4017 + 0.0000i -0.1849 + 0.0000i 0.0149 + 0.0000i 0.1964 + 0.0000i 0.3579 + 0.0000i 0.4977 + 0.0000i 0.6144 + 0.0000i
Columns 43 through 49
0.7073 + 0.0000i 0.7760 + 0.0000i 0.8209 + 0.0000i 0.8430 + 0.0000i 0.8442 + 0.0000i 0.8268 + 0.0000i 0.7939 + 0.0000i
Columns 50 through 56
0.7490 + 0.0000i 0.6958 + 0.0000i 0.6386 + 0.0000i 0.5813 + 0.0000i 0.5281 + 0.0000i 0.4828 + 0.0000i 0.4487 + 0.0000i
Columns 57 through 61
0.4288 + 0.0000i 0.4254 + 0.0000i 0.4400 + 0.0000i 0.4734 + 0.0000i 0.5257 + 0.0000i
4) t=0:0.5:2.5
t =
0 0.5000 1.0000 1.5000 2.0000 2.5000
t>=0&t<1
ans =
1 1 0 0 0 0
>> t>=1&t<2
ans =
0 0 1 1 0 0
>> t>=2&t<3
ans =
0 0 0 0 1 1
z41=(t.^2).*(t>=0&t<1)
z41 =
0 0.2500 0 0 0 0
>> z42=(t.^2-1).*(t>=1&t<2)
z42 =
0 0 0 1.2500 0 0
z43=(t.^2-2*t+1).*(t>=2&t<3)
z43 =
0 0 0 0 1.0000 2.2500
>> z4=z41+z42+z43
z4 =
0 0.2500 0 1.2500 1.0000 2.2500
实验2:
代码: a=[12,34,-4;34,7,87;3,65,7]
a =
12 34 -4
34 7 87
3 65 7
>> b=[1,3,-1;2,0,3;3,-2,7]
b =
1 3 -1
2 0 3
3 -2 7
>> i=eye(3)
i =
1 0 0
0 1 0
0 0 1
>> a+6*b
ans =
18 52 -10
46 7 105
21 53 49
>> a-b+i
ans =
12 31 -3
32 8 84
0 67 1
>> a*b
ans =
68 44 62
309 -72 596
154 -5 241
>> a.*b
ans =
12 102 4
68 0 261
9 -130 49
>> a^3
ans =
37226 233824 48604
247370 149188 600766
78688 454142 118820
>> a.^3
ans =
1728 39304 -64
39304 343 658503
27 274625 343
>> a/b
ans =
16.4000 -13.6000 7.6000
35.8000 -76.2000 50.2000
67.0000 -134.0000 68.0000
>> b/a
ans =
0.1027 -0.0062 -0.0069
0.0617 0.0403 -0.0366
0.0205 0.0855 -0.0507
>> [a,b]
ans =
12 34 -4 1 3 -1
34 7 87 2 0 3
3 65 7 3 -2 7
>> [a([1,3],:);b^2]
ans =
12 34 -4
3 65 7
4 5 1
11 0 19
20 -5 40
实验3:
代码 a=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25]
a =
1 2 3 4 5
6 7 8 9 10
11 12 13 14 15
16 17 18 19 20
21 22 23 24 25
>> b=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11]
b =
3 0 16
17 -6 9
0 23 -4
9 7 0
4 13 11
c=a*b
c =
93 150 77
258 335 237
423 520 397
588 705 557
753 890 717
d=c(3:5,2:3)
d =
520 397
705 557
890 717
实验4:
代码: a=100:1:999
1)
a =
Columns 1 through 23
100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122
b=find(mod(a,21)==0)
b =
Columns 1 through 23
6 27 48 69 90 111 132 153 174 195 216 237 258 279 300 321 342 363 384 405 426 447 468
Columns 24 through 43
489 510 531 552 573 594 615 636 657 678 699 720 741 762 783 804 825 846 867 888
>> length(b)
ans =
43
2) ch=abc123d4e56fg9
ch =
abc123d4e56fg9
>> k=find(ch>=a&ch<=z);
>> ch(k)=[]
ch =
1234569
实验结论及心得:
将式子转化为MATLAB语言很重要
15
长沙理工大学(学生所在学院和专业) 学生姓名
展开阅读全文