Matlab基础语法篇（下）

Matlab基础语法（下）

一、逻辑基础
二、结构基础
- （一）条件结构
- - （1）if-elseif-else-end
  - （2）switch-case-otherwise-end
- （二）循环结构
- [（三）小结 + 练习](#（三）小结 + 练习)
三、函数基础

一、逻辑基础

（一）逻辑运算符

MATLAB 推荐大家直接使用运算符进行计算，因此&、和~这三个符号的功能大家要牢记。

以"&"为例:

针对矩阵运算时:"逻辑与&"是对A和B进行计算的，计算时会比较A和 B对应位置的元素。需要符合算术运算中介绍的五种兼容式即可。
"逻辑与&"不仅可以作用在逻辑值0和1上，还可以用于普通的数值上，这时候，MATLAB 会将非零数值视为逻辑1，将数值零视为逻辑0进行运算。
"逻辑与&"也可以进行连续运算，例如1&2&3。
为避免运算优先级的问题，直接使用小括号来指定计算顺序。

matlab 复制代码

A = randi([-3,3],2,4)
B = randi([-3,3],1,4)
A & B
A | B
~A
xor(3,4) % 亦或 相当于 xor(1,1)
xor(A,B)

另：

支持连续运算

matlab 复制代码

A = randi([-3,3],2,4)
B = randi([-3,3],1,4)
C = randi([-3,3],1,4)
D = A & B & C

关系运算符的优先级高于逻辑运算符

matlab 复制代码

3>4 & 2>-1 % 相当于 0 & 1

辨析：&&和|| vs &和|
(1)&&和|| 只能对标量（只有一个元素） 进行逻辑运算，不能对有多个元素的向量或者矩阵 进行运算，而&和|可以
(2)&&和||有短路功能，可以提高运行效率

（二）all、any、find函数

all(A,dim) 与 & 类似

matlab 复制代码

B = randi([1,100],2,5)
B(3) = 0 % 线性索引改值
all(B,1)
all(B,2)

any(A,dim) 与 | 类似

matlab 复制代码

B = randi([1,100],2,5)
B(3) = 0 % 线性索引改值
any(B,1)
any(B,2)

find(X,n,direction)

matlab 复制代码

A = randi([0,2],2,3)
ind1 = find(A) %返回其中所有 非0元素 的 线性索引
ind2 = find(A,2) %返回其中 前两个非0元素 的 线性索引
ind3 = find(A,2,"last") %返回其中 后两个非0元素 的 线性索引
[row,col] = find(A) %按行列分开 生成 两个向量
[row,col,v] = find(A) % 多生成了一个A中非0元素的列向量

（三）练习

这里以5个同学为例

matlab 复制代码

score = randi([50,100],5,3)
% 先将分数矩阵转化为逻辑矩阵
any(score<60,2) % 关注同学 --> 生成列向量
all(score>=60,1) % 关注科目 --> 生成行向量

这里也以5个同学为例

matlab 复制代码

score = randi([50,100],5,3)
find(sum(score < 60,2) == 2) 
find(sum(score,2) > 260)

二、结构基础

（一）条件结构

（1）if-elseif-else-end

matlab 复制代码

if [1,2;0,2] % 矩阵中的元素都非 0，表达式才会被当所逻辑值 1
    res = 0;
else
    res = 10;
end

res

matlab 复制代码

% 若想要至少有一个元素非 0 就执行 if 语句
A = [1,2;0,2]
A(:)
% A(:) 先用线性索引把 A 拉成一个列向量，不能直接 any(A)，这样生成的是向量
if any(A(:)) 
    res = 0;
else
    res = 10;
end

res

matlab 复制代码

score = 87;
if (score >= 90) && (score <= 100)
    res = 1;
elseif (score >= 80) && (score < 90)  % elseif之间没有空格
    res = 2;
elseif (score >= 60) && (score < 80)
    res = 3;
elseif (score >= 0) && (score < 60)
    res = 4;
else
    res = 0;
end % 不要忘了写！

res

matlab 复制代码

% 找最大值
a = 10;
b = 20;
c = 15;
if a > b
    if a > c
        max = a;
    else
        max = c;
    end
else
    if b > c
        max = b;
    else
        max = c;
    end
end

max

（2）switch-case-otherwise-end

matlab 复制代码

% 注：开关表达式的计算结果必须是一个 数值标量或者是一个字符向量/字符串，不能是向量或矩阵
season = randi([1,4])
switch season
    case 1
        disp("春季");
    case 2
        disp("夏季");    
    case 3
        disp("秋季");
    otherwise
        disp("冬季"); % 加不加分号不影响输出结果
end
% 或
if season == 1
    disp("春季");
elseif season == 2
    disp("夏季");
elseif season == 3
    disp("秋季");
else
    disp("冬季");
end

（二）循环结构

（1）for-end

循环变量 ：用于迭代的变量名，它会在每次循环迭代中从向量或矩阵中取出列的值。数值向量或者矩阵则表示了循环变量可以取值的范围，一旦循环变量遍历完数值向量或者矩阵中的所有值，循环就会结束

matlab 复制代码

% for后跟的可以是一个向量，也可以是一个矩阵
A = randi([-3,3],2,3)
for i = A % 每次都会取出一 列 的值
    i
end

matlab 复制代码

x = 1:6
res_sum = 0;
for i = x
    res_sum = res_sum + i;
end

res_sum

matlab 复制代码

leap_year_num = 0;
for i = 1 : 9999
    if ((mod(i,4) == 0) && (mod(i,100) ~= 0)) || (mod(i,400) == 0)
        leap_year_num = leap_year_num + 1;
    end
end

leap_year_num

（2）while-end

matlab 复制代码

f(1) = 1;
f(2) = 1;
n = 2;
while f(n) <= 99999
    n = n + 1;
    f(n) = f(n-1) + f(n-2);
end
n
f(n)

一个无限循环，在命令行窗口中按下快捷键 Ctrl + C来中断程序的运行。

matlab 复制代码

while 1
    disp("死循环");
end

while 后面表达式的计算结果不一定非得是逻辑值1或0。如果表达式的计算结果是一个数值常数，则只有当这个常数为非零值时循环才会进行

matlab 复制代码

i = 5;
while i
    i;
    i = i - 1;
end
i

（3）break、continue

matlab 复制代码

for i = 1:10
    if mod(i,2) == 0
        continue;
    end
    i
end

matlab 复制代码

n = 9;
is_prime = 1;
for i = 2:n-1
    if mod(n,i) == 0
        is_prime = 0;
        break;
    end
end
is_prime

（三）小结 + 练习

matlab 复制代码

a = 6;
b = 10;
eps = 1e-8; %误差阈值
while 1
    c = (a+b)/2;
    fc = c^3 - 8*c^2 +c - 5;
    if abs(fc)<eps
        break;
    end
    fa = a^3 - 8*a^2 +a - 5;
    if fa * fc < 0
        b = c;
    else
        a = c;
    end
end

x0 = c;
x0

三、函数基础

（一）自定义函数

（1）m文件定义函数

无需关注具体的实现过程，需要编写额外的.m文件 ，定义完的函数保存为同函数名的m文件，放置在Matlab当前路径下

函数的基本结构（end可以省略）
函数的调用

matlab 复制代码

% code1.m 文件中
X = 1:20;
[max,min] = max_min__values(X)

matlab 复制代码

% max_min__values.m文件中
function[max,min] = max_min_values(X)
max = subfun1(X);
min = subfun2(X);

    function r = subfun1(X)
        x1 = sort(X,'descend');
        r = x1(1);
    end

    function r = subfun2(X)
        x1 = sort(X);
        r = x1(1);
    end
end

（2）匿名函数

matlab 复制代码

f1 = @(x,y)x.^2 + y.^2;
f1(2,3)

x = 1:5;
y = 0.1:0.1:0.5; %从 0.1，间隔 0.1，到 0.5
f1(x,y)

%双重匿名函数
f2 = @(a,b)@(x)a*x+b; % 外层变量：a,b  内层变量：x
f2(2,3) % 输入两个变量，返回一个单变量函数

matlab 复制代码

f = @(a)@(x)exp(x)+x^a+x^(sqrt(x))-100; % 根据 a 的函数，获得 x 的函数
A = 0:0.1:2;
% fzero(fun,x0) x0:预测初始值
x = @(a)fzero(f(a),4); % 将 x 变为 a 的函数
X = @(A)arrayfun(@(a)x(a),A);  % 将函数应用于每个数组元素
Y = X(A)
% 相当于 Solution = @(A)(@(a)fzero(f(a),4),A)

（二）特殊函数

贝塞尔函数、beta函数、gamma函数......

有需求时去MATLAB帮助中心查看

（三）常用函数

（1）数值运算、三角函数

matlab 复制代码

sqrt(1:9)
format long g % 取消掉科学计数法
sqrt(1:9)
sqrt(-4)
format short

exp(1:9)

（2）isempty函数

matlab 复制代码

% 法一：
A = [1:9];
x = 10;
% A == x 元素全为 0 的逻辑向量
~isempty(find(A == x))
% 法二：
A = [1:9];
x = 10;
ismember(x,A)
% 推广：
A = [1:9];
x = 5:10;
ismember(x,A) %找 x 是否在 A 中存在

（3）meshgrid函数

matlab 复制代码

% 基于向量 x 和 y 中包含的坐标来 返回 二维网络坐标
x = 0:4;
y = 0:3;
[X,Y] = meshgrid(x,y)
% 若写成[X,Y] = meshgrid(x)，相当于 [X,Y] = meshgrid(x,x)
Z = X.^2 + Y.^2

（4）rng函数

matlab 复制代码

% 设置随机数种子，这样能生成可重复的随机数
seed = 3;
rng(seed);
randi(10,3,3)
rng('shuffle'); %根据当前的时间来设置随机数种子

完