数学建模--MATLAB基本使用

1.线性方程组

这个是一个线性方程组(属于线性代数的范畴),Ax=b类型的方程,如果使用MATLAB进行求解,就需要分别表示A矩阵(线性方程组未知数前面的系数),b矩阵(表示等式右边的结果),inv是这个软件里面的一个函数,用来进行求解A的逆矩阵,因为Ax=b,那么x=A-1次方乘上b

另外,我们也可以使用这个A\b来进行求解,这个里面需要注意是A\b,而不是我们熟悉的A/b,他们所代表的含义是不一样的,我们在脚本里面敲代码的时候,需要在前面的2行加上结尾处的引号,最后一行不用加引号,否则我们在命令行窗口里面无法看到最后的输出结果;

2.二维图形的绘制

二维图形的绘制要是用这个plot函数,我们通过这个例子来认识一下函数的使用方法,以sin1/x为例吧:

(1)这个地方我们只绘制了-1~1范围里面的函数图像,这样便于观察函数的效果,学过高等数学的同学们对于这个函数就非常的熟悉,因为他是特例,这个函数是有界限的,但没有极限存在;

(2)下面我们来聊一聊这个里面每一行代码的作用,首先,第一行是生成行向量,我们可以简单地理解为,通过这一行,我们要告诉MATLAB我们要绘制图像的定义域是-1~1之间的,100000是在这个-1~1这个区间里面生成100000个采样点,

(3)这个里面使用了点运算符,只要是涉及到矩阵里面的元素运算,我们在大部分的情况下都要是用这个运算符,后续还会涉及到,这个运算符在矩阵的运算里面很常见,我们都知道,对于两个矩阵的乘法运算,他们并不是对应位置的元素相乘,但是如果我们使用点乘运算符,就可以让矩阵里面对应位置的元素进行相乘的运算;

**************************************************************************************************************

接下来我们举一个在一个坐标平面里面绘制2个图形的例子

我们在一个坐标平面里面绘制y=e^x*cos和y=10*e^-0.5x*sin(2pi x)的函数图像(0~pi/2区间里面)

这个里面的plot中的单引号里面表示的是一些标记符号;

b是蓝色,-表示实线;合在一起就是蓝色实线;

k表示黑色,:表示虚线;合在一起就是黑色虚线;

下面是效果图:

3.空间曲线的绘制

绘制空间曲线x^2+y^2+z^2=16,x+y=0(这2个方程是在一个括号里面进行联立的)

首先要转化成为参数方程,也就是x=2根号2sint;y=-2根号2cost;z=4cost;

下面是是否添加grid on的区别,读者可以自行感受:

这个是没有grid on的,可见就没有曲线的格子;

这个是添加了grid on的,差别就显而易见了;

4.函数极限实验

(1)符号运算,matlab使用sym创建符号对象,一次只能够创建一个符号对象,我们也可以使用syms一次性创建多个符号对象;

(2)极限运算

MATLAB提供limit函数求极限

limit(f,x,a)表示求f(x)在a点处的极限值;

limit(f,x,a,'right')是求这点处的右极限,如果我们把right换成left就是求左极限;

当然,也可以进行求解多元函数的极限:这个时候就需要嵌套limit;

z=limit(limit(f,x,x0),y,y0)就是求在x0,y0处的极限值;

f1里面就是我们要求的函数的表达式,在x趋近于1的极限值;

我们都知道sin(1/x)的极限是不存在的,所以程序的输出结果就是NaN;

5.MATLAB求导运算

MATLAB里面可以调用diff函数求函数的导函数;

例如diff(f,x,5)就是求符号函数f在x这个点处的5阶导数;

接下来我们使用MATLAB求一下y=x^2*exp(-x)的3阶导数;

6.MATLAB实现函数的泰勒幂级数展开

taylor里面就是我们要展开的函数表达式;

6.MATLAB积分实验

我们使用int(s,x)就是表示对s表达式里面的x进行积分;

顺便提一句:无论是积分求解还是极限的求解,我们都是使用inf代表无穷;

下面是例子:

(需要了解的是这个不定积分的求解是没有常数项的(实际上数学里面的求解是由常数的))

当然,MATLAB也是可以计算定积分的:

对于积分的上下限,我们在int(s,x,a,b)里面添加2个参数,a是积分下限,b是积分的上限

7.MATLAB实现微分方程

Dny表示y的n阶导,D2y就是2阶导;我们可以调用dsolve函数实现常微分方程的求解;

这个板块不同的版本软件的代码语言有所差别,这个是以我的2023为例

(1)我们去求dy/dx+y/x=sinx/x的通解:

diff就是求导数的运算;

(2)求解dy/dx+y=y^2*exp(-x)满足条件y(0)=-2的特解:

提前说明一下,笔者使用的是旧版本的教材,但是用的是新版MATLAB,2023里面可能会因为不支持原来的某些符号而报警告,但是我们依然是能够得到结果的;

D2y表示的是2阶导,D3y表示的是3阶导,以此类推;

(3)y''-6y'+13y=0的通解;

8.空间曲面(马鞍面&&双曲抛物面)

(1)同样是马鞍面,也叫做双曲抛物面,我们可以使用不同的函数进行绘制,不同的函数的展示的效果是不一样的,我们上面的是使用mesh函数绘制,这个函数主要是连接成网格点,下面的这个有颜色的马鞍面是用这个surf进行绘制的,代码的其他部分不变,只需要把第4行的函数的名字改变一下就可以了,通过这个案例我们也可以发现这两者之间的差别,前者线条有颜色,线条之间的补面没有颜色的填充,后者的曲面图线条是黑色的,之间使用颜色进行了填充;

(2)MATLAB为我们提供了grid函数决定是否显示分格线,他的调用格式为

grid od:显示分格子线

grid off:不显示分格子线

(3)相信你一定注意到了这个里面的点运算,点运算是MATLAB里面很有特色的一种运算,我们都知道,矩阵在进行乘法操作的时候,并不是对应的元素进行相乘,但是如果我们对于矩阵里面的元素使用点乘运算,就可以实现对应位置的元素相乘,在绘图里面这个运算有十分广泛的运用;

(4)这个里面的meshgrid也是一个函数,这个函数的功能就是生成网格数据,当X轴和Y轴的数据点个数相等的时候,我们就可以使用[x,y]=meshgrid(x,y);

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

如若内容造成侵权/违法违规/事实不符,请联系瞎胡闹网进行投诉反馈,一经查实,立即删除!

相关文章

体系班第十七节(经典递归)

1汉诺塔 从左移到最右,圆盘必须满足小压大原则 写一个大方法,大方法包括两步:第一步将最后一个圆盘上面的所有的放到第二个塔上面,然后将最后一个圆盘放到最后塔上面,再把第二个塔上面圆盘全放在第三个塔上面 #incl…

IO流(5)——转换流

不同编码读取出现乱码的问题 解决方法 字符输入转换流(InputStreamReader)

泰克Tektronix TDP1000高压差分探头

181/2461/8938产品概述&#xff1a; TDP1000、TDP0500和P6251高压差分探头提供当今开关电源&#xff08;SMPS&#xff09;、CAN/LIN总线和高速数字系统设计所需的出色高速电气和机械性能。 关键性能规格 1 GHz和500 MHz探头带宽<1 pF differential input capacitance1mω…

Postman-Installation has failed

如图&#xff1a; 解决方法&#xff1a; 打开文件夹 Postman-win64-Setup 点击Postman.exe 即可

解压即用,2024最简单好用AI开源换脸应用,整合包已备好

软件整合包&#xff1a;点击下载 关键词&#xff1a;#AI换脸 #开源应用 #可视化界面 #实时换脸 #高清修复 #多个模型 #人脸遮挡处理 #模糊修复 #性能优化 #操作简单 总结&#xff1a;本软件是一款2024年最强大、最易用的AI换脸开源应用。该应用界面经过汉化&#xff0c;操作简…

3D Occupancy 预测冠军方案:FB-OCC

文章结尾有视频和连接 背景知识 Occupancy 更像是一个语义分割任务&#xff0c;但是它是 3D 空间的语义分割它的我们对 Occupancy 分自己的期望是它能够具有通用的这种目标建模的能力&#xff0c;才能够不是不受制于这种目标框这种几何的矩形的这种约束而能够建模任意形状的这…

零信任SDP是什么,有什么作用

物理边界曾经是可信网络和不可信网络之间的有效分割。防火墙通常部署于网络的边缘&#xff0c;基于静态策略来限制网路流量。位于防火墙内部的用户会被授予较高信任等级来访问企业的敏感资源&#xff0c;因为他们被默认是可信的。 然而随着云计算、移动互联、物联网、人工智能…

JavaSE-----认识异常【详解】

目录 一.异常的概念与体系结构&#xff1a; 1.1异常的概念&#xff1a; 1.2一些常见的异常&#xff1a; 1.3异常的体系结构&#xff1a; 1.4异常的分类&#xff1a; 二.异常的处理机制&#xff1a; 2.1 抛出异常&#xff1a; 2.2异常的捕获&#xff1a; 2.3try-catch-&…

【模拟string函数的实现】

提示&#xff1a;文章写完后&#xff0c;目录可以自动生成&#xff0c;如何生成可参考右边的帮助文档 目录 前言 模拟string函数的实现 浅拷贝 深拷贝 vs和g下string结构的说明 总结 前言 模拟string函数的实现 浅拷贝 深拷贝 总结 前言 世上有两种耀眼的光芒&#…

如何理解Linux文件IO?

一、文件IO的概述 1、什么是文件&#xff1f; Linux下一切皆文件。普通文件、目录文件、管道文件、套接字文件、链接文件、字符设备文件、块设备文件。 2、什么是IO&#xff1f; input output&#xff1a;输入输出 3、什么是文件IO&#xff1f; 对文件的输入输出&#xff0c;把…

奇怪的比赛(Python,递归,状态压缩动态规划dp)

目录 前言&#xff1a;题目&#xff1a;思路&#xff1a;递归&#xff1a;代码及详细注释&#xff1a; 状态压缩dp&#xff1a;代码及详细注释&#xff1a; 总结&#xff1a; 前言&#xff1a; 这道题原本是蓝桥上的题&#xff0c;现在搜不到了&#xff0c;网上关于此题的讲解…

JWT令牌的使用

1、什么是jwt JWT (JSON Web Token) 是一种基于 JSON 的轻量级的开放标准&#xff08;RFC 7519&#xff09;&#xff0c;用于在不同系统之间安全地传输信息。JWT 由三部分组成&#xff0c;分别是头部&#xff08;Header&#xff09;、载荷&#xff08;Payload&#xff09;和签…

计算机通识——01.进制转换

前言 学习资料来自 C训练以及CSDN各博主的博客整合而来&#xff0c;内容涵盖计算机通识内容&#xff1a;进制转换、信息单位、数据校验、多媒体基础参数、HTTP \ HTTPS协议、OSI七层模型、IP基础 \ IPv6、网络拓扑机构、域名解析、常用网络命令和端口、数据结构常识等内容&…

【C++】了解一下编码

个人主页 &#xff1a; zxctscl 如有转载请先通知 文章目录 1. 前言2. ASCII编码3. unicode4. GBK5. 类型转换 1. 前言 看到string里面还有Template instantiations&#xff1a; string其实是basic_string<char>&#xff0c;它还是一个模板。 再看看wstring&#xff1…

win下 VirtualBox 自动启动脚本脚本

文章目录 一、找到VBoxManage二、测试脚本1、打开cmd2、输入命令 (直接把上面找到的VBoxManage.exe 拖入到cmd中&#xff0c;这样就不用输入路径了)3、效果展示 比如虚拟机中的系统名称叫“centos-mini” 三、设置自动启动脚本1、复制刚才测试好的命令到新建文本中2、修改文本名…

7. 字符串和集合(重点)

常见API API &#xff08;全称 Application Programming Interface&#xff1a;应用程序编程接口&#xff09;就是别人写好的一些程序&#xff0c;给咱们程序员直接拿去调用即可解决问题的。 1 包 1.1 什么是包&#xff1f; 包是用来分门别类的管理各种不同程序的&#xff…

python知识点总结(一)

这里写目录标题 一、什么是WSGI,uwsgi,uWSGI1、WSGI2、uWSGI3、uwsgi 二、python中为什么没有函数重载&#xff1f;三、Python中如何跨模块共享全局变量?四、内存泄露是什么?如何避免?五、谈谈lambda函数作用?六、写一个函数实现字符串反转&#xff0c;尽可能写出你知道的所…

说下你对TCP以及TCP三次握手四次挥手的理解?

参考自简单理解TCP三次握手四次挥手 什么是TCP协议&#xff1f; TCP( Transmission control protocol )即传输控制协议&#xff0c;是一种面向连接、可靠的数据传输协议&#xff0c;它是为了在不可靠的互联网上提供可靠的端到端字节流而专门设计的一个传输协议。 面向连接&a…

算法——前缀和之除自身以外数组的乘积、和为K的子数组、和可被K整除的子数组、连续数组、矩阵区域和

这几道题对于我们前面讲过的一维、二维前缀和进行了运用,包含了面对特殊情况的反操作 目录 4.除自身以外数组的乘积 4.1解析 4.2题解 5.和为K的子数组 5.1解析 5.2题解 6.和可被K整除的子数组 6.1解析 6.2题解 7.连续数组 7.1题解 7.2题解 8.矩阵区域和 8.1解析 …

Java基础入门day09

day09 万年历综合案例 说明&#xff1a;1900年的1月1日是礼拜一&#xff0c;所有后面的任何一天到底是礼拜几&#xff0c;一定是一个固定值 所有的日历都会从1900年1月1日是礼拜一开始算起 整体思路&#xff1a; 我们可以计算用户输入年份和月份距离1900年1月1日总共有多少天&…