0.618法matlab实验报告
（最新版）MATLAB实验报告
（最新版）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”窗口。
0.618黄金分割法
实验二黄金分割法
一、实验目的
1、在已知优化问题的搜索区间基础上，求解一维优化问题，加深对算法的理解。
2、验证黄金分割法的算法流程
二、实验原理
黄金分割法是建立在区间消去法原理基础上的试探方法，即在搜索区间[a,b]内适当插入两点a1、a2，并计算其函数值。
a1、a2将区间分成三段，应用函数的单峰性质，通过函数值大小的比较，删去其中一段，使搜索区间得以缩短。
然后再在保留下来的区间上做同样的处置。
如此迭代下去，使搜索区间无限缩小，从而得到极小点的数值近似解。
黄金分割法的计算步骤如下：
三、主要仪器及耗材
计算机，matlab 2010a开发平台。
四、实验内容和步骤
1、在matlab 2010a平台上按照黄金分割法算法狂徒编写算法程序
2、编译运行程序，输入参数，得到运行结果
3、分析运行结果，并判断程序的正确性
五、数据处理与分析
六、实验注意事项
注意matlab的语言规则，养成良好的变成习惯。
黄金分割法参考程序
运行过程。
MATLAB实验报告最终定稿
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 -=。
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/23．两个函数之间的操作求和(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. Taylor逼近的直观演示用Taylor 多项式逼近y = sin x．已知正弦函数的Taylor 逼近式为∑=----=≈nkk kkxxPx1121!)12()1()(sin．实验目的：利用Taylor多项式逼近y = sin x，并用图形直观的演示。
实验结果报告（含基本步骤、主要程序清单、运行结果及异常情况记录等）：1.将k从1取到5，得到相应的P = x-1/6*x^3+1/120*x^5-1/5040*x^7+1/362880*x^9;2.用MATLAB进行Taylor逼近，取x的范围是（-3.2，3.2）；程序清单如下：syms x; y = sin(x); p = x - (x^3)/6 + (x^5)/120 - (x^7)/5040 + (x^9)/362880 x1 = -3.2:0.01:3.2;ya = sin(x1);y1 = subs(p,x,x1);plot(x1,ya,'-',x1,y1)4.程序运行正常。
思考与深入：取y = sin x 的Taylor 多项式为P 的逼近效果很良好，基本接近y = sin x 的图像，不过随着k 的取值变多，逼近的效果会越来越好。
实验内容： 2. 数据插值在(,)[8,8][8,8]x y =-⨯-区域内绘制下面曲面的图形：2222sin()x y z x y+=+并比较线性、立方及样条插值的结果。
．实验目的：学会用MATLAB对函数进行线性、立方及样条插值，并比较结果。
实验结果报告（含基本步骤、主要程序清单、运行结果及异常情况记录等）：1.用MATLAB一次进行对函数的线性、立方、及样条插值，并进行算法误差分析。
2.主要程序如下：[x,y] = meshgrid(-8:1:8,-8:1:8);z = sin((x.^2 + y.^2).^0.5)./((x.^2 + y.^2).^0.5);surf(x,y,z),axis([-8,8,-8,8,-2,3])title('z的曲面图形'); %画出z的曲面图形%选较密的插值点，用默认的线性插值算法进行插值figure;[x1,y1] = meshgrid(-8:0.4:8,-8:0.4:8);z1=interp2(x,y,z,x1,y1);surf(x1,y1,z1),axis([-8,8,-8,8,-2,3])title('线性插值');%立方和样条插值figure;z1=interp2(x,y,z,x1,y1,'cubic');z2=interp2(x,y,z,x1,y1,'spline');surf(x1,y1,z1),axis([-8,8,-8,8,-2,3])title('立方插值');figure;surf(x1,y1,z2),axis([-8,8,-8,8,-2,3])title('样条插值');%算法误差的比较z = sin((x1.^2 + y1.^2).^0.5)./((x1.^2 + y1.^2).^0.5);figure;title('误差分析1');figure;surf(x1,y1,abs(z-z2)),axis([-8,8,-8,8,0,0.025]) title('误差分析2');3.结果如下：4.程序运行正常。
MATLAB实验报告
MATLAB实验报告MATLAB实验报告姓名：专业：学号：实验⼀MATLAB环境的熟悉与基本运算⼀、实验⽬的：1．熟悉MATLAB开发环境2．掌握矩阵、变量、表达式的各种基本运算⼆、实验基本知识：1.熟悉MATLAB环境:MATLAB桌⾯和命令窗⼝、命令历史窗⼝、帮助信息浏览器、⼯作空间浏览器⽂件和搜索路径浏览器。
2.掌握MATLAB常⽤命令3.MATLAB变量与运算符变量命名规则如下：（1）变量名可以由英语字母、数字和下划线组成（2）变量名应以英⽂字母开头（3）长度不⼤于31个（4）区分⼤⼩写MATLAB中设置了⼀些特殊的变量与常量，列于下表。
MATLAB运算符，通过下⾯⼏个表来说明MATLAB的各种常⽤运算符表2MATLAB算术运算符表3MATLAB关系运算符表4MATLAB逻辑运算符表5MATLAB特殊运算4.MATLAB的⼀维、⼆维数组的寻访表6⼦数组访问与赋值常⽤的相关指令格式5.MATLAB的基本运算表7两种运算指令形式和实质内涵的异同表6.MATLAB的常⽤函数表8标准数组⽣成函数表9数组操作函数三、实验内容1、学习安装MATLAB软件。
2、学习使⽤help命令，例如在命令窗⼝输⼊helpeye，然后根据帮助说明，学习使⽤指令eye（其它不会⽤的指令，依照此⽅法类推）3、学习使⽤clc、clear，观察commandwindow、commandhistory和workspace等窗⼝的变化结果。
4、初步程序的编写练习，新建M-file，保存（⾃⼰设定⽂件名，例如exerc1、exerc2、exerc3……），学习使⽤MATLAB的基本运算符、数组寻访指令、标准数组⽣成函数和数组操作函数。
注意：每⼀次M-file的修改后，都要存盘。
四、实验结果练习A：（1）helprand，然后随机⽣成⼀个2×6的数组，观察commandwindow、commandhistory和workspace等窗⼝的变化结果。
最优化方法实验
《最优化方法》实验报告实验序号：01 实验项目名称：线性规划及MATLAB应用《最优化方法》实验报告实验序号：02 实验项目名称：0.618黄金分割法的应用结果分析：根据以上结果可知，在区间[0,3]上，函数g(x)=x^3-2*x+1的最小值点在x=0.9271处，此时最小值为0。
第二题：P50 例题3.1程序：function [t,f]=golden3(a,b) %黄金分割函数的m文件t2=a+0.382*(b-a);f2=2*(t2)^2-(t2)-1;t1=a+0.618*(b-a); %按照黄金分割点赋值，更准确可直接算f1=2*(t1)^2-(t1)-1;while abs(t1-t2)>0.16; %判定是否满足精度if f1<f2a=t2;t2=t1;f2=f1;t1=a+0.618*(b-a);f1=2*(t1)^2-(t1)-1;elseb=t1;t1=t2;f1=f2;t2=a+0.382*(b-a);f2=2*(t2)^2-(t2)-1;endendt=(t1+t2)/2; %满足条件取区间中间值输出第四题：P64 T3程序：function [t,d]=newtow2(t0)t0=2.5;t=t0-(4*(t0)^3-12*(t0)^2-12*(t0)-16)/(12*(t0)^2-24*(t0)-12);k=1;T(1)=t;while abs(t-t0)>0.000005t0=t;t=t0-(4*(t0)^3-12*(t0)^2-12*(t0)-16)/(12*(t0)^2-24*(t0)-12); k=k+1;T(k)=t;endt1=t0;d=(t1)^4-4*(t1)^3-6*(t1)^2-16*(t1)+4;kTend运行结果：当x(0)=2.5当x(0)=3四．实验小结：1.通过这次实验，加深了对0.618法的理解。
2.在学习0.618法的过程中，又巩固了倒数、求解函数值等相关知识。
matlab实验报告(实验1)
开课学院及实验室：机电学院 2012年11 月16日
（2）在E盘新建一个文件夹并把其加入搜索路径。
截屏相关过程结果并保存。
（3）逐步按照以下步骤进行并记录各个步骤的结果：
1）在指令窗中输入clear指令；
5）再次在指令窗中输入clear指令后，输入
8）把cow.m暂时移动到另一个无关的文件夹，如“我的文档”。
（2）保存全部变量为数据文件Mydata.mat；在对应文件夹中找到该文件，截屏后保存结果；
（4）删除全部内存变量；
Editor/Debugger和脚本编写初步
把以下指令编写为M脚本文件，并运行：
t=0:pi/50:4*pi;
y0=exp(-t/3); y=exp(-t/3).*sin(3*t); plot(t,y,'-r',t,y0,':b',t,-y0,':b') grid
什么？此方法也是查询函数注释的有效方法，请牢记。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
function [x,k]=GSe(f,a,b,delta) % 0.618法求解单峰函数极小点 f=inline(f); N=10000; for k=1:N
m=a+0.382*(b-a); n=a+0.618*(b-a); if f(m)>f(n)
a=m; m=n; else b=n; n=m; end if abs(b-a)<delta x=0.5*(b+a);break; end end
实验报告
实验题目: 学生姓名: 学 号: 实验时间:
0.618 法的 MATLAB 实现
一．实验名称:
0.618 法求解单峰函数极小点
二．实验目的及要求:
1. 了解并熟悉 0.618 法的方法原理, 以及它的 MATLAB 实现. 2. 运用 0.618 法解单峰函数的极小点.
三．实验内容：
1. 0.618 法方法原理: 定 理 : 设 f 是 区 间 [a, b] 上 的 单 峰 函 数 , x(1) , x(2) [a, b] , 且 x(1) x(2) . 如 果
1 a1 0.382 (b1 a1), 1 a1 0.618(b1 a1). 令k 1.
② 若 bk ak L , 则 停 止 计 算 . 否 则 , 当 f (k ) f (k ) 时 , 转 步 骤 ③; 当
f (k ) f (k ) 时, 转步骤④.
bk ak L
N
f (k ) f (k )
Y
x 0.5(ak bk )
Y
置 ak1 k , bk1 bk , k1 k k1 ak1 0.618 (bk1 ak1) 计算函数值 f (k1)
置 k : k 1
2. 代码及数值算例:
0.618 法计算试探点的公式:
2. 0.618 法的算法步骤:
k ak 0.382(bk ak ), k ak 0.618(bk ak ).
①置初始区间 [a1,b1] 及精度要求 L 0 , 计算试探点 1 和 1 , 计算函数值 f (1) 和
f (1) . 计算公式是
(2) 数值算例:
min f (x)def 2x2 x 1 ,
初始区间[a1,b1] [1,1] , 精度 L 0.16 .
(i) 键入命令： >> [x,k]=GSe('2*x^2-x-1',-1,1,0.16) (ii) 运行结果： x=
0.2229
k=
6
五．总结:
0.618 法(黄金分割法)适用于单峰函数, 故应先确定目标函数的单峰区间, 方可进行迭代 计算. 但单峰区间不是很明显就能确定, 故可用进退法寻找并确定单峰区间.
③ 置 ak1 k , bk1 bk , k1 k , k1 ak1 0.618 (bk1 ak1) , 计 算 函 数 值
f (k1) , 转步骤⑤.
④ 置 ak1 ak , bk1 k , k1 k , k1 ak1 0.382 (bk1 ak1) , 计 算 函 数 值 f (k1) , 转步骤⑤.
如果 f (x(1) ) f (x(2) ) , 则 x [x(1) ,b] ;
如果 f (x(1) ) f (x(2) ) , 则 x [a，据上述定理, 通过取试探点使包含极小点的区间(不确定区 间)不断缩短, 当区间长度小到一定程度时, 区间上各点的函数值均接近极小值, 因此任意 一点都可作为极小点的近似.
f (x(1) ) f (x(2) ) , 则对每一个 x [a, x(1) ] , 有 f (x) f (x(2) ) ; 如果 f (x(1) ) f (x(2) ) ,
则对每一个 x [x(2), b] , 有 f (x) f (x(1) ) .
根据上述定理, 只需选择两个试探点, 就可将包含极小点的区间缩短. 事实上, 必有
⑤置 k : k 1, 返回步骤②.
四．实验流程图及其 MATLAB 实现:
1. 流程图：
开始
置初始区间[a1,b1] 及精度要求 L 0
计算试探点 1 和 1 计算函数值 f (1) 和 f (1)
N
置 ak1 ak , bk1 k , k1 k k1 ak1 0.382 (bk1 ak1) 计算函数值 f (k1)