扬州大学MATLAB实验7.1.3
- 格式:doc
- 大小:67.00 KB
- 文档页数:11
MATLAB实验报告一、实验目的本次 MATLAB 实验旨在深入了解和掌握 MATLAB 软件的基本操作和应用,通过实际编程和数据处理,提高解决问题的能力,培养编程思维和逻辑分析能力。
二、实验环境本次实验使用的是 MATLAB R2020a 版本,运行在 Windows 10 操作系统上。
计算机配置为英特尔酷睿 i5 处理器,8GB 内存。
三、实验内容(一)矩阵运算1、矩阵的创建使用直接输入、函数生成和从外部文件导入等方式创建矩阵。
例如,通过`1 2 3; 4 5 6; 7 8 9` 直接输入创建一个 3 行 3 列的矩阵;使用`ones(3,3)`函数创建一个 3 行 3 列元素全为 1 的矩阵。
2、矩阵的基本运算包括矩阵的加减乘除、求逆、转置等。
例如,对于两个相同维度的矩阵`A` 和`B` ,可以进行加法运算`C = A + B` 。
3、矩阵的特征值和特征向量计算通过`eig` 函数计算矩阵的特征值和特征向量,加深对线性代数知识的理解和应用。
(二)函数编写1、自定义函数使用`function` 关键字定义自己的函数,例如编写一个计算两个数之和的函数`function s = add(a,b) s = a + b; end` 。
2、函数的调用在主程序中调用自定义函数,并传递参数进行计算。
3、函数的参数传递了解值传递和引用传递的区别,以及如何根据实际需求选择合适的参数传递方式。
(三)绘图功能1、二维图形绘制使用`plot` 函数绘制简单的折线图、曲线等,如`x = 0:01:2pi; y = sin(x); plot(x,y)`绘制正弦曲线。
2、图形的修饰通过设置坐标轴范围、标题、标签、线条颜色和样式等属性,使图形更加清晰和美观。
3、三维图形绘制尝试使用`mesh` 、`surf` 等函数绘制三维图形,如绘制一个球面`x,y,z = sphere(50); surf(x,y,z)`。
(四)数据处理与分析1、数据的读取和写入使用`load` 和`save` 函数从外部文件读取数据和将数据保存到文件中。
(最新版)MATLAB实验报告实验一典型环节的MATLAB仿真一、实验目的1.熟悉MATLAB桌面和命令窗口,初步了解SIMULINK功能模块的使用方法。
2.通过观察典型环节在单位阶跃信号作用下的动态特性,加深对各典型环节响应曲线的理解。
3.定性了解各参数变化对典型环节动态特性的影响。
二、SIMULINK的使用MATLAB中SIMULINK是一个用来对动态系统进行建模、仿真和分析的软件包。
利用SIMULINK功能模块可以快速的建立控制系统的模型,进行仿真和调试。
1.运行MATLAB软件,在命令窗口栏“>>”提示符下键入simulink命令,按Enter 键或在工具栏单击按钮,即可进入如图1-1所示的SIMULINK仿真环境下。
2.选择File菜单下New下的Model命令,新建一个simulink 仿真环境常规模板。
3.在simulink仿真环境下,创建所需要的系统。
以图1-2所示的系统为例,说明基本设计步骤如下:1)进入线性系统模块库,构建传递函数。
点击simulink下的“Continuous”,再将右边窗口中“Transfer Fen”的图标用左键拖至新建的“untitled”窗口。
2)改变模块参数。
在simulink仿真环境“untitled”窗口中双击该图标,即可改变传递函数。
其中方括号内的数字分别为传递函数的分子、分母各次幂由高到低的系数,数字之间用空格隔开;设置完成后,选择OK,即完成该模块的设置。
3)建立其它传递函数模块。
按照上述方法,在不同的simulink 的模块库中,建立系统所需的传递函数模块。
例:比例环节用“Math”右边窗口“Gain”的图标。
4)选取阶跃信号输入函数。
用鼠标点击simulink下的“Source”,将右边窗口中“Step”图标用左键拖至新建的“untitled”窗口,形成一个阶跃函数输入模块。
5)选择输出方式。
用鼠标点击simulink下的“Sinks”,就进入输出方式模块库,通常选用“Scope”的示波器图标,将其用左键拖至新建的“untitled”窗口。
实验总要求1、封面必须注明实验名称、实验时间和实验地点,实验人员班级、学号(全号)和姓名等。
2、内容方面:注明实验所用设备、仪器及实验步骤方法;记录清楚实验所得的原始数据和图像,并按实验要求绘制相关图表、曲线或计算相关数据;认真分析所得实验结果,得出明确实验结论。
3、图形可以打印出来并剪贴上去,文字必须用标准试验纸手写。
实验一MATLAB绘图基础一、实验目的了解MATLAB常用命令和常见的内建函数使用。
熟悉矩阵基本运算以及点运算。
掌握MATLAB绘图的基本操作:向量初始化、向量基本运算、绘图命令plot,plot3,mesh,surf 使用、绘制多个图形的方法。
二、实验内容建立并执行M文件multi_plot.m,使之画出如图的曲线。
三、实验方法(参考程序)024681012Plot of y=sin(2x) and its derivative四、实验要求1. 分析给出的MA TLAB 参考程序,理解MA TLAB 程序设计的思维方法及其结构。
2. 添加或更改程序中的指令和参数,预想其效果并验证,并对各语句做出详细注释。
对不熟悉的指令可通过HELP 查看帮助文件了解其使用方法。
达到熟悉MA TLAB 画图操作的目的。
3. 总结MATLAB 中常用指令的作用及其调用格式。
五、实验思考1、实现同时画出多图还有其它方法,请思考怎样实现,并给出一种实现方法。
(参考程序如下)t=0:pi/100:4*pi;y1=sin(2*t);y2=2*cos(2*t);plot(t,y1,'-b');hold on; %保持原图plot(t,y2,'-g');grid onaxis([0 4*pi -2 2])title('Plot of y=sin(2x) and its derivative')Plot of y=sin(2x) and its deriv ativ e024681012024681012-2-1012xyPlot of y=sin(2x)024681012-2-1012xyPlot derivative of y=sin(2x);y=2cos(2x)t=0:pi/100:4*pi; y1=sin(2*t); y2=2*cos(2*t);024681012-2-1.5-1-0.500.511.52Plot of y=sin(2x) and its deriv ativ et=0:pi/100:4*pi; y1=sin(2*t); y2=2*cos(2*t); plot(t,y1,'r--'); hold on ;plot(t,y2,'-b'); grid onaxis([0 4*pi -2 2])title('Plot of y=sin(2x) and its derivative')2468101214Plot of y=sin(2x)xyPlot of y=sin(2x) and its deriv ativ exyt=0:pi/100:4*pi; y1=sin(2*t); y2=2*cos(2*t); plot(t,y1,'r--');title('Plot of y=sin(2x)'); xlabel('x'),ylabel('y'); figure(2) plot(t,y2,'-b');title('Plot of y=sin(2x) and its derivative') xlabel('x'),ylabel('y'); grid onaxis([0 4*pi -2 2])2、思考三维曲线(plot3)与曲面(mesh, surf)的用法,(1)绘制参数方程233,)3cos(,)3sin()(t z e t t y e t t t x t t ===--的三维曲线;t=0:pi/30:10*pi;plot3(t.^3.*sin(3.*t).*exp(-t),t.^3.*cos(3.*t).*exp(-t),t.^2);2(2)绘制二元函数xyy xe x x y xf z ----==22)2(),(2,在XOY 平面内选择一个区域(-3:0.1:3,-2:0.1:2),然后绘制出其三维表面图形。
《MATLAB与控制系统仿真》实验报告一、实验目的本实验旨在通过MATLAB软件进行控制系统的仿真,并通过仿真结果分析控制系统的性能。
二、实验器材1.计算机2.MATLAB软件三、实验内容1.搭建控制系统模型在MATLAB软件中,通过使用控制系统工具箱,我们可以搭建不同类型的控制系统模型。
本实验中我们选择了一个简单的比例控制系统模型。
2.设定输入信号我们需要为控制系统提供输入信号进行仿真。
在MATLAB中,我们可以使用信号工具箱来产生不同类型的信号。
本实验中,我们选择了一个阶跃信号作为输入信号。
3.运行仿真通过设置模型参数、输入信号以及仿真时间等相关参数后,我们可以运行仿真。
MATLAB会根据系统模型和输入信号产生输出信号,并显示在仿真界面上。
4.分析控制系统性能根据仿真结果,我们可以对控制系统的性能进行分析。
常见的性能指标包括系统的稳态误差、超调量、响应时间等。
四、实验步骤1. 打开MATLAB软件,并在命令窗口中输入“controlSystemDesigner”命令,打开控制系统工具箱。
2.在控制系统工具箱中选择比例控制器模型,并设置相应的增益参数。
3.在信号工具箱中选择阶跃信号,并设置相应的幅值和起始时间。
4.在仿真界面中设置仿真时间,并点击运行按钮,开始仿真。
5.根据仿真结果,分析控制系统的性能指标,并记录下相应的数值,并根据数值进行分析和讨论。
五、实验结果与分析根据运行仿真获得的结果,我们可以得到控制系统的输出信号曲线。
通过观察输出信号的稳态值、超调量、响应时间等性能指标,我们可以对控制系统的性能进行分析和评价。
六、实验总结通过本次实验,我们学习了如何使用MATLAB软件进行控制系统仿真,并提取控制系统的性能指标。
通过实验,我们可以更加直观地理解控制系统的工作原理,为控制系统设计和分析提供了重要的工具和思路。
七、实验心得通过本次实验,我深刻理解了控制系统仿真的重要性和必要性。
MATLAB软件提供了强大的仿真工具和功能,能够帮助我们更好地理解和分析控制系统的性能。
第一次练习教学要求:熟练掌握Matlab 软件的基本命令和操作,会作二维、三维几何图形,能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。
补充命令vpa(x,n) 显示x 的n 位有效数字,教材102页fplot(‘f(x)’,[a,b]) 函数作图命令,画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位(1-9班)或4位(10班以上) 1.1 计算30sin limx mx mx x →-与3sin limx mx mxx →∞- syms xlimit((902*x-sin(902*x))/x^3) ans =366935404/3limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos1000xmxy e =,求''y syms xdiff(exp(x)*cos(902*x/1000),2) ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算221100xy e dxdy +⎰⎰dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.13941.4 计算4224x dx m x +⎰ syms xint(x^4/(902^2+4*x^2)) ans =(91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求syms xdiff(exp(x)*cos(902*x),10) ans =-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.6 0x =的泰勒展式(最高次幂为4).syms xtaylor(sqrt(902/1000+x),5,x) ans =-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,)n n n x x x n --=+=用循环语句编程给出该数列的前20项(要求将结果用向量的形式给出)。
信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。
2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。
4.学会运用MATLAB进行连续信号时移、反折和尺度变换。
5.学会运用MATLAB进行连续时间微分、积分运算。
6.学会运用MATLAB进行连续信号相加、相乘运算。
7.学会运用MATLAB进行连续信号的奇偶分解。
二、实验任务将实验书中的例题和解析看懂,并在MATLAB软件中练习例题,最终将作业完成。
三、实验内容1.MATLAB软件基本运算入门。
1). MATLAB软件的数值计算:算数运算向量运算:1.向量元素要用”[ ]”括起来,元素之间可用空格、逗号分隔生成行向量,用分号分隔生成列向量。
2.x=x0:step:xn.其中x0位初始值,step表示步长或者增量,xn为结束值。
矩阵运算:1.矩阵”[ ]”括起来;矩阵每一行的各个元素必须用”,”或者空格分开;矩阵的不同行之间必须用分号”;”或者ENTER分开。
2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。
3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例:计算一个函数并绘制出在对应区间上对应的值。
2).MATLAB软件的符号运算:定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语:>>plot(x,y)2).输出多个图像表顺序:例如m和n表示在一个窗口中显示m行n列个图像,p表示第p个区域,表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围:axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母:xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前:title(‘……’)6).输出:grid on举例1:举例2:3.matlab程序流程控制1).for循环:for循环变量=初值:增量:终值循环体End2).while循环结构:while 逻辑表达式循环体End3).If分支:(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号:y=k*exp(a*t)举例:2).正弦信号:y=k*sin(w*t+phi)3).复指数信号:举例:4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号:y=tripuls(t,width,skew)(skew的取值在-1~+1之间,若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算:语句格式:d iff(function,’variable’,n)Function:需要进行求导运算的函数,variable:求导运算的独立变量,n:求导阶数2).连续时间信号的积分运算:语句格式:int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件,并且能够输入简单的语句并输出相应的结果和波形图,也在一定程度上巩固了c语言的一些语法。
MATLAB实验指导书(共5篇)第一篇:MATLAB实验指导书MATLAB 实验指导书皖西学院信息工程学院实验一 MATLAB编程环境及简单命令的执行一、实验目的1.熟悉MATLAB编程环境二、实验环境1.计算机2.MATLAB7.0集成环境三、实验说明1.首先应熟悉MATLAB7.0运行环境,正确操作2.实验学时:2学时四、实验内容和步骤1.实验内容(1)命令窗口的使用。
(2)工作空间窗口的使用。
(3)工作目录、搜索路径的设置。
(4)命令历史记录窗口的使用。
(5)帮助系统的使用。
(6)了解各菜单的功能。
2.实验步骤(1)启动MATLAB,熟悉MATLAB的桌面。
(2)进入MATLAB7.0集成环境。
(3)在命令窗口执行命令完成以下运算,观察workspace的变化,记录运算结果。
1)(365-52⨯2-70)÷3 2)>>area=pi*2.5^2 3)已知x=3,y=4,在MATLAB中求z:x2y3 z=2(x-y)4)将下面的矩阵赋值给变量m1,在workspace中察看m1在内存中占用的字节数。
⎡162313⎤⎢511108⎥⎥m1=⎢⎢97612⎥⎢⎥414151⎣⎦执行以下命令>>m1(2 , 3)>>m1(11)>>m1(: , 3)>>m1(2 : 3 , 1 : 3)>>m1(1 ,4)+ m1(2 ,3)+ m1(3 ,2)+ m1(4 ,1)5)执行命令>>helpabs 查看函数abs的用法及用途,计算abs(3 + 4i)6)执行命令>>x=0:0.1:6*pi;>>y=5*sin(x);>>plot(x,y)7)运行MATLAB的演示程序,>>demo,以便对MATLAB有一个总体了解。
五、思考题1、以下变量名是否合法?为什么?(1)x2(2)3col(3)_row (4)for2、求以下变量的值,并在MATLAB中验证。
matlab操作实验报告一、实验目的1、学会使用matlab建立.m文件。
2、学会二机五节点的潮流计算计算原理。
3、学会使用matlab命令窗编写程序。
4、学会潮流计算的matlab的程序。
5、学会matlab中simulink模块库的模型用法。
二、实验原理潮流分布(1)描述电力系统运行状态的技术术语,它表明电力系统在某一确定运行方式和接线方式下,系统从电源经网络到负荷各处的电压、电流、功率的大小和方向的分布情况。
(2)电力系统潮流分布主要取决于负荷的分布、电力网参数、以及和供电电源之间的关系。
潮流计算的方法1)建立描述电力系统运行状态的数学模型;2)确定解算数学模型的方法;3)制定程序框图,编写计算机计算程序,并进行计算;4)对计算结果进行分析。
对图1所示电力系统拓扑图进行潮流计算该拓扑图为2机5节点的系统应用Matpower计算潮流技巧的核心在于输入好三个矩阵和部分参数,清晰的知道输入参数、矩阵中每一个元素的含义。
参数一、MATPOWER CaseFormat:Version2mpc.version='2';解释:目前普遍采用2形式的算法。
参数二、system MVA base mpc.baseMVA=100;解释:采用有铭值图1mpc.baseMVA=100;(Matpower只能计算有铭值得网络)矩阵一%%bus data %bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin解释:bus data母线参数也就是我们所说的节点参数,下面逐条注释:1bus number(positive integer):第一列表示节点的编号(括号里面注释正整数); 2bus type:第二列表示节点的类型,一般只用得到1、2、3三种节点类型,4类型的节点目前没有接触到。
PQ bus=1PV bu=2reference bus=3isolated bu=43Pd,(MW):表示负荷所需要的有功功率(所有数据都是正数)(有铭值)。
matlab实验报告Matlab实验报告实验题目:利用Matlab进行数据处理与分析实验目的:通过使用Matlab进行数据处理与分析的实践,掌握Matlab的基本操作和数据处理的方法。
实验内容:1. 数据读取:从文件中读取原始数据,并进行初步的观察和分析。
2. 数据清洗:对原始数据进行清洗,包括去除空值、异常值等。
3. 数据可视化:利用Matlab的绘图函数,将数据可视化呈现,以方便对数据进行进一步的分析和理解。
4. 数据分析:对清洗后的数据进行统计分析,包括求平均值、方差、相关系数等。
5. 模型建立:根据数据分析的结果,建立合适的数学模型,并使用Matlab进行参数估计和模型验证。
6. 结果验证:利用实验数据和模型进行实验结果的对比,验证模型的准确性和可靠性。
实验步骤:1. 准备实验数据:从实验样本中获得原始数据,并将其存储为文本文件。
2. 使用Matlab导入数据:使用Matlab的数据导入函数,将文本文件中的数据导入到Matlab的工作空间中。
3. 数据分析与处理:使用Matlab的数据处理函数,对导入的数据进行清洗和处理,去除异常值和空值,并进行初步的数据观察和分析。
4. 数据可视化:利用Matlab的绘图函数,绘制数据的直方图、散点图、折线图等,以展示数据的分布和趋势。
5. 数据统计分析:使用Matlab的统计分析函数,对处理后的数据进行统计分析,包括计算平均值、方差、相关系数等。
6. 模型建立与验证:根据数据分析的结果,建立合适的数学模型,并使用Matlab进行参数估计和模型验证。
7. 结果对比和讨论:将实验结果与模型预测结果进行对比,并进行结果的讨论和分析。
8. 实验结论:总结实验结果并给出结论。
实验结果:根据实验数据的分析和处理,得出如下结论:1. 数据呈现正态分布,符合正态性假设。
2. 数据之间存在显著的正相关关系,相关系数为0.8,结果具有统计学意义。
3. 建立的数学模型与实验数据拟合良好,模型预测结果与实验结果吻合度高。
Matlab 实验任务书一、 matlab语言基础(一)实验目的1、学习了解MA TLAB的实验环境。
2、学习MA TLAB命令的基本操作。
3、学习M文件的编辑及基本操作(二)实验步骤1、学习了解MA TLAB的实验环境通过启动MA TLAB,来学习其实验环境,通过以下命令来了解matlab的使用,并且要求用中文作简要的记录。
help control;help step ;help impulse;save及load 命令的使用。
2、学习MA TLAB命令的基本操作。
要求掌握矩阵的基本操作:输入、转置、求逆、行列式、矩阵提取、矩阵大小等以及+,—,*(乘),/(左除),\(右除),^(幂)和点运算的操作,对于每一项内容,自己给出例子进行实验,并记录内容和结果。
3、学习M文件的编辑及基本操作分别就下列内容编写一个程序和一个函数程序,比较两者的异同:(1)编写一个累加程序,给定累加起始值、中止值和累加的步长。
记录程序内容及调用运行的结果。
(2)编写一个带参数的累加程序,返回参数为累加结果,输入参数为累加起始值、中止值和累加的步长。
记录程序内容及调用运行的结果。
(三)实验报告要求按照上述内容编写实验程序完成实验,并记录实验过程、程序及实验结果,完成实验报告二、控制系统的数学描述(一)实验目的1、学习系统数学模型的表示方法和转换的方法。
2、学习控制系统各类标准型的matlab求解方法。
3、学习控制系统稳定性的分析方法。
(二)实验内容与步骤1、通过学习系统数学模型的表示方法和转换的方法,按照下列内容完成实验并记录实验过程及实验结果。
给定系统Num=[1 1.3 2 2.5];den=[1 0.3 1.2 1]分别求出其状态空间表达式、零极点表达式及部分分式表达式,以及他们之间的转换函数。
2、给定系统分别利用matlab函数求出这两个系统进行串联、并联、负反馈连接后的系统状态方程及传递函数。
3、利用rmodel函数(用法利用帮助获得)构造一个三阶系统,求出该系统的特征向量矩阵和特征值标准型,并利用特征方程验证特征向量的正确性。
1.实验目的(1)学习MATLAB语言的基本符号运算;(2)学习MATLAB语言的矩阵符号运算;(3)MAPLE函数调用。
2.实验内容在下面的实验操作中,认真理解和查询每项操作的作用和目的。
(1)基本符号运算1)符号微分,积分。
>> syms t>> f1=sin(2*t);>> df1=diff(f1)df1 =2*cos(2*t)>> if1=int(f1)if1 =-1/2*cos(2*t)2)泰勒级数展开。
>> tf1=taylor(f1,8)tf1 =2*t-4/3*t^3+4/15*t^5-8/315*t^7>> tf1=taylor(f1,8)tf1 =2*t-4/3*t^3+4/15*t^5-8/315*t^73)符号代数方程求解。
>> syms a b c x>> f=a*x^2+b*x+c;>> ef=sole(f)>> ef=solve(f)ef =[ 1/2/a*(-b+(b^2-4*a*c)^(1/2))][ 1/2/a*(-b-(b^2-4*a*c)^(1/2))]4)符号微分方程求解>> f='D2x+2*Dx+10*x=0';g='Dx(0)=1,x(0)=0';>> dfg=dsolve(f,g)dfg =1/3*exp(-t)*sin(3*t)5)积分变换>> syms t>> f1=exp(-2*t)*sin(5*t); >> F1=laplace(f1)F1 =5/((s+2)^2+25)>> syms s>> F2=1/(s+2)^2;>> f2=ilaplace(F2)f2 =t*exp(-2*t)>> syms t>> f1=sin(t);>> Fz1=ztrans(f1)Fz1 =z*sin(1)/(z^2-2*z*cos(1)+1)>> syms z>> Fz2=z/(z-2)+z/(z-3); >> fz2=iztrans(Fz2)fz2 =2^n+3^n6)计算精度。
>> g=sym('[pi/2 2;3 exp(1)]')g =[ pi/2, 2][ 3, exp(1)]>> g1=vpa(g)g1 =[ 1.5707963267948966192313216916398, 2.][ 3., 2.7182818284590452353602874713527]>> digitsDigits = 32>> g1=vpa(g,5)g1 =[ 1.5708, 2.][ 3., 2.7183](2)符号矩阵运算1)创建与修改符号矩阵。
>> G1=sym('[1/(s+1),s/(s+1)/(s+2);1/(s+1)/(s+2),s/(s+2)]')G1 =[ 1/(s+1), s/(s+1)/(s+2)][ 1/(s+1)/(s+2), s/(s+2)]>> G2=subs(G1,G1(2,2),'0')G2 =[ 1/(s+1), s/(s+1)/(s+2)][ 1/(s+1)/(s+2), 0]>> G3=G1(1,1)G3 =1/(s+1)2)符号线性代数inv 符号矩阵求逆det 符号矩阵行列式eig 符号矩阵特征值transpose 符号矩阵转置charpoly 符号特征多项式3)常规符号运算。
>> syms s>> d1=1/(s+1);d2=1/(s+2);d=d1*d2d =1/(s+1)/(s+2)>> ad=sym('[s+1 s;0 s+2]');G=d*adG =[ 1/(s+2), s/(s+1)/(s+2)][ 0, 1/(s+1)]>> n1=[1 2 3 4 5];n2=[1 2 3];>> p1=poly2sym(n1);p2=poly2sym(n2); >> p=p1+p2p =x^4+2*x^3+4*x^2+6*x+8>> pn=sym2poly(p)pn =1 2 4 6 8>> a=[0 1;-2 -3];>> gcha=poly2sym(poly(eig(a)))gcha =x^2+3*x+2>> sym sans =s>> mg=s*eye(2)-amg =[ s, -1][ 2, s+3]>> geig=transpose(eig(mg))geig =[ s+2, s+1]>> gdet=det(mg)gdet =s^2+3*s+2>> groot=transpose(solve(gdet))groot =[ -2, -1]>> G=inv(mg)G =[ (s+3)/(s^2+3*s+2), 1/(s^2+3*s+2)][ -2/(s^2+3*s+2), s/(s^2+3*s+2)]>> g21=ilaplace(G(1,1));>> g22=ilaplace(G(2,2));>> g12=ilaplace(G(1,2));>> g21=ilaplace(G(2,1));>> g=[g11 g12;g21 g22]g =[ -exp(-2*t)+2*exp(-t), exp(-t)-exp(-2*t)][ -2*exp(-t)+2*exp(-2*t), 2*exp(-2*t)-exp(-t)](3)MAPLE函数调用maple 调用maple内核函数mhelp maple库函数帮助命令1)MAPLE下的泰勒级数展开。
>> maple('readlib(mtaylor);');>> tmap=maple('mtaylor(sin(2*t),t)')tmap =2*t-4/3*t^3+4/15*t^5>> syms t;f=sin(2*t);>> tmat=taylor(f)tmat =2*t-4/3*t^3+4/15*t^5>> mhelp mtaylormtaylor - multivariate Taylor series expansionCalling Sequence:mtaylor(f, v)mtaylor(f, v, n)mtaylor(f, v, n, w)Parameters:f - an algebraic expressionv - a list or set of names or equationsn - (optional) non-negative integerw - (optional) list of positive integersDescription:- The mtaylor function computes a multivariate Taylor series expansion of the input expression f, with respect to the variables v, to order n, using the variable weights w.- The result type has changed in Maple V to be simply a polynomial in the variables with no order term.- The variables v can be a list or set of names or equations. If v[i] is an equation, then the left-hand side of v[i] is the variable, and the right-hand side is the point of expansion. If v[i] is a name, then v[i] = 0 is assumed as the point of expansion.- If the third argument n is present then it specifies the ``truncation order''of the series. The concept of ``truncation order'' used is ``total degree'' in the variables. If n is not present, the truncation order used is the value of the global variable Order, which is 6 by default.- If the fourth argument w is present it specifies the variable weights to beused (by default all 1). A weight of 2 will halve the order in thecorresponding variable to which the series is computed.- Note: mtaylor restricts its domain to ``pure'' Taylor series, those serieswith non-negative powers in the variables.- This function must be loaded using readlib(mtaylor).Examples:> readlib(mtaylor):> mtaylor(sin(x^2+y^2), [x,y]);2 2x + y> mtaylor(sin(x^2+y^2), [x,y], 8);2 2 6 2 4 4 2 6x + y - 1/6 x - 1/2 y x - 1/2 y x - 1/6 y> mtaylor(sin(x^2+y^2), [x,y], 8, [2,1]);2 2 6y + x - 1/6 y> mtaylor(sin(x^2+y^2), [x=1,y], 3);2 2sin(1) + 2 cos(1) (x - 1) + (-2 sin(1) + cos(1)) (x - 1) + cos(1) y> mtaylor(f(x,y), [x,y], 3);2f(0, 0) + D[1](f)(0, 0) x + D[2](f)(0, 0) y + 1/2 D[1, 1](f)(0, 0) x2+ x D[1, 2](f)(0, 0) y + 1/2 D[2, 2](f)(0, 0) ySee Also:taylor, series, coeftayl2)MAPLE库函数与MA TLAB函数。