matlab上机实验报告
MatlabSimulink课程上机实验汇报
模糊控制
模糊控制概述
模糊控制是以模糊集合论，模糊语言变量及模糊逻辑推理为基础的计算 机智能控制。在传统的控制领域里，控制系统动态模式的精确与否是影 响控制优劣的最主要关键，系统动态的信息越详细，则越能达到精确控 制的目的。传统的控制理论对于明确系统有强而有力的控制能力，但对 于过于复杂或难以精确描述的系统，则显得无能为力。因此便尝试着以 模糊数学来处理这些控制问题。
PID控制器参数的工程整定方法，主要有临界比例法、反应曲线法和衰减法。三种方法各有其特点， 其共同点都是通过试验，然后按照工程经验公式对控制器参数进行整定。但无论采用哪一种方法所得到的 控制器参数，都需要在实际运行中进行最后调整与完善。现在一般采用的是临界比例法。
PID控制
1.2 PID控制器的参数整定
模糊控制
模糊集合
经典集合论中任意一个元素与任意一个集合之间的关系，只是“属于”或“不属于”，两者必居其一 而且只居其一。它描述的是有明确分界线的元素组合。经典集合用0或1简单地表示“属于”或“不属于” 的分类对比，而模糊集合则是把它扩展成用0~1之间的连续变化值来描述元素的属于程度。经典集合中元 素必须符合特征函数，而模糊集合实际上是将经典集合论中的特征函数表示扩展到用隶属度函数来表示。
隶属度函数实质上反映的是事物的渐变性，要遵守一些基本规则。 1．表示隶属度函数的模糊集合必须是凸模糊集合 2．变量所取隶属度函数通常是对称和平衡的 3．隶属度函数要遵从语意顺序和避免不恰当重叠 4．隶属度函数的选择需要考虑重叠指数
模糊控制
模糊集合
隶属度函数举例
MATLAB实验报告
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基础运用实验报告1
（黑体3号）实验报告（黑体3号）课程：Matlab基础运用班级：姓名：学号成绩：指导教师：实验日期及时间：2009.9.16座位号：必修/选修：实验序号：实验名称：（标题均为黑体小4号）一、实验目的与要求：实验目的：Matlab是功能强大的科学及工程计算软件。
是以矩阵计算为基础的计算工具，本次实验通过学生对Matlab运行环境的熟悉和一些简单操作，掌握帮助命令、绘图命令以及矩阵操作等使用方法。
实验要求：1．上机结束时须经指导老师检查实验结果后方可离开。
2．课代表在实验后三天内将实验报告收齐交指导老师。
二、实验原理、内容和步骤（一）Matlab界面和矩阵运算基础1，建立自己的工作目录，将该目录设为当前目录，并将其设置成MATLAB的搜索路径。
2，熟悉Matlab的帮助功能。
采用两种方式获得关于函数“exp”的解释。
（1）利用Help Browser.（2）在Command Window下通过输入命令“help exp”。
（3）利用Help navigator3，从屏幕上输入矩阵A=[1 2 3；4 5 6；7 8 9]或A=[1，2，3；4，5，6；7，8，9]观察输出结果；试用回车代替分号，观察输出结果；4，选择File|new菜单中的M-file，输入B=[9 ，8，7；6，5，4；3，2，1]，保存为tmp.m文件，退出编辑环境。
此时在工作环境中使用tmp命令就可调出B矩阵；5，输入who和whos观察结果，了解其作用；6，用至少两种方式创建以下矩阵。
A=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡987654321>>a=[1 2 3;4 5 6;7 8 9]a =1 2 34 5 67 8 9或使用以下命令调用>> tmpa =1 2 34 5 67 8 9>>a=[1:3;4:6;7:9]a =1 2 34 5 67 8 9（二）基本绘图1， 绘制 sin (1/t), -1<t<1,标上坐标轴和标题.。
Matlab simulink 上机实验报告 简单版
201006113 11002 Matlab上机实验报告
◆实验一: Smulink动态仿真集成环境
➢ 1.目的要求
➢熟悉simulink环境, 掌握simulink的仿真方法。
➢ 2.掌握要点
➢熟悉simulink环境, 掌握simulink的仿真方法。
➢ 3.实验内容
➢熟悉simulink环境；
➢熟悉基本的模块库以及功能模块
➢搭建简单的电路进行仿真；
➢对分析参数对结果的影响；
1.建立如图所示的仿真系统.
完成过程:
********* ***** 结果如下:
◆ 2.建立如图所示的仿真系统.
◆将红色区域部分创建并封装装成子系统
完成过程:
没有设置子系统时:
没有设置子系统时的结果如下:
以下开始设置子系统并封装: 修改变量后:
最终如下图所示:
开始封装设置过程: 设置子系统各个参数
设置完成后如下图所示:
双击设置好的封装并分别输入与变量对应的参数如下:
运行结果如下:。
Matlab上机报告
《MATLAB 课程》课程实习报告（一） 实习日期：2016-9-13学院：数学与统计学院学院 专业：数学与应用数学 姓名：王子璇 学号：201413230061， 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况，并保存全部变量。
（1）21185sin 2ez +=代码：>> z1=2*sin(85*pi/180)/(1+exp(2))截图：（3）()0.3,9.2,8.2,...,8.29.2,0.3,233.0ln 3.0sin 23.03.03---=+++-=-a a a e e z a a 提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。
代码：>> a=-3.0:0.1:3.0;>> z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)截图：（4）5.2:5.0:0其中,32,1221,110,2224=⎪⎩⎪⎨⎧<≤+-<≤-<≤=t t t t t t t t z提示：用逻辑表达式求分段函数值。
代码：>> t=0:0.5:2.5;>> z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3)截图：2， 设有矩阵A 和B⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=25242322212019181716151413121110987654321A ，⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--=11134079423096171603B （1） 求它们的乘积C 。
（2） 将矩阵C 的右下角23⨯子矩阵赋给D 。
（3） 查看MATLAB 工作空间的使用情况。
代码：（1）>> A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25]>> B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]>> C=A*B（2）>> F=size(C);>> D=C(F(1)-2:F(1),F(2)-1:F(2));>> whos截图：（1）（2）3，下面是一个线性方程组：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡52.067.095.061141514111121321x x x（1） 求方程的解。
MATLAB技术上机实践报告
目录上机实践一MATLAB在“高等数学”中的应用练习（4学时）（1）上机练习——MATLAB解常微分方程（SJ101.m）（2）上机练习——MATLAB二元函数作图（SJ102.m）（3）上机编程——微分运算和积分运算（SJ103.m）（4）上机编程——线性方程组的数值解法（SJ104.m）上机实践二MATLAB在“力学”中的应用练习（4学时）（1）上机练习——质点直线运动的位置、速度和加速度（SJ201.m）（2）上机练习——轻质杆的斜抛运动（SJ202.m）（3）上机编程——质点的抛体运动（SJ203.m）（4）上机编程——物体在恒力作用下的直线运动（SJ204.m）上机实践三MATLAB在“电磁学”中的应用练习（4学时）（1）上机练习——电磁波的发射方法一（SJ301.m）（2）上机练习——电磁波的发射方法二（SJ302.m）（3）上机编程——静电场场强分布计算（SJ303.m）（4）上机编程——带电粒子在电场中的运动（SJ304.m）上机实践四MATLAB在“振动和波动”中的应用练习（4学时）（1）上机练习——旋转矢量法（SJ401.m）（2）上机练习——用惠更斯作图法确定波阵面（SJ402.m）（3）上机编程——同方向简谐振动的合成（SJ403.m）（4）上机编程——简谐波横波的演示（SJ404.m）上机实践五MATLAB在“量子物理”中的应用练习（4学时）（1）上机练习——黑体辐射随波长的变化规律（SJ501.m）（2）上机练习——势垒和隧道效应（SJ502.m）（3）上机编程——双缝干涉图样的动画模拟（SJ503.m）上机实践六MATLAB在“电子信息技术”中的应用练习（4学时）（1）上机练习——含受控源的电阻电路计算（SJ601.m）（2）上机练习——方波分解为多次正弦波之和（SJ602.m）（3）上机编程——电阻电路的计算（SJ603.m）（4）上机编程——连续信号的MATLAB描述（SJ604.m）上机实践一 MATLAB 在“高等数学”中的应用练习（4学时）【上机实践目的】使用MATLAB 软件编写相应的M 文件进行计算模拟高等数学中的一些基本问题。
matlab上机实验报告
《MATLAB 语言及其应用实验》上机实验报告1姓名： 学号：请各位同学将命令和结果由MATLAB 中复制（或截屏）后保存至word 文档中。
1. 计算矩阵 与 的和、乘积、各元素的乘积以及Kronecker 积。
2.对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
3.对矩阵912563827⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦进行特征值分解，给出特征值与特征向量。
4.已知向量[4＋5i, 6, 7-1.5j, 5+9j]，以其为对角线元素生成对角矩阵。
5. 生成含有8个元素的高斯向量a，要求元素的均值、方差分别为5、4，再生成含有7个元素的随机向量，要求元素服从上的均匀分布。
求两者的卷积。
6. 求解多项式x3-7x2+2x+40的根。
7. 计算多项式乘法(x2+2x+2)(x2+5x+4)。
8. 对下式进行部分分式展开：9. 已知y=sin(x)，x从0到2 ，间隔 x=0.02 ，求y的最大值、最小值、均值和标准差。
10. 用符号函数法求解方程at2+b*t+c=0。
11. 用符号计算验证三角等式：sin( 1)cos( 2)-cos( 1)sin( 2) =sin( 1- 2)12. 因式分解：13. ⎥⎥⎦⎤⎢⎢⎣⎡=)sin()log(12x x e x x a f ax ，用符号微分求df/dx ，并用符号定积分来验证。
14. 求解非线性方程组1211223sin 2cos x x x x x x -=⎧⎨+=⎩。
>> syms x1 x2>> [x1,x2]=solve(x1-3*x2-sin(x1),2*x1+x2-cos(x2))x1 =0.49662797440907460178544085171994x2 =0.0067214622395756734146654770697884。
(完整word)Matlab实验报告
实验一:Matlab操作环境熟悉一、实验目的1．初步了解Matlab操作环境.2．学习使用图形函数计算器命令funtool及其环境。
二、实验内容熟悉Matlab操作环境,认识命令窗口、内存工作区窗口、历史命令窗口;学会使用format 命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息;学会使用图形函数计算器funtool，并进行下列计算：1．单函数运算操作。
求下列函数的符号导数(1)y=sin（x)；（2） y=(1+x)^3*（2-x）；求下列函数的符号积分(1)y=cos(x）;(2)y=1/（1+x^2）；（3）y=1/sqrt（1—x^2)；(4）y=(x1)/(x+1）/(x+2)求反函数(1)y=（x-1)/(2＊x+3）; （2） y=exp(x）；（3） y=log（x+sqrt（1+x^2));代数式的化简(1）(x+1)*（x-1）*(x-2）/（x-3）/（x—4)；(2）sin（x）^2+cos(x）^2；（3）x+sin(x）+2*x—3*cos(x）+4＊x*sin（x）;2．函数与参数的运算操作。
从y=x^2通过参数的选择去观察下列函数的图形变化(1)y1=(x+1）^2(2）y2=（x+2)^2(3) y3=2＊x^2 (4） y4=x^2+2 (5） y5=x^4 (6） y6=x^2/2 3．两个函数之间的操作求和(1)sin（x)+cos（x) (2） 1+x+x^2+x^3+x^4+x^5乘积(1)exp(—x）＊sin（x） (2） sin(x)＊x商(1)sin(x）/cos（x); （2) x/（1+x^2); （3) 1/(x—1)/(x—2）; 求复合函数(1)y=exp（u) u=sin（x） (2） y=sqrt（u) u=1+exp(x^2）(3） y=sin（u） u=asin（x） (4) y=sinh(u) u=-x实验二：MATLAB基本操作与用法一、实验目的1.掌握用MATLAB命令窗口进行简单数学运算。
MATLAB第一次上机实验报告
第一章上机实验一.实验内容（一）熟悉Matlab的基本命令（1）工作空间管理：①Who和whos：显示当前工作空间的变量情况，例如想查询A和B两个变量的空间占用情况，则可以用命令whos A B。
如果只给出whos命令，则将显示工作空间中全部变量的变量名称、大小、内存占用情况及变量类型。
Who命令给出的结果是简单的，只显示当前工作空间中的变量名称。
②clear：删除工作空间中的变量。
Clear A B表示删除A和B两个变量；若只给出clear命令，则将删除工作空间内的所有变量。
③which：显示出某个MATLAB文件所在目录，如which xxure将给出xxure.m 文件所在目录。
④quit：退出⑤workspace：该命令将打开工作空间浏览器窗口。
（2）路径编辑：①path：显示MATLAB包含路径。
②cd：返回上一级目录。
③addpath：将一个新的字符串表示的子目录名添加到MATLAB包含的路径下。
④rmpath：从当前MATLAB包含的路径下清除某个字符串表示的子目录。
（3）显示格式：①format long：双精度数表示小数点后14位有效数字，单精度数表示小数点后7位有效数字。
②short e：5位科学计数表示。
③short g：从format short和format short e中自动选择最佳计数方式（默认设置）④long g：从format long和format long e中自动选择最佳计数方式⑤hex：十六进制表示⑥rat：近似有理数表示（4）窗口清理：①clc：命令窗口清屏。
②cla：可以清除当前坐标系下的内容。
③clf：图形窗口清屏。
④close：用于关闭所有的图形窗口。
（5）操作系统命令：①copyfile：复制文件②pwd(Identify current directory)：显示当前的工作目录名称。
（6）熟悉Matlab的联机帮助命令：①help：可以达到帮助信息目录。
matlab上机实验报告ppt
matlab上机实验报告pptMATLAB上机实验报告摘要：本实验报告利用MATLAB软件进行了一系列实验，包括数据处理、图像处理、信号处理等内容。
通过实验，我们掌握了MATLAB软件的基本操作和应用技巧，提高了数据分析和处理的能力。
1. 实验背景MATLAB是一种用于数学计算、数据分析和可视化的高级技术计算语言和交互式环境。
它是工程师和科学家们进行算法开发、数据分析、数据可视化和数值计算的首选工具。
本次实验旨在通过实际操作，掌握MATLAB的基本操作和应用技巧。
2. 实验内容本次实验主要包括以下内容：（1）数据处理：利用MATLAB对一组实验数据进行处理，包括数据的导入、清洗、分析和可视化。
（2）图像处理：利用MATLAB对一幅图像进行处理，包括图像的读取、处理和保存。
（3）信号处理：利用MATLAB对一组信号进行处理，包括信号的生成、滤波和频谱分析。
3. 实验过程（1）数据处理：首先，我们利用MATLAB将实验数据导入到工作空间中，然后对数据进行清洗和分析，最后利用MATLAB绘制出数据的可视化图表。
（2）图像处理：我们利用MATLAB读取一幅图像，并对图像进行处理，比如调整图像的亮度、对比度等参数，最后保存处理后的图像。
（3）信号处理：我们利用MATLAB生成一组信号，并对信号进行滤波处理，然后利用MATLAB进行信号的频谱分析。
4. 实验结果通过本次实验，我们成功地利用MATLAB对实验数据进行了处理和分析，得到了清晰的数据可视化图表；对一幅图像进行了处理，并保存了处理后的图像；对一组信号进行了滤波处理，并进行了频谱分析。
实验结果表明，MATLAB是一款功能强大、灵活多样的工程计算软件，能够满足工程师和科学家们的各种需求。
5. 实验结论本次实验通过MATLAB软件的实际操作，使我们掌握了MATLAB的基本操作和应用技巧，提高了我们的数据分析和处理能力。
同时，也加深了我们对MATLAB软件的理解和认识，为今后的工程计算和科学研究打下了坚实的基础。
matlab上机实验心得
MATLAB上机实验心得1. 引言在学习MATLAB课程期间，我们进行了一系列的上机实验。
通过这些实验，我深刻体会到了MATLAB在数学建模和数据分析方面的强大功能。
本文将详细介绍我在实验中的学习心得和体会，并分享一些使用MATLAB进行数据处理和可视化的技巧。
2. 实验一：MATLAB基础在第一次实验中，我们掌握了MATLAB的基本操作和语法。
通过编写简单的脚本，我学会了如何定义变量、进行算术运算、使用条件语句和循环结构等。
我还学会了如何使用MATLAB自带的函数库来解决常见的数学问题。
这次实验让我对MATLAB有了初步的认识，并为后续实验打下了坚实的基础。
3. 实验二：数据处理与可视化在第二次实验中，我们探索了MATLAB在数据处理和可视化方面的能力。
我们使用了一些常见的数据处理函数，如读取文件、筛选数据、计算统计量等。
我们还学习了如何使用plot函数绘制线图、scatter函数绘制散点图以及histogram函数绘制直方图等。
通过这次实验，我意识到MATLAB在数据处理和可视化方面的高效和便捷。
使用MATLAB，我们可以快速地对大量数据进行处理和分析，并通过可视化方式直观地展示数据的特征和规律。
4. 实验三：数学建模第三次实验是最具挑战性的一次，我们需要运用MATLAB解决实际问题并进行数学建模。
在实验中，我们学习了如何将实际问题转化为数学模型，并使用MATLAB求解。
我们通过编写脚本来解决最优化问题、微分方程求解等。
这次实验让我深刻理解了数学建模的重要性，并提高了我的问题解决能力。
MATLAB 的强大计算能力和丰富的函数库为数学建模提供了极大的便利。
5. 实验四：图像处理在第四次实验中，我们学习了MATLAB在图像处理方面的应用。
我们掌握了如何读取、显示、修改和保存图像。
我们还学会了一些常见的图像处理算法，如灰度变换、直方图均衡化、滤波器等。
这次实验让我对图像处理有了初步的认识，并意识到MATLAB在该领域有着广泛应用。
matlab上机实验报告
matlab上机实验报告《Matlab上机实验报告：探索数据分析与图像处理的奇妙世界》摘要：本实验报告基于Matlab平台，探索了数据分析与图像处理的应用。
通过实验，我们发现Matlab具有强大的数据处理和图像处理能力，能够帮助我们快速高效地分析数据和处理图像。
本文将详细介绍实验过程和结果，以及对Matlab在数据分析与图像处理中的应用前景进行探讨。
1. 实验目的本次实验旨在通过Matlab平台进行数据分析与图像处理实验，探索Matlab在这两个领域的应用能力，为学习者提供实践机会，加深对数据分析与图像处理的理解，并为未来的科研工作和工程实践提供技术支持。
2. 实验内容2.1 数据分析我们首先使用Matlab对一组实验数据进行分析，包括数据的可视化、统计特征提取、数据拟合等操作。
通过Matlab强大的数据处理工具，我们可以快速准确地分析数据，提取出其中的规律和特征，为后续的决策和应用提供支持。
2.2 图像处理其次，我们利用Matlab对一组图像进行处理，包括图像的增强、滤波、边缘检测等操作。
通过Matlab丰富的图像处理工具，我们可以对图像进行各种操作，改善图像质量，提取图像特征，甚至实现图像识别和分析。
3. 实验结果在数据分析方面，我们利用Matlab对实验数据进行了可视化展示，并提取了数据的统计特征，得到了数据的分布规律和趋势。
在图像处理方面，我们利用Matlab对图像进行了增强和边缘检测，得到了处理后的图像，观察到了图像质量的改善和边缘特征的提取。
4. 实验结论通过本次实验，我们深刻认识到Matlab在数据分析与图像处理领域的强大能力，能够帮助我们快速高效地分析数据和处理图像，为科研工作和工程实践提供技术支持。
我们相信，在未来的发展中，Matlab将在数据分析与图像处理领域发挥越来越重要的作用。
5. 展望基于本次实验的结果，我们对Matlab在数据分析与图像处理领域的应用前景进行了展望。
我们相信，随着科技的不断发展，Matlab将在数据分析与图像处理领域发挥越来越重要的作用，为人类社会的进步和发展做出更大的贡献。
MATLAB上机实验报告_2
4.已知某系统的开环传递函数为, 计算该系统的相位裕度和幅值裕度。
5.已知某系统的闭环传递函数为, 设ζ＝0.707,ωn＝1:1:5, 画出该系统的单位阶跃响应曲线。
四. 实验设计方案（实验步骤或开发过程）
七. 附录（源程序清单）
1.clear all
closeall
clc
t=0:0.001:3;
u=sin(300*t)+2*cos(200*t);
y=fft(u);
yy=abs(y);
plot(t,yy)
2.d=[1 10 20 30 41 54];
r=roots(d)
3.clear all
closeall
num=1;
den=[12 7 1 0];
bode(num,den)
figure(2)
nyquist(num,den)Байду номын сангаас
4.num=5;
den=conv([1,0],conv([1,4],[1,5]));
sys=tf(num,den);
w=logspace(0,4,50);
bode(sys,w);
grid;
六. 实验结果
五. 实验中存在问题及解决办法
六. 实验结果
六. 实验结果
1. 5、
2.r =
-7.8739
-1.4201 + 0.9942i
-1.4201 - 0.9942i
0.3570 + 1.4679i
0.3570 - 1.4679i由此看出系统不稳定
3.
4.Gm = 36（幅值裕度） Pm =83.5835（相角裕度）
MATLAB上机实验报告
MATLAB上机报告专业:姓名:学号:一、向量运算1.1向量的生成1.直接输入法--------------------------------------------------- 32.冒号表达式法 ------------------------------------------------- 33.函数法 ------------------------------------------------------- 3 2. 1向量的加减和数乘运算-------------------------------------------- 3向量的点、叉积运算--------------------------------------------- 31.点积运算------------------------------------------------------ 32.向量的叉积运算 --------------------------------------------- 3二、矩阵运算2.1矩阵元素的表示及相关操作1.元素的下标表示法2.矩阵元素的赋值2.2矩阵的创建1.直接输入法52.抽取法3.拼接法 ----------------------------------------------- 64.函数法-------------------------------------------------------- 75.拼接函数和变形函数法----------------------------------------- 86.加载法7.M文件法--------------------------------------------- 10 2.3矩阵的代数运算1.求矩阵行列式的值2.矩阵的加减、数乘与乘法--------------------- 103.求矩阵的逆矩阵4.矩阵的除法5.矩阵的秩------------------------ 116.求矩阵的特征值与特征向量7.矩阵的乘嘉与开方运算------------- 128.矩阵的指数与对数9,矩阵转置 ----------------------------------- 1310.矩阵的提取与翻转--------------------------------------------- 14三、数组运算3.1多维数组的创建1.下标赋值法2.工具阵函数法-------------------------------------- 153 .拼接和变形函数法 ---------------------------------------------- 16 3.2数组的代数运算1・数组的加减，数乘与乘法---------------------------------------- 172.数组的除法----------------------------------------------------- 183.数组的乘募与开方4.数组的指数与对数--------------------------- 195.数组或矩阵的单纯转置 ------------------------------------------ 20 3.3数组的关系与逻辑运算1.数组的关系运算------------------------------------------------ 20203 .与逻辑运算相关的函数 ----------------------------------------- 213.4数组和矩阵函数的通用形式------------------------------------------ 22四、字符串运算4.1字符串变量与一维字符数组1.给字符串变量赋值2.一维字符数组的生成------------------------- 22 4.2对字符串的多项操作1.求字符串长度 -------------------------------------------------- 232.字符串与一维数组数值的相互转换------------------------------- 233•比字^^串--------------------------------------------------------- 234•查找字符串5.显示字符串--------------------------------------- 244.3二维字符数组----------------------------------------------- 24五、习题部分1 .单项选择题------------------------------------------- 242 •判断题----------------------------------------------- 24六、实验二部分一、实验目的------------------------------------------------ 26二、实验内容------------------------------------------------ 26三、实验步骤------------------------------------------------ 263.1向量的生成和运算1）向量的生成直接输入法、冒号表达式法----------------------------------- 26函数法、练习题 -------------------------------------------- 272）向量的运算------------------------------------------------- 283.2矩阵的创建、引用和运算1）矩阵的创建和引用直接输入法、抽取法--------------------------------------- 29拼接法，函数法------------------------------------------- 30拼接函数和变形函数法、练习题----------------------------- 31 2）矩阵的运算矩阵的加减、数乘与乘法、求矩阵的逆矩阵------------------- 34矩阵的除法、练习题---------------------------------------- 353. 3多维数组的创建及运算1）多维数组的创建--------------------------------------------- 362）多维数组的运算--------------------------------------------- 37 练习题---------------------------------------------------- 383.4字符串的操作1）字符串的创建---------------------------------------------- 382）求字符串长度---------------------------------------------- 393）字符串与一维数值数组的相互转换------------------------- 39练习题--------------------------------------------------- 40相关知识点及练习结果一、向量运算1.直接输入法» A=[2,3,4,5,6],B=[1;2;3;4;5]/C=[4 5 6 7 8 9];A= 2 3 4 5 6 B= 123452.冒号表达式法» A=l:2:10,b=l:10/C=10:-l:lA= 1 3 5 7 9b= 1 2 3 4 5 6 7 8 9 10C= 10 9 8 7 6 5 4 3 2 13.函数法Linspace()是线性等分函数，logspace ()是对数等分函数» A=linspace(l,30/10),B=logspace(0/4/5)A = 1.0000 4.2222 7.4444 10.6667 13.8889 17.1111 20.333323.5556 26.7778 30.0000B= 1 10 100 1000 100002.1向量的加减和数乘运算只有维数相同的行向量之间可以相加减，维数相同的列向量也可以相加减，标量数值可以与向量直接相乘除。
MATLAB上机实验实验报告
MATLAB上机实验一一、实验目的初步熟悉MATLAB 工作环境,熟悉命令窗口,学会使用帮助窗口查找帮助信息; 命令窗口二、实验内容1 熟悉MATLAB 平台的工作环境;2 熟悉MATLAB 的5 个工作窗口;3 MATLAB 的优先搜索顺序;三、实验步骤1. 熟悉MATLAB 的5 个基本窗口①Command Window 命令窗口②Workspace 工作空间窗口③Command History 命令历史记录窗口④Current Directory 当前目录窗口⑤Help Window 帮助窗口1 命令窗口Command Window;在命令窗口中依次输入以下命令：>>x=1>> y=1 2 34 5 67 8 9;>> z1=1:10,z2=1:2:5;>> w=linspace1,10,10;>> t1=ones3,t2=ones1,3,t3=ones3,1>> t4=ones3,t4=eye4x =1z1 =1 2 3 4 5 6 7 8 9 10t1 =1 1 11 1 11 1t2 =1 1 1t3 =111t4 =1 1 11 1 11 1 1t4 =1 0 0 00 1 0 00 0 1 00 0 0 1思考题：①变量如何声明,变量名须遵守什么规则、是否区分大小写;答：1变量声明1.局部变量每个函数都有自己的局部变量,这些变量只能在定义它的函数内部使用;当函数运行时,局部变量保存在函数的工作空间中,一旦函数退出,这些局部变量将不复存在;脚本没有输入输出参数,由一系列MATLAB命令组成的M文件没有单独的工作空间,只能共享调用者的工作空间;当从命令行调用,脚本变量存在基本工作空间中；当从函数调用,脚本变量存在函数空间中;2.全局变量在函数或基本工作空间内,用global声明的变量为全局变量;例如声明a为全局变量：>>global a声明了全局变量的函数或基本工作空间,共享该全局变量,都可以给它曲赋值;如果函数的子函数也要使用全局变量,也必须用global声明;3.永久变量永久变量用persistent声明,只能在M文件函数中定义和使用,只允许声明它的函数存取;当声明它的函数退出时,MATLAB不会从内存中清除它,例如声明a为永久变量：>>persistent a2变量命名规则如下：始于字母,由字母、数字或下划线组成；区分大小写；可任意长,但使用前N个字符;N与硬件有关,由函数namelengthmax返回,一般N=63；不能使用关键字作为变量名关键字在后面给出；避免使用函数名作为变量名;如果变量采用函数名,该函数失效;3区分大小写②试说明分号、逗号、冒号的用法;答：分号：加上分号“;”其作用是将计算机结果存入内存,但不显示在屏幕上,反之,语句结尾若不加“;”,则表示在语句执行后,在将计算结果存入内存的同时,还将运算结果显示出来;逗号：分割列表冒号：从什么到什么,比如1:10意思是——1到10③linspace 称为“线性等分”函数,说明它的用法;可使用help命令,格式如下：>>help linspace④说明函数ones 、zeros 、eye 的用法;答；1ones函数：全部元素都为1的常数矩阵； 2zeros函数：全部元素都为0的矩阵； 3eye函数：单位矩阵；4linspace函数：如a=linspacen1,n2,n3,表示在线性空间上,行矢量的值从n1到n22 工作空间窗口Workspace;单击工作空间窗口右上角的按钮,将其从MATLAB 主界面分离出来;①在工作空间查看各个变量,或在命令窗口用who, whos注意大小写查看各个变量;②在工作空间双击变量,弹出Array Editor 窗口数组编辑器窗口,即可修改变量;③使用save 命令把工作空间的全部变量保存为文件;>>save④输入下列命令：>>clear all %清除工作空间的所有变量观察工作空间的变量是否被清空;使用load 命令把刚才保存的变量载入工作空间;>>load⑤清除命令窗口命令：>>clc3 历史命令窗口Command History;打开历史命令窗口,可以看到每次运行MATLAB 的时间和曾在命令窗口输入过的命令,练习以下几种利用历史命令窗口重复执行输入过的命令的方法;①在历史命令窗口中选中要重复执行的一行或几行命令,右击,出现快捷菜单,选择Copy,然后再Paste 到命令窗口;②在历史命令窗口中双击要执行的一行命令,或者选中要重复执行的一行或几行命令后,用鼠标将其拖动到命令窗口中执行;③在历史命令窗口中选中要重复执行的一行或几行命令,右击,出现快捷菜单,选择Evaluate Selection,也可以执行;④或者在命令窗口使用方向键的上下键得到以前输入的命令;例如,按方向键“↑”一次,就重新将用户最后一次输入的命令调到MATLAB 提示符下;重复地按方向上键“↑”,就会在每次按下的时候调用再往前一次输入的命令;类似地,按方向键“↓”的时候,就往后调用一次输入的命令;按方向键“←”或者方向键“→”就会在提示符的命令中左右移动光标,这样用户就可以用类似于在字处理软件中编辑文本的方法编辑这些命令;4 当前目录命令窗口Current Directory;MATLAB 的当前目录即是系统默认的实施打开、装载、编辑和保存文件等操作时的文件夹;打开当前目录窗口后,可以看到用“save”命令所保存的文件是保存在目录C:\MATLAB6p5\work 下;5 帮助窗口Help Window;单击工具栏的图标,或选择菜单View|Help,或选择菜单Help|MATLAB Help 都能启动帮助窗口;①通过Index 选项卡查找log2函数的用法,在Search index for 栏中输入需要查找的词汇“log2”,在左下侧就列出与之最匹配的词汇条目,选择“log21”,右侧的窗口就会显示相应的内容;②也可以通过Search 选项卡查找log2 函数的用法;Search 选项卡与Index 选项卡不同,Index 只在专用术语表中查找,而Search 搜索的是整个HTML 帮助文件;2. MATLAB 的数值显示格式设置屏幕显示方式有紧凑Compact和松散Loose两种,其中Loose 为默认方式;>>a=ones1,30>>format compact>>a数字显示格式有short、long、short e、long e 等,请参照教材的列表练习一遍;>>format long>>pi>>format short>>pi>>format long>>pi>>format +>>pi>>-pi3. 变量的搜索顺序在命令窗口中输入以下指令：>>pi>>sinpi;>>exist'pi'>>pi=0;>>exist'pi'>>pi>>clear pi>>exist'pi'>>pi思考题：① 3 次执行exist'pi'的结果一样吗如果不一样,试解释为什么答：不一样,pi原来是库存函数,但是如果被赋值则系统默认被赋予的值为pi后来的值,但是当执行clear pi之后所赋的值被清空,因此pi的值又成为②圆周率pi 是系统的默认常量,为什么会被改变为0答：pi原来是库存函数,但是如果被赋值则系统默认被赋予的值为pi后来的值,但是当执行clear pi之后所赋的值被清空,因此pi的值又成为实验二 MATLAB语言基础一、实验目的基本掌握MATLAB向量、矩阵、数组的生成及其基本运算区分数组运算和矩阵运算、常用的数学函数;交接字符串的操作;二、实验内容1 向量的生成和运算;2 矩阵的创建、引用和运算;3 多维数组的创建及运算;4 字符串的操作;三、实验步骤1.向量的生成和运算1 向量的生成直接输入法：A =2 3 4 5 6>> B=1;2;3;4;5B =12345冒号生成发：>> A=1:2:10 ,B=1:10 ,C=10:-1:1A =1 3 5 7 9B =1 2 3 4 5 6 7 8 9 10C =10 9 8 7 6 5 4 3 2 1函数法：Linspace 是线性等分函数,logspace 是对数等分函数;>> A=linspace1,10 ,B=linspace1,30,10A =Columns 1 through 9Columns 10 through 18Columns 19 through 27Columns 28 through 36Columns 37 through 45Columns 46 through 54Columns 55 through 63Columns 64 through 72Columns 73 through 81Columns 82 through 90Columns 91 through 99Column 100B =Columns 1 through 9Column 10>> A=logspace0,4,5A =1 10 100 1000 10000练习：使用logspace 创建1-4 的有10个元素的行向量;答案：>> A=logspace1,10,4piA =+010Columns 1 through 9Columns 10 through 122 向量的运算维数相同的行向量只见可以相加减,维数相同的列向量也可以相加减,标量可以与向量直接相乘除;>> A=1 2 3 4 5, B=3:7,A =1 2 3 4 5B =3 4 5 6 7>> AT=A', BT=B',AT =12345BT =34567>> E1=A+B, E2=A-BE1 =4 6 8 10 12E2 =-2 -2 -2 -2 -2 >> F=AT-BT,F =-2-2-2-2-2>> G1=3A, G2=B/3,G1 =3 6 9 12 15G2 =向量的点积与叉积运算;>> A=ones1,10;B=1:10; BT=B';>> E1=dotA,BE1 =55>> E2=ABTE2 =55>> clear>> A=1:3,B=3:5,A =1 2 3B =3 4 5>> E=crossA,BE =-2 4 -22.矩阵的创建、引用和运算1 矩阵的创建和引用矩阵是由nm 元素构成的矩阵结构,行向量和列向量是矩阵的特殊形式;直接输入法：>> A=1 2 3;4 5 6A =1 2 3 4 5 6 >> B= 1 4 72 5 83 6 9B =1 4 72 5 83 6 9 >> A1ans =1>> A4:endans =5 36 >> B:,1ans =123>> B:ans =123456789>> B5ans =5抽取法>> clear>> A=1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16A =1 2 3 45 6 7 89 10 11 1213 14 15 16>> B=A1:3,2:3B =2 36 710 11>> C=A1 3,2 4C =2 410 12>> A1 3;2 4ans =1 95 13函数法：>> A=ones3,4A =1 1 1 11 1 1 11 1 1 1>> B=zero3Undefined function or method 'zero' for input arguments of type 'double'. >> B=zeros3B =0 0 00 0 00 0 0>> C=eyes3,2Undefined function or method 'eyes' for input arguments of type 'double'. >> C=eye3,2C =1 00 10 0>> D=magic3D =8 1 63 5 74 9 2拼接法>> clear>> A=ones3,4A =1 1 1 11 1 1 1>> B=zeros3B =0 0 00 0 00 0 0>> C=eye4C =1 0 0 00 1 0 00 0 1 00 0 0 1>> D=A BD =1 1 1 1 0 0 0 1 1 1 1 0 0 0 1 1 1 1 0 0 0 >> F=A;CF =1 1 1 11 1 1 11 1 1 10 1 0 0 0 0 1 0 0 0 0 1 拼接函数和变形函数法：>> clear>> A=0 1;1 1A =0 11 1>> B=2ones2B =2 22 2>> cat 1,A,B,Aans =0 11 12 22 20 11 1>> cat 2 A,B,Acat 2 A,B,AError: Unexpected MATLAB expression.>> cat2,A,B,Aans =0 1 2 2 0 11 12 2 1 1>> repmatA,2,2ans =0 1 0 11 1 1 10 1 0 11 1 1 1>> repmatA,2ans =0 1 0 11 1 1 10 1 0 11 1 1 1练习：使用函数法、拼接法、拼接函数法和变形函数法,按照要求创建以下矩阵：A为43⨯的3⨯的单位矩阵、D为33⨯的0矩阵/C为33⨯的全1矩阵、B为3魔方阵、E由C 和D纵向拼接而成,F抽取E的2---5行元素生成、G由Ｆ经变形为46⨯的大3⨯的矩阵而得、以Ｇ为子矩阵用复制函数repmat生成8矩阵H;答案：>> A=ones3,4A =1 1 1 11 1 1 11 1 1 1>> B=zeros3,3B =0 0 00 0 00 0 0>> C=eye3C =1 0 00 1 00 0 1>> D=magic3D =8 1 63 5 74 9 2>> E=C;DE =1 0 00 1 00 0 18 1 63 5 74 9 2>> F=2:5,:F=2:5,:Error: Expression or statement is incorrect--possibly unbalanced , {, or . >> F=E2:5,:F =0 1 00 0 18 1 63 5 7>> G=respaceE,3,4Undefined command/function 'respace'.>> G=respaceF,3,4Undefined command/function 'respace'.>> G=reshapeF,3,4G =0 3 1 10 1 5 68 0 0 7>> H=repmatG,2,2H =0 3 1 1 0 3 1 1 0 1 5 6 0 1 5 6 8 0 0 7 8 0 0 7 0 3 1 1 0 3 1 1 0 1 5 6 0 1 5 6 8 0 0 7 8 0 0 7 2)矩阵的运算矩阵的加减、数乘与乘法已知矩阵：>> A=1 23 -1,A =1 23 -1>> B=-1 01 2B =-1 01 2>> A+Bans =0 24 1>> 2Aans =2 46 -2>> 2A-3Bans =5 43 -8>> ABans =1 4-4 -2矩阵的逆矩阵>> format rat;A=1 0 1;2 1 2;0 4 61 0 12 1 2 0 4 6 >> A1=invAA1 =-1/3 2/3 -1/6 -2 1 0 4/3 -2/3 1/6 >> AA1ans =1 0 0 0 1 0 0 0 1 矩阵的除法>> a=1 2 1;3 1 4;2 2 1,b=1 1 2,d=b'a =1 2 1 3 1 42 2 1b =1 1 2112>> c1=binva,c2=b/ac1 =6/7 3/7 -4/7 c2 =6/7 3/7 -4/7 >> c3=invad, c4=a\bc3 =12/7-4/7Error using ==> mldivideMatrix dimensions must agree.>> c3=invad, c4=a\dc3 =12/7-4/7c4 =12/7-4/7练习：按下列要求求出各种的矩阵运算的值求矩阵的秩、特征值和特征向量、矩阵的乘幂与开方；矩阵的指数与对数矩阵的提取与翻转答案：>> A=6 3 4 3-2 5 7 -48 -1 -3 -7A =6 3 4 3-2 5 7 -48 -1 -3 -7>> B=rankAB =3>> rb=rankArb =3>> X,Lamda=eigAError using ==> eigMatrix must be square.>> X,Lamda=eigsAError using ==> eigsA must be a square matrix or a function which computes Ax. >> C=6 3 4-2 5 78 -1 -3C =6 3 4-2 5 78 -1 -3>> X,Lamda=eigsCX =Lamda =0 00 00 0>> X,Lamda=eigCX =Lamda =0 0 0 0 0 0 >> X,Lamda=eigCX =Lamda =0 0 0 0 0 0 >> D=A^2Error using ==> mpower Matrix must be square. >> D=C^2D =34 12 626 22 34>> E=sqrtmCE =+ - -+ - -- + +>> F=expmCF =+004>> G=logmCWarning: Principal matrix logarithm is not defined for A with nonpositive real eigenvalues. A non-principal matrix logarithm is returned.> In funm at 153In logm at 27G =+ - -- + + >> H=fliplrCH =4 3 6 75 -2 -3 -1 8 >> I=triuCI =6 3 4 0 57 0 0 -3 >> J=trilCJ =6 0 0 -2 5 0 8 -1 -3 >> K=diagCK =65-33.多维数组的创建及运算1）多维数组的创建>> A1=1,2,3;4 5 6;7,8,9;A2=reshape10:18,3,3 A2 =10 13 1611 14 1712 15 18>> T1:,:,1=ones3;T1:,:,2=zeros3T1:,:,1 =1 1 11 1 11 1 1T1:,:,2 =0 0 00 0 00 0 0>> T2=ones3,3,2T2:,:,1 =1 1 11 1 11 1 1T2:,:,2 =1 1 11 1 11 1 1>> T3=cat3,A1,A2,T4=repmatA1,1,1,2T3:,:,1 =1 2 34 5 67 8 9T3:,:,2 =10 13 1611 14 1712 15 18T4:,:,1 =1 2 34 5 67 8 9T4:,:,2 =1 2 34 5 67 8 92）多维数组的创建数组运算用小圆点加在运算符的前面表示,以区分矩阵的运算;特点是两个数组相对应的元素进行运算;>> A=1:6;B=ones1,6;>> C1=A+B,C2=A-BC1 =2 3 4 5 6 7C2 =0 1 2 3 4 5 >> C3=A.B,C4=B./A,C5=A.\BC3 =1 2 3 4 5 6C4 =C5 =关系运算或逻辑运算的结果都是逻辑值;>> I=A>3,C6=AII =0 0 0 1 1 1C6 =4 5 6>> A1=A-3,I2=A1&AA1 =-2 -1 0 1 2 3 I2 =1 1 0 1 1 1>> I3=~II3 =1 1 1 0 0 0 4.字符串的操作1 字符串的创建>> S1='Ilike MATLAB'S1 =Ilike MATLAB>> S2='I''m a stuent.'S2 =I'm a stuent.>> S3=S2,'and',S1S3 =I'm a MATLAB2）求字符串长度>> lengthS1ans =12>> sizeS1ans =1 123）字符串与一维数值数组的相互转换>> CS1=absS1CS1 =73 108 105 107 101 32 77 65 84 76 65 66>> CS2=doubleS1CS2 =73 108 105 107 101 32 77 65 84 76 65 66>> charCS2ans =Ilike MATLAB>> setstrCS2ans =Ilike MATLAB练习：用char 和向量生成的方法创建如下字符串AaBbCc......XxYyZz.>> S1=65:90;S2=97:122;>> C=S1;S2;>> C=C:';>> S3=doubleC;charS3ans =AaBbCcDdEeFfGgHhIiJjKkLlMmNnOoPpRrSsTtUuVvWwXxYyZz实验三 MALTAB数值运算一、实验目的掌握MATLAB的数值运算及其运算中所用到的函数,掌握结构数组和细胞数组的操作;二、实验内容1）多项式运算;2）多项式插值和拟合;3）数值为积分;4）结构数组和细胞数组;三、实验步骤1、多项式运算1）多项式表示;在MATLAB中,多项式表示成向量的形式;如：432S S S+-+在MATLAB中表示为359>> s=1 3 -5 0 92）多项式的加减法相当于向量的加减法,但必须注意阶次要相同;如不同,低阶次的要补0;如多项式23547+-++相加;s s s s++与多项式432239s s>> s1=0 0 2 3 11>> s2=1 2 -5 4 7>> s3=s1+s2答；s1 =0 0 2 3 11s2 =1 2 -5 4 7s3 =1 2 -3 7 183）多项式的乘、除法分别用函数conv和的deconv实现; >> s1=2 3 11>> s2=1 3 -5 4 7>> s3=convs1,s2>> s4=deconvs3,s1答；s1 =2 3 11s2 =1 3 -5 4 7s3 =2 9 10 26 -29 65 77s4 =1 3 -5 4 74）多项式求根用函数roots;>> s1=2 4 2>> rootss1答；s1 =2 4 2ans =-1-15）多项式求值用函数polyval>> s1=2 4 1 -3>> polyvals1,3>> x=1:10>> y=polyvals1,x答；s1 =2 4 1 -3ans =90x =1 2 3 4 5 6 7 8 9 10y =Columns 1 through 84 31 90 193 352 579 886 1285 Columns 9 through 10练习：求23(1)(3)(1)21s s ss s+++++的“商”及余数;>> s1=1 0 1;s2=1 3;s3=1 1;>> s4=1 0 2 1;>> q,r=deconvconvconvs1,s2,s3,s4答；q =1 4r =0 0 2 -5 -12、多项式插值和拟合有一组实验数据如附表1-1所示;附表1-1分别用拟合二阶至三阶和插值线性和三次样条的方法来估算X=时Y的值;以下是实现一阶拟合的语句;>> x=1:10>> y=16 32 70 142 260 436 682 1010 1342 1960>> p1=polyfitx,y,1 %一阶拟合>> y1=polyvalp1, %计算多项式P1在x=的值答；x =1 2 3 4 5 6 7 8 9 10y =Columns 1 through 816 32 70 142 260 436 682 1010Columns 9 through 101342 1960p1 =y1 =+0033、数值微积分1）差分使用diff函数实现;>> x=1:2:9>> diffx答；x =1 3 5 7 9ans =2 2 2 22）可以用因变量和自变量差分的结果相处得到数值微分;>> x=linspace0,2pi,100;>> y=sinx;>> plotx,y>> y1=diffy./diffx;>> plotx1:end-1,y1答；3）cumsum函数求累计积分,trapz函数用梯形法求定积分,即曲线的面积;>> x=ones1,10>> cumsumx>> x=linspace0,pi,100;>> y=sinx;>> trapzx,y>> p=cumsumy;>> p100pi/100-1答；x =1 1 1 1 1 1 1 1 1 1ans =1 2 3 4 5 6 7 8 9 10ans =ans =练习:图A1是瑞士地图,为了算出其国土面积,首先对地图作如下测量：以由西向东方向为X轴,由南向北方向为Y 轴,选择方便的原点,并将从最西边界点到最东边界点在X轴上的区间适当划分为若干段,在每个分点的Y方向测出南边界点和北边界点的Y坐标Y1和Y2,这样就得到表1,根据地图比例尺回到18mm相当于40Km,试由测量数据计算瑞士国km比较;地图的数据见附表1-2单位mm;土近似面积,与其精确值412282附表1-2X 7 13 34 48 56 61 91Y1 44 45 47 50 50 38 30 30 34 36 34 41 45 46续表提示：由高等数学的知识,一条曲线的定积分是它与x轴所围成的面积,那么两条曲线所围成的面积可由两条曲线的定积分相减得到;4、结构数组与细胞数组1）机构数组的创建;>> ='Jack';>> student2.n>> ='Lucy';>> student或者用struct函数创建;>> student=struct'number',{'001','002'},'name',{'Jack','Lucy'};答；student =1x2 struct array with fields:numbername2）机构数组的操作;>> student1.subject= %添加subject域并赋予空值>> student1.score=>> filednamesstudeng>>fieldnamesstudeng>>fieldnamesstudent>> getfieldstudent,{2},'name'>> student=rmfieldstudent,'subject' %删除subject域>> student=setfieldstudent,{1},'score',90;>> student2.score=88; %比较和上一条语句是否效果一样答；student =1x2 struct array with fields:numbernamesubjectstudent =1x2 struct array with fields:numbernamesubjectscoreUndefined function or variable 'studeng'.练习：创建一结构数组stusorce,其域为：No,Name,English,Math,Chinese,Total,Average;结构数组的大小为2×2;3）细胞数组的创建;>> A={'How are you',ones3;1 2;3 4,{'cell'}}; %直接创建>> B1,1={'Hello world'}; %由各个细胞元素创建>> B1,2={magic3};>> B2,1={1 2 3 4};答或者用cell函数先创建空的细胞数组,然后再给各个元素赋值c=cell1,2;>> c1,1={'Hello world'};>> c1,2={magic3};>> c1,3={1 2 3 4};4）细胞数组的操作;>> ans1=A1,1>> ans2=A1,1>> whos ans1 ans2>> celldispA>> a1=A{2,1}1,2>> a2 a3=dealA{1:2}答；ans1 ='How are you'ans2 ='How are you'Name Size Bytes Class Attributesans1 1x1 84 cellans2 1x1 84 cellA{1,1} =How are youA{2,1} =1 23 4A{1,2} =1 1 11 1 11 1 1A{2,2}{1} =cella1 =2a2 =How are youa3 =1 23 4实验四 MALTAB符号运算一、实验目的掌握符号变量和符号表达式的创建,掌握MALTAB的symbol工具箱的一些基本运用;二、实验内容1）符号变量、表达式、方程及函数的表示; 2）符号微积分运算;3）符号表达式的操作和转换;4）符号微分方程求解;三、实验步骤1、符号运算的引入在数值运算中如果求x xxπsinlim→,则可以不断让的让x趋近0,一球的表达式趋近什么数,但终究不能令x=0,因为在数值运算中0不是能作除数的;MATLAB的符号运算能解决这内问题;输入如下命令：>> f=sym'sinpix/x'>> limitf,'x',0答；f =sinpix/xans =pi2、符号常量、符号变量、符号表达式的创建1）使用sym创建输入以下命令,观察Workspace中A、B、f是什么内性的数据,占用多少字节的内存空间;>> A=sym'1' %符号常量>> B=sym'x' %符号变量>> f=sym'2x^2+3x-1' %符号表达式>> clear>> f1=sym'1+2' %有单引号,表示字符串>> f2=sym1+2 %无单引号>> f2=sym1+2>> f4=sym'2x+3' %为什么出错>> x=1>> f4=sym2x+3答；A =1B =x2x^2+3x-1f1 =1+2f2 =3f2 =3f4 =2x+3x =1f4 =5通过看MATLAB的帮助可知,sym的参数可以使字符串或是数值类型,无论是哪种类型都会生成符号类型数据; 2）使用syms创建>> clear>> syms x y z>> x,y,z>> f1=x^2+2x+1>> f2=expy+expz^2>> f3=f1+f2答；x =xy =yz =zf1 =x^2+2x+1f2 =expy+expz^2f3 =x^2+2x+1+expy+expz^23、符号矩阵创建>> syms a1 a2 a3 a4>> A=a1 a2;a3 a4>> A1,A3答；A =a1, a2a3, a4ans =a1a24、符号算术运算1）符号向量相乘、相除符号量相成和数值量相乘一样,分成矩阵乘和数组乘;>> a=sym5;b=sym7;>> c1=ab>> c2=a/b>> a=sym5;B=sym3 4 5;>> C1=aB,C2=a\B>> syms a b>> A=5 a;b 3;B=2a b;2b a;>> C1=AB,C2=A.B>> C3=A\B,C4=A./B答；c1 =35c2 =5/7C1 =15, 20, 25C2 =3/5, 4/5, 1C1 =10a+2ab, 5b+a^22ab+6b, b^2+3aC2 =10a, ab2b^2, 3aC3 =2ab-3/-15+ab, a^2-3b/-15+ab2ba-5/-15+ab, -5a-b^2/-15+abC4 =5/2/a, a/b1/2, 3/a2）符号数值任意精度控制和运算任意精度的VPA运算可以使用命令digits设定默认的精度和vpa对指定对象以新的精度进行计算来实现; >> a1=sym'2sqrt5+pi'>> a=sym'2sqrt5+pi'>> b=sym2sqrt5+pi>> digits>> vpaa>> digits15>> c1=vpaa,56>> c2=vpab,56答a1 =2sqrt5+pia =2sqrt5+pib =Digits = 32ans =ans =c1 =c2 =注意观察c1和c2的数值类型,c1和c2是否相等;3） 符号类型与数值类型的转换使用命令sym 可以把数值型对象转换成有理数性符号对象,命令vpa 可以讲数值型对象转换为任意精度的VPA 型符号对象;使用double,numeric 函数可以将有理数型和VPA 型符号对象转换成数值对象.>> clear>> a1=sym'2sqrt5+pi'>> b1=doublea1 %符号转数值>> b2=isnumericb1 %判断是否转换成了数值>> a2=vpaa1,70 %数值转符号答；a1 =2sqrt5+pib1 =b2 =1a2 =85、 符号表达式的操作和转换1独立变量的确定原则独立变量的确定原则：在符号表达式中默认变量是惟一的;MATLAB 会对单个英文小写字母除i 、j 外进行搜索,且以x 为首选独立变量;如果表达式中字母不唯一,且无x,就选在字母表最接近x 的字母;如果有相连的字母,则选择在字母表中较后的那一个;例如：'*3'z y +中,y 是默认独立变量;)'*sin('b t a +,t 是默认独立变量;输入以下命令,观察并分析结果;>> clear>> f=sym'a+b+i+j+x+y+xz'>> findsymf>> findsymf,1>> findsymf,2>> findsymf,3>> findsymf,4>> findsymf,5>> findsymf,6答；f =a+b+i+j+x+y+xzans =a, b, j, x, xz, yans =xans =x,xzans =x,xz,yans =x,xz,y,jans =x,xz,y,j,bans =x,xz,y,j,b,a2符号表达式的化简符号表达式化简主要包括表达式美化pretty 、合并同类项collcet 、多项式展开expand 、因式分解factor 、化简simple 或simplify 等函数;①合并同类项collect;分别按x 的同幂项和e 指数同幂项合并表达式： 2(1)()t t x xe x e --+++;>> syms x t;>> f=x^2+x^exp-t+1x+exp-t;>> f1=collectf>> f2=collectf,'exp-t'答；f1 =x^3+exp-tx^2+x^exp-t+1x+x^exp-t+1exp-tf2 =x^2+x^exp-t+1exp-t+x^2+x^exp-t+1x②对显示格式加以美化pretty;针对上例,用格式美化函数可以使显示的格式更符合数学书写习惯;>> prettyf1>> prettyf2答；f1 =x^3+exp-tx^2+x^exp-t+1x+x^exp-t+1exp-tf2 =x^2+x^exp-t+1exp-t+x^2+x^exp-t+1x>> prettyf1prettyf23 2 exp-t exp-tx + exp-t x + x + 1 x + x + 1 exp-t2 exp-t 2 exp-tx + x + 1 exp-t + x + x + 1 x注意与直接输出的f1和f2对比;③多项式展开expand;展开12x-成x不同次幂的多项式.(1)>> syms x>> f=x-1^12;>> expandf>> prettyexpandf答；ans =1+x^12-12x^11+66x^10-220x^9+495x^8-792x^7+924x^6-792x^5+495x^4-220x^3+66x^2-12x 12 11 10 9 8 7 6 51 + x - 12 x + 66 x - 220 x + 495 x - 792 x + 924 x - 792 x4 3 2+ 495 x - 220 x + 66 x - 12 x④因式分解factor;将表达式121x-做因式分解;>> syms x;f=x^12-1;>> prettyfactorf答；ans =1+x^12-12x^11+66x^10-220x^9+495x^8-792x^7+924x^6-792x^5+495x^4-220x^3+66x^2-12x 12 11 10 9 8 7 6 51 + x - 12 x + 66 x - 220 x + 495 x - 792 x + 924 x - 792 x4 3 2+ 495 x - 220 x + 66 x - 12 x>> syms x;f=x^12-1;prettyfactorf2 2 2 4 2x - 1 1 + x + x 1 + x 1 - x + x 1 + x x - x + 1⑤化简simple或simplify;将函数f=;>> clear>> syms x;f=1/x^3+6/x^2+12/x+8^1/3;>> g1=simplef>> g2=simplifyf答；g1 =2x+1/xg2 =2x+1^3/x^3^1/36、符号表达式的变量替换subs函数可以对符号表达式中的符号变量进行替换>> clear>> f=sym'x+y^2+4x+10'>> f1=subsf,'x','s' %使用s替换x>> f2=subsf,'x+y','z'答；f =x+y^2+4x+10f1 =s+y^2+4s+10f2 =z^2+4x+107、符号极限、符号积分与微分1）求极限函数的调用格式limitF,x,a %返回符号对象F当x→a时的极限limitF,a %返回符号对象F当独立变量→a时的极限limitF %返回符号对象F当独立变量→0a=0时的极限limitF,x,a,’right’ %返回符号对象F当x→a时的右极限limitF,x,a,’left’ %返回符号对象F当x→a时的左极限例一：>> clear>> f=sym'sinx/x+ax'>> limitf,'x',0 %以x为自变量求极限>> limitf,'a',0 %以a为自变量求极限>> limitf %在默认情况下以x为自变量求极限>> findsymf %得到变量并且按字母表顺序排列答‘f =sinx/x+axans =1ans =sinx/xans =1ans =a, x例二：>> clear>> f=sym'sqrt1+1/n;>> limitf,n,inf %求n趋于无穷大时的极限2）求积分函数的调用格式intF %求符号对象F关于默认变量的不定积分intF,v %求符号对象F关于指定变量v的不定积分intF,a,b %求符号对象F关于默认变量的从a到b的定积分intF,v,a,b %求符号对象F关于指定变量的从a到b的定积分3）求微分方程的调用格式diffF %求符号对象F关于默认变量的微分diffF,v %求符号对象F关于指定变量v的微分diffF,n %求符号对象F关于默认变量的n阶微分,n为自然数1、2、3……diffF,v,n %求符号对象F关于指定变量v的n阶微分8、符号方程求解1常规方程求解函数的调用格式g=solveeq %求方程或表达式或字串eq关于默认变量的解g=solveeq,var %求方程或表达式或字串eq关于指定变量var的解g=solveeq1,eq2,….,eqn,var1,var2,…,varn %求方程或表达式或字串eq1,eq2,eq3,……eqn关于指定变量组var1,var2,……,varn的解求一元二次方程20++=的解;其求解方法有多种形式：ax bx c①seq=solve'ax^2+bx+c'②seq=solve'ax^2+bx+c=0'③eq=’ax^2+bx+c’;④eq=’ax^2+bx+c=0’;seq=solveeq⑤sym x a b ceq=ax^2+bx+cseq=solveeq2常微分方程求解求解常微分方程的函数是dsolve;应用此函数可以求得常微分方程组的通解,以及给定边界条件或初始条件后的特解;常微分方程求解函数的调用格式：r=dsolve’eq1,eq2,…’,’cond1,cond2,…’,’v’r=dsolve’eq1’,’eq2’,…,’cond1’,’cond2’,…,’v’说明：①以上两式均可给出方程eq1,、qeq2对应初始条件cond1、cond2之下的一v 作为解变量的各微分方程的解; ②常微分方程解的默认变量为t;③第二式中最多可接受的输入式是12个;④微分方程的表达方法;在用MATLAB 求解常微分方程时,用大写字母Dy 表示dxdy ,用D2y 表示22d y dx ,依此类推; 边界条件以类似于ya=b 给出;其中y 为因变量,a 、b 为常数.如果初始条件给的不够,求出的解为含有C1、C2等待定常数的通解;例一 求微分方程x dx dy 2=的通解. 练习：1求222132lim x x x x →--+; 2求函数()cos 2sin 2f x x x =-的积分；求函数()g x =3计算定积分60(sin 2)x dx π+⎰ 4求下列线性方程组的解5求解但y0=2,在z0=7时,微分方程组的解;实验五 MATLAB 程序设计一、实验目的掌握MATLAB 程序设计的主要方法,熟练编写MATLAB 函数.二、 实验内容1M 文件的编辑;。
MATLAB上机实验实验报告
MATLAB上机实验实验报告实验名称：用MATLAB实现多项式拟合及插值一、实验目的：通过使用MATLAB实现多项式拟合及插值的方法，掌握MATLAB软件的基本操作和函数应用，进一步了解多项式拟合及插值的原理和实现过程。
二、实验原理：多项式拟合及插值是一种常见的数值分析方法，通过对已知数据点集合的拟合或插值，构造出一个多项式函数，用于近似表示原始数据。
1.多项式拟合：通过最小二乘法原理，选择一个合适的多项式函数，使得拟合出的多项式与已知数据点之间的误差最小。
拟合函数可以是一次、二次或高阶多项式。
2.多项式插值：通过已知数据点的横纵坐标值，构造一个满足这些点的多项式函数。
插值函数可以是一次、二次或高阶多项式。
插值函数经过每个已知数据点。
三、实验步骤：1.数据准备：选择一组已知数据，包含横纵坐标值。
数据点的个数可以根据具体情况自行确定。
2.多项式拟合：使用MATLAB中的polyfit函数，根据已知数据点进行多项式拟合。
根据拟合结果，获取拟合的多项式系数。
3.多项式插值：使用MATLAB中的polyfit函数，根据已知数据点进行多项式插值。
通过plot函数绘制原始数据点的散点图和插值多项式的曲线图。
可以尝试不同阶数的多项式插值。
4.结果分析：根据实验结果，分析拟合与插值的效果。
对比拟合结果与原始数据的误差大小，评估拟合的准确性。
对比插值结果与原始数据的差异，评估插值的精确度。
五、实验总结：通过这次实验，我熟练掌握了使用MATLAB实现多项式拟合及插值的方法。
在实验中，我了解了多项式拟合的原理，以及如何利用最小二乘法求取多项式拟合的系数。
同时，我也学会了如何使用MATLAB中的polyfit函数实现多项式拟合和插值。
通过实验结果的分析，我对拟合和插值的实际应用和效果有了更加深入的认识。
[1]MATLAB官方文档[2]高等数值分析教程以上为MATLAB上机实验实验报告，共计1200字。
MATLAB技术上机实践报告
目录上机实践一 MATLAB在“高等数学”中的应用练习（4学时）(1)上机练习—-MATLAB解常微分方程(SJ101.m）（2）上机练习-—MATLAB二元函数作图（SJ102.m）(3）上机编程——微分运算和积分运算（SJ103.m）（4)上机编程——线性方程组的数值解法(SJ104。
m）上机实践二 MATLAB在“力学"中的应用练习（4学时）(1）上机练习—-质点直线运动的位置、速度和加速度(SJ201。
m）(2)上机练习-—轻质杆的斜抛运动(SJ202。
m）（3）上机编程—-质点的抛体运动(SJ203。
m）(4）上机编程-—物体在恒力作用下的直线运动（SJ204。
m）上机实践三 MATLAB在“电磁学”中的应用练习（4学时）（1）上机练习——电磁波的发射方法一（SJ301.m）(2）上机练习—-电磁波的发射方法二(SJ302。
m）（3）上机编程——静电场场强分布计算（SJ303.m）（4）上机编程——带电粒子在电场中的运动（SJ304。
m）上机实践四 MATLAB在“振动和波动”中的应用练习（4学时）（1）上机练习——旋转矢量法(SJ401。
m）（2)上机练习—-用惠更斯作图法确定波阵面（SJ402.m）（3)上机编程——同方向简谐振动的合成（SJ403.m）（4)上机编程——简谐波横波的演示（SJ404。
m）上机实践五 MATLAB在“量子物理"中的应用练习（4学时）（1)上机练习-—黑体辐射随波长的变化规律(SJ501。
m）（2）上机练习——势垒和隧道效应(SJ502。
m）（3)上机编程——双缝干涉图样的动画模拟（SJ503.m）上机实践六 MATLAB在“电子信息技术"中的应用练习(4学时）（1)上机练习--含受控源的电阻电路计算（SJ601.m）（2）上机练习——方波分解为多次正弦波之和（SJ602。
m）（3)上机编程——电阻电路的计算（SJ603。
matlab上机实验报告
matlab上机实验报告实验1MATLAB初步1. 举例说明format命令的功能；Format命令⽤来调整数据的显⽰格式。
>> format long %15位定点表⽰，若改成short则为5位定点表⽰>> 2^ans =2. 说明clc,clf和clear命令的区别；clc即clear command(清屏)clear清理内存所有变量clf即clear figure(清理图形窗⼝3. 举例说明help命令的功能；在help命令后加⼊所要查询的函数或类名，将会返回所查询的对象的使⽤说明。
例如:>> help cosCOS Cosine.COS(X) is the cosine of the elements of X.See also acos, cosd.Overloaded functions or methods (ones with the same name in other directories)help sym/Reference page in Help browserdoc cos4. 如何观察和调⽤历史命令；在MTALAB界⾯的左下⾓有个Command History窗⼝，⾥⾯可以看到历史命令，如果双击该条历史命令即可实现调⽤。
5. 如何设置当前路径。
在MATLAB界⾯上⽅的⼯具栏有个Current Folder(当前⽬录)，更改该⽬录就可以设置当前路径。
实验2矩阵⽣成与元素标识1. 矩阵的创建=32943253036120031A ，=757065605550454035302520151051B （尽可能⽤最简单⽅法）；创建矩阵A ： >> A=[1 3 0 0;2 61 3 0;3 5 2 3;4 9 2 3]创建矩阵B ： >>B=[1 5 10 15;20 25 30 35;40 45 50 55;60 65 70 75]2. ⽣成⼀个单位矩阵C ，其结构与B 相同；>> eye(size(B))3.找出矩阵A 中所有⼤于8的元素位置（⾏号和列号）；>> C=find(B>35 & B<65)>> B(C)=B(C)+95. 提取矩阵B 的主对⾓组元素。
MATLAB上机实验报告
Matlab上机实验报告实验二读入MATLAB下自带图像pout.tif1）利用亮度变换函数，调整图像亮度。
a)调整范围设定[0 1]，[1 0]，观测显示效果；b)调整范围设定[0.5 0.75]，[1 0]，观测显示效果。
解：a, I=imread('pout.tif');colormap;imshow(I);j=imadjust(I,[0 1],[1 0],1.5);figure;subimage(j);b, >> I=imread('pout.tif');colormap;imshow(I);j=imadjust(I,[0.5 0.75],[1 0],1.5);figure;subimage(j);2）利用对比度拉伸函数，压缩高值灰度（c值自行设定）。
解：I=imread('pout.tif');colormap;subplot(1,2,1);imshow(I);xlabel('a)原始图像');J=double(I);J=100*log(J+1);I=uint8(J);subplot(1,2,2);subimage(J);xlabel('b)非线性变换');3）利用直方图函数，生成并绘制图像直方图。
解：I=imread('pout.tif');subplot(1,2,1);imshow(I);title('原始图像');subplot(1,2,2);imhist(I);4）利用直方图修正函数，生成均衡化后的图像直方图（n值自定设定）。
解：I=imread('pout.tif');figure(1);subplot(1,2,1);imshow(I);xlabel('a)原始图像');J=histeq(I);figure(1);subplot(1,2,2);imshow(J);xlabel('b)直方图均衡');figure(2);imhist(I,100);figure(3);imhist(J,100);实验三1.运行例3、4，显示并分析输出结果，说明逆滤波和维纳滤波的区别。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
姓名：班级：学号：台号：第一部分第二题（1）>> A=[1,2,3,4;4,3,2,1;2,3,4,1;3,2,4,1]A =1 2 3 44 3 2 12 3 4 13 24 1（2）>> B=[1+4j,2+3j,3+2j,4+1j;4+1j,3+2j,2+3j,1+4j;2+3j,3+2j,4+1j,1+4j;3+2j,2+3j,4+1j,1+4j]B =1.0000 + 4.0000i2.0000 +3.0000i 3.0000 + 2.0000i4.0000 + 1.0000i4.0000 + 1.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 1.0000 + 4.0000i2.0000 +3.0000i 3.0000 + 2.0000i4.0000 + 1.0000i 1.0000 + 4.0000i3.0000 + 2.0000i 2.0000 + 3.0000i4.0000 + 1.0000i 1.0000 + 4.0000i第三题>> A=magic(8);>> B=A(2:2:end,:)B =9 55 54 12 13 51 50 1640 26 27 37 36 30 31 3341 23 22 44 45 19 18 488 58 59 5 4 62 63 1第四题>> i=0:63;s=sum(2.^i)s =1.8447e+019第五题（1）>> t=[-1:0.001:1];>> y=sin(1./t);Warning: Divide by zero.>> plot(t,y)（2）t=[-pi:0.05:-1.8,-1.799:0.001:-1.2,-1.2:0.05:1.2,1.201:0.001:1.8,1.81:0.05:pi]; >> y=sin(tan(t))-tan(sin(t));>> plot(t,y)第六题>> xx=[-2:0.1:-1.2,-1.1:0.02:-0.9,-0.8:0.1:0.8,0.9:0.02:1.1,1.2:0.1:2]; >> yy=[-1:0.1:-0.2,-0.1:0.02:0.1,0.2:0.1:1];[x,y]=meshgrid(xx,yy); >> z=1./(sqrt((1-x).^2+y.^2))+1./(sqrt((1+x).^2+y.^2)); Warning: Divide by zero.Warning: Divide by zero.>> subplot(224),surf(x,y,z)>> subplot(221),surf(x,y,z),view(0,90)>> subplot(222),surf(x,y,z),view(90,0)>> subplot(223),surf(x,y,z),view(0,0)第七题（1）>> syms x;f=(3^x+9^x)^(1/x);l=limit(f,x,inf)l =9（2）>> syms x y;f=x*y/(sqrt(x*y+1)-1);limit(limit(f,x,0),y,0)ans =2（3）>> syms x y;f=(1-cos(x^2+y^2))*exp(x^2+y^2)/(x^2+y^2);limit(limit(f,x,0),y,0) ans =第八题>> syms t;x=log(cos(t));y=cos(t)-t*sin(t);diff(y,t)/diff(x,t)ans =-(-2*sin(t)-t*cos(t))/sin(t)*cos(t)>> f=diff(y,t,2)/diff(x,t,2);subs(f,t,sym(pi)/3)ans =3/8-1/24*pi*3^(1/2)第九题>> syms x y t>> s=int(exp(-t^2),t,0,x*y);>> x/y*diff(f,x,2)-2*diff(diff(f,x),y)+diff(f,y,2)ans =2*x^2*y^2*exp(-x^2*y^2)-2*exp(-x^2*y^2)-2*x^3*y*exp(-x^2*y^2)第十题（1）>> syms k n;symsum(1/((2*k)^2-1),k,1,inf)ans =1/2>> limit(symsum(1/((2*k)^2-1),k,1,n),n,inf)ans =1/2（2）>> limit(n*symsum(1/(n^2+k*pi),k,1,n),n,inf)ans =1第十一题（1）>> syms a t;x=a*(cos(t)+t*sin(t));y=a*(sin(t)-t*cos(t));>> f=x^2+y^2;I=int(f*sqrt(diff(x,t)^2+diff(y,t)^2),t,0,2*pi)I =2*a^2*pi^2*(a^2)^(1/2)+4*a^2*pi^4*(a^2)^(1/2)（2）>> syms x y a b c t;x=c*cos(t)/a;y=c*sin(t)/b;>> P=y*x^3+exp(y);Q=x*y^3+x*exp(y)-2*y;>> ds=[diff(x,t);diff(y,t)];I=int([P Q]*ds,t,0,pi)I =-2/15*c*(-2*c^4+15*b^4)/b^4/a第十二题>> syms a b c d e;A=vander([a b c d e])A =[ a^4, a^3, a^2, a, 1][ b^4, b^3, b^2, b, 1][ c^4, c^3, c^2, c, 1][ d^4, d^3, d^2, d, 1][ e^4, e^3, e^2, e, 1]>> det(A),simple(ans)ans =a^4*b^3*c^2*d-a^4*b^3*c^2*e-a^4*b^3*d^2*c+a^4*b^3*d^2*e+a^4*b^3*e^2*c-a^4*b^3*e^2* d-a^4*c^3*b^2*d+a^4*c^3*b^2*e+a^4*c^3*d^2*b-a^4*c^3*d^2*e-a^4*c^3*e^2*b+a^4*c^3*e^ 2*d+a^4*d^3*b^2*c-a^4*d^3*b^2*e-a^4*d^3*c^2*b+a^4*d^3*c^2*e+a^4*d^3*e^2*b-a^4*d^3* e^2*c-a^4*e^3*b^2*c+a^4*e^3*b^2*d+a^4*e^3*c^2*b-a^4*e^3*c^2*d-a^4*e^3*d^2*b+a^4*e^ 3*d^2*c-b^4*a^3*c^2*d+b^4*a^3*c^2*e+b^4*a^3*d^2*c-b^4*a^3*d^2*e-b^4*a^3*e^2*c+b^4* a^3*e^2*d+b^4*c^3*a^2*d-b^4*c^3*a^2*e-b^4*c^3*d^2*a+b^4*c^3*d^2*e+b^4*c^3*e^2*a-b^ 4*c^3*e^2*d-b^4*d^3*a^2*c+b^4*d^3*a^2*e+b^4*d^3*c^2*a-b^4*d^3*c^2*e-b^4*d^3*e^2*a+ b^4*d^3*e^2*c+b^4*e^3*a^2*c-b^4*e^3*a^2*d-b^4*e^3*c^2*a+b^4*e^3*c^2*d+b^4*e^3*d^2 *a-b^4*e^3*d^2*c+c^4*a^3*b^2*d-c^4*a^3*b^2*e-c^4*a^3*d^2*b+c^4*a^3*d^2*e+c^4*a^3*e ^2*b-c^4*a^3*e^2*d-c^4*b^3*a^2*d+c^4*b^3*a^2*e+c^4*b^3*d^2*a-c^4*b^3*d^2*e-c^4*b^3 *e^2*a+c^4*b^3*e^2*d+c^4*d^3*a^2*b-c^4*d^3*a^2*e-c^4*d^3*b^2*a+c^4*d^3*b^2*e+c^4* d^3*e^2*a-c^4*d^3*e^2*b-c^4*e^3*a^2*b+c^4*e^3*a^2*d+c^4*e^3*b^2*a-c^4*e^3*b^2*d-c^4 *e^3*d^2*a+c^4*e^3*d^2*b-d^4*a^3*b^2*c+d^4*a^3*b^2*e+d^4*a^3*c^2*b-d^4*a^3*c^2*e-d ^4*a^3*e^2*b+d^4*a^3*e^2*c+d^4*b^3*a^2*c-d^4*b^3*a^2*e-d^4*b^3*c^2*a+d^4*b^3*c^2* e+d^4*b^3*e^2*a-d^4*b^3*e^2*c-d^4*c^3*a^2*b+d^4*c^3*a^2*e+d^4*c^3*b^2*a-d^4*c^3*b^ 2*e-d^4*c^3*e^2*a+d^4*c^3*e^2*b+d^4*e^3*a^2*b-d^4*e^3*a^2*c-d^4*e^3*b^2*a+d^4*e^3* b^2*c+d^4*e^3*c^2*a-d^4*e^3*c^2*b+e^4*a^3*b^2*c-e^4*a^3*b^2*d-e^4*a^3*c^2*b+e^4*a^ 3*c^2*d+e^4*a^3*d^2*b-e^4*a^3*d^2*c-e^4*b^3*a^2*c+e^4*b^3*a^2*d+e^4*b^3*c^2*a-e^4* b^3*c^2*d-e^4*b^3*d^2*a+e^4*b^3*d^2*c+e^4*c^3*a^2*b-e^4*c^3*a^2*d-e^4*c^3*b^2*a+e^ 4*c^3*b^2*d+e^4*c^3*d^2*a-e^4*c^3*d^2*b-e^4*d^3*a^2*b+e^4*d^3*a^2*c+e^4*d^3*b^2*a-e^4*d^3*b^2*c-e^4*d^3*c^2*a+e^4*d^3*c^2*bsimplify:a^4*b^3*c^2*d-a^4*b^3*c^2*e-a^4*b^3*d^2*c+a^4*b^3*d^2*e+a^4*b^3*e^2*c-a^4*b^3*e^2* d-a^4*c^3*b^2*d+a^4*c^3*b^2*e+a^4*c^3*d^2*b-a^4*c^3*d^2*e-a^4*c^3*e^2*b+a^4*c^3*e^ 2*d+a^4*d^3*b^2*c-a^4*d^3*b^2*e-a^4*d^3*c^2*b+a^4*d^3*c^2*e+a^4*d^3*e^2*b-a^4*d^3* e^2*c-a^4*e^3*b^2*c+a^4*e^3*b^2*d+a^4*e^3*c^2*b-a^4*e^3*c^2*d-a^4*e^3*d^2*b+a^4*e^ 3*d^2*c-b^4*a^3*c^2*d+b^4*a^3*c^2*e+b^4*a^3*d^2*c-b^4*a^3*d^2*e-b^4*a^3*e^2*c+b^4* a^3*e^2*d+b^4*c^3*a^2*d-b^4*c^3*a^2*e-b^4*c^3*d^2*a+b^4*c^3*d^2*e+b^4*c^3*e^2*a-b^ 4*c^3*e^2*d-b^4*d^3*a^2*c+b^4*d^3*a^2*e+b^4*d^3*c^2*a-b^4*d^3*c^2*e-b^4*d^3*e^2*a+ b^4*d^3*e^2*c+b^4*e^3*a^2*c-b^4*e^3*a^2*d-b^4*e^3*c^2*a+b^4*e^3*c^2*d+b^4*e^3*d^2 *a-b^4*e^3*d^2*c+c^4*a^3*b^2*d-c^4*a^3*b^2*e-c^4*a^3*d^2*b+c^4*a^3*d^2*e+c^4*a^3*e ^2*b-c^4*a^3*e^2*d-c^4*b^3*a^2*d+c^4*b^3*a^2*e+c^4*b^3*d^2*a-c^4*b^3*d^2*e-c^4*b^3 *e^2*a+c^4*b^3*e^2*d+c^4*d^3*a^2*b-c^4*d^3*a^2*e-c^4*d^3*b^2*a+c^4*d^3*b^2*e+c^4* d^3*e^2*a-c^4*d^3*e^2*b-c^4*e^3*a^2*b+c^4*e^3*a^2*d+c^4*e^3*b^2*a-c^4*e^3*b^2*d-c^4 *e^3*d^2*a+c^4*e^3*d^2*b-d^4*a^3*b^2*c+d^4*a^3*b^2*e+d^4*a^3*c^2*b-d^4*a^3*c^2*e-d ^4*a^3*e^2*b+d^4*a^3*e^2*c+d^4*b^3*a^2*c-d^4*b^3*a^2*e-d^4*b^3*c^2*a+d^4*b^3*c^2* e+d^4*b^3*e^2*a-d^4*b^3*e^2*c-d^4*c^3*a^2*b+d^4*c^3*a^2*e+d^4*c^3*b^2*a-d^4*c^3*b^ 2*e-d^4*c^3*e^2*a+d^4*c^3*e^2*b+d^4*e^3*a^2*b-d^4*e^3*a^2*c-d^4*e^3*b^2*a+d^4*e^3* b^2*c+d^4*e^3*c^2*a-d^4*e^3*c^2*b+e^4*a^3*b^2*c-e^4*a^3*b^2*d-e^4*a^3*c^2*b+e^4*a^ 3*c^2*d+e^4*a^3*d^2*b-e^4*a^3*d^2*c-e^4*b^3*a^2*c+e^4*b^3*a^2*d+e^4*b^3*c^2*a-e^4* b^3*c^2*d-e^4*b^3*d^2*a+e^4*b^3*d^2*c+e^4*c^3*a^2*b-e^4*c^3*a^2*d-e^4*c^3*b^2*a+e^ 4*c^3*b^2*d+e^4*c^3*d^2*a-e^4*c^3*d^2*b-e^4*d^3*a^2*b+e^4*d^3*a^2*c+e^4*d^3*b^2*a-e^4*d^3*b^2*c-e^4*d^3*c^2*a+e^4*d^3*c^2*bradsimp:a^4*b^3*c^2*d-a^4*b^3*c^2*e-a^4*b^3*d^2*c+a^4*b^3*d^2*e+a^4*b^3*e^2*c-a^4*b^3*e^2* d-a^4*c^3*b^2*d+a^4*c^3*b^2*e+a^4*c^3*d^2*b-a^4*c^3*d^2*e-a^4*c^3*e^2*b+a^4*c^3*e^ 2*d+a^4*d^3*b^2*c-a^4*d^3*b^2*e-a^4*d^3*c^2*b+a^4*d^3*c^2*e+a^4*d^3*e^2*b-a^4*d^3* e^2*c-a^4*e^3*b^2*c+a^4*e^3*b^2*d+a^4*e^3*c^2*b-a^4*e^3*c^2*d-a^4*e^3*d^2*b+a^4*e^ 3*d^2*c-b^4*a^3*c^2*d+b^4*a^3*c^2*e+b^4*a^3*d^2*c-b^4*a^3*d^2*e-b^4*a^3*e^2*c+b^4* a^3*e^2*d+b^4*c^3*a^2*d-b^4*c^3*a^2*e-b^4*c^3*d^2*a+b^4*c^3*d^2*e+b^4*c^3*e^2*a-b^ 4*c^3*e^2*d-b^4*d^3*a^2*c+b^4*d^3*a^2*e+b^4*d^3*c^2*a-b^4*d^3*c^2*e-b^4*d^3*e^2*a+ b^4*d^3*e^2*c+b^4*e^3*a^2*c-b^4*e^3*a^2*d-b^4*e^3*c^2*a+b^4*e^3*c^2*d+b^4*e^3*d^2 *a-b^4*e^3*d^2*c+c^4*a^3*b^2*d-c^4*a^3*b^2*e-c^4*a^3*d^2*b+c^4*a^3*d^2*e+c^4*a^3*e ^2*b-c^4*a^3*e^2*d-c^4*b^3*a^2*d+c^4*b^3*a^2*e+c^4*b^3*d^2*a-c^4*b^3*d^2*e-c^4*b^3 *e^2*a+c^4*b^3*e^2*d+c^4*d^3*a^2*b-c^4*d^3*a^2*e-c^4*d^3*b^2*a+c^4*d^3*b^2*e+c^4* d^3*e^2*a-c^4*d^3*e^2*b-c^4*e^3*a^2*b+c^4*e^3*a^2*d+c^4*e^3*b^2*a-c^4*e^3*b^2*d-c^4 *e^3*d^2*a+c^4*e^3*d^2*b-d^4*a^3*b^2*c+d^4*a^3*b^2*e+d^4*a^3*c^2*b-d^4*a^3*c^2*e-d^4*a^3*e^2*b+d^4*a^3*e^2*c+d^4*b^3*a^2*c-d^4*b^3*a^2*e-d^4*b^3*c^2*a+d^4*b^3*c^2* e+d^4*b^3*e^2*a-d^4*b^3*e^2*c-d^4*c^3*a^2*b+d^4*c^3*a^2*e+d^4*c^3*b^2*a-d^4*c^3*b^ 2*e-d^4*c^3*e^2*a+d^4*c^3*e^2*b+d^4*e^3*a^2*b-d^4*e^3*a^2*c-d^4*e^3*b^2*a+d^4*e^3* b^2*c+d^4*e^3*c^2*a-d^4*e^3*c^2*b+e^4*a^3*b^2*c-e^4*a^3*b^2*d-e^4*a^3*c^2*b+e^4*a^ 3*c^2*d+e^4*a^3*d^2*b-e^4*a^3*d^2*c-e^4*b^3*a^2*c+e^4*b^3*a^2*d+e^4*b^3*c^2*a-e^4* b^3*c^2*d-e^4*b^3*d^2*a+e^4*b^3*d^2*c+e^4*c^3*a^2*b-e^4*c^3*a^2*d-e^4*c^3*b^2*a+e^ 4*c^3*b^2*d+e^4*c^3*d^2*a-e^4*c^3*d^2*b-e^4*d^3*a^2*b+e^4*d^3*a^2*c+e^4*d^3*b^2*a-e^4*d^3*b^2*c-e^4*d^3*c^2*a+e^4*d^3*c^2*bcombine(trig):a^4*b^3*c^2*d-a^4*b^3*c^2*e-a^4*b^3*d^2*c+a^4*b^3*d^2*e+a^4*b^3*e^2*c-a^4*b^3*e^2* d-a^4*c^3*b^2*d+a^4*c^3*b^2*e+a^4*c^3*d^2*b-a^4*c^3*d^2*e-a^4*c^3*e^2*b+a^4*c^3*e^ 2*d+a^4*d^3*b^2*c-a^4*d^3*b^2*e-a^4*d^3*c^2*b+a^4*d^3*c^2*e+a^4*d^3*e^2*b-a^4*d^3* e^2*c-a^4*e^3*b^2*c+a^4*e^3*b^2*d+a^4*e^3*c^2*b-a^4*e^3*c^2*d-a^4*e^3*d^2*b+a^4*e^ 3*d^2*c-b^4*a^3*c^2*d+b^4*a^3*c^2*e+b^4*a^3*d^2*c-b^4*a^3*d^2*e-b^4*a^3*e^2*c+b^4* a^3*e^2*d+b^4*c^3*a^2*d-b^4*c^3*a^2*e-b^4*c^3*d^2*a+b^4*c^3*d^2*e+b^4*c^3*e^2*a-b^ 4*c^3*e^2*d-b^4*d^3*a^2*c+b^4*d^3*a^2*e+b^4*d^3*c^2*a-b^4*d^3*c^2*e-b^4*d^3*e^2*a+ b^4*d^3*e^2*c+b^4*e^3*a^2*c-b^4*e^3*a^2*d-b^4*e^3*c^2*a+b^4*e^3*c^2*d+b^4*e^3*d^2 *a-b^4*e^3*d^2*c+c^4*a^3*b^2*d-c^4*a^3*b^2*e-c^4*a^3*d^2*b+c^4*a^3*d^2*e+c^4*a^3*e ^2*b-c^4*a^3*e^2*d-c^4*b^3*a^2*d+c^4*b^3*a^2*e+c^4*b^3*d^2*a-c^4*b^3*d^2*e-c^4*b^3 *e^2*a+c^4*b^3*e^2*d+c^4*d^3*a^2*b-c^4*d^3*a^2*e-c^4*d^3*b^2*a+c^4*d^3*b^2*e+c^4* d^3*e^2*a-c^4*d^3*e^2*b-c^4*e^3*a^2*b+c^4*e^3*a^2*d+c^4*e^3*b^2*a-c^4*e^3*b^2*d-c^4 *e^3*d^2*a+c^4*e^3*d^2*b-d^4*a^3*b^2*c+d^4*a^3*b^2*e+d^4*a^3*c^2*b-d^4*a^3*c^2*e-d ^4*a^3*e^2*b+d^4*a^3*e^2*c+d^4*b^3*a^2*c-d^4*b^3*a^2*e-d^4*b^3*c^2*a+d^4*b^3*c^2* e+d^4*b^3*e^2*a-d^4*b^3*e^2*c-d^4*c^3*a^2*b+d^4*c^3*a^2*e+d^4*c^3*b^2*a-d^4*c^3*b^ 2*e-d^4*c^3*e^2*a+d^4*c^3*e^2*b+d^4*e^3*a^2*b-d^4*e^3*a^2*c-d^4*e^3*b^2*a+d^4*e^3* b^2*c+d^4*e^3*c^2*a-d^4*e^3*c^2*b+e^4*a^3*b^2*c-e^4*a^3*b^2*d-e^4*a^3*c^2*b+e^4*a^ 3*c^2*d+e^4*a^3*d^2*b-e^4*a^3*d^2*c-e^4*b^3*a^2*c+e^4*b^3*a^2*d+e^4*b^3*c^2*a-e^4* b^3*c^2*d-e^4*b^3*d^2*a+e^4*b^3*d^2*c+e^4*c^3*a^2*b-e^4*c^3*a^2*d-e^4*c^3*b^2*a+e^ 4*c^3*b^2*d+e^4*c^3*d^2*a-e^4*c^3*d^2*b-e^4*d^3*a^2*b+e^4*d^3*a^2*c+e^4*d^3*b^2*a-e^4*d^3*b^2*c-e^4*d^3*c^2*a+e^4*d^3*c^2*bfactor:(c-d)*(b-d)*(b-c)*(a-d)*(a-c)*(a-b)*(-d+e)*(e-c)*(e-b)*(e-a)expand:a^4*b^3*c^2*d-a^4*b^3*c^2*e-a^4*b^3*d^2*c+a^4*b^3*d^2*e+a^4*b^3*e^2*c-a^4*b^3*e^2* d-a^4*c^3*b^2*d+a^4*c^3*b^2*e+a^4*c^3*d^2*b-a^4*c^3*d^2*e-a^4*c^3*e^2*b+a^4*c^3*e^2*d+a^4*d^3*b^2*c-a^4*d^3*b^2*e-a^4*d^3*c^2*b+a^4*d^3*c^2*e+a^4*d^3*e^2*b-a^4*d^3* e^2*c-a^4*e^3*b^2*c+a^4*e^3*b^2*d+a^4*e^3*c^2*b-a^4*e^3*c^2*d-a^4*e^3*d^2*b+a^4*e^ 3*d^2*c-b^4*a^3*c^2*d+b^4*a^3*c^2*e+b^4*a^3*d^2*c-b^4*a^3*d^2*e-b^4*a^3*e^2*c+b^4* a^3*e^2*d+b^4*c^3*a^2*d-b^4*c^3*a^2*e-b^4*c^3*d^2*a+b^4*c^3*d^2*e+b^4*c^3*e^2*a-b^ 4*c^3*e^2*d-b^4*d^3*a^2*c+b^4*d^3*a^2*e+b^4*d^3*c^2*a-b^4*d^3*c^2*e-b^4*d^3*e^2*a+ b^4*d^3*e^2*c+b^4*e^3*a^2*c-b^4*e^3*a^2*d-b^4*e^3*c^2*a+b^4*e^3*c^2*d+b^4*e^3*d^2 *a-b^4*e^3*d^2*c+c^4*a^3*b^2*d-c^4*a^3*b^2*e-c^4*a^3*d^2*b+c^4*a^3*d^2*e+c^4*a^3*e ^2*b-c^4*a^3*e^2*d-c^4*b^3*a^2*d+c^4*b^3*a^2*e+c^4*b^3*d^2*a-c^4*b^3*d^2*e-c^4*b^3 *e^2*a+c^4*b^3*e^2*d+c^4*d^3*a^2*b-c^4*d^3*a^2*e-c^4*d^3*b^2*a+c^4*d^3*b^2*e+c^4* d^3*e^2*a-c^4*d^3*e^2*b-c^4*e^3*a^2*b+c^4*e^3*a^2*d+c^4*e^3*b^2*a-c^4*e^3*b^2*d-c^4 *e^3*d^2*a+c^4*e^3*d^2*b-d^4*a^3*b^2*c+d^4*a^3*b^2*e+d^4*a^3*c^2*b-d^4*a^3*c^2*e-d ^4*a^3*e^2*b+d^4*a^3*e^2*c+d^4*b^3*a^2*c-d^4*b^3*a^2*e-d^4*b^3*c^2*a+d^4*b^3*c^2* e+d^4*b^3*e^2*a-d^4*b^3*e^2*c-d^4*c^3*a^2*b+d^4*c^3*a^2*e+d^4*c^3*b^2*a-d^4*c^3*b^ 2*e-d^4*c^3*e^2*a+d^4*c^3*e^2*b+d^4*e^3*a^2*b-d^4*e^3*a^2*c-d^4*e^3*b^2*a+d^4*e^3* b^2*c+d^4*e^3*c^2*a-d^4*e^3*c^2*b+e^4*a^3*b^2*c-e^4*a^3*b^2*d-e^4*a^3*c^2*b+e^4*a^ 3*c^2*d+e^4*a^3*d^2*b-e^4*a^3*d^2*c-e^4*b^3*a^2*c+e^4*b^3*a^2*d+e^4*b^3*c^2*a-e^4* b^3*c^2*d-e^4*b^3*d^2*a+e^4*b^3*d^2*c+e^4*c^3*a^2*b-e^4*c^3*a^2*d-e^4*c^3*b^2*a+e^ 4*c^3*b^2*d+e^4*c^3*d^2*a-e^4*c^3*d^2*b-e^4*d^3*a^2*b+e^4*d^3*a^2*c+e^4*d^3*b^2*a-e^4*d^3*b^2*c-e^4*d^3*c^2*a+e^4*d^3*c^2*bcombine:a^4*b^3*c^2*d-a^4*b^3*c^2*e-a^4*b^3*d^2*c+a^4*b^3*d^2*e+a^4*b^3*e^2*c-a^4*b^3*e^2* d-a^4*c^3*b^2*d+a^4*c^3*b^2*e+a^4*c^3*d^2*b-a^4*c^3*d^2*e-a^4*c^3*e^2*b+a^4*c^3*e^ 2*d+a^4*d^3*b^2*c-a^4*d^3*b^2*e-a^4*d^3*c^2*b+a^4*d^3*c^2*e+a^4*d^3*e^2*b-a^4*d^3* e^2*c-a^4*e^3*b^2*c+a^4*e^3*b^2*d+a^4*e^3*c^2*b-a^4*e^3*c^2*d-a^4*e^3*d^2*b+a^4*e^ 3*d^2*c-b^4*a^3*c^2*d+b^4*a^3*c^2*e+b^4*a^3*d^2*c-b^4*a^3*d^2*e-b^4*a^3*e^2*c+b^4* a^3*e^2*d+b^4*c^3*a^2*d-b^4*c^3*a^2*e-b^4*c^3*d^2*a+b^4*c^3*d^2*e+b^4*c^3*e^2*a-b^ 4*c^3*e^2*d-b^4*d^3*a^2*c+b^4*d^3*a^2*e+b^4*d^3*c^2*a-b^4*d^3*c^2*e-b^4*d^3*e^2*a+ b^4*d^3*e^2*c+b^4*e^3*a^2*c-b^4*e^3*a^2*d-b^4*e^3*c^2*a+b^4*e^3*c^2*d+b^4*e^3*d^2 *a-b^4*e^3*d^2*c+c^4*a^3*b^2*d-c^4*a^3*b^2*e-c^4*a^3*d^2*b+c^4*a^3*d^2*e+c^4*a^3*e ^2*b-c^4*a^3*e^2*d-c^4*b^3*a^2*d+c^4*b^3*a^2*e+c^4*b^3*d^2*a-c^4*b^3*d^2*e-c^4*b^3 *e^2*a+c^4*b^3*e^2*d+c^4*d^3*a^2*b-c^4*d^3*a^2*e-c^4*d^3*b^2*a+c^4*d^3*b^2*e+c^4* d^3*e^2*a-c^4*d^3*e^2*b-c^4*e^3*a^2*b+c^4*e^3*a^2*d+c^4*e^3*b^2*a-c^4*e^3*b^2*d-c^4 *e^3*d^2*a+c^4*e^3*d^2*b-d^4*a^3*b^2*c+d^4*a^3*b^2*e+d^4*a^3*c^2*b-d^4*a^3*c^2*e-d ^4*a^3*e^2*b+d^4*a^3*e^2*c+d^4*b^3*a^2*c-d^4*b^3*a^2*e-d^4*b^3*c^2*a+d^4*b^3*c^2* e+d^4*b^3*e^2*a-d^4*b^3*e^2*c-d^4*c^3*a^2*b+d^4*c^3*a^2*e+d^4*c^3*b^2*a-d^4*c^3*b^ 2*e-d^4*c^3*e^2*a+d^4*c^3*e^2*b+d^4*e^3*a^2*b-d^4*e^3*a^2*c-d^4*e^3*b^2*a+d^4*e^3* b^2*c+d^4*e^3*c^2*a-d^4*e^3*c^2*b+e^4*a^3*b^2*c-e^4*a^3*b^2*d-e^4*a^3*c^2*b+e^4*a^ 3*c^2*d+e^4*a^3*d^2*b-e^4*a^3*d^2*c-e^4*b^3*a^2*c+e^4*b^3*a^2*d+e^4*b^3*c^2*a-e^4* b^3*c^2*d-e^4*b^3*d^2*a+e^4*b^3*d^2*c+e^4*c^3*a^2*b-e^4*c^3*a^2*d-e^4*c^3*b^2*a+e^ 4*c^3*b^2*d+e^4*c^3*d^2*a-e^4*c^3*d^2*b-e^4*d^3*a^2*b+e^4*d^3*a^2*c+e^4*d^3*b^2*a-e^4*d^3*b^2*c-e^4*d^3*c^2*a+e^4*d^3*c^2*bconvert(exp):a^4*b^3*c^2*d-a^4*b^3*c^2*e-a^4*b^3*d^2*c+a^4*b^3*d^2*e+a^4*b^3*e^2*c-a^4*b^3*e^2* d-a^4*c^3*b^2*d+a^4*c^3*b^2*e+a^4*c^3*d^2*b-a^4*c^3*d^2*e-a^4*c^3*e^2*b+a^4*c^3*e^ 2*d+a^4*d^3*b^2*c-a^4*d^3*b^2*e-a^4*d^3*c^2*b+a^4*d^3*c^2*e+a^4*d^3*e^2*b-a^4*d^3* e^2*c-a^4*e^3*b^2*c+a^4*e^3*b^2*d+a^4*e^3*c^2*b-a^4*e^3*c^2*d-a^4*e^3*d^2*b+a^4*e^ 3*d^2*c-b^4*a^3*c^2*d+b^4*a^3*c^2*e+b^4*a^3*d^2*c-b^4*a^3*d^2*e-b^4*a^3*e^2*c+b^4* a^3*e^2*d+b^4*c^3*a^2*d-b^4*c^3*a^2*e-b^4*c^3*d^2*a+b^4*c^3*d^2*e+b^4*c^3*e^2*a-b^ 4*c^3*e^2*d-b^4*d^3*a^2*c+b^4*d^3*a^2*e+b^4*d^3*c^2*a-b^4*d^3*c^2*e-b^4*d^3*e^2*a+ b^4*d^3*e^2*c+b^4*e^3*a^2*c-b^4*e^3*a^2*d-b^4*e^3*c^2*a+b^4*e^3*c^2*d+b^4*e^3*d^2 *a-b^4*e^3*d^2*c+c^4*a^3*b^2*d-c^4*a^3*b^2*e-c^4*a^3*d^2*b+c^4*a^3*d^2*e+c^4*a^3*e ^2*b-c^4*a^3*e^2*d-c^4*b^3*a^2*d+c^4*b^3*a^2*e+c^4*b^3*d^2*a-c^4*b^3*d^2*e-c^4*b^3 *e^2*a+c^4*b^3*e^2*d+c^4*d^3*a^2*b-c^4*d^3*a^2*e-c^4*d^3*b^2*a+c^4*d^3*b^2*e+c^4* d^3*e^2*a-c^4*d^3*e^2*b-c^4*e^3*a^2*b+c^4*e^3*a^2*d+c^4*e^3*b^2*a-c^4*e^3*b^2*d-c^4 *e^3*d^2*a+c^4*e^3*d^2*b-d^4*a^3*b^2*c+d^4*a^3*b^2*e+d^4*a^3*c^2*b-d^4*a^3*c^2*e-d ^4*a^3*e^2*b+d^4*a^3*e^2*c+d^4*b^3*a^2*c-d^4*b^3*a^2*e-d^4*b^3*c^2*a+d^4*b^3*c^2* e+d^4*b^3*e^2*a-d^4*b^3*e^2*c-d^4*c^3*a^2*b+d^4*c^3*a^2*e+d^4*c^3*b^2*a-d^4*c^3*b^ 2*e-d^4*c^3*e^2*a+d^4*c^3*e^2*b+d^4*e^3*a^2*b-d^4*e^3*a^2*c-d^4*e^3*b^2*a+d^4*e^3* b^2*c+d^4*e^3*c^2*a-d^4*e^3*c^2*b+e^4*a^3*b^2*c-e^4*a^3*b^2*d-e^4*a^3*c^2*b+e^4*a^ 3*c^2*d+e^4*a^3*d^2*b-e^4*a^3*d^2*c-e^4*b^3*a^2*c+e^4*b^3*a^2*d+e^4*b^3*c^2*a-e^4* b^3*c^2*d-e^4*b^3*d^2*a+e^4*b^3*d^2*c+e^4*c^3*a^2*b-e^4*c^3*a^2*d-e^4*c^3*b^2*a+e^ 4*c^3*b^2*d+e^4*c^3*d^2*a-e^4*c^3*d^2*b-e^4*d^3*a^2*b+e^4*d^3*a^2*c+e^4*d^3*b^2*a-e^4*d^3*b^2*c-e^4*d^3*c^2*a+e^4*d^3*c^2*bconvert(sincos):a^4*b^3*c^2*d-a^4*b^3*c^2*e-a^4*b^3*d^2*c+a^4*b^3*d^2*e+a^4*b^3*e^2*c-a^4*b^3*e^2* d-a^4*c^3*b^2*d+a^4*c^3*b^2*e+a^4*c^3*d^2*b-a^4*c^3*d^2*e-a^4*c^3*e^2*b+a^4*c^3*e^ 2*d+a^4*d^3*b^2*c-a^4*d^3*b^2*e-a^4*d^3*c^2*b+a^4*d^3*c^2*e+a^4*d^3*e^2*b-a^4*d^3* e^2*c-a^4*e^3*b^2*c+a^4*e^3*b^2*d+a^4*e^3*c^2*b-a^4*e^3*c^2*d-a^4*e^3*d^2*b+a^4*e^ 3*d^2*c-b^4*a^3*c^2*d+b^4*a^3*c^2*e+b^4*a^3*d^2*c-b^4*a^3*d^2*e-b^4*a^3*e^2*c+b^4* a^3*e^2*d+b^4*c^3*a^2*d-b^4*c^3*a^2*e-b^4*c^3*d^2*a+b^4*c^3*d^2*e+b^4*c^3*e^2*a-b^ 4*c^3*e^2*d-b^4*d^3*a^2*c+b^4*d^3*a^2*e+b^4*d^3*c^2*a-b^4*d^3*c^2*e-b^4*d^3*e^2*a+ b^4*d^3*e^2*c+b^4*e^3*a^2*c-b^4*e^3*a^2*d-b^4*e^3*c^2*a+b^4*e^3*c^2*d+b^4*e^3*d^2 *a-b^4*e^3*d^2*c+c^4*a^3*b^2*d-c^4*a^3*b^2*e-c^4*a^3*d^2*b+c^4*a^3*d^2*e+c^4*a^3*e ^2*b-c^4*a^3*e^2*d-c^4*b^3*a^2*d+c^4*b^3*a^2*e+c^4*b^3*d^2*a-c^4*b^3*d^2*e-c^4*b^3 *e^2*a+c^4*b^3*e^2*d+c^4*d^3*a^2*b-c^4*d^3*a^2*e-c^4*d^3*b^2*a+c^4*d^3*b^2*e+c^4* d^3*e^2*a-c^4*d^3*e^2*b-c^4*e^3*a^2*b+c^4*e^3*a^2*d+c^4*e^3*b^2*a-c^4*e^3*b^2*d-c^4 *e^3*d^2*a+c^4*e^3*d^2*b-d^4*a^3*b^2*c+d^4*a^3*b^2*e+d^4*a^3*c^2*b-d^4*a^3*c^2*e-d ^4*a^3*e^2*b+d^4*a^3*e^2*c+d^4*b^3*a^2*c-d^4*b^3*a^2*e-d^4*b^3*c^2*a+d^4*b^3*c^2* e+d^4*b^3*e^2*a-d^4*b^3*e^2*c-d^4*c^3*a^2*b+d^4*c^3*a^2*e+d^4*c^3*b^2*a-d^4*c^3*b^2*e-d^4*c^3*e^2*a+d^4*c^3*e^2*b+d^4*e^3*a^2*b-d^4*e^3*a^2*c-d^4*e^3*b^2*a+d^4*e^3* b^2*c+d^4*e^3*c^2*a-d^4*e^3*c^2*b+e^4*a^3*b^2*c-e^4*a^3*b^2*d-e^4*a^3*c^2*b+e^4*a^ 3*c^2*d+e^4*a^3*d^2*b-e^4*a^3*d^2*c-e^4*b^3*a^2*c+e^4*b^3*a^2*d+e^4*b^3*c^2*a-e^4* b^3*c^2*d-e^4*b^3*d^2*a+e^4*b^3*d^2*c+e^4*c^3*a^2*b-e^4*c^3*a^2*d-e^4*c^3*b^2*a+e^ 4*c^3*b^2*d+e^4*c^3*d^2*a-e^4*c^3*d^2*b-e^4*d^3*a^2*b+e^4*d^3*a^2*c+e^4*d^3*b^2*a-e^4*d^3*b^2*c-e^4*d^3*c^2*a+e^4*d^3*c^2*bconvert(tan):a^4*b^3*c^2*d-a^4*b^3*c^2*e-a^4*b^3*d^2*c+a^4*b^3*d^2*e+a^4*b^3*e^2*c-a^4*b^3*e^2* d-a^4*c^3*b^2*d+a^4*c^3*b^2*e+a^4*c^3*d^2*b-a^4*c^3*d^2*e-a^4*c^3*e^2*b+a^4*c^3*e^ 2*d+a^4*d^3*b^2*c-a^4*d^3*b^2*e-a^4*d^3*c^2*b+a^4*d^3*c^2*e+a^4*d^3*e^2*b-a^4*d^3* e^2*c-a^4*e^3*b^2*c+a^4*e^3*b^2*d+a^4*e^3*c^2*b-a^4*e^3*c^2*d-a^4*e^3*d^2*b+a^4*e^ 3*d^2*c-b^4*a^3*c^2*d+b^4*a^3*c^2*e+b^4*a^3*d^2*c-b^4*a^3*d^2*e-b^4*a^3*e^2*c+b^4* a^3*e^2*d+b^4*c^3*a^2*d-b^4*c^3*a^2*e-b^4*c^3*d^2*a+b^4*c^3*d^2*e+b^4*c^3*e^2*a-b^ 4*c^3*e^2*d-b^4*d^3*a^2*c+b^4*d^3*a^2*e+b^4*d^3*c^2*a-b^4*d^3*c^2*e-b^4*d^3*e^2*a+ b^4*d^3*e^2*c+b^4*e^3*a^2*c-b^4*e^3*a^2*d-b^4*e^3*c^2*a+b^4*e^3*c^2*d+b^4*e^3*d^2 *a-b^4*e^3*d^2*c+c^4*a^3*b^2*d-c^4*a^3*b^2*e-c^4*a^3*d^2*b+c^4*a^3*d^2*e+c^4*a^3*e ^2*b-c^4*a^3*e^2*d-c^4*b^3*a^2*d+c^4*b^3*a^2*e+c^4*b^3*d^2*a-c^4*b^3*d^2*e-c^4*b^3 *e^2*a+c^4*b^3*e^2*d+c^4*d^3*a^2*b-c^4*d^3*a^2*e-c^4*d^3*b^2*a+c^4*d^3*b^2*e+c^4* d^3*e^2*a-c^4*d^3*e^2*b-c^4*e^3*a^2*b+c^4*e^3*a^2*d+c^4*e^3*b^2*a-c^4*e^3*b^2*d-c^4 *e^3*d^2*a+c^4*e^3*d^2*b-d^4*a^3*b^2*c+d^4*a^3*b^2*e+d^4*a^3*c^2*b-d^4*a^3*c^2*e-d ^4*a^3*e^2*b+d^4*a^3*e^2*c+d^4*b^3*a^2*c-d^4*b^3*a^2*e-d^4*b^3*c^2*a+d^4*b^3*c^2* e+d^4*b^3*e^2*a-d^4*b^3*e^2*c-d^4*c^3*a^2*b+d^4*c^3*a^2*e+d^4*c^3*b^2*a-d^4*c^3*b^ 2*e-d^4*c^3*e^2*a+d^4*c^3*e^2*b+d^4*e^3*a^2*b-d^4*e^3*a^2*c-d^4*e^3*b^2*a+d^4*e^3* b^2*c+d^4*e^3*c^2*a-d^4*e^3*c^2*b+e^4*a^3*b^2*c-e^4*a^3*b^2*d-e^4*a^3*c^2*b+e^4*a^ 3*c^2*d+e^4*a^3*d^2*b-e^4*a^3*d^2*c-e^4*b^3*a^2*c+e^4*b^3*a^2*d+e^4*b^3*c^2*a-e^4* b^3*c^2*d-e^4*b^3*d^2*a+e^4*b^3*d^2*c+e^4*c^3*a^2*b-e^4*c^3*a^2*d-e^4*c^3*b^2*a+e^ 4*c^3*b^2*d+e^4*c^3*d^2*a-e^4*c^3*d^2*b-e^4*d^3*a^2*b+e^4*d^3*a^2*c+e^4*d^3*b^2*a-e^4*d^3*b^2*c-e^4*d^3*c^2*a+e^4*d^3*c^2*bcollect(e):a^4*b^3*c^2*d-a^4*b^3*d^2*c-a^4*c^3*b^2*d+a^4*c^3*d^2*b+a^4*d^3*b^2*c-a^4*d^3*c^2* b-b^4*a^3*c^2*d+b^4*a^3*d^2*c+b^4*c^3*a^2*d-b^4*c^3*d^2*a-b^4*d^3*a^2*c+b^4*d^3*c^ 2*a+c^4*a^3*b^2*d-c^4*a^3*d^2*b-c^4*b^3*a^2*d+c^4*b^3*d^2*a+c^4*d^3*a^2*b-c^4*d^3* b^2*a-d^4*a^3*b^2*c+d^4*a^3*c^2*b+d^4*b^3*a^2*c-d^4*b^3*c^2*a-d^4*c^3*a^2*b+d^4*c^ 3*b^2*a+(a^3*b^2*c-a^3*b^2*d-a^3*c^2*b+a^3*c^2*d+a^3*d^2*b-a^3*d^2*c-b^3*a^2*c+b^3 *a^2*d+b^3*c^2*a-b^3*c^2*d-b^3*d^2*a+b^3*d^2*c+c^3*a^2*b-c^3*a^2*d-c^3*b^2*a+c^3*b ^2*d+c^3*d^2*a-c^3*d^2*b-d^3*a^2*b+d^3*a^2*c+d^3*b^2*a-d^3*b^2*c-d^3*c^2*a+d^3*c^2 *b)*e^4+(-a^4*b^2*c-c^4*d*b^2-c^4*d^2*a+c^4*d^2*b+a^4*b^2*d+a^4*c^2*b-a^4*c^2*d-c^4 *a^2*b+c^4*a^2*d+c^4*b^2*a+b^4*a^2*c-b^4*a^2*d-b^4*c^2*a+c^2*d*b^4+b^4*d^2*a-d^2*c*b^4-a^4*d^2*b+a^4*d^2*c+d^4*a^2*b-d^4*a^2*c-d^4*b^2*a+c*d^4*b^2+d^4*c^2*a-c^2*d^ 4*b)*e^3+(b^4*a^3*d-c^4*b^3*a-a^4*c^3*b-c*d^4*b^3+c^4*d*b^3-a^4*b^3*d-b^4*a^3*c+a^4 *c^3*d+a^4*b^3*c-a^4*d^3*c-c^3*d*b^4+d^4*b^3*a+a^4*d^3*b+b^4*c^3*a-c^4*a^3*d-d^4*a ^3*b+c^4*d^3*a+c^4*a^3*b-c^4*d^3*b-b^4*d^3*a+c*d^3*b^4+d^4*a^3*c-d^4*c^3*a+c^3*d^ 4*b)*e^2+(a^4*b^3*d^2+c^4*a^3*d^2-b^4*c^3*a^2-a^4*b^3*c^2+a^4*d^3*c^2+d^4*c^3*a^2-c^4*d^2*b^3-a^4*c^3*d^2+c^4*d^3*b^2+c^4*b^3*a^2-c^4*d^3*a^2-a^4*d^3*b^2-b^4*a^3*d^ 2+b^4*d^3*a^2+b^4*a^3*c^2-c^3*d^4*b^2-d^4*a^3*c^2+c^2*d^4*b^3+a^4*c^3*b^2+c^3*d^2 *b^4-c^2*d^3*b^4-c^4*a^3*b^2+d^4*a^3*b^2-d^4*b^3*a^2)*emwcos2sin:a^4*b^3*c^2*d-a^4*b^3*c^2*e-a^4*b^3*d^2*c+a^4*b^3*d^2*e+a^4*b^3*e^2*c-a^4*b^3*e^2* d-a^4*c^3*b^2*d+a^4*c^3*b^2*e+a^4*c^3*d^2*b-a^4*c^3*d^2*e-a^4*c^3*e^2*b+a^4*c^3*e^ 2*d+a^4*d^3*b^2*c-a^4*d^3*b^2*e-a^4*d^3*c^2*b+a^4*d^3*c^2*e+a^4*d^3*e^2*b-a^4*d^3* e^2*c-a^4*e^3*b^2*c+a^4*e^3*b^2*d+a^4*e^3*c^2*b-a^4*e^3*c^2*d-a^4*e^3*d^2*b+a^4*e^ 3*d^2*c-b^4*a^3*c^2*d+b^4*a^3*c^2*e+b^4*a^3*d^2*c-b^4*a^3*d^2*e-b^4*a^3*e^2*c+b^4* a^3*e^2*d+b^4*c^3*a^2*d-b^4*c^3*a^2*e-b^4*c^3*d^2*a+b^4*c^3*d^2*e+b^4*c^3*e^2*a-b^ 4*c^3*e^2*d-b^4*d^3*a^2*c+b^4*d^3*a^2*e+b^4*d^3*c^2*a-b^4*d^3*c^2*e-b^4*d^3*e^2*a+ b^4*d^3*e^2*c+b^4*e^3*a^2*c-b^4*e^3*a^2*d-b^4*e^3*c^2*a+b^4*e^3*c^2*d+b^4*e^3*d^2 *a-b^4*e^3*d^2*c+c^4*a^3*b^2*d-c^4*a^3*b^2*e-c^4*a^3*d^2*b+c^4*a^3*d^2*e+c^4*a^3*e ^2*b-c^4*a^3*e^2*d-c^4*b^3*a^2*d+c^4*b^3*a^2*e+c^4*b^3*d^2*a-c^4*b^3*d^2*e-c^4*b^3 *e^2*a+c^4*b^3*e^2*d+c^4*d^3*a^2*b-c^4*d^3*a^2*e-c^4*d^3*b^2*a+c^4*d^3*b^2*e+c^4* d^3*e^2*a-c^4*d^3*e^2*b-c^4*e^3*a^2*b+c^4*e^3*a^2*d+c^4*e^3*b^2*a-c^4*e^3*b^2*d-c^4 *e^3*d^2*a+c^4*e^3*d^2*b-d^4*a^3*b^2*c+d^4*a^3*b^2*e+d^4*a^3*c^2*b-d^4*a^3*c^2*e-d ^4*a^3*e^2*b+d^4*a^3*e^2*c+d^4*b^3*a^2*c-d^4*b^3*a^2*e-d^4*b^3*c^2*a+d^4*b^3*c^2* e+d^4*b^3*e^2*a-d^4*b^3*e^2*c-d^4*c^3*a^2*b+d^4*c^3*a^2*e+d^4*c^3*b^2*a-d^4*c^3*b^ 2*e-d^4*c^3*e^2*a+d^4*c^3*e^2*b+d^4*e^3*a^2*b-d^4*e^3*a^2*c-d^4*e^3*b^2*a+d^4*e^3* b^2*c+d^4*e^3*c^2*a-d^4*e^3*c^2*b+e^4*a^3*b^2*c-e^4*a^3*b^2*d-e^4*a^3*c^2*b+e^4*a^ 3*c^2*d+e^4*a^3*d^2*b-e^4*a^3*d^2*c-e^4*b^3*a^2*c+e^4*b^3*a^2*d+e^4*b^3*c^2*a-e^4* b^3*c^2*d-e^4*b^3*d^2*a+e^4*b^3*d^2*c+e^4*c^3*a^2*b-e^4*c^3*a^2*d-e^4*c^3*b^2*a+e^ 4*c^3*b^2*d+e^4*c^3*d^2*a-e^4*c^3*d^2*b-e^4*d^3*a^2*b+e^4*d^3*a^2*c+e^4*d^3*b^2*a-e^4*d^3*b^2*c-e^4*d^3*c^2*a+e^4*d^3*c^2*bans =(c-d)*(b-d)*(b-c)*(a-d)*(a-c)*(a-b)*(-d+e)*(e-c)*(e-b)*(e-a)>>第十三题>> A=[-2,0.5,-0.5,0.5;0,-1.5,0.5,-0.5;2,0.5,-4.5,0.5;2,1,-2,-2];[V J]=jordan(sym(A))V =[ 0, 1/2, 1/2, -1/4][ 0, 0, 1/2, 1][ 1/4, 1/2, 1/2, -1/4][ 1/4, 1/2, 1, -1/4]J =[ -4, 0, 0, 0][ 0, -2, 1, 0][ 0, 0, -2, 1][ 0, 0, 0, -2]第十四题数值方法>> A=[3,-6,-4,0,5;1,4,2,-2,4;-6,3,-6,7,3;-13,10,0,-11,0;0,4,0,3,4];>> B=[3,-2,1;-2,9,2;-2,-1,9];>> C=[-2,1,-1;4,1,2;5,-6,1;6,-4,-4;-6,6,-3];>> X=lyap(A,B,C)X =-2.3192 -0.4678 0.1505-3.6284 0.1579 0.06295.4246 -1.0516 -0.5090-0.5718 2.5848 -0.36493.0417 -0.6265 0.1580>> norm(A*X+X*B+C)ans =3.8830e-014解析方法>> editfunction X=lyap(A,B,C)if nargin==2,C=B;B=A';end[nr,nc]=size(C);A0=kron(A,eye(nc))+kron(eye(nr),B'); tryC1=C';x0=-inv(A0)*C1(:);X=reshape(x0,nc,nr)'; catch,error('singular matrix found.'),end>> A=[3,-6,-4,0,5;1,4,2,-2,4;-6,3,-6,7,3;-13,10,0,-11,0;0,4,0,3,4];>> B=[3,-2,1;-2,-9,2;-2,-1,9];>> C=[-2,1,-1;4,1,2;5,-6,1;6,-4,-4;-6,6,-3];X=lyap(sym(A),B,C)X =[ -434641749950/107136516451, -4664546747350/321409549353, 503105815912/321409549353][ 3809507498/107136516451, 8059112319373/321409549353, -880921527508/321409549353][ 1016580400173/107136516451, 8334897743767/321409549353, -1419901706449/321409549353][ 288938859984/107136516451, 6956912657222/321409549353, -927293592476/321409549353][ 827401644798/107136516451, 10256166034813/321409549353, -1209595497577/321409549353]>> A*X+X*B+Cans =[ 0, 0, 0][ 0, 0, 0][ 0, 0, 0][ 0, 0, 0][ 0, 0, 0]第十五题（1）>> A=[-4.5,0,0.5,-1.5;-0.5,-4,0.5,-0.5;1.5,1,-2.5,1.5;0,-1,-1,-3];>> A=sym(A);syms t;>> expm(A*t)ans =[ 1/2*exp(-3*t)-1/2*t*exp(-3*t)+1/2*exp(-5*t)+1/2*t^2*exp(-3*t), 1/2*exp(-5*t)-1/2*exp(-3*t)+t*exp(-3*t),1/2*t*exp(-3*t)+1/2*t^2*exp(-3*t),1/2*exp(-5*t)-1/2*exp(-3*t)-1/2*t*exp(-3*t)+1/2*t^2*exp(-3*t)][ 1/2*t*exp(-3*t)+1/2*exp(-5*t)-1/2*exp(-3*t), 1/2*exp(-3*t)+1/2*exp(-5*t),1/2*t*exp(-3*t), 1/2*t*exp(-3*t)+1/2*exp(-5*t)-1/2*exp(-3*t)][ 1/2*t*exp(-3*t)-1/2*exp(-5*t)+1/2*exp(-3*t), -1/2*exp(-5*t)+1/2*exp(-3*t),exp(-3*t)+1/2*t*exp(-3*t), 1/2*t*exp(-3*t)-1/2*exp(-5*t)+1/2*exp(-3*t)][ -1/2*t^2*exp(-3*t), -t*exp(-3*t), -1/2*t^2*exp(-3*t)-t*exp(-3*t), exp(-3*t)-1/2*t^2*exp(-3*t)](2)>> A=[-4.5,0,0.5,-1.5;-0.5,-4,0.5,-0.5;1.5,1,-2.5,1.5;0,-1,-1,-3];>> A=sym(A);syms t;>> sin(A*t)ans =[ -sin(9/2*t), 0, sin(1/2*t), -sin(3/2*t)][ -sin(1/2*t), -sin(4*t), sin(1/2*t), -sin(1/2*t)][ sin(3/2*t), sin(t), -sin(5/2*t), sin(3/2*t)][ 0, -sin(t), -sin(t), -sin(3*t)](3)第二部分第一题（1）>> syms a t;f=sin(a*t)/t;laplace(f)ans =atan(a/s)（2）>> syms t a;f=t^5*sin(a*t);laplace(f)ans =60*i*(-1/(s-i*a)^6+1/(s+i*a)^6)（3）>> syms t a;f=t^8*cos(a*t);laplace(f)ans =20160/(s-i*a)^9+20160/(s+i*a)^9第二题（1）>> syms s a b;F=1/(s^2*(s^2-a^2)*(a+b));ilaplace(F)ans =1/(a+b)*(-1/a^2*t+1/a^3*sinh(a*t))（2）>> syms s a b;F=sqrt(s-a)-sqrt(s-b);ilaplace(F)ans =1/2/t^(3/2)/pi^(1/2)*(exp(b*t)-exp(a*t))（3）>> syms s a b;F=log((s-a)/(s-b));ilaplace(F)ans =1/t*(exp(b*t)-exp(a*t))第三题（1）>> syms x;f=x^2*(3*sym(pi)-2*abs(x));F=fourier(f)F =-6*(4+pi^2*dirac(2,w)*w^4)/w^4>> ifourier(F)ans =x^2*(-4*x*heaviside(x)+3*pi+2*x)（2）>> syms t;f=t^2*(t-2*sym(pi))^2;F=fourier(f)F =2*pi*(-4*pi^2*dirac(2,w)+4*i*pi*dirac(3,w)+dirac(4,w)) >> ifourier(F)ans =x^2*(-2*pi+x)^2第四题（1）>> syms k a T;f=cos(k*a*T);F=ztrans(f)F =(z-cos(a*T))*z/(z^2-2*z*cos(a*T)+1)>> f1=iztrans(F)f1 =cos(a*T*n)（2）>> syms k T a;f=(k*T)^2*exp(-a*k*T);F=ztrans(f)F =T^2*z*exp(-a*T)*(z+exp(-a*T))/(z-exp(-a*T))^3>> f1=iztrans(F)f1 =T^2*(1/exp(a*T))^n*n^2（3）>> syms a k T;f=(a*k*T-1+exp(-a*k*T))/a;F=ztrans(f) F =1/a*(a*T*z/(z-1)^2-z/(z-1)+z/exp(-a*T)/(z/exp(-a*T)-1)) >> iztrans(F)ans =((1/exp(a*T))^n+a*T*n-1)/a第五题（1）>> syms x;x1=solve('exp(-(x+1)^2+pi/2)*sin(5*x+2)') x1 =-2/5>> subs('exp(-(x+1)^2+pi/2)*sin(5*x+2)',x,x1)ans =>> f=inline('exp(-(x+1).^2+pi/2).*sin(5*x+2)','x');>> x2=fsolve(f,0)Optimization terminated: first-order optimality is less than options.TolFun. x2 =0.2283>> subs('exp(-(x+1)^2+pi/2)*sin(5*x+2)',x,x2)ans =4.7509e-008（2）>> syms x;y1=solve('(x^2+y^2+x*y)*exp(-x^2-y^2-x*y)=0','y')y1 =(-1/2+1/2*i*3^(1/2))*x(-1/2-1/2*i*3^(1/2))*x>> y2=simple(subs('(x^2+y^2+x*y)*exp(-x^2-y*2-x*y)','y',y1))y2 =第六题>> syms x c;y=int((exp(x)-c*x)^2,x,0,1)y =1/2*exp(1)^2+1/3*c^2-1/2-2*c>> editfunction y=f(c)y=1/2*exp(1)^2+1/3*c-1/2-2*c>> x=fminsearch('f',0)x =3.1691e+026>> ezplot(y,[0,5])第七题>> editfunction y=f1(x)y=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1) function [c,ce]=f2(x);ce=[];c=[x(1)+x(2);x(1)*x(2)-x(1)-x(2)+1.5;-10-x(1)*x(2)]; >> edit>> A=[];B=[];Aeq=[];Beq=[];xm=[-10;-10];xM=[10;10];>> x0=(xm+xM)/2;>> ff=optimset;ff.TolX=1e-10;ff.TolFun=1e-20;>> x=fmincon('f1',x0,A,B,Aeq,Beq,xm,xM,'f2',ff)Maximum number of function evaluations exceeded;increase OPTIONS.MaxFunEvals.x =0.41950.4195>> i=1;x=x0;>> while (1)[x,a,b]=fmincon('f1',x,A,B,Aeq,Beq,xm,xM,'f2',ff);if b>0,break;endi=i+1;end>> x,ix =1.1825-1.7398i =5第八题function y=f3 (x)y=-(592*x(1)+381*x(2)+273*x(3)+55*x(4)+48*x(5)+37*x(6)+23*x(7)); function [c,ceq]=f4(x)c=[];ceq=[];>> A=[3534,2356,1767,589,528,451,304];B=119567;intlist=ones(7,1);Aeq=[];Beq=[];xm=zeros(7,1);xM=inf*ones(7,1);x0=zeros(7,1);ff=optimset;ff.TolFun=1e-6;ff.TolX=1e-6;ff.TolCon=1e-20;ff.MaxSQPIter=1000;settings=[0];ix=(intlist==1);xm(ix)=ceil(xm(ix));xM(ix)=floor(xM(ix)); [errmsg,f,x]=bnb20('f3',x0,intlist,xm,xM,A,B,Aeq,Beq,f4, settings,ff)if length(errmsg)==0,x=round(x),endx =3221。