中国石油大学MATLAB实验报告一

matlab实习心得matlab实习心得通过实习，让我们发现自己的不足，认识自己的缺点，在工作中慢慢改进，积极进取，去实现自己的目标，为自己的将来奋斗。

下面带来matlab实习心得，欢迎阅览!matlab实习心得【1】Matlab是一个强大的数学工具,它的应用广泛,涉及到各个领域.它使用起来十分方便,不用麻烦去定义变量.它的绘图能力很强,甚至可以模拟出三维视图.矩阵是它应用的核心,许多工程繁琐的运算都需要靠矩阵来化简,这正是它的生命力所在.但是,他的函数很多,开始学时记的比较痛苦,我已经深深感觉到了.不过看多了也就熟了,感觉和学五笔差不多.它的语法简单,像我学过C语言的学起来还是蛮容易的.它的数组定义十分符合自然,是从1开始的,数组元素的调用也很接近数学的表达.此外,函数的名字也很符合英文规则,反正我用得很开心就是了.通过学习matlab,我又一次锻炼了自己的思维.它学起来得心应手也让我明白了学习一门语言(c语言)对学习其他语言的帮助指导作用.同时,它也加强了我理论联系实际的能力.这是一个专业课的基础工具,学好它是必要的.matlab实习心得【2】从中国石油大学来到大庆油田勘探开发研究院实习，刚刚来到这里时我内心充满了激动和喜悦，对这个有着辉煌成绩的科研单位充满了期待。

到了这里没多久我就参加培训中心组织的入院教育，对研究院的辉煌过去有了更深刻的了解，那时我就下定决心在研究院好好实习，多学习一些现场经验，为以后的工作打下良好的基础。

一个月的入院培训结束后，经过自己申请和中心安排我被分到了地震资料处理一室，实习期间地震一室的副主任陈志德是我的现场导师，陈主任有着丰富的理论知识和现场经验，能够成为他的学生我感到非常的高兴。

最初他把我分到了高分辨率处理解释项目组，工作内容主要是跟着项目组成员用landmark做一些解释性的工作，期间我跟着张晶师傅学习到了如何加载地震数据和测井数据;如何应用测井数据做了合成地震记录，大学生创业贷款，另外我还做了一些地层解释方面的工作。

目录一、基础题 (2)二、绘图题 (3)三、句柄图形和图形用户界面制作 (9)四、Simulink仿真题 (13)五、专题实验总结 (17)六、参考文献 (17)一．基础题实验目的：1、掌握数组的创建和寻访；2、掌握MA TLAB 数组的运算；3、熟悉MA TLAB 关系操作和逻辑操作；4、掌握函数的调用。

实验要求：1、熟练掌握数组运算；2、熟悉MA TLAB 的各种数据类型及函数调用；3、初步了解M 文件及其编写、调试与运行。

实验内容：某公司投资2000万元建成一条生产线。

投产后，在时刻t 的追加成本和追加收益分别为G(t)= (百万元/年), H(t)= (百万元/年)。

试确定该生产线在合适何时停产可获最大利润？最大利润是多少？解：构造函数f(t)=H(t)-G(t)=13-t-3t 2/3=0 ;令t 1/3=x,则f(t)=-t 3-3t 2+13 可得矩阵P=[-1,-3,0,13]求最佳生产时间的源程序如下：p=[-1,-3,0,13]; x=roots(p); t=x.^3运行结果如下： t =3.6768 +21.4316i 3.6768 -21.4316i4.6465考虑到实际情况，显然两个虚数根应该舍掉。

所以将t=4.6425带入，求积分。

代码：t=4.6465; x=0:0.01:t;y=13-x-3*x.^(2/3); trapz(x,y)运行结果： ans =26.32083/225tt ++3/218t-结论：比较以上三组数据，可知最佳生产时间t=4.6465年，可获得的最大收益为 26.3208（百万元/年）。

减去20（百万元）投资，可得最终利润为6.3208（百万元）。

二．作图题实验目的：1、进一步熟悉M 文件调试过程；2、熟练掌握MA TLAB 二维、三维图形的绘制；3、掌握图形的修饰； 实验要求：1、进一步熟悉和掌握MA TLAB 的编程及调试；2、掌握二维、三维图形的绘制；3、掌握图形交互指令的使用；实验内容：1、二维绘图： （1）函数)sin(2x ey x-=,求：a. 绘制[0 ,8]区间内的图形，加注x,y 轴及图形名称；b. 最大值和最小值；c. 零值；求解：a、命令：ezplot('2*exp(-x)*sin(x)',[0,8]) %绘图title(f) , xlabel('x') ,ylabel('y') %加轴名及图名图像:b、先求最小值：从图像上可以看出，在[3,4]之间有最小值，所以命令：[xmin,fmin]=fminbnd(f,3,4)Hold onPlot(xmin,fmin)运行结果：xmin =3.9270 fmin =-0.0279增加一句标注： text(xmin,fmin,'(3.9270,-0.0279)')求最小值求最大值分析：要求最大值，即为求y=-f(x)的最小值，从图像上可以看出，在[0,2]之间有最大值命令：f=inline('-2*exp(-x)*sin(x)') ezplot(f,[0,8][xmin,fmin]=fminbnd(f,0,2) hold onplot(xmin,fmin,'*') Hold onPlot(xmin,fmin)运行结果： xmin = 0.7854 fmin = -0.6448增加标注：text(xmin,fmin,'(0.7854,-0.6448)')（如上图）综上可知：fmax= 0.6448 fmin= -0.0279c 、为了便于观察，首先增加一条直线y=0 命令：hold onezplot('0',[0,8])图像：0123456780.10.20.30.40.50.60.7x2 exp(-x) sin(x)由图像可以看出，在x=0,x=3,x=6附近有零点，所以命令及运行结果如下：fzero(f,0) %求在x=0附近的零点 ans = 0fzero(f,3) %求在x=3附近的零点 ans =3.1416fzero(f,6) %求在x=6附近的零点 ans =6.2832所以零点有三个分别为：x 1=0 x 2=3.1416 x 3=6.2832（2）在同一个画面上建立几个坐标系, 用subplot(m,n,p)命令；把一个画面分成m×n 个图形区域, p 代表当前的区域号，在每个区域中分别画一个图（函数为：;cos sin 2;cos ;sin x x u x z x y ===xxv cos sin =）； 命令：>> subplot(2,2,1)>> ezplot('sin(x)',[-2*pi,2*pi]) >> subplot(2,2,2)>> ezplot('cos(x)',[-2*pi,2*pi]) >> subplot(2,2,3)>> ezplot('2*sin(x)*cos(x)',[-2*pi,2*pi]) >> subplot(2,2,4)>> ezplot('sin(x)/cos(x)',[-2*pi,2*pi])图像：2、三维绘图：（1）绘制[-8 8]区间内函数2222)sin(),(yx y x y x f ++=的三维网格曲面图，给该图加上标题；（2）绘制[-8 8]区间内单叶双曲面 125169222=-+z y x 的三维网格曲面图，给该图加上标题；(3) 求函数 f(x,y) = 3x 2+10y 2+3xy-3x +2y 在原点附近的一个极小值点和极小值；作函数 f(x,y)在|x|<2, |y|<1内的图；解：（1）命令：ezsurf('sin(sqrt(x^2+y^2))/sqrt(x^2+y^2)',[-8,8])title('f=sin(sqrt(x^2+y^2))/sqrt(x^2+y^2)')图像：（2）命令：xa=-8:0.2:8; ya=xa;[x,y]=meshgrid(xa,ya); a=3;b=4;c=5;z1=sqrt(c.^2*(x.^2/(a.^2)+y.^2/(b.^2)-1)); z2=-sqrt(c.^2*(x.^2/(a.^2)+y.^2/(b.^2)-1)); surf(x,y,real(z1)) hold on ;surf(x,y,real(z2)) shading interp图像：(3)求极值点命令：>> f=inline('3*x(1)^2+10*x(2)^2+3*x(1)*x(2)-3*x(1)+2*x(2)')f =Inline function:f(x) = 3*x(1)^2+10*x(2)^2+3*x(1)*x(2)-3*x(1)+2*x(2)>> [x,fval]=fminsearch(f,[0,0])运行结果：x = 0.5946 -0.1892fval = -1.0811所以原点附近的极小值点为（0.5946，-0.1892），极小值为：-1.0811 绘图命令：>> x=-2:0.1:2;>> y=-1:0.1:1;>> [xx,yy]=meshgrid(x,y);>> zz=3.*xx.^2+10.*yy.^2+3.*xx.*yy-3.*xx+2.*yy>> surf(xx,yy,zz)图像：三、句柄图形和图形用户界面制作实验目的：1、熟悉句柄图形体系的对象树结构；2、熟练掌握句柄图形体系的对象属性的创建、设置、查询；3、熟练掌握句柄的获取；4、熟练掌握图形用户界面（GUI）的制作；实验要求：1、会获取和显示图形对象的句柄；2、会设置菜单和子菜单；3、会设置用户控件；实验内容：利用底层绘图指令绘制一条余弦曲线。

实验报告课程名称：可视化计算机语言实验项目名称：matlab程序设计学院：信息工程学院专业：电子信息工程指导教师：报告人：学号：班级：实验时间：实验报告提交时间：教务部制一、实验目的与要求：1、熟练掌握matlab的程序流程控制结构..2、熟悉掌握M文件的结构和函数调用..3、掌握内联函数和函数句柄的使用..4、了解程序性能剖析窗口..二、内容和步骤：MATLAB的语法规则简洁;编程效率高;作为一个完整的程序语言;MATLAB也有各种程序流程控制;文件格式和函数调用的规则;通过对函数的调用就能够组成庞大的程序;完成复杂的功能..1.使用程序流程控制Fibonacci数列的各元素为：1、1、2、3、5、8、...满足一下关系F1=1 F2=1 F n=F n-1+F n-2用M函数文件实现;数列的元素个数为输入变量..（1）按M函数文件格式创建文件开头..function f=shiyan0501n%SHIYAN0501%Fibonacci数列%n 元素个数%f 构成Fibonacci数列向量%%copyright 2015-11-25（2）用while循环实现程序功能f1=1;f2=1;i=2;while i<=nfi+1=fi-1+fi;i=i+1;End运行后;输入参数10;运行结果如下所示：（3）使用for循环实验for i=2:nfi+1=fi-1+fi;end运行结果如下所示：（4）当某个元素大于50时;退出循环结构;程序修改如下：for i=2:nif fi>10breakelsefi+1=fi-1+fi;endend当某个元素大于50;程序便退出循环结构当n=10;最后一个元素刚好是大于50的;正好退出循环结构当n=12时;第十个元素对应的fi大于50;于是退出循环结构;不再运行f11;f12了..（5）将该.m文件生成P码文件>> pcode shiyan0501将shiyan0501.m删除;重新运行该文件夹;结果如下所示：练习：将该M文件函数改为M脚本文件;将数列元素个数通过键盘输入;程序应该如何修改①将M文件函数另存为M脚本文件；②在函数前加n=input'please input a num:'便可通过键盘键入文件f=mf;endk=f1/2^2n1f^22n1+1;end调用程序shiyan0502;运行结果如下所示：>> y=shiyan05020.7y =0.7754运行结果与使用子函数factorial的运行结果一致..（4）使用程序性能剖析..选择菜单‘view’->‘profile’命令；或使用在命令窗口输入‘profile viewer’命令都可以打开程序性能剖析窗口..在程序性能剖析窗口的‘命令输入栏’中输入需要剖析的命令;‘y=shiyan05020.7’;然后单击‘start profiling’按钮;查看剖析报告..（5）程序的调试..当有多个函数调用时;由于函数变量的工作空间是独立的;被调用的函数执行结束后变量消失;因此调试时要使用matlab调试器查看运行过程中的变量值..①设置断点..在需要查看的程序的地方设置断点;>> shiyan05020.7K>> yy =0.7000K>> x=0.5x =0.5000在K>>输入x=0.5即可将x的值从0.7修改为0.5.去除断点;查看结果：ans =0.7236结果与x值为0.7时不一致;当x=0.7时;y =0.7754;由此可知;x的值已被修改..②单步运行>> shiyan05020.79 n=n+1;K>>按单步运行键：可看到箭头不断的移动：6使用函数句柄..在命令窗口使用函数句柄调用函数..>> h_shiyan0502=shiyan0502h_shiyan0502 =shiyan0502>> y=fevalh_shiyan0502;0.5y =0.5236（7）使用全局变量..Matlabe的编程不提倡使用全局变量;本例中的程序主要是为了查看全局变量的概念..将n 作为全局变量;子函数factorial不修改;子函数cal程序和主函数shiyan0502修改如下：function y=shiyan0502x%shiyan0502 arcsinxglobal n;n=1;if absx<1y=x;while caln>0.0001y=y+calnx^2n+1;n=n+1;endelsedisp'输入错误';y=0;returnendfunction k=caln1global nfor m=1:nk=factorial2n/2^2nfactorialn^222n+1;Endglobal为设置的全局变量;子函数没有输入变量;而用全局变量n传递..在全局变量前设置断点;然后运行程序;当程序运行到断点处停止;接着不断使用单步运行调试;结果如下所示：|||练习：使用单步运行调试;查看全局变量n的变化;并在工作空间查看n；3、利用泛函命令实现数值分析①创建函数shiyan0503实现上述表达式关系..function y=shiyan0503t%shiyan0503 y=sint.^2.expat-babsta=0.1;b=0.5;y=sint.^2.expat-babst;②查看该函数的输出波形;如下图所示：注：1、报告内的项目或内容设置;可根据实际情况加以调整和补充..2、教师批改学生实验报告时间应在学生提交实验报告时间后10日内..。

MATLAB实验报告机电工程学院姓名：李祖达学号： 5901111120专业：机械设计制造及其自动化目录实验一熟悉MATLAB环境认识MATLAB (3)实验二 MATLAB运算础 (6)实验三 MATLAB矩阵分析和理 (8)实验四求余弦的积分并绘出像 (11)实验五、六用matlab求解多项式并用plot绘制函数图象（常微分程） (13)实验七函数件 (16)实验八、九 MATLAB程序设计（循环、择） (19)实验十采用SIMULINK的系统真 (24)实验十一菜单设计 (27)实验一熟悉MATLAB 环境认识MATLAB一、 实验目的熟悉matlab 的安装与启动；熟悉matlab 用户界面；熟悉matlab 功能、建模元素；熟悉matlab 优化建模过程。

二、 实验设备与仪器1.微机2.matlab 仿真软件 三、 实验步骤1. 了解matlab 的硬件和软件必备环境；2. 启动matlab ；3. 熟悉标题栏，菜单栏，工具栏，元素选择窗口，状态栏，控制栏以及系统布局区；4. 学习优化建模过程。

四、 实验报告要求1. 写出matlab 系统界面的各个构成；以及系统布局区的组成；以及每一部分的功能；2. 优化建模过程应用举例 五、实验内容（一）、Matlab 操作界面1. 命令窗口（command window ）2. 命令历史窗口（command history ）3. 工作空间管理窗口（workspace ）4. 当前路径窗口（current directory ）（二）、实现下列优化建模过程1、简单矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A 的输入步骤。

指令：A=[1 2 3 ;4 5 6 ;7 8 9] A =1 2 3 4 5 6 7 8 92、矩阵的分行输入。

A=[1,2,3 4,5,6 7,8,9] 指令： >> A=[1,2,3 4,5,67,8,9]A =1 2 3 4 5 6 7 8 9 3、指令的续行输入S=1-1/2+1/3-1/4+ ... 1/5-1/6+1/7-1/8 指令：S=1-1/2+1/3-1/4+ ... 1/5-1/6+1/7-1/8 S = 0.63454、画出衰减振荡曲线t e y t 3sin 3-=及其它的包络线30t e y -=。