【MATLAB】常见语法结构和实用操作
创始人
2024-02-23 11:35:43
0

目录

语法结构

循环(遍历)结构

函数结构

画图操作

绘制函数图像

ezplot

figure多图绘制

创建多项式 

按照系数创建多项式

 利用多项式零点创建多项式

解方程组


语法结构

循环(遍历)结构

与其它编程语言类似,Matlab也有循环或者说是遍历的语法。下面以输出1到10为例:

for i = 1:10disp(i);
end

可以看出它的结构主要有三个部分:

第一部分是关键字for和结束标志end:

第二部分是循环初始值和结束位置,这里注意末尾的10也会进入循环,类似于C语言中的for循环的边界 i<=n ,当 i 为 n 是也会进入循环:

第三部分就是循环内容,也就是循环内部要实现什么目标了,这里我就单纯输出1到10所以数字即可:

函数结构

在我们看matlab中的函数之前,首先看一下C语言中的函数结构:

void func(int x, int y)
{return;
}

我们很明显可以看出C语言中函数结构主要包含下面几个部分:函数名、返回值、参数列表、函数内容。

之后再让我们看看matlab中的函数,matlab中的函数也由相似的几部分构成:

第一部分是function函数关键字和返回返回值关键字return:

 第二部分是一个类似函数调用的部分,主要由返回值c、函数名func(注意这里的函数名也是函数文件名)以及参数a构成:

 最后一部分激素函数内容部分了:

画图操作

绘制函数图像

matlab中画图的方式有很多种,在看如何画图之前,我们先来看几个基本的指令:hold onhold off grid on。其中grid on比较好理解,就是在我们画的图上面多了一些网格;而hold on 和hold off有点不太同,hold on可以让多个图并存,即如果此后我们更新图像那么原本的图也会发生改变,即牵一发动全身;hold off则不同,它的每一个图都是独立的,并且相互覆盖

grid on

syms y x;
ezplot(y==x,[-10,10]);
grid on;

hold on 和hold off的区别

hold on:

syms y x;
ezplot(y==x^2,[-6,6]);
hold on;
ezplot(y==x,[-10,10]);

 hold off:

syms y x;
ezplot(y==x^2,[-6,6]);
hold off;
ezplot(p==0,[-10,10]);
hold on
hold off

我们可以看到明显区别hold on的图保持了原来的两幅图,而hold off的图则覆盖掉了原来的第一副图。了解了上面的这些命令后我们先来看看其中比较常见的几个函数:

ezplot

这个函数是比较常见的绘图函数,它包括几种常用的用法、形式:

syms x a y;
y = sin(x);
ezplot('sin(x)');
ezplot(y);

其中包括字符串形式和关于函数的形式,其中ezplot(y),y=sin(x)就是关于函数的形式,y是一个x的函数;还有一种就是字符串形式,例如'sin(x)',其中单引号里面包括的就是一个字符串。 

但是往往并不是所以的函数都是显函数,有时候我们也会拿到一些隐函数,例如:y^{3}-6y+\sqrt{5-x}=0,此时我们就需要利用一个函数f来表示,然后令函数f为0得到的方程画出来的图就是隐函数的图形。

syms y x;
f = y^3-6*y-sqrt(5-x);
ezplot(f==0,[-6,6]);

figure多图绘制

我们知道,ezplot绘制的是一个图形,那么如果我们想要绘制一个框里面多个图像应该怎么做呢?这时候就需要用到figure的多图绘制。首先我们要绘制一个图,需要将x细分为一堆等长区间:

x = 0:0.0001:2; 

这里第一个0表示x的起始值,第三个2表示x的末尾值,而中间的0.001表示将x细分的步长。所以我们这里定义出来的x是一个矩阵,即如果在定义x为自变量的函数时,应该注意点了。再准备一个函数y,令y = x^{2}-6

y = x.*x - 6;

然后我们就可以开始准备画图了,首先准备一个figure框,也就是我们画图常用的哪个框框,并且使用hold on指令保持住这个框,方便接下来画图。然后操作分框函数subplot,这个函数用于将当前的figure框分为几个部分,也就是我们画几个图,它参数的三个数字,前面表示划分的行数,第二个表示划分的列数,第三个数字表示当前图像是第几个框(横向数)。以6个框为例:

subplot(321);
subplot(322);
subplot(323);
subplot(324);
subplot(325);
subplot(326);

然后就是对图像进行绘制(填充)了,我们使用fill函数,第一个参数为自变量,第二个参数为因变量,第三个参数为图像的颜色(英文字符或英文字符缩写均可)。最后的title函数里面跟的就是图像名,这样一个函数图像就画出来了!

x = 0:0.0001:2; 
y = x.*x - 6;
figure;
hold on;
subplot(223);
fill(x,y,'green');
title('y-x');

创建多项式 

下面我们看看在matlab中如何创建多项式,这是我看了其他几个很棒的博主的文章后学习到的,链接我放在了文章结尾。首先是创建多项式的思路,分为两种:按照系数来创建多项式;按照根来创建多项式,下面我们一个个看:

按照系数创建多项式

s = [9,8,7,6,5,4,3];
y1 = poly2sym(s);
disp(y1);

这里使用的是函数poly2sym来创建多项式,我们要创建一个多项式,需要向这个函数传递多项式各项的系数。以数组s为例,我们数组里有7个数字,将数组传递进去后我们得到的多项式的系数从小到大分别为数组最大位到数组最小位,例如上面的s数组创建出来的多项式为:

9x^{6}+8x^{5}+7x^{4}+6x^{3}+5x^{2}+4x^{1}+3

其实我们还可以利用一个与ploy2sym类似的函数ploy2str来创建多项式:

s = [9,8,7,6,5,4,3];
y2 = poly2str(s,'x');
disp(y2)

 不过我们可以利用whos看到这两个方法还是有区别的,它们创建出来的多项式的数据类型是不同的:

 利用多项式零点创建多项式

r = [0,1,3];
y = poly2sym(poly(r));
disp(y);

除了刚才利用系数来创建多项式的方法,我们还可以利用多项式的零点来创建多项式,例如我们得到多项式一组零点r=[0,1,3]。那么我们也可以先使用ploy函数处理零点,再利用poly2sym函数来创建多项式。

解方程组

前面说了一些语法的结构,下面我们来看一个更加实际的作用:求解线性方程组的解。要知道,之前我们求简单方程的解使用的是solve这个函数,但是我们有时候往往需要去求一些方程组,这时候就需要寻找别的方法了。

回顾一下线性代数,我们求解方程组还有一个方法就是使用矩阵求解,例如方程组:

\left\{\begin{matrix} 3x+5y-6z=9 \\x-8y+7z=6 \\2x+6y+3z=3 \end{matrix}\right.

我们使用矩阵求解的大致步骤就是求出系数矩阵,然后将两边乘系数矩阵逆矩阵消去系数的解就是X:

\begin{pmatrix} 3& 5 &-6 \\ 1& -8 &7 \\ 2&6 &3 \end{pmatrix}X = \begin{pmatrix} 9\\6 \\3 \end{pmatrix} =>X = \begin{pmatrix} 3& 5 &-6 \\ 1& -8 &7 \\ 2&6 &3 \end{pmatrix}^{-1}\begin{pmatrix} 9\\6 \\3 \end{pmatrix}

所以我们使用matlab求解线性方程组的方法也类似,首先定义一个系数矩阵A,然后再定义一个矩阵B,X的结果就是A^{^-1}B。对应着,matlab表现A^{^-1}的方法是A后面跟一个反斜杠 \

A = [3 5 -6;1 -8 7;2  6 3];
B = [9;6;3];
X = A\B;
disp(X);

参考资料:

https://blog.csdn.net/qq_38431572/category_9417566.html

matlab多项式构造/求根/某处的值_执念斩长河的博客-CSDN博客_matlab多项式构造

matlab中 hold on 与 hold off,figure作用_江江江在奋斗的博客-CSDN博客_figure hold on

matlab中的ezplot函数详解_「已注销」的博客-CSDN博客_ezplot函数

相关内容

热门资讯

涉嫌内幕交易!千亿锂矿巨头被移... 锂矿龙头之一的赣锋锂业突发公告。 12月29日晚间,江西赣锋锂业集团股份有限公司(以下简称“赣锋锂业...
中央商场控股子公司泗阳雨润被江... 观点网讯:12月30日,南京中央商场(集团)股份有限公司发布公告,披露其控股子公司泗阳雨润中央购物广...
法治阳光照高墙:湖南芙蓉律所律... 为提升服刑人员的法治观念和回归社会的适应能力,帮助刑满释放人员顺利融入社会,近日,湖南芙蓉(常德)律...
强化执业纪律意识,筑牢合规执业... 为进一步强化律师执业纪律意识,持续提升律所规范化、专业化管理水平,2025年12月29日下午,湖南芙...
海口美兰区海府街道多部门联合化... “没想到这么快就能解决问题,现在上下楼方便多了,感谢你们为老百姓办实事!”12月29日,邮政西院小区...
重磅!国补最新政策落地 国家发展改革委 财政部关于2026年实施大规模设备更新和消费品以旧换新政策的通知 各省、自治区、直...
扫描“主播”|真人互动擦边影游... 视频编辑:张兆亿(01:24) “这能看?”“这种游戏是能播的吗?”“这几个女生衣服不能好好穿吗?...
明年“两新”政策部署来了! 每经记者|周逸斐 每经编辑|贾运可 作为扩内需的关键抓手,“两新”政策(大规模设备更新和消费品以旧...
提升政策温度,缩小经济温差 作为宏观记者,这一年笔者穿梭于各种新闻发布会、政策吹风会,解读从中央到各个部委发布的经济数据、政策文...
上海11处历史文化名镇:实行历... 新民晚报周馨摄影 赓续历史文脉,激发发展活力,实现经济与人文相互促进、共同繁荣。今天下午,市十六届人...