Matlab基础语法篇（下）

最新推荐文章于 2024-07-19 22:16:04 发布

wanchen_Gabby

最新推荐文章于 2024-07-19 22:16:04 发布

阅读量522

点赞数 10

分类专栏： MATLAB 数学建模文章标签： matlab android 开发语言

本文链接：https://blog.csdn.net/wanchen_Gabby/article/details/140295208

版权

数学建模同时被 2 个专栏收录

6 篇文章 0 订阅

订阅专栏

MATLAB

3 篇文章 0 订阅

订阅专栏

一、逻辑基础

（一）逻辑运算符

在这里插入图片描述
MATLAB 推荐大家直接使用运算符进行计算，因此&、和~这三个符号的功能大家要牢记。
以“&”为例:

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

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)

在这里插入图片描述
另：

支持连续运算

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

在这里插入图片描述

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

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

在这里插入图片描述

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

（二）all、any、find函数

在这里插入图片描述

all(A,dim) 与 & 类似

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

在这里插入图片描述

any(A,dim) 与 | 类似

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

在这里插入图片描述

find(X,n,direction)

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个同学为例

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

在这里插入图片描述

这里也以5个同学为例

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

在这里插入图片描述

二、结构基础

（一）条件结构

（1）if-elseif-else-end

在这里插入图片描述

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

res

在这里插入图片描述

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

res

在这里插入图片描述

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

在这里插入图片描述

% 找最大值
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

在这里插入图片描述

% 注：开关表达式的计算结果必须是一个 数值标量或者是一个字符向量/字符串，不能是向量或矩阵
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

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

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

在这里插入图片描述

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

res_sum

在这里插入图片描述

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

在这里插入图片描述

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来中断程序的运行。

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

在这里插入图片描述

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

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

在这里插入图片描述

（3）break、continue

在这里插入图片描述

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

在这里插入图片描述

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

在这里插入图片描述

（三）小结 + 练习

在这里插入图片描述

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可以省略）
函数的调用

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

% 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）匿名函数

在这里插入图片描述

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) % 输入两个变量，返回一个单变量函数

在这里插入图片描述

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）数值运算、三角函数

在这里插入图片描述

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

exp(1:9)

在这里插入图片描述

（2）isempty函数

在这里插入图片描述

% 法一：
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函数

在这里插入图片描述

% 基于向量 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函数

在这里插入图片描述

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

在这里插入图片描述

完

wanchen_Gabby

关注

10
点赞
踩
9

收藏

觉得还不错? 一键收藏
0
评论
Matlab基础语法篇（下）

包含MATLAB的逻辑基础（逻辑运算符、all、any、find函数）、结构基础（条件结构语句、循环结构语句）、函数基础（两种自定义函数、特殊函数、常用函数——数值类、isempty、meshgrid、rng）的详细介绍
复制链接

扫一扫