数学建模Matlab之基础操作

作者由于后续课程也要学习Matlab,并且之前也进行了一些数学建模的练习(虽然是论文手),所以花了几天零碎时间学习Matlab的基础操作,特此整理。


基本运算

a=5+5 %加法,同理减法
b=2^3  %立方
c=5*2 %乘法
x = 1; %分号用于不显示,这样在命令行窗口就没有任何输出
y = x + 5
eps %浮点数的相对误差
1i^2 %虚数
pi 

1.原来;不是可有可无的,只要写上;就代表在命令行窗口中不产生输出

2.虚数也可以是ly


变量和向量

%% 变量
% 输入,如果输入的不是正确格式,则继续接收,直到格式正确
x=input('请输入数字:');
y=x*2% 数据显示格式
% 修改数据显示格式
format long ,pi % 可以修改为小数点后16位
format  short ,pi % 默认short formal,四个小数位
% 定义时指定数据显示格式
format short %直接声明为短型
a=3.1415926a=5
b=2
a>b %大于
a==b %等于
a<b %小于
% 如果为真,则返回1,否则返回0
%% 向量
% 冒号创建,开始:步长:结尾
x=2:2:10% linspace函数创建线性间隔变量,(开始,结束,元素个数)
y=linspace(1,20,10)a=[1 2 3 4 5]; % 一行矩阵,可以看成是向量的一部分
b=[6 7 8 9 10];
c=a.*b % 加上. 就表示元素之间的运算
sum(c) % 将c向量的所有元素求和d=dot(a,b) % 使用dot函数对两个向量进行点积,就是x1y1+x2y2……+xnyn这个和的结果a=[1 2 3 ];
b=[6 7 8 ];
c=cross(a,b) % 三维向量的交叉积,右手法则从a到b就是c的方向pause %用于暂停指令,下面的指令暂时不会执行

 其中,y=linspace()是一个较为常见的创建线性变量的方法。


矩阵

基础的不讲,首先是矩阵。

%% 矩阵
A=[1 2 3; 4 5 2; 3 2 7]
A=[1 2 3 4 5 2 3 2 7]
B=A' %把A进行转置,行变列,列变行
C=A(:) %从左到右,从上到下,竖向展示矩阵的内容

多种矩阵定义中,作者感觉第二种较为常见,层次很分明。很少见到A(:)展示,A‘注意是转置

D=inv(A) %求矩阵的逆矩阵,注意,只有非奇异矩阵(行列式值不为0)才能求逆矩阵
A*D %A乘A的逆,就是单位矩阵E

注意:inv(A)是逆矩阵,非常常见。并且一定要注意的是,只有非奇异矩阵(即对应行列式不为0)才能有逆矩阵!那么,如何才能快速求得某一个方阵对应行列式值呢?方法如下:

det_A = det(A)
F=[1,2,3;4,5,6;7,8,9]
det_F = det(F) %在实际的数值计算中,由于浮点数的表示和计算误差,得到一个非常接近于0的行列式是可能的。

对于上文的A矩阵,可以得到det_A=-34,很明显不为0;对于新建的F矩阵,可以得到det_F为:

这……这已经很接近于0了,那它到底等不等于0啊?

原来,在实际的数值计算中,由于浮点数的表示和计算误差,得到一个非常接近于0的行列式是可能的。通常使用一些容差值来决定是否将计算出的行列式值视为0,本质上就相当于设置一段趋近于0的区间,只要行列式结果位于这个区间,就认为求得的行列式值为0,即对应方阵为奇异矩阵,没有逆矩阵。

tol = 1e-10;
if abs(det_F) < toldisp('F 是奇异矩阵');
elsedisp('F 不是奇异矩阵');
end

后续可以考虑封装为函数,通过这个流程,我们就可以对F矩阵是否能求逆矩阵进行判断:

然后是创建随机矩阵,作者目前这几天还没遇到过使用它的情况:

E=zeros(10,5,3) %创建一个10行5列3维的全0矩阵,原来3表示3维%%伪随机数(伪随机矩阵生成)
E(:,:,1) = rand(10,5) %rand(m,n)生成m行n列的均匀分布的伪随机数矩阵
%rand(m,n,‘double’)生成指定精度的均匀分布的伪随机数,参数还可以是’single’,默认情况下rand生成的全是整数randi(10) %开区间0~10,貌似只能生成整数E(:,:,2) = randi(5,10,5) %参数1界定最大范围。
E(:,:,2) = randi([1,6],10,5) %参数1也可以用两个数表示一个区间

另外补充一些矩阵的常见使用:

A = [1,2,3,4,5,6,5,4,6] % 没有;分割,就是一整行
B = 1:2:9 %第二个参数为步长,不可缺省,1 3 5 7 9 
B = 1:3:9 % 1 4 7
C = repmat(B,3,2) %以B矩阵为元素,生成3行2列的新的矩阵
D = ones(2,4) %生成一个2行4列的全1矩阵

其实通过B=1:9来获得行向量,还是较为常见的。作者在后续定义初始路径时,使用的就是这样的定义方法。zeros(2,3)或者ones(2,3)都是一种类型,使用固定的常量定义矩阵。

然后就是矩阵及其元素的四则运算,普通的*对应的就是矩阵的乘法,若使用.*或者./,那么涉及到的就是矩阵对应元素之间的运算了(别忘了矩阵相加减是矩阵元素相加减哦!):

% 矩阵及其元素的四则运算
A = [1 2 3 4; 5 6 7 8]
B = [1 1 2 2; 2 2 1 1]
C = A + B
D = A - B
E = A * B' %注意这个是‘转置
F = A .* B % .*表示对应元素相乘,不是矩阵的乘法
G = A / B % 相当于A*B的逆,逆矩阵是inv(B)
G = A \ B % 相当于B*A的逆
H = A ./ B % ./表示对应项相除

在实际操作过程中,矩阵转置对于初学者实在是不好区分,明明正确表述就是B^T,Matlab非要写成逆矩阵B'的形式……初学者们一定要注意哈!逆矩阵是inv()函数得到的,即inverse matrix。

% 矩阵的下标
A = magic(5)
% magic幻阵:数字是从1到n^2 的整数,并且每一行、每一列以及主对角线和副对角线上的数字之和都相等。
% 幻阵常用于密码学、设计与艺术、数学研究和解决某些优化和排列问题
B = A(2,3) %取矩阵二行三列元素值
C = A(3,:) % :为取全部,那么这条语句表示取第三行
D = A(:,4) %取第四列
[m,n] = find(A > 20) %找到大于20的序号值/矩阵
%取的是索引值,n是行,m是列

这里取元素值蛮重要,A(2,3)就是取矩阵A二行三列的元素,:表示全部,相当于占位符*。


元胞数组

Matlab的元胞数组(cell array)是一种特殊的数组类型,它可以容纳不同类型和大小的数据。与常规数组不同,每个元胞可以包含任何类型的数据。

{ }法

C = {'text', [1, 2, 3], 5}; %这里,C 是一个包含三个元胞的元胞数组。
% 第一个元胞包含一个字符串,第二个元胞包含一个向量,第三个元胞包含一个标量。C{1} % 返回 'text'
C{2} % 返回 [1, 2, 3]
C{3} % 返回 5
C{1} = 'new text'; % 修改第一个元胞的内容

通过{ }创建元胞数组,比较简单,推荐使用。调用元素时注意一定不是平常数组的()方法,而是同样需要{ }进行元素调用。可以对元胞数组进行简单的函数操作(简单看看就行):

% 使用 cellfun 函数来对每个元胞应用函数
result = cellfun(@length, C); % 返回每个元胞的长度:8,3,1% 可以使用 cellfun 和匿名函数来计算每个向量的和:
C = {[1, 2, 3], [4, 5], [6, 7, 8, 9]};
sums = cellfun(@(x) sum(x), C); %sums 是一个数组,其中包含 C 中每个向量的和。%使用循环来操作元胞数组的内容,显示 C 的每个元胞的内容。
for i = 1:length(C)disp(C{i});
end

cell( )法

%% 元胞数组2
A = cell(1,6) %这行代码创建了一个1行6列的空元胞数组A。
A{2} = eye(3) %将一个3x3的单位矩阵赋值给A的第2个元素。
A{5} = magic(5) %将一个5x5的幻方矩阵赋值给A的第5个元素。
B = A{5}

通过cell函数创建的元胞数组,如果数组中的元素比较少,建议使用{ },直接初始化。但是如果元素比较多,那作者认为使用cell先创建空元胞数组应该是更好地办法吧!


结构体

在MATLAB中,结构体(或称为“结构”)是一种数据类型,它允许将多个不同类型的数据组合在一起。结构体中的每个数据项都被称为字段,每个字段都有一个与之相关的字段名。

作者爱用的定义

student.name = '张三';
student.age = 20;
student.grade = '三年级';

就这么简单!定义都不定义,直接加.表示前者就是结构体,后者就是字段。运行结果如下:

对于多个同一类型的结构体(比如不同的班级),可以进行如下定义:

class(1).name = '张三';
class(1).id = '001';
class(1).scores = [95, 88, 76];class(2).name = '李四';
class(2).id = '002';
class(2).scores = [89, 92, 85];

哇去,简直是清晰明了啊!


长得像C语言的定义

books = struct('name',{{'Machine Learning','Data Mining'}},'price',[30,40]) %包含两个字段:name 和 price。
% 为结构体的某个字段分配一个元胞数组,你需要使用双大括号 {{}}。books.name %访问结构体的字段
books.name(1)
books.name{1}

这就有的讲了,上述代码设置books为结构体,其中有两个字段:name 和 price。以一个字段一个值这个顺序定义的,是不是不难理解?

读者可能会发现,这name后面,怎么跟了{{ }}?难道我在学vue在这传递变量了?其实,在结构体定义中,如果要为字段赋元胞数组类型的值,那么这个元胞数组就不能是原本的{ },必须要嵌套两层{{ }}才能表示这是个元胞数组,不信,你看上述代码执行的结果:

name(1)和name{1}有什么不同呢?不是说元胞数组调用元素使用的是{ }吗?不妨我们运行一下:

糟了,要长脑子了:

  1. books.name(1) 返回的是name字段中的第一个元素的单元格,而不是单元格中的实际内容。因此,返回的是一个包含'Machine Learning'1x1元胞数组

  2. books.name{1} 通过使用花括号 {} 来索引元胞数组,可以直接访问单元格中的实际内容。因此,这将返回字符串 'Machine Learning'


扩展

既然结构体中赋值元胞数组需要{{ }},那么作者想,使用{ }给name赋值会出现什么情况呢?

books_single = struct('name',{'Machine Learning','Data Mining'},'price',[30,40])
books_single(1).name
books_single(1).price
books_single(2).name
books_single(2).price

结果:

欸,变成1*2的结构体数组了!这是怎么回事呢?其实,这是matlab的规定,如果是单个{},则表示将{}中的所有元素分别分配给前面的字段name。在上面的例子中,{}的两个值分别赋给name,这就使得struct有两个name,导致整个的struct变成1*2的结构体数组。

可以通过book(1)这样调用第一个结构体,第一个结构体中name被分配第一个值,第二个结构体中name被分配第二个值,此时price的值在两个结构体中是相同的。

如果同时给price分配两个数(和name一样用{}分配),那结果和作者想的不同,作者以为是笛卡尔积一共是1*4结构体数组,没想到还是1*2:

books_single2 = struct('name',{'Machine Learning','Data Mining'},'price',{[30,40],[50,60]})
books_single2(1).name
books_single2(1).price
books_single2(2).name
books_single2(2).price

所以作者认为,单独的{}本质上就是为了方便创建结构体数组的,name参数1和price参数1一起,然后按照顺序,name参数2和price参数2一起构成下一个结构体。那如果再给name增加一个数会发生什么呢?

books_single3 = struct('name',{'Machine Learning','Data Mining','Matlab'},'price',{[30,40],[50,60]})

结果报错啦!name字段是一个长度为3的字符串数组,而price字段是一个长度为2的数组的数组。这两个字段的长度不匹配,导致了错误。只能让两者长度保持一致,或者其中之一为标量值,没有{},如上文的books_single。

% 保证结构体字段长度相同
books_single3 = struct('name',{'Machine Learning','Data Mining','Matlab'},'price',{30,40,50})
books_single3 = struct('name',{'Machine Learning','Data Mining','Matlab'},'price',{{[30,40]},{[50,60]},{[70,80]}})

与元胞数组的比较

在数学建模中的应用:

  1. 结构体:
    • 用于存储模型的参数和配置信息。
    • 例子: model.parameters, model.settings
  2. 元胞数组:
    • 用于存储不同类型或大小的数据集,例如多个数据集或多个模型。
    • 例子: 不同类型的输入数据或模型。 
  • 如果你需要存储具有命名字段的相关数据结构体可能更合适。
  • 如果你只是需要存储不同类型或大小的数据,并不关心字段名称,那么元胞数组可能是更好的选择。
  • 在数学建模中,这两种数据类型都很有用,具体的选择取决于你的具体需求和使用场景。

流程结构

%% 循环结构
sum=0;
for i=1:5sum=sum+1;
end
% 间隔也可以为负数
for a = 1.0: -0.1: 0.0disp(a)
end
% 也可以不是常见的循环结构,可以是数组,但是注意,如果之前已经赋过a的话,也需要这么写
for a = [2 3 4 5 6]disp(a)
end
% 如果a本来就有值,也不能写成 for a ... disp(a) ...end的格式,matlab会检测表达式无效% while循环
i=0;
sum=0;
while(i<10)sum=sum+i;i=i+1;
end
fprintf('累加和为: %d \n', sum);%% 分支结构
if sum==3'成立'
else'不成立'
endswitch sum>0case 1disp('sum大于零');otherwisedisp('sum小于等于0');
end

到时候看看具体算法,看多了就熟悉了。


基本绘图操作

代码来自于CSDN最火的两个matlab教程。

波形函数绘图

%% 基本绘图操作
%1.二维平面绘图
x = 0:0.01:2*pi %定义x的范围,第二个参数表示步长
y = sin(x)
figure %建立一个幕布
plot(x,y) %绘制当前二维平面图
title('y = sin(x)') %标题
xlabel('x') %x轴
ylabel('sin(x)') %y轴
xlim([0 2*pi]) % x坐标值的范围

x = 0:0.01:20;
y1 = 200*exp(-0.05*x).*sin(x);
y2 = 0.8*exp(-0.5*x).*sin(10*x);
figure
[AX,H1,H2] = plotyy(x,y1,x,y2,'plot'); %共用一个x的坐标系,在y上有不同的取值
%设置相应的标签
set(get(AX(1),'Ylabel'),'String','Slow Decay')
set(get(AX(2),'Ylabel'),'String','Fast Decay')
xlabel('Time(\musec)')
title('Multiple Decay Rates')
set(H1,'LineStyle','--')
set(H2,'LineStyle',':')

双峰函数绘图 

%2.三维立体绘图
t = 0: pi/50: 10*pi;
plot3(sin(t),cos(t),t)
xlabel('sin(t)')
ylabel('cos(t)')
zlabel('t')
%hold on
%hold off %不保留当前操作
grid on %把图片绘制出来,在图片中加一些网格线
axis square %使整个图(连同坐标系)呈方体[x,y,z] = peaks(30); %peaks命令用于产生双峰函数或者是用双峰函数绘图
mesh(x,y,z)
grid

散点图

%% 数据可视化——散点图
% x表示年龄,y1表示人体血压。绘制年龄与血压的关系,使用scatter散点图函数
x=[75;78;51;82;77;88;41;78;78;61;71;74;62;81;75;64;80;72;51;80;56;73];
y1=[208;146;168;149;208;102;130;155;163;154;145;147;143;161;145;120;153;158;123;163;177;148];scatter(x,y1)
scatter(x,y1,'r')  %填充颜色绘制为红色% 利用向量配置多彩颜色,每个点的颜色由c确定。
c = linspace(1,10,length(x)); %生成了一个线性间隔的向量c。c的长度与x相同,值从1到10。
scatter(x,y1,[],c) %这些值作为颜色数据传递给scatter函数,使散点具有从一种颜色渐变到另一种颜色的效果。% 填充并标记散点彩色图
sz = 30; % 设置散点的大小为30。
c = linspace(1,10,length(x)); %同样利用向量配置多彩颜色
scatter(x,y1,sz,c,'filled') %填充散点
多彩填充散点图

条形图

%% 数据可视化——条形图
x = {'计算机系统基础','汇编语言','信息安全','机器学习'}; %x个数一定要与y个数对应
y = [94,54,65,87];
%在调用 bar 函数时,MATLAB 不支持直接使用字符串数组作为 x 坐标轴标签。您可以使用类别数组 categorical 来解决这个问题。
x = categorical(x); % 将x转换为类别数组
% bar(x,y)
b = bar(categorical(x), y);
b.FaceColor = [0.5 0.7 1]; % 设置为淡蓝色 (RGB 色值)
xlabel('科目') %经过我的深思熟虑,我认为还是要保证整个图像的完整性,不要仅仅依赖于图题
ylabel('成绩')
title('成绩分布图')


函数

简单说明一下,在MATLAB中,不能在命令窗口中直接定义函数。函数必须保存在一个扩展名为.m的文件中,并且文件名必须与函数名相同。这文件的创建,左上角新建--函数,就会有对应函数模板的文件出现啦:

格式为function 输出形参表 = 函数名(输入形参表)。对于函数,要有以下保证:

1.确保文件已经正确保存,文件名与函数名匹配,且扩展名为.m(例如f.m)。
2.确保函数文件位于当前的MATLAB工作路径中,或者您在调用函数时提供了完整的文件路径。

给个例子:

% function 输出形参表 = 函数名(输入形参表)
function y = funcdemo1(x)
%UNTITLED 此处显示有关此函数的摘要
%   此处显示详细说明
%   x为一个整数
if x>0y=2*x;
elseif x==0y=0;
elsey=x*x;
end
disp(y)
end

到时候调用啊,直接调用funcdemo1(x)就行了:

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/125965.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

八、【快速选择工具组】

文章目录 对象选择工具快速选择工具魔棒工具 对象选择工具 当我们选择对象选择工具时&#xff0c;需要先注意上边有一个循环的圆&#xff0c;它会进行内容识别&#xff0c;当识别完成会停止旋转。这个时候我们按住n键&#xff0c;或者将鼠标放上对应的图形时会出现选中的颜色。…

windows 远程连接 ubuntu桌面xrdp

更新 sudo apt update安装组件 sudo apt-get install xorg sudo apt-get install xserver-xorg-core sudo apt-get install xorgxrdp sudo apt install xfce4 xfce4-goodies xorg dbus-x11 x11-xserver-utilsxrdp sudo apt install xrdp sudo systemctl status xrdp sudo …

Linux CentOS7 vim宏操作

vim的macro就是用来解决重复的问题。在vim寄存器的文章里面已经对macro有所涉及&#xff0c;macro的操作都是以文本的方式存放在寄存器中。 宏是一组命令的集合&#xff0c;应用极其广泛&#xff0c;包括MS Office中的word编辑器&#xff0c;excel编辑器和各种文本编辑器&…

使用Thrift实现跨语言RPC调用

&#x1f4cb; 个人简介 &#x1f496; 作者简介&#xff1a;大家好&#xff0c;我是阿牛&#xff0c;全栈领域优质创作者。&#x1f61c;&#x1f4dd; 个人主页&#xff1a;馆主阿牛&#x1f525;&#x1f389; 支持我&#xff1a;点赞&#x1f44d;收藏⭐️留言&#x1f4d…

基于SpringBoot的智能物流管理系统

目录 前言 一、技术栈 二、系统功能介绍 顾客信息管理 员工信息管理 员工信息管理 门店信息管理 门店信息管理 订单信息管理 三、核心代码 1、登录模块 2、文件上传模块 3、代码封装 前言 随着信息技术在管理上越来越深入而广泛的应用&#xff0c;管理信息系统的实施…

re学习(38)HGAME2020-re-Level-Week1-maze

题目描述 You won’t figure out anything if you give in to fear. 学习资料: https://ctf-wiki.github.io/ctf-wiki/reverse/maze/maze-zh/ 附加说明&#xff1a;请走最短路线 题解 分析题目 一看题目&#xff1a;maze 可以确定是一个迷宫题 void __fastcall __noreturn…

wisemodel 始智AI - 小记

文章目录 关于 wisemodel 始智AI 关于 wisemodel 始智AI https://www.wisemodel.cn/home 旨在打造中国版 “HuggingFace” 该社区汇聚了清华 / 智谱 chatglm2-6B、Stable Diffusion V1.5、alphafold2、seamless m4t large 等模型&#xff0c;以及 shareGPT、ultrachat、moss-…

Pytorch之shuffleNet图像分类

&#x1f482; 个人主页:风间琉璃&#x1f91f; 版权: 本文由【风间琉璃】原创、在CSDN首发、需要转载请联系博主&#x1f4ac; 如果文章对你有帮助、欢迎关注、点赞、收藏(一键三连)和订阅专栏哦 前言 ShuffleNet是Face&#xff08;旷视&#xff09;在2017年发布的一个高效率…

黑马JVM总结(二十五)

&#xff08;1&#xff09;字节码指令-cinit 构造方法可以分为两类&#xff0c;一类是cinit 一类init cinit是整个类的构造方法 putstatic&#xff1a;进行static变量的赋值&#xff0c;是到常量池里找到名字一个叫做i的变量 &#xff08;2&#xff09;字节码指令-init in…

Waves 14混音特效插件合集mac/win

Waves14是一款音频处理软件&#xff0c;主要用于音频编辑、混音和母带处理。该软件提供了各种插件&#xff0c;包括EQ、压缩、混响、延迟、失真等&#xff0c;以及一些专业的音频处理工具&#xff0c;如L2限幅器、Linear Phase EQ和多频道扬声器管理。 Mac软件下载&#xff1a;…

信创之国产浪潮电脑+统信UOS操作系统体验2:安装visual studio code和cmake搭建C++开发环镜

一、引言 在linux以及unix下开发C程序很多人都做过&#xff0c;笔者以前是通过vi编辑器编辑或者在windows下编辑甚至编译再将源文件ftp保存到服务器上&#xff0c;再在服务器端使用gcc编译&#xff0c;用gdb调试&#xff0c;但现在使用的国产电脑&#xff0c;客户端就是linux&…

迄今为止丨ChatGPT最强指令,一个可以让机器人生成机器人的Prompt,价值百万!

原文&#xff1a; 【ChatGPT调教】ChatGPT最强指令、让机器人为你生成机器人&#xff01;-CSDN博客 说明&#xff1a;最好看原文 昨天&#xff0c;发现了一条可能是迄今为止&#xff0c;我见过最牛的&#xff0c;商业价值最高的ChatGPT指令。 通过这条指令&#xff0c;可以…