【数学建模】Matlab基础知识

MATLAB基础知识

一、基本常识

注：matlab中使用的是脚本文件，后缀是.m

1.注释

%% 独占一行的注释（有上下横线分割–可以用来分区）
% 普通注释

多行注释：选中，快捷键Ctrl+R
取消多行注释：选中，快捷键Ctrl+T

2.分号

(1)

a = 3;
a = 5

语句后面加上分号（;），表示不显示运行结果，否则会有输出，如上例。

(2) 分号也可以用于区分行。

clear;clc

3.clear和clc

clear：
clear 关键字 —— 删除某个变量
clear —— 删除工作区所有变量

clc ：
清除命令行窗口中所有文本，让屏幕干净。

这两条通常一起使用，用于”初始化“，防止之前的结果影响后续操作/新脚本文件

4.输入输出（disp 和 input）

输出：

a = 100
disp(a)
disp('这是字符串')
b = [1 2 3]
disp(b);
%注：disp函数后面有没有分号都可以！

输入:

% 一般我们会将输入的数、向量、矩阵、字符串等赋给一个变量，这里我们赋给A
A = input('请输入A：');%加不加分号有所区别，加分号会有一个输出
B = input('请输入B：')%

5.判断语句

if所在的行不需要冒号，语句的最后一定要以end结尾 ；中间的语句要注意缩进。

a = input('请输入考试分数:')
if a >= 85  
    disp('成绩优秀')
elseif a >= 60 
    disp('成绩合格')
else
    disp('成绩挂科')
end %表示输入完整了，不然会要求继续输入

二、数据类型

1.数字

可以进行普通运算 加减乘除平方开方指数对数，具体可以在右上角搜索框中查找对应运算符。
如：
3+3
5-1
6*6
8/8

2.字符与字符串

使用 ’ ’ （单引号）或 " " （双引号）表示

str = '字符串‘

关于字符串的函数

(1)合并两个字符串

a. 函数strcat(str1,str2,str3…strn)
eg：strcat('字符串1','字符串2')

b.[’字符串1‘,‘字符串2’]（用逗号或者空格隔开）

(2)将数字转换为字符串

函数 ： num2str()

c = 100
num2str(c)
% 拼接字符串
disp(['c的取值为' num2str(c)])
disp(strcat('c的取值为', num2str(c)))

3.向量（行/列为1的矩阵）

行向量

a = [1 2 3]
%或者 (同一行中用空格或者逗号分割）
a = [1,2,3]

列向量

%分号可以用来分隔每一行的元素
a = [1;2;3] 

4.矩阵

E = [1,2;3,4;5,6]
%同行之间用逗号和空格隔开，不同行之间用分号隔开

4.1sum求和

向量求和

无论行向量还是列向量，都是直接求和

E = [1,2,3]%行向量
sum(E)
E = [1;2;3]%列向量
sum(E)

矩阵求和

矩阵求和需要按行和列的方向做区分

E = [1 2;3 4;5 6]

% a = sum(x) （未经指定，按列求和，得到一个行向量）
  a = sum(E)

% dim = 1 表示按照列求和（得到行向量）
% dim = 2 表示按照行求和（得到列向量）
  a = sum(E,1)
  a = sum(E,2)

a = sum(x(:) % 对整个矩阵求和
E(:)  % 按列向量提取矩阵
a = sum(E(:))

%另一种方法
a = sum(sum(E))

a = sum(E( : ) )

4.2 提取矩阵中指定位置元素

（1）取指定行和列的一个元素（输出一个值）

A=[1 1 4 1/3 3;1 1 4 1/3 3;1/4 1/4 1 1/3 1/2;3 3 3 1 3;1/3 1/3 2 1/3 1]; 
% 矩阵不用非得写成一行！
A
% 上面的选中不方便，这样可以方便输出整个矩阵
A(2,1)
A(3,2)

（2）取指定的某一行的全部元素（输出的是一个行向量）

A
A(2,:)% :表示取所有元素
A(5,:)

（3）取指定的某一列的全部元素（输出的是一个列向量）

A
A(:,1)
A(:,3)

（4）取指定的某些行的全部元素（输出的是一个矩阵）

A
A([2,5],:)      % 只取第二行和第五行（一共2行）
A(2:5,:)        % 取第二行到第五行（一共4行）
A(2:2:5,:)     % 取第二行和第四行 （从2开始，每次递增2个单位，到5结束）

%注 等差数列的生成
1 : 3 : 10
10:-1:1 %(从10开始，每次-1，到1结束）
1:10  %没说的话，等差的差是1

A(2:end,:)      % 取第二行到最后一行
A(2:end-1,:)    % 取第二行到倒数第二行

取全部元素(按列拼接的，最终输出的是一个列向量)

A
A(:)

4.3 size函数

size函数的用法

size(A) % （1）用于求矩阵大小，返回一个行向量，第一个元素是矩阵的行数，第二个元素为列数
[r,c] = size(A) % （2）赋值 将A的行数赋给r A的列数赋给c
r = size(A,1)  % （3）返回行数
c = size(A,2   %      返回列数

4.4 repmat函数

B = repmat(A,m,n):将矩阵A复制m×n块，即把A作为B的元素，B由m×n个A平铺而成。

A = [1,2,3;4,5,6]
B = repmat(A,3,2)

4.5 矩阵运算

矩阵加减乘除

A + B 
A - B
A * B
A / B

求矩阵的逆

inv(B)

A = [1,2,3,4]
inv(A)
A/B = A*inv(B)

两个形状相同的矩阵对应元素之间的乘除法

需要使用 “ .*" 和 " ./"

A = [1,2;3,4]
B = [1,0;1,1]
A .* B
A ./ B %MATLAB中 inf表示无穷

每个元素同时和常数相乘或相除

四种运算符都可以用

A = [1,2;3,4]
A * 2 % 和下面的效果一样
A .* 2
A / 2 % 和下面的效果一样
A ./ 2

每个元素同时乘方

A = [1,2;3,4]
A .^ 2 
A ^ 2 % =A*A
A * A

4.6 求特征值和特征向量

% 在Matlab中，计算矩阵A的特征值和特征向量的函数是eig(A),其中最常用的两个用法：
A = [1 2 3 ;2 2 1;2 0 3]
% （1）E=eig(A)：求矩阵A的全部特征值，构成向量E。%每个特征值对应一组特征向量
E=eig(A)
% （2）[V,D]=eig(A)：求矩阵A的全部特征值，构成对角阵D，并求A的特征向量构成V的列向量。（V的每一列都是D中与之相同列的特征值的特征向量）
[V,D]=eig(A)

4.7 find函数的基本用法

用来返回向量或者矩阵中不为0的元素的位置索引

X = [1 0 4 -3 0 0 0 8 6]
ind = find(X)
% 其有多种用法，比如返回前2个不为0的元素的位置：
ind = find(X,2)


% 如果X是一个矩阵（二维，有行和列）
X = [1 -3 0;0 0 8;4 0 6] 
%返回 1 3 4 8 9，因为看成了一个列向量
ind = find(X)

% 如果需要按照行列的顺序输出
[r,c] = find(X)
[r,c] = find(X,1)%只找第一个非0元素

4.8 矩阵与常数的大小判断运算

X = [1 -3 0;0 0 8;4 0 6]
X > 0
X == 4