实验1matlab运算基础,实验1MATLAB运算基础实验报告模板.doc

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

长沙理工大学(学生所在学院和专业) 学生姓名

展开阅读全文