MATLAB实验参考2018

合集下载

MATLAB语言及其应用实验报告.doc

MATLAB语言及其应用实验报告.doc一、实验目的1、熟悉MATLAB语言及其基本操作；2、掌握利用MATLAB进行数据分析和可视化的方法；3、掌握MATLAB应用于科学计算的基本方法。

二、实验环境MATLAB R2018a三、实验内容及方法3.1 实验内容1、打开MATLAB环境；2、读入数据文件，并对数据进行分析和处理；3、通过MATLAB绘制图表，对数据进行可视化。

四、实验细节及流程4.1 读入数据文件本实验使用的数据文件为sales.xlsx，其中包括2018年各个季度的销售数据。

首先，通过以下命令导入数据：data=xlsread('sales.xlsx');4.2 数据分析和处理1、计算各季度销售总额通过以下代码计算每个季度的销售额并求和，得到每年的销售总额：Q1=data(:,2);TotalSales=sum([Q1 Q2 Q3 Q4],2);2、计算增长率根据每年的销售总额，计算出每年的增长率。

具体代码如下：GrowthRate(1)=0;for i=2:length(TotalSales)GrowthRate(i)=((TotalSales(i)-TotalSales(i-1))/TotalSales(i-1))*100; endSalesGrowth=[TotalSales GrowthRate];3、计算每个季度的均值和标准差meanQ1=mean(Q1);stdQ1=std(Q1);4、计算出每年第一个季度的销售额所占比例首先，我们将第一个季度的销售额单独提出来，具体代码如下：4.3 数据可视化1、柱形图对于销售总额，使用柱形图进行可视化，具体代码如下：结果如下图所示：图1 销售总额2、线性图3、箱形图boxplot([Q1 Q2 Q3 Q4],{'Q1','Q2','Q3','Q4'});图3 每个季度的销售额4、饼图pie(FirstQSalesRatio(:,1));五、结论本实验通过对销售数据的分析和可视化，得出以下结论：1、2018年销售总额呈逐年上升趋势，其中2017年到2018年的增长率最高；2、每年第一个季度的销售额所占比例在40%至45%之间，与其他季度相比，显著高于其他季度；3、2018年第二季度的销售额经过调整后，表现出了相对较高的波动。

平行平面腔自再现模FoxLi数值迭代解法及MATLAB实现

激光原理课程设计题目：方形镜平行平面腔自再现模Fox -Li数值迭代解法及MATLAB实现院系理学院专业班级0910101学生姓名指导教师提交日期2018 年 4 月15 日目录、设计目地1b5E2RGbCAP、设计要求和设计指标I plEanqFDPw、设计内容l DXDiTa9E3d3.1Fox-Li 3.2matlab3.2.13.2.23.2.3 平行平面腔地迭代解法1RTCrpUDGiT 实现35PCzVD7HxA迭代解法地过程3jLBHrnAILg程序实现 4X HAQX74J0X自再现模形成地判断6LDAYtRyKfE3.3 GUI 界面地制作8Zzz6ZB2Ltk四、本设计改进建议9dvzfvkwMI1五、设计感想9rqyn14ZNXI六、主要参考文献9EmxvxOtOco附录一、设计目地为了加深对激光原理中Fox-Li平行平面腔地迭代解法地理解，学习matlab地使用,锻炼运用数值方法解决专业问题地能力.SixE2yXPq5二、设计要求和设计指标在matlab用Fox-Li平行平面腔地迭代解法求得激光器腔镜面上地光场分布求出距离镜面中点为x处地光场地振幅A和相位P,并作出二维图像.在得出地一维图像地基础上作出镜面上光强地二维分布.6ewMyirQFL二、设计内容3.1 Fox-Li平行平面腔地迭代解法谐振腔是激光器必备条件之一，它使激光反复通过增益物质，从而实现光地自激振荡.在激光地发展史上最早提出地是平行平面腔,又称为F—P腔,它由两块平行平面反射镜组成，第一台红宝石激光器地谐振腔就是用它来做成地.kavU42VRUs对于开放式光腔，镜面上稳态场分布地形成可以看成是光在两个界面间往返传播地结果.因此，两个界面上地场必然是互相关联地：一个镜面上地场可以视为由另一个镜面上地场所产生，于是求解镜面上稳态场地分布问题就归结为求解一个积分方程.y6v3ALoS89 考虑在开腔中往返传播地一列波.设初始时刻在镜I上有某一个场分布u1,则当波在腔中经第一次渡越而到达镜II时,将在镜II上形成一个新地场分布u2,场u2经第二次渡越后又将在镜I上形成一个新地场分布u3.每次渡越时，波都将因为衍射损失一部分能量,并引起能量分布变化，如此重复下去……由于衍射主要是发生在镜地边缘附近，因此在传播过程中，镜边缘附近地场将衰落得更快，经多次衍射后所形成地场分布，其边缘振幅往往都很小＜与中心处比较），具有这种特征地场分布受衍射地影响也将比较小. 可以预期：在经过足够多次渡越之后，能形成这样一种稳态场：分布不再受衍射地影响在腔内往返一次后能够“再现”出发时地场分布，即实现了模地“自再现”，具体过程图 1 所示：M2ub6vSTnP图1开腔中自再现模地形成光学中地惠更斯一菲涅尔原理是从理论上分析衍射问题地基础，该原理地严格数学表示是菲涅尔一基尔霍夫衍射积分.设已知空间任意曲面S 上光波场地振幅和相位分布函数为u(x ； y ),由它所要考察地空间任一点P 处场分布为u(x, y),二者之间有以下关系式：OYujCfmUCw_Lk ■? e s u(x',y'—(1cos^dSS:- 式中,'为(x , y )与(x, y)连线地长度，B 为S 面上点(x , y )处地法线和上述连线之间地夹角，ds •为S 面上地面积元,k 为波矢地模•而对于方形镜平行平面镜eUts8ZQVRdu(x,y) = ikj(x-x)2(y-y)2L2将'按(x-X')/L, (y-y')/L 地幕级数展开,当满足a2/L'L (L/a)2和b2/l_L (L/b)2时才W x；x')2+2(y p _ikL 呦(x；x')2+(y；L')2]e =e 2L2L =ee方形镜中a 二b .3.2 matlab 实现3.2.1迭代解法地过程本文采用Fox — Li 数值迭代法得到了了镜面上自再现模在 x 方向地分布并推广到整个镜面,最终动态显示每次渡越镜面上光场分布.sQsAEJkW5T虽然是复数积分，但其和实数积分实现方法相同，即取一定步长，用矩形面积地和代替函数地定积分.F 面是程序框图:从而得到v(x, y)2 2i，kL a b ik[—L]二(~T)e ~」v(x',y')e2L 2Ldx'dy''L将上式分离变量.令v(x, y)二 v(x)v(y)得到v(x)二 x aK x (x,x')vdx..av(y)二 y .；K y (y,y')vdy ^K x (x,x') (X _x')2 ~2L-K y (y,y')(y _y ')2~2L~ike322程序实现源程序：clear,clc global steps L k a lamda=input('波长lamda='＞。

matlab2018部分中文部分英文

matlab2018部分中文部分英文如果您在MATLAB 2018中遇到部分中文部分英文的情况，可能是由于某些因素导致中文字符和英文字符混合在一起。

以下是一些可能的原因和解决方法：1. 字体问题：确保您使用的字体支持中文字符。

如果您的编辑器或系统默认字体不支持中文，您需要更改字体设置。

在MATLAB编辑器中，选择“工具”>“选项”>“字体”，选择支持中文字符的字体，例如“宋体”或“微软雅黑”。

2. 语言环境设置：检查您的MATLAB语言环境设置。

在MATLAB命令窗口中输入以下命令：```matlabsetenv('MATLAB_LANG','zh_CN')```这将设置MATLAB语言环境为中文。

如果您想设置为英文，可以将'zh_CN'替换为'en_US'。

3. 文件编码：确保您的MATLAB脚本或代码文件使用正确的编码格式。

如果文件编码不正确，可能会导致中文字符和英文字符混合在一起。

使用文本编辑器（如Notepad++）打开文件，并检查并更改文件编码为UTF-8。

4. MATLAB版本问题：确保您使用的是支持中文的MATLAB版本。

MATLAB的不同版本可能对中文支持有所不同。

如果您的MATLAB版本不支持中文，您可以考虑升级到较新的版本。

5. 操作系统语言环境：确保您的操作系统语言环境设置为中文。

如果您的操作系统语言环境设置为英文，可能会导致MATLAB中的中文字符显示为乱码。

您可以进入操作系统设置，将语言环境更改为中文。

尝试以上方法，并检查是否可以解决您的问题。

如果问题仍然存在，请考虑寻求技术支持或与MATLAB社区联系以获得帮助。

我的MATLAB（R2018b）界面布局

MATLAB的布局可以根据个人的偏好来调整，并可以将调整好的Layout保存。
可以用鼠标拖动各个窗口来调整位置，还可以通过上图中的Layout选项来选择部分窗口是否显示、调整一些工具栏的位置，点击SaveLayout可以保存调整好的这是我调整好的布局，MATLAB界面的主要部分是编辑器，右侧是workspace、command window的位置，上方是HOME、PLOTS等选项卡。
可以用鼠标拖动各个窗口来调整位置还可以通过上图中的layout选项来选择部分窗口是否显示调整一些工具栏的位置点击savelayout可以保存调整好的这是我调整好的布局matlab界面的主要部分是编辑器右侧是workspacecommandwindow的位置上方是homeplots等选项卡
我的 MATLAB（ R2018b）界面布局
---2021年7月10日添加以下内容：后来又将背景换成了绿色，效果如下
பைடு நூலகம்

MATLAB基础教程(第四版) 第1章 MATLAB R2018b简介

2019年9月27日
第6页
MATLAB R2018b 基础教程
MATLAB 概述（3/3）
MATLAB有以下几个特点。
高级科学计算语言。代码、数据文件的集成管理环境。算法设计开发的交互式工具。用于线性代数、统计、傅立叶分析、滤波器设计、优化和数值计算的基本数学
函数。 2-D和3-D数据可视化。创建自定义工程师图形界面的工具。与第三方算法开发工具—— C/C++、FORTRAN、Java、COM、Microsoft
Real-Time Windows Target
2019年9月27日
第10页
MATLAB R2018b 基础教程
启动 MATLAB R2018b （1/1）
选择“开始”|“所有程序”| MATLAB| R2018b | MATLAB R2018b命令，或者直接双击桌面上的MATLAB图标，启动 MATLAB R2018b。
2019年9月27日
第5页
MATLAB R2018b 基础教程
MATLAB 概述（2/3）
MATLAB可以被广泛地应用于不同领域，如信号与图像处理、控制系统设计与仿真、通信系统设计与仿真、测量测试与数据采集、金融数理分析以及生物科学等。在MATLAB中内嵌了丰富的数学、统计和工程计算函数，使用这些函数进行问题的分析解答，无论是问题的提出还是结果的表达都采用工程师习惯的数学描述方法，这一特点使MATLAB成为了数学分析、算法开发及应用程序开发的良好环境。MATLAB是MathWorks产品家族中所有产品的基础。附加的工具箱扩展MATLAB基本环境用于解决特定领域的工程问题。
Excel—— 集成开发基于 MATLAB的算法。

2018年高教社杯数学建模c题matlab

2018年高教社杯数学建模C题Matlab一、背景介绍2018年高教社杯全国大学生数学建模竞赛是由我国高等教育学会主办的一项全国性学科竞赛。

该竞赛以促进大学生数学建模能力的培养和提高为宗旨，得到了全国各地高校的广泛参与和支持。

其中，C题是该竞赛中的一个重要环节，涉及使用Matlab进行数学建模，要求参赛选手具备一定的编程和数学建模能力。

二、竞赛题目2018年高教社杯数学建模C题的具体内容是基于某一物理问题，通过建立数学模型并利用Matlab编程分析解决相应问题。

题目要求参赛选手熟练掌握Matlab的基本语法和数学建模方法，能够灵活应用各种算法和技巧解决实际问题。

三、题目分析该竞赛题目所涉及的物理问题可能涉及到动力学、流体力学、热传导等多个领域，因此参赛选手需要具备较强的物理基础知识和数学建模能力。

利用Matlab进行编程要求选手具备一定的计算机编程基础和数值计算能力。

四、Matlab应用在解决数学建模问题时，Matlab是一种非常优秀的数学建模工具，它具有强大的数学计算能力和丰富的绘图函数，能够有效地帮助选手分析和解决复杂的数学问题。

Matlab还支持各种算法的实现和优化，能够帮助选手提高数学建模的效率和精度。

五、比赛经验共享参加2018年高教社杯数学建模C题的竞赛选手，可以共享自己在数学建模和Matlab编程过程中的经验和收获。

他们可以讲述自己在解题过程中所遇到的困难和挑战，以及如何克服这些困难，找到合适的解题方法。

他们可以展示自己对数学建模和Matlab编程的独特见解和理解，这对其他竞赛选手也是一种宝贵的学习和借鉴。

六、总结2018年高教社杯数学建模C题的竞赛不仅考察了选手在数学建模和Matlab编程方面的能力，还促进了选手之间的学习和交流。

通过共享各自的经验和思考，可以帮助选手更好地提高自己在数学建模和Matlab应用方面的能力，促进全国范围内大学生数学建模水平的提高。

这也为相关领域的研究和教育提供了一个重要的交流评台，促进了学术研究和创新成果的产生。

matlabr2018教程习题

matlabr2018教程习题matlabr2018是由微软公司提供的一套针对数据库的管理系统。

它主要用于对数据库进行维护、添加、删除和管理。

在应用 MATLAB软件进行数据库维护时，需通过 MATLAB来进行数据库管理。

该系统主要是为数据中心而设计，能够提供包括文件管理、统计分析、可视化等功能。

通过将系统和数据库连接起来，使数据库管理更简单也更强大。

......一、选项中选项是指可以作为变量使用的，但在实际使用过程中，不能单独作为一个变量使用。

选项选项是不能单独作为一个变量使用的。

当项目数据为()时，系统就会自动地将该选项项包含到项目中进行设置。

在出现这种情况时，如果是由于数据需要进行查找时，我们就会将该选项项包含到项目中。

如果不存在上述情况的话，系统就会将该选项项包含到项目中进行显示。

其中包含选项项作为变量使用时，如果项目出现多个选项变量时也会自动地将该选项项包含到项目中。

二、“在一个时钟周期内，数据会发生多少次改变？答：通常情况下，一个时钟周期内只会发生一次变化。

但有时，为了便于比较，我们会把这种情况称为“非同步”。

例如，在时钟周期发生了两次变化时，我们可以将这两次变化称为“同步”，但同时会把两次变化称为“非同步”。

因此，两次变数之间要存在一定的时间间隔。

”但如果需要的话，那么这就需要通过其他方法来实现了。

三、如果不能将多个变量的变量值加在一起，就会形成多个连续变量。

”4、变量名称时，应使用什么变量命名规则？①变量名称的使用规则：可以分为：可替换变量和重复变量三种：前两种命名规则较为简单、适用于单一变量情况，后一种命名规则较为复杂、适用于多个变量情况。

②变量名称的使用规则：变量名称一般使用名称与其所在的字符串相匹配。

③变量名称与字符串相同时，要使用一个组合变量名称；如果有多个组合变量要同时使用两个组合变量名称时（通常以第一个为限),必须使用组合变量的第一个字符串名称；如果需要两个组合变量名称时（通常以第一个为限）可以使用一个组合单元格作为变量组或按以下规则命名：a.变量名称与其字符串具有相同含义。

2018实验1认识matlab环境

X X X X 大学（自动化学院）课内实验报告实验名称：认识matlab编程环境专业名称：班级：学生姓名：学号：指导教师：实验日期：实验1 认识matlab编程环境一、实验目的（1）熟悉启动和退出matlab软件方法（2）掌握主要窗口的使用方法二、实验设备及条件( 1).计算机（2）Matlab7.0集成环境三、实验内容及要求（1）启动matlab,熟悉matlab界面（2）练习下面指令：clear ,clc, help, doc, who, whos 熟悉各命令的作用并写出总结.（3）建立自己的当前工作目录mywork建立自己搜索路径：点 FILE 菜单下的setpath,如下图行进设置：在add folder 下添加搜索路径，然后保存。

(4) 在命令窗口中执行命令完成以下运算，观察workspace变化，记录运算结果。

1）计算(365-52*2-70)/32）计算area=pi*2.5^2h = exp(log(realmax))toobig = pi*h3）已知x=3,y=4,z=x+y在Matlab中求z4）执行命令>>help abs 查看函数abs的用法及用途，计算abs(3+4i)5）字符串赋值：>>a=’This is my book‘实验结果四、心得体会通过本次实验，使我们了解并熟悉了MATLAB的运行环境，以及启动及退出的方法，并且掌握主要窗口的使用方法。

因为是第一次接触到matlab软件，对于matlab这门技术并不是很了解，对它并不能很熟悉的运用，有很多地方需要多次尝试与练习，但是随着学习的深入我对matlab的了解也相应深入，了解到Matlab有很多的优点，计算方便、快捷、准确、操作方便、流程简单。

虽然所以可能会有些困难，但我会好好学习这门课程。

matlab机械臂模拟正弦运动规律

一、概述机械臂是一种能够模仿人类手臂运动的装置，广泛应用于工业生产线、医疗设备和科学研究等领域。

在机械臂的控制过程中，正弦运动是一种常见的运动模式，它具有周期性和规律性的特点，在很多实际应用中都有重要的作用。

本文将介绍如何使用Matlab软件对机械臂进行正弦运动的模拟与控制。

二、正弦运动的基本原理1. 正弦运动的定义正弦运动是指在运动过程中，物体的位移、速度或加速度随时间变化呈现出正弦曲线的规律。

其数学表达式为：y = A*sin(ωt + φ)。

其中，A表示振幅，ω表示角速度，φ表示相位差，t表示时间。

2. 正弦运动的特点正弦运动具有周期性、规律性和简谐性的特点，可以用于描述很多实际运动的规律。

在机械臂控制中，通过对正弦运动的模拟与控制，可以实现机械臂的精确定位、柔性操作等功能。

三、Matlab软件在机械臂正弦运动模拟中的应用1. Matlab工具箱的使用Matlab提供了丰富的工具箱，包括控制系统工具箱、仿真工具箱等，这些工具箱可以很好地支持机械臂正弦运动的模拟与控制。

在进行机械臂正弦运动模拟时，可以充分利用Matlab提供的工具箱进行快速开发和验证。

2. 机械臂建模在进行机械臂正弦运动模拟之前，首先需要对机械臂进行建模。

利用Matlab的建模工具，可以将机械臂的结构、动力学特性等进行建模，并得到相应的数学模型。

3. 正弦运动控制算法在机械臂正弦运动的控制过程中，需要设计相应的控制算法。

Matlab提供了丰富的控制系统设计工具，可以根据机械臂的数学模型，设计出适合的正弦运动控制算法。

4. 仿真与验证利用Matlab进行机械臂正弦运动的仿真与验证是非常重要的一步。

通过对模型进行仿真，可以验证控制算法的有效性，确保机械臂能够稳定、精确地进行正弦运动。

四、实例分析：利用Matlab对机械臂进行正弦运动模拟以某型号机械臂为例，介绍如何利用Matlab对机械臂进行正弦运动模拟的具体步骤。

1. 机械臂建模利用Matlab的建模工具对该型号机械臂进行建模，包括机械结构、电机特性、传感器等。

MATLAB作业1参考答案(2018)

MATLAB 作业1 参考答案1 、编写一个矩阵相加函数mat_add() ，使其具体的调用格式为A=mat_add(A1 ,A2 ,A3 ,…) ，要求该函数能接受任意多个矩阵进行加法运算。

（注：varargin 变量的应用）【求解】可以编写下面的函数，用varargin 变量来表示可变输入变量function A=mat_add(varargin)A=0;for i=1:length(varargin), A=A+varargin{i}; end如果想得到合适的错误显示，则可以试用try, catch 结构。

function A=mat_add(varargin)tryA=0;for i=1:length(varargin), A=A+varargin{i}; endcatch, error(lasterr); end2、请绘制出sin()xy 的三维图和等高线。

【求解】给出下面命令即可得出的图形。

>> [x,y]=meshgrid(-pi:.1:pi);surf(x,y,sin(x.*y)), figure; contour(x,y,sin(x.*y),30)3、试求出如下极限。

222222001cos()lim ()x y x y x y x y e +→→-++【求解】极限问题可以由下面语句直接求解。

>> fc=(1-cos(x^2+y^2))*exp(x^2+y^2)/(x^2+y^2);limit(limit(fc,x,0),y,0)ans =4、试求出下面函数的导数。

22atan ln()y x y x=+ 【求解】为隐函数，故需要用隐函数求导公式得出导数。

>> syms x,y; f=atan(y/x)-log(x^2+y^2);f1=simple(-diff(f,x)/diff(f,y))f1 =(y+2*x)/(x-2*y)5、假设20(,)xyt f x y e dt -=⎰，试求222222x f f f y x x y y ∂∂∂-+∂∂∂∂。

控制系统仿真实验报告

控制系统仿真实验报告班级：测控1402班姓名：王玮学号：072018年01月实验一经典的连续系统仿真建模方法一实验目的:1 了解和掌握利用仿真技术对控制系统进行分析的原理和步骤。

2 掌握机理分析建模方法。

3 深入理解阶常微分方程组数值积分解法的原理和程序结构，学习用Matlab编写数值积分法仿真程序。

4 掌握和理解四阶Runge-Kutta法，加深理解仿真步长与算法稳定性的关系。

二实验内容:1. 编写四阶 Runge_Kutta 公式的计算程序，对非线性模型（3）式进行仿真。

（1）将阀位u 增大10％和减小10％，观察响应曲线的形状;（2）研究仿真步长对稳定性的影响，仿真步长取多大时RK4 算法变得不稳定（3）利用 MATLAB 中的ode45()函数进行求解，比较与（1）中的仿真结果有何区别。

2. 编写四阶 Runge_Kutta 公式的计算程序，对线性状态方程（18）式进行仿真（1）将阀位增大10％和减小10％，观察响应曲线的形状;（2）研究仿真步长对稳定性的影响，仿真步长取多大时RK4 算法变得不稳定（4）阀位增大10％和减小10％，利用MATLAB 中的ode45()函数进行求解阶跃响应，比较与（1）中的仿真结果有何区别。

三程序代码:龙格库塔:%RK4文件clccloseH=[,]';u=; h=1;TT=[];XX=[];for i=1:h:200k1=f(H,u);k2=f(H+h*k1/2,u);k3=f(H+h*k2/2,u);k4=f(H+h*k3,u);H=H+h*(k1+2*k2+2*k3+k4)/6;TT=[TT i];XX=[XX H];end;hold onplot(TT,XX(1,:),'--',TT,XX(2,:)); xlabel('time')ylabel('H')gtext('H1')gtext('H2')hold on水箱模型:function dH=f(H,u)k=;u=;Qd=;A=2;a1=;a2=;dH=zeros(2,1);dH(1)=1/A*(k*u+Qd-a1*sqrt(H(1)));dH(2)=1/A*(a1*sqrt(H(1))-a2*sqrt(H(2)));2编写四阶 Runge_Kutta 公式的计算程序，对线性状态方程（18）式进行仿真：1 阀值u对仿真结果的影响U=;h=1; U=;h=1;U=;h=1;2 步长h对仿真结果的影响:U=;h=5; U=;h=20;U=;h=39 U=;h=50由以上结果知,仿真步长越大,仿真结果越不稳定。

[精品]最速下降路线的确定--实验报告

[精品]最速下降路线的确定--实验报告实验目的：通过编程实现最速下降法，确定多元函数的最优解实验环境：MATLAB R2018b实验原理：最速下降法是求解函数极小值的一种方法，该方法通过在当前点沿着当前点梯度方向移动到下一个点，使得函数值下降最快，直到达到极小值。

具体实现过程如下：1. 给定初始点x0和判断停止迭代的精度tol2. 计算x0处的梯度g03. 如果梯度的范数小于tol，则停止迭代4. 将x0移动到x1 = x0 - α*g0，其中α是步长，通常通过线性搜索确定5. 重复2-4步，直到达到精度要求或者达到最大迭代次数实验步骤：1. 定义多元函数function f = fun(x)f = 4*x(1)^2 + 2*x(1)*x(2) + 3*x(2)^2;3. 最速下降法实现function x = steepestDescent(x0, tol, max_iter)i = 1;x(:,i) = x0;g = grad(x0);alpha = 0.1; % 步长while norm(g) >= tol && i <= max_iterp = -g; % 下降方向x(:,i+1) = x(:,i) + alpha*p;g = grad(x(:,i+1));i = i+1;endend4. 运行程序并输出结果x = steepestDescent([3;3],1e-6,1000);disp(['极小值点：(', num2str(x(1,end)), ',', num2str(x(2,end)), ')']);disp(['极小值：', num2str(fun(x(:,end)))]);disp(['迭代次数：', num2str(size(x,2)-1)]);运行结果：极小值点：(-0.3333333331214,-0.1111111112736)极小值：0.4444444445835迭代次数：12实验结论：通过最速下降法，可以准确地求解多元函数的极小值点和极小值。

《控制系统matlab仿真实训》参考文献

《控制系统matlab仿真实训》参考文献
以下是《控制系统matlab仿真实训》参考文献：
1. 陈红. MATLAB在控制系统仿真中的应用[J]. 科技资讯, 2019(36): 159-160.
2. 朱宏伟, 李玉峰. 基于MATLAB的控制系统仿真实验教学研究[J]. 实验技术与管理, 2018(1): 19-21.
3. 王晓峰, 赵丽华. MATLAB在控制系统仿真实验中的应用[J]. 实验技术与管理, 2017(1): 24-26.
4. 吴兆华, 杨志宏. MATLAB在控制系统仿真中的应用[J]. 实验技术与管理, 2016(1): 28-30.
5. 赵宇, 刘洋. 基于MATLAB的控制系统仿真实验教学研究[J]. 实验技术与管理, 2015(1): 25-27.
6. 李晓峰, 张晓芳. MATLAB在控制系统仿真中的应用[J]. 实验技术与管理, 2014(1): 28-30.
7. 陈晓明, 王建华. 基于MATLAB的控制系统仿真实验教学研究[J]. 实验技术
与管理, 2013(1): 25-27.
8. 王宏伟, 李明. MATLAB在控制系统仿真中的应用[J]. 实验技术与管理, 2012(1): 28-30.
9. 陈红, 王晓峰. 基于MATLAB的控制系统仿真实验教学研究[J]. 实验技术与管理, 2011(1): 25-27.
10. 赵宇, 刘洋. MATLAB在控制系统仿真中的应用[J]. 实验技术与管理, 2010(1): 28-30.
以上是《控制系统matlab仿真实训》参考文献，供您参考。

matlab证明涡旋光束自聚焦

matlab证明涡旋光束自聚焦matlab证明涡旋光束自聚焦引言光束的自聚焦现象一直是光学研究领域的一个重要课题。

特别是涡旋光束的自聚焦现象，因其独特的光学性质和广泛的应用前景而备受关注。

本文将使用MATLAB仿真工具来证明涡旋光束的自聚焦现象，并探讨其涡旋数和光束参数对自聚焦过程的影响。

1. 涡旋光束的特性涡旋光束是一种相位呈螺旋状分布的光束，具有自旋角动量和轨道角动量等特性。

它的相位分布可以表示为exp(ilφ)，其中l是涡旋数，φ是光束的极角。

涡旋光束在传播过程中会产生自聚焦现象，即光束的横向尺寸会自动调整，使得整个光束逐渐变细，并最终形成一个强度极大的光束核心。

2. 自聚焦方程模型涡旋光束的自聚焦过程可以通过非线性薛定谔方程来描述。

在本文中，我们使用一个经典的非线性薛定谔方程模型来模拟涡旋光束的自聚焦现象：i ∂A/∂z + Δ⊥A + β|A|^2 A = 0 (1)其中A是复振幅，z是传播距离，Δ⊥是横向Laplacian算子，β是非线性系数，|A|^2 表示复振幅的模的平方。

3. MATLAB仿真模拟为了证明涡旋光束的自聚焦现象，我们使用MATLAB进行数值仿真。

我们需要设定涡旋数、光束起始尺寸和非线性系数等参数。

我们使用有限差分法或薛定谔方程传播算法对方程(1)进行数值求解。

在MATLAB中，我们可以定义一个自定义的函数来代表方程(1)。

通过使用该函数，我们可以通过迭代法计算出光束在传播过程中的强度和相位分布。

4. 涡旋光束的自聚焦过程在MATLAB中，我们可以观察到涡旋光束的自聚焦过程。

随着传播距离的增加，光束的横向尺寸逐渐减小，同时光束的强度也增强。

当传播距离足够大时，光束核心会形成一个明显的峰值，而光束边缘则变得非常模糊。

5. 涡旋数和光束参数的影响通过改变涡旋数和光束参数，我们可以观察到涡旋光束自聚焦过程的差异。

当涡旋数增加时，光束的自聚焦速度加快，并且光束核心的强度增加。

MATLAB作业2参考答案（2018年）

MA TLAB 作业二参考答案1、试求下面线性代数方程的解析解与数值解，并检验解的正确性。

2932114010110503848246303356684953X -----⎡⎤⎡⎤⎢⎥⎢⎥----⎢⎥⎢⎥=⎢⎥⎢⎥---⎢⎥⎢⎥------⎣⎦⎣⎦【求解】求出A , [A;B ] 两个矩阵的秩，可见二者相同，所以方程不是矛盾方程，应该有无穷多解。

>> A=[2,-9,3,-2,-1; 10,-1,10,5,0; 8,-2,-4,-6,3; -5,-6,-6,-8,-4]; B=[-1,-4,0; -3,-8,-4; 0,3,3; 9,-5,3]; [rank(A), rank([A B])] ans =4 4用下面的语句可以求出方程的解析解，并可以验证该解没有误差。

>> x0=null(sym(A));x_analytical=sym(A)\B; syms a; x=a*[x0 x0 x0]+x_analyticalx =[ a+967/1535, a-943/1535, a-159/1535][ -1535/1524*a, -1535/1524*a, -1535/1524*a][ -3659/1524*a-1807/1535,-3659/1524*a-257/1535,-3659/1524*a-141/1535] [ 1321/508*a+759/1535, 1321/508*a-56/1535, 1321/508*a-628/1535] [ -170/127*a-694/307, -170/127*a+719/307, -170/127*a+103/307] >> A*x-B ans = [ 0, 0, 0][ 0, 0, 0] [ 0, 0, 0] [ 0, 0, 0]用数值解方法也可以求出方程的解，但会存在误差，且与任意常数a 的值有关。

>> x0=null(A); x_numerical=A\B; syms a; x=a*[x0 x0 x0]+x_numerical; vpa(x,10)ans =[ .2474402553*a+.1396556436, .2474402553*a-.6840666849, .2474402553*a-.1418420333][-.2492262414*a+.4938507789,-.2492262414*a+.7023776988e-1,-.2492262414*a+.3853511888e-1][ -.5940839201*a, -.5940839201*a, -.5940839201*a][ .6434420813*a-.7805411315, .6434420813*a-.2178190763,.6434420813*a-.5086089095][-.3312192394*a-1.604263460, -.3312192394*a+2.435364854, -.3312192394*a+.3867176824] >> A*x-B[ 1/18014398509481984*a, 1/18014398509481984*a, 1/18014398509481984*a] [ -5/4503599627370496*a, -5/4503599627370496*a, -5/4503599627370496*a][ -25/18014398509481984*a, -25/18014398509481984*a, -25/18014398509481984*a] [ 13/18014398509481984*a, 13/18014398509481984*a, 13/18014398509481984*a]2、求解下面的联立方程，并检验得出的高精度数值解（准解析解）的精度。

用MATLAB求二元泰勒展开式

第九章用MATLAB 求二元泰勒展开式1级数求和命令S=symsum(u,t,a,b)的功能是计算级数和S=bt a u =∑。

其中u 是包含符号变量t 的表达式，是待求和级数的通项。

当u 的表达式中只含一个变量时，参数t 可省略。

例9.11.1判断下列级数是否收敛，如收敛则求其和：11n n ∞=∑，211n n ∞=∑，220(1)n n x x ∞=+∑ 解创建符号变量n 和x,用symsum 命令计算各级数的和：syms n x ↙symsum(1/n,1,inf) ↙ans=inf 知级数11n n ∞=∑发散至无穷大。

symsum(1/n^2,1,inf) ↙ans=1/6*pi^2 知级数211n n ∞=∑收敛，且其和为26π 对级数220(1)n n x x ∞=+∑，由于其通项中包含两个变量x 和n ，故使用symsum 命令时须指定求和变量是n ：un=x^2/(1+x^2)^n; ↙symsun(un,n,0,inf) ↙ans=1+x^2得和函数为1+x^2对有些级数，symsum 命令不能求得其和，从而也无法得知其敛散性。

此时，可使用MATLAB 的数值计算功能进行处理。

例9.11.2试求级数211ln(1)n n ∞=+∑的和解用symsum 命令求解：syms n ↙symsum(log(1+1/n^2),1,inf) ↙ans=sum(log(1+1/n^2),n=1..inf)此结果表示symsum 命令不能求得其和。

我们转而采用数值方法计算部分和211ln(1)n n k S k ==+∑。

将下面的程序存入一个m 文件中： clear alln=9000; %部分和的项数Sn=0;for k=1:nSn= Sn+log(1+1/k^2);endfprintf(‘Sn=%f,(n=%d)’, Sn,n)执行该程序，显示结果为Sn=1.301735，（n=9000）再对程序中的变量n 分别赋值n=9000 ，n=900000 ，n=9000000并执行程序，得计算结果为：Sn=1.301835，（n=9000 ）Sn=1.301845，（n=900000）Sn=1.301846，（n=9000000）由此看出，随着n 增大，Sn 趋于1.30185。

如何在Matlab2018中配置VLfeat？

如何在Matlab2018中配置VLfeat？最近需要⽤Matlab运⾏SIFT代码，在⽹上⾟⾟苦苦找到了SIFT的代码，但是⼀直在报错报错原因如下：Matlab\vlfeat\toolbox\mex\mexw64\vl_sift.mexw64'⽆效: 缺少依赖共享库:；vl_sift.mexw64'需要 'VCOMP100.DLL'。

刚开始以为是缺少DLL⽂件的原因，跟着⽹上的解决⽅案⼀步步进⾏了操作，但是还是没有效果，后⾯定位到⽂件夹，发现VCOMP100.DLL 存在，后来发现是没有配置VLfeat的原因，⽹上也有很多配置的过程，有些试了没什么反应，现在将我⾃⼰配置VLfeat的过程进⾏⼀个记录，也是为了以后有需要可以有个参考；1.在官⽹下载/install-matlab.html，如果不想⾃⼰编译，⼀定要下载右边对应的，⽐如我是matlab程序就下载下图这个：点击，可以看到官⽹其实已经将具体的步骤写出来了，点击Download 跳转页⾯下载安装过程要注意找到toolbox，如果下载的是编译好了的⽂件vlfeat-0.9.21-bin.tar.gz进⾏解压，先确认⾥⾯有mex⽂件夹，然后直接将vlfeat-0.9.20-bin\复制到matlab根⽬录下的toolbox⾥。

2. 在matlab命令⾏直接输⼊：run('D:\matlabb\toolbox\vlfeat-0.9.21\toolbox/vl_setup')即可。

3.为了判断是否安装成功，在Matlab命令⾏中输⼊：vl_version；如果安装成功会显⽰安装的vlfeat的版本号4.到这个步骤位置，只是进⾏⼀次安装，⼀旦重新打开matlab，就需要重新配置vlfeat，。