Matlab实验报告三
Matlab实验报告(三)-MATLAB绘图
实验目的1．掌握MATLAB的基本绘图命令。
2．掌握运用MATLAB绘制一维、二维、三维图形的方法.3．给图形加以修饰。
一、预备知识1．基本绘图命令plotplot绘图命令一共有三种形式：⑴plot（y)是plot命令中最为简单的形式，当y为向量时,以y的元素为纵坐标，元素相应的序列号为横坐标,绘制出连线；若y为实矩阵，则按照列绘出每列元素和其序列号的对应关系，曲线数等于矩阵的列数;当y为复矩阵时，则按列以每列元素的实部为横坐标，以虚部为纵坐标，绘出曲线，曲线数等于列数。
⑵ plot（x,y,［linspec］）其中linspec是可选的，用它来说明线型。
当x和y为同维向量时，以x为横坐标，y为纵坐标绘制曲线;当x是向量，y是每行元素数目和x维数相同的矩阵时,将绘出以x为横坐标，以y中每行元素为纵坐标的多条曲线，曲线数等于矩阵行数；当x为矩阵，y为相应向量时，使用该命令也能绘出相应图形。
⑶ plot（x1,y1,x2,y2，x3,y3……）能够绘制多条曲线,每条曲线分别以x和y为横纵坐标，各条曲线互不影响。
线型和颜色MATLAB可以对线型和颜色进行设定，线型和颜色种类如下：线：—实线：点线 -.虚点线——折线点：.圆点 +加号 *星号 x x型 o 空心小圆颜色：y 黄 r 红 g 绿 b 蓝 w 白 k 黑 m 紫 c 青特殊的二维图形函数表5 特殊2维绘图函数[1］ 直方图在实际中,常会遇到离散数据,当需要比较数据、分析数据在总量中的比例时，直方图就是一种理想的选择，但要注意该方法适用于数据较少的情况。
直方图的绘图函数有以下两种基本形式。
·bar(x,y) 绘制m*n 矩阵的直方图.其中y 为m ＊n 矩阵或向量,x 必须单向递增。
·bar(y） 绘制y 向量的直方图，x 向量默认为x=1：m close all; ％关闭所有的图形视窗。
x=1:10；y=rand （size(x ）); bar(x,y ）; ％绘制直方图.123456789100.51Bar(）函数还有barh （）和errorbar （）两种形式，barh()用来绘制水平方向的直方图，其参数与bar()相同，当知道资料的误差值时，可用errorbar （)绘制出误差范围,其一般语法形式为：errorbar （x,y,l,u)其中x,y 是其绘制曲线的坐标,l ，u 是曲线误差的最小值和最大值，制图时，l 向量在曲线下方,u 向量在曲线上方。
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的实验报告
基于matlab的实验报告实验报告：基于MATLAB 的实验一、实验目的通过使用MATLAB 软件，掌握如何进行数据分析、图像处理、算法实现等一系列实验操作，提高实验者的实践能力和动手能力。
二、实验原理MATLAB 是一种在科学计算和技术开发领域广泛应用的计算机软件。
它能进行矩阵计算、绘制函数和数据图像、实现算法以及进行数据分析等。
通过掌握MATLAB 的使用，能够快速、高效地解决各种科学和工程问题。
三、实验内容1. 数据分析：使用MATLAB 的数据分析工具进行数据的导入、处理和分析。
2. 图像处理：利用MATLAB 的图像处理工具包对图像进行滤波、增强、分割等操作。
3. 算法实现：使用MATLAB 实现常用的算法，如排序、搜索、图像压缩等。
四、实验步骤1. 数据分析：（1）使用MATLAB 的读取数据函数将数据导入MATLAB 环境中。
（2）利用MATLAB 的数据处理函数进行数据清洗和预处理。
（3）使用MATLAB 的统计工具进行数据分析，如求平均值、标准差等。
（4）利用MATLAB 的绘图函数将分析结果可视化。
2. 图像处理：（1）使用MATLAB 的读取图像函数将图像导入MATLAB 环境中。
（2）利用MATLAB 的图像处理工具包进行滤波操作，如均值滤波、中值滤波等。
（3）使用MATLAB 的图像增强函数对图像进行锐化、变换等操作。
（4）利用MATLAB 的图像分割算法对图像进行分割。
3. 算法实现：（1）使用MATLAB 编写排序算法，如冒泡排序、快速排序等。
（2）使用MATLAB 编写搜索算法，如二分查找、线性搜索等。
（3）使用MATLAB 实现图像压缩算法，如离散余弦变换（DCT）。
五、实验结果实验中，我们使用MATLAB 完成了数据分析、图像处理和算法实现的一系列实验操作。
通过数据分析，我们成功导入了数据并对其进行了清洗和预处理，最后得到了数据的统计结果。
在图像处理方面，我们对图像进行了滤波、增强和分割等操作，最终得到了处理后的图像。
MATLAB实验报告(1-4)
信号与系统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语言的一些语法。
实验3 Matlab程序设计1实验报告
实验3 Matlab程序设计1实验报告I am going to write the report for the Matlab program design experiment 3. In this experiment, I was required to use Matlab to design a program to solve a specific problem.I had to write the program, test it, and then write areport on the results.First, I started by analyzing the problem and breaking it down into smaller, more manageable parts. This is a crucial step in programming, as it allows me to understand the problem fully and come up with an effective solution. Once I had a clear understanding of the problem, I began writing the program in Matlab.I encountered some challenges along the way, such as debugging errors and optimizing the code for efficiency. However, with some perseverance and problem-solving skills, I was able to overcome these challenges and successfully write the program.After writing the program, I tested it with different inputs to ensure that it produced the correct outputs inall cases. Testing is an essential part of programming, asit helps to identify and fix any errors or bugs in the code.Once the program was tested and working correctly, I wrote a report on the results. In the report, I explainedthe problem, my approach to solving it, and the results of the program. I also included any insights or observations that I gained from working on the program.Overall, this experiment was a valuable learning experience for me. It allowed me to practice my programming skills and problem-solving abilities, and it also gave me a better understanding of how to use Matlab for real-world applications.中文回答：我要写实验3的Matlab程序设计实验报告。
MATLAB实验报告(1-4)
信号与系统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语言及其应用》实验报告(第三章--绘制二维和三维图形)
——绘制二维和三维图形
姓名：学号：
实验名称
绘制二维和三维图形
实验目标
利用Matlab常见函数完成二维图形的绘制和图形的标注；实现三维曲线和曲面图形的绘制。
实验要求
熟悉Matlab基本绘图函数、图形处理函数，了解三维曲线和曲面图形的绘制方法。
实验步骤
1、用Matlab基本绘图函数绘制二维图形：根据已知数据，用plot函数画出正弦函数曲线，并进行相应标注。
enon
实验内容
1．二维曲线绘图
例：精细指令实例
2.三维曲线绘图
【例】三维曲线绘图基本指令演示一：plot3
t=(0:0.02:2)*pi;x=sin(t);y=cos(t);z=cos(2*t);
plot3(x,y,z,'b-',x,y,z, 'rd')三维曲线绘图（蓝实线和红菱形）
box on
legend('链','宝石')在右上角建立图例
subplot(121);
surf(x1,y1,z1);
subplot(122);
[x2,y2,z2]=sphere (30);
surf(x2,y2,z2);
clear;clf;
z=peaks;
subplot(1,2,1);mesh(z);% 透视
hidden off
subplot(1,2,2);mesh(z);%不透视
2、用三维曲线绘图基本指令plot 3绘制三维曲线图：t=0~2pi;x=sin(t);y=cos(t);z=cos(2*t);用plot3函数画出关于x,y,z的三维曲线图，并适当加标注。
MATLAB实验三
效果图：2. 当k = 1,2,3时，在同一幅图用同时画出函数sin()y kx 的图形。
MATLAB 代码：x = 0:pi/100:2*pi; y1 = sin(x); y2 = sin(2*x); y3 = sin(3*x); plot(x,y1,x,y2,x,y3) axis([0 2*pi -1 1])legend('y = sin(x)','y = sin(2x)','y = sin(3x)')注：legend 的作用是对图中的不同曲线做标注效果图：以上画图的效果也可以用hold on实现：x = 0:pi/100:2*pi;y1 = sin(x);y2 = sin(2*x);y3 = sin(3*x);plot(x,y1);hold on;plot(x,y2);plot(x,y3);hold offsubplot(2,2,4); plot(x,y)title('y = exp(x)*sin(4x)')注：title 的作用是给图形写上标题 效果图：4. 作隐函数3330x y xy +-=（笛卡尔儿叶形线）的图形 MATLAB 代码：ezplot('x^3+y^3-3*x*y',[-2 2 -2 2])axis equal注： axis equal 的作用是把纵横坐标的刻度比例一致 效果图：5. 在极坐标系下画出心形线2(1cos )r θ=+ MATLAB 代码：theta = 0:pi/100:2*pi; r = 2*(1+cos(theta)); polar(theta,r,'r')注：参数’r’是表示用红色作图，参数与颜色之间的关系如下表所示：b 蓝色g 绿色r 红色c 青绿色m 洋红色y 黄色k 黑色效果图：5. 画出空间螺旋线cos(),sin(),x t y t z t ===的图形MATLAB 代码：t = 0:pi/100:10*pi;x = cos(t);y = sin(t);z = t;plot3(x,y,z)title('x = cos(t),y = sin(t) z = t')效果图：7. 画出MATLAB测试函数peaks的图形MATLAB代码：[X,Y,Z] = peaks;surf(X,Y,Z)shading interpcolormap(cool)注：shading可使用的参数及含义如下：faceted 网络修饰，这是缺省的方式flat 去掉黑色线条，根据小方格的值确定颜色interp 颜色整体改变，根据小方块四角的值差补过度点的值确定颜色colormap为色图设定函数，其参数及含义如下：hsv ——饱和值色图gray ——线性灰度色图hot ——暖色色图cool ——冷色色图bone ——兰色调灰色图copper ——铜色色图pink ——粉红色图prism ——光谱色图jet ——饱和值色图IIflag ——红、白、蓝交替色图六、课堂练习1. 在同一坐标系中画出函数1sin(2*3)y x=+，2sin(32)y x=+的图形，且对于不以同的曲线使用不同的线型；>> x=0:pi/100:2*pi;y1=sin(2*x+3);y2=sin(3*x+2);plot(x,y1,'r-',x,y2,'g+')legend('y1=sin(2*x+3)','y2=sin(3*x+2)')>>2. 当k 分别取1, 2, 3, 4的时候，在同一个窗口的四个子图中分别画出函数*sin(*)y t k t ，每个子图加上适当的标题；>> t = 0:pi/100:2*pi;% 画第一幅子图y = t.*sin(1*t);subplot(2,2,1);3. 画出三个两两相切的圆（假设这个三个圆的方程分别是：22221,(3)4x y x y +=-+=，22(4)9x y +-=）；>> ezplot('x^2+y^2=1',[-7 7 -7 7]);hold onezplot('(x-3)^2+y^2=4',[-7 7 -7 7]);hold onezplot('x^2+(y-4)^2=9',[-7 7 -7 7]);4. 在极坐标系下画出三叶形曲线3*sin(3*)r θ=，并尝试不同的线形和颜色； theta = 0:pi/100:2*pi;>> r=3*sin(3*theta);>> polar(theta,r,'r')theta = 0:pi/100:2*pi; r=3*sin(3*theta); polar(theta,r,'g+')5. 画出星形线332*cos ,2*sin x t y t ==的图形； >> t = 0:pi/100:10*pi;x=2*cos(t).^3;>> y=2*sin(t).^3;>> plot(x,y)6. 画出双纽线22222()3*()x y x y +=-的图形； ezplot('(x^2+y^2)^2=3*(x^2-y^2)',[-3 3 -3 3])7. 画出三维曲线图：*cos ,*sin ,x t t y t t z t ===>> t=0:pi/100:2*pi;>> x=t.*cos(t);>> y=t.*sin(t);>> z=t;>> plot3(x,y,z)8. 分别用mesh和surf函数画出马鞍面2246x yz=-+9x=-pi:pi/50:pi;y=x;[X,Y]=meshgrid(x,y); z=-X.^2./4+Y.^2./6.*9; mesh(X,Y,z)>> x=-pi:pi/50:pi;y=x;[X,Y]=meshgrid(x,y); z=-X.^2./4+Y.^2./6.*9; >> surf(X,Y,z)9. 画出球心在原点的球。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《Matlab》课内实验任务书
学生姓名:安海旭2014309010104
及学号：
学院:理学院
班级:数学141
课程名称：Matlab
题目:Matlab的优化函数和微分方程求解
指导教师
王冠讲师
姓名及职称：
朱振菊实验师
2015年10月13日
《Matlab 》课内实验
实验三 Matlab 的优化函数和微分方程求解
一、实验目的
1. 掌 求解最优化问题方法；
2. 掌握利用Matlab 求解非线性最小二乘问题的方法；
3. 掌握利用Matlab 求解非线性方程组的方法；
4. 掌握利用Matlab 求解微分方程的方法。
二、实验学时和类型
本次实验为2学时、设计性实验。
三、实验内容
1.求解无约束优化问题
)12424(min 2212
2211++++x x x x x e x ，取初值(-1,1)。
函数式文件：
Function y=mfw(x)
Y=4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1 命令窗口：
[x,f,e]=fminsearch(@mfw,[-1,1])
2.求函数22122),(y xy x y x f ++=在}25
4|),{(2
2≤+=y x y x D 上的最大值。
函数式文件：
Function y=mfw(x)
Y=2*x(1)^2+12*x(1)*x(2)+x(2)^2 Function [f,g]=mfw1(x) F=x(1)^2+4*x(2)^2-25 G=[] 命令窗口：
[x,f,e]=fmincon(@mfw,[1,1],[],[],[],[],[],[],@mfw1);
3.求解非线性方程组⎪⎩⎪⎨⎧=----=--0
1)5.0()2(0
12
221221x x x x 。
函数式文件：
Function y=mfw(x)
Y=[x(1)^2-x(2)-1;
(x(1)-2)^2-(x(2)-0.5)^2-1] 命令窗口：
[x,f,e]=fsolve(@mfw,[1,1]); 4. 变量t 和y 满足函数关系))exp()(exp(111
t k kt k
k ak y ----=
，通过下面的数据确定参数a ，k 和1k 。
函数式文件：
Function y=mfw (x,)
Y=x(1)*x(2)*(exp(-x(3)*t)-exp(-x(2)*t))/(x(2)-x(3)) T=[0.083,0.167,0.25,0.5,0.75,1,1.5,2.25,3,4,6,8,10,12]
Y=[10.9,21.1,27.3,36.4,35.5,38.4,34.8,24.2,23.6,15.7,8.2,8.3,2.2,1.8] 命令窗口：
[x,f,e] = lsqcurvefit(@mfw ,[1,-1,1],t,y) ; 求下述微分方程的数值解 5．⎩⎨
⎧=-==+=0)0('0
)0()exp(*'y x
y y x t y t x
函数式文件：
function dy=md(t,x) dy=[t*x(2)+exp(t);
x(2)-x(1)];
命令窗口：
[X,Y]=ode45(@md,[0,1],[0,0])；
00.20.40.60.81
6．1)0('0)0()cos(sin '''==+-=y y t y yy y 令X=y ，Y=y ’=X ’则y ’’=Y ’则Y ’=X*Y-sinX+cost 函数式文件：
function dy=md(t,x) dy=[x(2);
x(1)*x(2)-sin(x(1))+cos(t)];
命令窗口：
[X,Y]=ode45(@md,[0,2],[0,1]);
00.51 1.52
四、实验体会或遇到的问题：
1.将代码写在.m文件里，这样只需要该文件里的就可以而不必每次都去敲命令
2.要尽量多的设想你的问题的复杂性，当然，要一步一步复杂，这样才能保证你的程序的适用性很强。