一、MATLAB基础知识整理

最新推荐文章于 2024-08-28 14:17:10 发布

渡鸦摸鱼哥

最新推荐文章于 2024-08-28 14:17:10 发布

阅读量3.8k

点赞数 4

分类专栏： matlab 文章标签： matlab 矩阵算法傅立叶分析

本文链接：https://blog.csdn.net/weixin_44583303/article/details/108097329

版权

matlab 专栏收录该内容

1 篇文章 1 订阅

订阅专栏

一、MATLAB基础知识

1.1 matlab系统环境

1.续行符 …

2.如何设置当前文件夹？

(1)在当前文件夹工具栏或当前文件夹窗口中选择某文件夹为当前文件夹

(2)使用cd命令

***一定要先建立文件夹，再将这个文件夹设为当前文件夹

3.matlab搜索路径

变量 -> 内部函数 -> 程序文件 -> 当前文件夹下的程序文件

-> 文件搜索路径文件夹中的程序文件

用path命令设置文件搜索路径。例如：

>>path(path,'d:\work')

1.2数值数据类型

1.数值数据类型的分类

*整型

*浮点型

*复数型

(1)整型

无符号整数：无符号8位整数、无符号16位整数、无符号32位整数、无符号64位整数
无符号8位整数数据范围：00000000~11111111 (0~2^8-1) (0~255)
带符号整数：带符号8位整数、带符号16位整数、带符号32位整数、带符号64位整数
带符号8位整数数据范围：10000000~01111111 (-2^7~27-1) (-128~127)

MATLAB提供了数据类型转换函数来实现数据类型的转换，例如:

uint8函数将数值数据转换为无符号8位函数，int8函数将数值数据转换为带符号8位整数

实例：

>> x=int8(129)		% 调用int8函数，把129转换为带符号整数

>> x=127							%由于带符号8位整数最大值为127，所以转换后取最大值127

>> x=uint(129)					%调用uint8函数，把129转化为无符号整数

>> x=129							%无符号8位整数最大值为255，所以可以转换为129

(2)浮点型

浮点型数据分为单精度型和双精度型，单精度型实数在内存中占用4个字节，而双精度型实数在内存中占用8个字节，所以双精度型的数据精度更高。
single函数：将其它类型的数据转换为单精度型。
double函数：将其它类型的数据转换为双精度型。

实例：

>> class(4)

>> ans=double			%调用class函数得到4的数据类型

>> class(single(4))		%得到4的数据类型为single

>> ans=single

(3)复型

复型数据包括实部和虚部两个部分，实部和虚部默认为双精度型，虚数单位用i或j来表示。
real函数：求复数的实部
imag函数：求复数的虚部

实例;

>> 6+5i

>> ans=6.0000 + 5.0000i

>> 6+5j

>> ans=6.0000 + 5.0000i					%表示同一个数据

2.数值数据的输出格式

format命令的格式：format格式符

>> format long			%输出格式设置为long格式
>> 50/3					%输出50/3的数值
ans =
16.666666666666668
>> format				%返回默认格式
>> 50/3
ans =
16.6667

***注意：***format命令只影响数据输出格式，而不影响数据的计算和存储。

3.常用数学函数

(1)函数的调用格式为

函数名(函数自变量的值)

函数的自变量规定为矩阵变量，当然也可以是标量，标量本身是矩阵的一种特例。

函数在运算时是将函数逐项作用于矩阵的每个元素上，所以最后运算的结果就是一个自变量同型的矩阵。

>> A=[4,2;3,6]			%定义两行两列的矩阵
A =
4 2
3 6
>> B=exp(A)				%调用exp()函数
B =
54.5982 7.3891
20.0855 403.4288		%求自然指数所得

(2)常用函数的应用

① 三角函数有以弧度为单位的函数和以角度为单位的函数，如果是以角度为单位的函数就在函数名后面加“d”，以示区别。

>> sin(pi/2)			%求二分之π弧度的正弦函数值
ans =
1
>> sind(90)				%求90度的正弦函数
ans =	
1

② abs函数可以求实数的绝对值、复数的模、字符串的ASCII码值。

>> abs(-4)				%求-4的绝对值
ans =
4
>> abs(3+4i)			%求复数的模
ans =
5
>> abs('a')				%求字符串的ASCII值
ans =
97

③用于取整的函数有fix、floor、ceil、round。

>> round(4.7)			%round函数是按照四舍五入的规则来取整
ans =
5
>> fix(-3.2)			%ceil是向上取整，取大于等于这个数的第一个整数。
ans =
-3
>> floor(3.6)			%floor是向下取整，取小于等于这个数的第一个整数。
ans =
3
>> ceil(-3.8)			%fix是固定取靠近0的那个整数，也就是舍去小数取整。
ans =
-3

④函数应用举例。

分别求一个三位正整数的个位数字、十位数字和百位数字。

>> m=345;
>> m1=rem(m,10)				%调用rem函数，求m的余数
m1 =
5
>> m2=rem(fix(m/10),10)		%fix取整
m2 =
4
>> m3=fix(m/100)
m3 =
3

求[1，100]区间的所有素数。

>> x=1:100;							%生成1-100内的全部整数组成的向量x
>> k=isprime(x);					%调用isprime函数生成k向量，k向量中要么为1要么为0，为1还是为0取决于x向量对应位置的元素是否为素数
>> k1=find(k);						%调用find函数生成k1向量，k1向量中的元素为k向量中非0向量的序号也就是x中素数元素的序号
>> p=x(k1)							%输出x中全部素数
p =
1 至 13 列
2 3 5 7 11 13 17 19 23 29 31 37 41
14 至 25 列
43 47 53 59 61 67 71 73 79 83 89 97

**matlab函数大全参见：**https://wenku.baidu.com/view/387dbb85960590c69ec3768d.html

1.3变量及其操作

变量与赋值语句

变量本质上讲的是内存单元的一个抽象
在matlab中，变量名是以字母开头，后接字母、数字或下划线的字符串序列，最多63个字符。
变量名区分字母的大小写。
标准函数名以及命令必须用小写字母。

赋值语句的两种格式：

变量=表达式 %赋值给左边的变量
表达式 %赋值给预定义变量ans
如果在表达式后边加上分号，就仅仅执行赋值操作，不会再显示运算后变量的结果

例计算表达式 5+cos 47° 1+ 𝑥−𝑦 的值，并将结果赋给变量z，然后显示计算结果。其中，𝑥 = √7 − 2𝑖, 𝑦 = 𝑒 ^(π/2)。

>> x=sqrt(7)-2i;			%调用求平方根函数sqrt
>> y=exp(pi/2);				%调用自然指数函数exp
>> z=(5+cosd(47))/(1+abs(x-y))
z =
1.4395

预定义变量

预定义变量是由系统本身定义的变量

ans是默认赋值变量
i和j代表虚数单位
pi代表圆周率
NaN代表非数

变量的管理

（1）内存变量的删除与修改

在工作区直接右击进行修改
命令窗口使用who命令与whos命令查看

（2）内存变量文件

用于保存matlab工作区变量的文件叫做内存变量文件，其扩展名为.mat,也叫MAT文件。

生成变量文件：
- save命令：创建内存变量文件。
- load命令：装入内存变量文件。
- ```
>> save mydata a x
>> load mydata
```

1.4 matlab矩阵的表示

1.矩阵的建立

（1）利用直接输入法建立矩阵:将矩阵的元素用中括号括起来，按矩阵行的顺序输入各元素，同一行的各元素之间用逗号或空格分隔，不同行的元素之间用分号分隔。

>> A=[1,2,3;4,5,6;7,8,9]
A =
1 2 3
4 5 6
7 8 9

（2）利用已建好的矩阵建立更大的矩阵: 一个大矩阵可以由已经建立好的小矩阵拼接而成

>> A=[1,2,3;4,5,6;7,8,9];
>> B=[-1,-2,-3;-4,-5,-6;-7,-8,-9];
>> C=[A,B;B,A]
C =
1 2 3 -1 -2 -3
4 5 6 -4 -5 -6
7 8 9 -7 -8 -9
-1 -2 -3 1 2 3
-4 -5 -6 4 5 6
-7 -8 -9 7 8 9

（3）可以用实部矩阵和虚部矩阵构成复数矩阵

>>B=[1,2,3;4,5,6];
>>C=[6,7,8;9,10,11];
>>A=B+i*C
A =
1.0000 + 6.0000i 2.0000 + 7.0000i 3.0000 + 8.0000i
4.0000 + 9.0000i 5.0000 +10.0000i 6.0000 +11.0000i

2.冒号表达式

（1）格式：e1:e2:e3 e1为初始值，e2为步长，e3为终止值

>> t=0:1:5			%初始值为0，步长为1，终止值为5
t =
0 1 2 3 4 5

省略步长e2，则步长为1。例如， t=0:5与t=0:1:5等价。

（2）格式：linspace(a,b,n) a为第一个元素，b为最后一个元素，n为元素总数

>> x=linspace(0,pi,6)			%跪一个元素为0，最后一个元素为π，元素总数为6
x =
0 0.6283 1.2566 1.8850 2.5133 3.1416

当n省略时，自动产生100个元素。

3.结构矩阵

（1）结构矩阵

由结构数据构成的矩阵就是结构矩阵，结构矩阵里的每个元素就是结构数据类型。

格式为：结构矩阵元素.成员名=表达式

>> a(1).x1=10; a(1).x2='liu'; a(1).x3=[11,21;34,78];
>> a(2).x1=12; a(2).x2='wang'; a(2).x3=[34,191;27,578];
>> a(3).x1=14; a(3).x2='cai'; a(3).x3=[13,890;67,231];

（2）单元矩阵

建立单元矩阵和一般矩阵相似，直接输入就可以了，只是单元矩阵元素用大括号括起来。

>> b= {10,'liu',[11,21;34,78];12,'wang',[34,191;27,578];14,'cai',...[13,890;67,231]}
b =
[10] 'liu' [2x2 double]
[12] 'wang' [2x2 double]
[14] 'cai' [2x2 double]

1.5矩阵元素的引用

1.矩阵元素的引用方式

（1）通过下标来引用矩阵的元素

A(3,2)表示A矩阵第3行第2列的元素。

>> A(3,2)=200

>> A=[1,2,3;4,5,6];
>> A(4,5)=10
A =
1 2 3 0 0
4 5 6 0 0
0 0 0 0 0
0 0 0 0 10

**注意：**如果给出的行下标或列下标大于原来矩阵的行数和列数，那么MATLAB将自动扩展原来的矩阵，并将扩展后没有赋值的矩阵元素置为0。

（2）通过序号来引用

在MATLAB中，矩阵元素按列存储，即首先存储矩阵的第一列元素，然后存储第二列元素，…，一直到矩阵的最后一列元素。
矩阵元素的序号就是矩阵元素在内存中的排列顺序。

>> A=[1,2,3;4,5,6]
A =
1 2 3
4 5 6
>> A(3)
ans =
2

序号与下标是一一对应的，以m×n矩阵A为例，矩阵元素A(i,j)的序号为(j-1)×m+i

*矩阵元素的序号与下标可以利用sub2ind和ind2sub函数实现相互转换。

*sub2ind函数：将矩阵中指定元素的行、列下标转换成存储的序号。调用格式为：

D=sub2ind(S,I,J) S为行数和列数组成的向量，I为转换矩阵元素的行下标，J为转换矩阵元素的列下标

>> A=[1:3;4:6]
A =
1 2 3
4 5 6
>> D=sub2ind(size(A),[1,2;2,2],[1,1;3,2])		%所对应的序号为D中所示
D =
1 2
6 4

*ind2sub函数：将把矩阵元素的序号转换成对应的下标，调用格式为：

[I,J]=ind2sub(S,D) I为行下标，J为列下表，S为行数和列数组成的向量，D为序号

>> [I,J]=ind2sub([3,3],[1,3,5])			%三行三列矩阵的序号1，序号3，序号5所对应的下标
I =
1 3 2
J =
1 1 2

2．利用冒号表达式获得子矩阵

子矩阵是指由矩阵中的一部分元素构成的矩阵。

A(i,:) 第i行的全部元素
A(:,j) 第j列的全部元素
A(i:i+m,k:k+m) 第i～i+m行内且在第k～k+m列中的所有元素
A(i:i+m,:) 第i～i+m行的全部元素

>> A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15]
A =
1 2 3 4 5
6 7 8 9 10
11 12 13 14 15
>> A(1:2,:)				%第一行和第二行所有元素
ans =
1 2 3 4 5
6 7 8 9 10
>> A(2:3,1:2:5)			%第二行到第三行中第一列第三列和第五列的元素
ans =
6 8 10
11 13 15

end运算符：表示某一维的末尾元素下标。

>>A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20];
>>A(end,:)				%最后一行所有元素
ans =
16 17 18 19 20
>> A([1,4],3:end)		%第一行和第四行中第三列到最后一列的所有元素
ans =
3 4 5
18 19 20

3．利用空矩阵删除矩阵的元素

空矩阵是指没有任何元素的矩阵。

>> x=[]
x =
[]				%X是一个空矩阵。

>> A=[1,2,3,0,0;7,0,9,2,6;1,4,-1,1,8]
A =
1 2 3 0 0
7 0 9 2 6
1 4 -1 1 8
>> A(:,[2,4])=[]			%第2列和第四列赋空矩阵，也就是删除第2列和第4列
A =
1 3 0
7 9 6
1 -1 8

4．改变矩阵的形状

reshape(A,m,n)：在矩阵总元素保持不变的前提下，将矩阵A重新排成m×n的二维矩阵。

**注意：**reshape函数只是改变原矩阵的行数和列数，但并不改变原矩阵元素个数及其存储顺序。

>> x=[23,45,65,34,65,34,98,45,78,65,43,76];
>> y=reshape(x,3,4)
y
23 34 98 65
45 65 45 43
65 34 78 76

**A(😃 ：**将矩阵A的每一列元素堆叠起来，成为一个列向量。

>> A =[-45,65,71;27,35,91]
A =
-45 65 71
27 35 91
>> B=A(:)					%A(:)等价于reshape(A,6,1)。
B =
-45
27
65
35
71
91

1.6matlab基本运算

matlab的运算都是针对矩阵而言的

1.算数运算

（1）基本算数运算

基本算术运算符：＋(加)、－(减)、*(乘)、/ (右除)、 \ (左除)、 ^(乘方)。

MATLAB的算术运算是在矩阵意义下进行的。
单个数据的算术运算只是矩阵运算的一种特例。

**注意：**MATLAB的运算是在矩阵意义下进行的，单个数据的算术运算只是矩阵运算的一种特例。

加减运算

若两矩阵同型，则运算时两矩阵的相应元素相加减。
若两矩阵不同型，则MATLAB将给出错误信息。
一个标量也可以和矩阵进行加减运算，这时把标量和矩阵的每一个元素进行加减运算。

乘法运算

矩阵A和B进行乘法运算，要求A的列数与B的行数相等，此时则称A、B矩阵是可乘的，或称A和B两矩阵维数和大小相容。
如果两者的维数或大小不相容，则将给出错误信息，提示用户两个矩阵是不可乘的。

除法运算

在MATLAB中，有两种矩阵除法运算：右除/和左除\。如果A矩阵是非奇异方阵，则B/A等效于B * inv(A)，A\B等效于inv(A) * B。

>> A=[1,2,3;4,2,6;7,4,9];
>> B=[4,3,2;7,5,1;12,7,92];
>> C1=B/A
C1 =
-0.1667 -3.3333 2.5000
-0.8333 -7.6667 5.5000
12.8333 63.6667 -36.5000
>> C2=A\B
C2 =
0.5000 -0.5000 44.5000
1.0000 0.0000 46.0000
0.5000 1.1667 -44.8333
>> 3/4
ans =
0.7500
>> 4\3
ans =
0.7500
>> a=[10.5,25]
a =
10.5000 25.0000
>> a/5
ans =
2.1000 5.0000
>> 5\a
ans =
2.1000 5.0000

乘方运算

一个矩阵的乘方运算可以表示成A^x，要求A为方阵，x为标量。

>> A=[1,2,3;4,5,6;7,8,0];
>> A^2
ans =
30 36 15
66 81 42
39 54 69

点运算

点运算符：.* 、./ 、.\和.^。
两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵同型。

>> A=[1,2,3;4,5,6;7,8,9];
>> B=[-1,0,1;1,-1,0;0,1,1];
>> C=A.*B
C =
-1 0 3
4 -5 0
0 8 9
>> D=A*B
D =
1 1 4
1 1 10
1 1 16

当x=0.1、0.4、0.7、1时，分别求y=sin x cos x的值。

>> x=0.1:0.3:1;
>> y=sin(x).*cos(x);

2．关系运算

关系运算符：<(小于)、<=(小于或等于)、>(大于)、>=(大于或等于)、==(等于)、～=(不等于)。
当两个比较量是标量时，直接比较两数的大小。若关系成立，关系表达式结果为1，否则为0。

>> 3>4
ans =
0
>> x=5
x =
5
>> x==5
ans =
1

当参与比较的量是两个同型的矩阵时，比较是对两矩阵相同位置的元素按标量关系运算规则逐个进行，最终的关系运算的结果是一个与原矩阵同型的矩阵，它的元素由0或1组成。
当参与比较的一个是标量，而另一个是矩阵时，则把标量与矩阵的每一个元素按标量关系运算规则逐个比较，最终的关系运算的结果是一个与原矩阵同型的矩阵，它的元素由0或1组成。

建立3阶方阵A，判断A的元素是否为偶数。

>> A =[24,35,13;22,63,23;39,47,80]
A =
24 35 13
22 63 23
39 47 80
>> P=rem(A,2)==0
P =
1 0 0
1 0 0
0 0 1

3．逻辑运算

逻辑运算符：&(与)、|(或)和～(非)。设参与逻辑运算的是两个标量a和b，那么运算规则为：

a&b a、b全为非零时，运算结果为1，否则为0。
a|b a、b中只要有一个为非零时，运算结果为1 。
~a 当a为零时，运算结果为1；当a为非零时，运算结果为0。

>> 3<4 & 6>5
ans =
1
>> ~(9==1)
ans =
1
>> ~9==1
ans =
0

运算符优先级

优先级	运算符
1	圆括号 ( )
2	矩阵转置和乘方：转置（.’）、共轭转置（’）、乘方（.^{）、矩阵乘方（}）
3	一元加法（+）、—元减法（-）、取反（~）
4	乘法（.）、矩阵乘法（）、右除（./）、左除（.\）、矩阵右除（/）、矩阵左除（\）
5	加法（+）、减法（-）、逻辑非（~）
6	冒号运算符（:）
7	小于（<）、小于等于（<=）、大于（>）、大于等于（>=）、等于（==）、不等于（~=）
8	逐元素逻辑与（&）
9	逐元素逻辑或（
10	避绕式逻辑与，或者捷径逻辑与（&&）
11	避绕式逻辑或，或者捷径逻辑或（

若参与逻辑运算的是两个同型矩阵，那么将对矩阵相同位置上的元素按标量规则逐个进行运算，最终运算结果是一个与原矩阵同型的矩阵，其元素由1或0组成。
若参与逻辑运算的一个是标量，一个是矩阵，那么将在标量与矩阵中的每个元素之间按标量规则逐个进行运算，最终运算结果是一个与矩阵同型的矩阵，其元素由1或0组成。

水仙花数是指各位数字的立方之和等于该数本身的三位正整数。求全部水仙花数。

>> m=100:999;
>> m1=rem(m,10);
>> m2=rem(fix(m/10),10);
>> m3=fix(m/100);
>> k=find(m==m1.*m1.*m1+m2.*m2.*m2+m3.*m3.*m3)
k =
54 271 272 308
>> s=m(k)
s =
153 370 371 407

1.7字符串处理

1．字符串的表示

在MATLAB中，字符串是用单引号括起来的字符序列。

>> xm='Central South University'
>> xm(1:3)
ans =
Cen

若字符串中的字符含有单引号，则该单引号字符要用两个单引号来表示。

>> 'I''m a teacher.'
ans =
I'm a teacher.

建立多行字符串，形成字符串矩阵。

>> ch=['abcdef';'123456'];		%多行字符串构成了矩阵，所以要求每一行的字符串长度相同
>> ch(2,3)
ans =
3

例题：
- 建立一个字符串向量，然后对该向量做如下处理：
- ① 取第1～5个字符组成的子字符串。
- ② 将字符串倒过来重新排列。
- ③ 将字符串中的小写字母变成相应的大写字母，其余字符不变。
- ④ 统计字符串中小写字母的个数。

>> ch='ABc123d4e56Fg9';
>> subch=ch(1:5)
subch =
ABc12
>> revch=ch(end:-1:1)		%冒号表达式，初始值为最后一个元素
revch =
9gF65e4d321cBA
>> k=find(ch>='a'&ch<='z')		%寻找小写字母的位置
k =
3 7 9 13
>> ch(k)=ch(k)-('a'-'A')		%根据ASCII计算
ch =
ABC123D4E56FG9
>> length(k)
ans =
4

2．字符串的操作

（1）字符串的执行

格式： eval(s)

>> t=pi;
>> m='[t,sin(t),cos(t)]';
>> y=eval(m)
y =
3.1416 0.0000 -1.0000

（2）字符串与数值之间的转换

abs和double函数都可以用来获取字符串矩阵所对应的ASCII码数值矩阵。
char函数可以把ASCII码矩阵转换为字符串矩阵。

>> s1='MATLAB';
>> a=abs(s1)
a =
77 65 84 76 65 66
>> char(a+32)
ans =
matlab

（3）字符串的比较

字符串的比较有两种方法：利用关系运算符或字符串比较函数。

关系运算符比较：两个字符串里的每个字符依次按ASCII值大小逐个进行比较，比较的结果是一个数值向量，向量中的元素要么是1 ，要么是0。

>> 'www0'>='W123'
ans =
1 1 1 0

字符串比较函数用于判断字符串是否相等，有4种比较方式，函数如下：
- strcmp(s1,s2)：用来比较字符串s1和s2是否相等，如果相等，返回结果为1，否则返回0。
- strncmp(s1,s2,n)：用来比较两个字符串前n个字符是否相等，如果相等，返回1，否则返回0。
- strcmpi(s1,s2)：在忽略字母大小写前提下，比较字符串s1和s2是否相等，如果相等，返回1，否则返回0。
- strncmpi(s1,s2,n)：在忽略字符串大小写前提下，比较两个字符串前n个字符是否相等，如果相等，返回1，否则返回0。

>> strcmp('www0', 'w123')
ans =
0
>> strncmpi('Www0','w123',1)
ans =
1

（4）字符串的查找与替换

findstr(s1,s2)：返回短字符串在长字符串中的开始位置。
strrep(s1,s2,s3)：将字符串s1中的所有子字符串s2替换为字符串s3。

>> p=findstr('This is a test!','is')
p =
3 6
>> p=findstr('is','This is a test!')
p =
3 6
>> result=strrep('This is a test!','test','class')
result =
This is a class!