matlab 二元函数 定义,实验五用matlab求二元函数及极值.doc
实验五用matlab求二元函数及极值
实验五?? 用matlab求二元函数的极值
?
1.计算二元函数的极值
对于二元函数的极值问题,根据二元函数极值的必要和充分条件,可分为以下几个步骤:
步骤1.定义二元函数.
步骤2.求解方程组,得到驻点.
步骤3.对于每一个驻点,求出二阶偏导数
步骤4. 对于每一个驻点,计算判别式,如果,则该驻点是极值点,当为极小值, 为极大值;如果,需进一步判断此驻点是否为极值点; 如果则该驻点不是极值点.
2.计算二元函数在区域D内的最大值和最小值
设函数在有界区域上连续,则在上必定有最大值和最小值。求在上的最大值和最小值的一般步骤为:
步骤1. 计算在内所有驻点处的函数值;
步骤2. 计算在的各个边界线上的最大值和最小值;
步骤3. 将上述各函数值进行比较,最终确定出在内的最大值和最小值。
3.函数求偏导数的MATLAB命令
MATLAB中主要用diff求函数的偏导数,用jacobian求Jacobian矩阵。
?
?
diff(f,x,n)? 求函数f关于自变量x的n阶导数。
jacobian(f,x) 求向量函数f关于自变量x(x也为向量)的jacobian矩阵。可以用help diff, help jacobian查阅有关这些命令的详细信息
例1? 求函数的极值点和极值.
首先用diff命令求z关于x,y的偏导数
>>clear;? syms x y;
>>z=x^4-8*x*y+2*y^2-3;
>>diff(z,x)
>>diff(z,y)
结果为
ans =4*x^3-8*y
??? ans =-8*x+4*y
即再求解方程,求得各驻点的坐标。一般方程组的符号解用solve命令,当方程组不存在符号解时,solve将给出数值解。求解方程的MATLAB代码为:
>>clear;
>>[x,y]=solve('4*x^3-8*y=0','-8*x+4*y=0','x','y')
结果有三个驻点,分别是P(-2,-4),Q(0,0),R(2,4).下面再求判别式中的二阶偏导数:
>>clear;? syms x y;
>>z=x^4-8*x*y+2*y^2-3;
>>A=diff(z,x,2)
>>B=diff(diff(z,x),y)
>>C=diff(z,y,2)
结果为
A=2*x^2
B =-8
??? C =4
由判别法可知和都是函数的极小值点,而点Q(0,0)不是极值点,实际上,和是函数的最小值点。当然,我们可以通过画函数图形来观测极值点与鞍点。
>>clear;
>>x=-5:0.2:5;? y=-5:0.2:5;
>>[X,Y]=meshgrid(x,y);
>>Z=X.^4-8*X.*Y+2*Y.^2-3;
>>mesh(X,Y,Z)
>>xlabel('x'),ylabel('y'),zlabel('z')
结果如图16.5.1
?
图16.5.1? 函数曲面图
可见在图6.1中不容易观测极值点,这是因为z的取值范围为[-500,100],是一幅远景图,局部信息丢失较多,观测不到图像细节.可以通过画等值线来观测极值.
>>contour(X,Y,Z, 600)
>>xlabel('x'),ylabel('y')
结果如图16.5.2
图16.5.2? 等值线图
由图16.5.2可见,随着图形灰度的逐渐变浅,函数值逐渐减小,图形中有两个明显的极小值点和.根据提梯度与等高线之间的关系,梯度的方向是等高线的法方向,且指向函数增加的方向.由此可知,极值点应该有等高线环绕,而点周围没有等高线环绕,不是极值点,是鞍点.
例2 求函数在条件下的极值..构造Lagrange函数
求Lagrange函数的自由极值.先求关于的一阶偏导数
>>clear; syms x y k
>>l=x*y+k*(x+y-1);
>>diff(l,x)
>>diff(l,y)
>>diff(l,k)
得再解方程
>>clear; syms x y k
>>[x,y,k]=solve('y+k=0','x+k=0','x+y-1=0','x','y','k')
得进过判断,此点为函数的极大值点,此时函数达到最大值.
?
例3 抛物面被平面截成一个椭圆,求这个椭圆到原点的最长与最短距离.
这个问题实际上就是求函数
在条件及下的最大值和最小值问题.构造Lagrange函数
求Lagrange函数的自由极值.先求关于的一阶偏导数
>>clear; syms x y z u v
>>l=x^2+y^2+z^2+u*(x^2+y^2-z)+v*(x+y+z-1);
>>diff(l,x)
>>diff(l,y)
>>diff(l,z)
>>diff(l,u)
>>diff(l,v)
得
再解方程
>>clear;
>>[x,y,z,u,v]=solv
matlab 二元函数 定义,实验五用matlab求二元函数及极值.doc相关推荐
- matlab 信号的频谱分析,实验五基于Matlab的信号频谱分析(复杂)
实验五基于Matlab的信号频谱分析(复杂) 本次实验注意:<实验五MALTAB基础知识(简单)> <实验五 基于Matlab的信号频谱分析(复杂)> 选作一个即可 实验五 基 ...
- matlab系统的根轨迹,实验五 利用MATLAB绘制系统根轨迹
<实验五 利用MATLAB绘制系统根轨迹>由会员分享,可在线阅读,更多相关<实验五 利用MATLAB绘制系统根轨迹(6页珍藏版)>请在人人文库网上搜索. 1.实验五 利用MAT ...
- 利用matlab实现卷积实验报告,实验五 使用matlab实现卷积的运算
实验五 使用matlab实现卷积的运算 一 实验目的 1. 2. 二 实验内容 学习MATLAB语言的编程方法及熟悉MATLAB指令: 深刻理解卷积运算,利用离散卷积实现连续卷积运算: 1. 完成f1 ...
- matlab复杂周期信号类建立,实验五 基于Matlab的信号频谱分析(复杂)
本次实验注意:<实验五MALTAB基础知识(简单)> <实验五 基于Matlab的信号频谱分析(复杂)> 选作一个即可 实验五 基于Matlab的信号频谱分析 (一) 实验目的 ...
- 函数定义涉及的三要素C语言,(函数定义与三要素答案.doc
(函数定义与三要素答案 函数定义及三要素 1.与函数y=10lg(2x-1)的图象相同的函数是( ) A.y= B.y=2x-1 C.y= D.y= [答案] B [解析] y=10lg(2x-1)= ...
- matlab实验5函数文件,实验五M文件和MATLAB设计技术总结.doc
实验五 M文件和MATLAB程序设计 一.实验目的matlab作为一种高级计算机语言,不仅可以命令行方式完成操作,也具有数据结构.控制流.输入输出等能力,本次实验通过熟悉和掌握m文件的建立与使用方法, ...
- 【通信原理】实验五 基于Matlab的2ASK和2FSK调制解调
目录 一.实验目的 二.实验器材 三.实验原理 1.二进制振幅键控(2ASK) 2.二进制频移键控(2FSK) 四.示例演示 1.2ASK的模拟调制程序如下: 2.2ASK的开关键控法,调制程序如下: ...
- matlab绘制系统函数频率特性,实验七用MATLAB进行系统频率特性曲线绘制
实验七用MATLAB进行系统频率特性曲线绘制 自动控制原理 课程验证性实验报告 实验名称 用MATLAB进行系统频率特性曲线绘制 实验时间 2013年 05月30日 学生姓名 实验地点 070312 ...
- matlab 2009a使用教程,实验一 安装MATLAB R2009a软件及其简单操作
实验一安装MATLAB R2009a软件及其简单操作 一.实验目的 1.掌握MATLAB R2009a软件的安装. 2.熟悉启动和退出matlab的方法. 3.熟悉matlab命令窗口的组成. 4.掌 ...
- 利用MATLAB进行系统时域分析,实验二 利用matlab进行系统的时域分析
实验二 利用matlab进行系统的时域分析 实验二 利用MATLAB进行系统的时域分析 1.实验目的 在理论学习的基础上,通过本实验熟悉LTI连续时间系统的时域分析方法, 熟悉系统的零输入响应.零状态 ...
最新文章
- POJ2299 Ultra-QuickSort
- JSP 和 Servlet 区别
- 5、ORDER BY:对查询结果排序
- PAT甲级 -- 1090 Highest Price in Supply Chain (25 分)
- element的分页怎么改中文配置_vue想改变element分页样式
- “约见”面试官系列之常见面试题之第四十五篇CSS优先级(建议收藏)
- 「北京」京东 JD.COM 招聘中/高级前端工程师
- 2_C语言中的数据类型 (四)整数与无符号数
- 【HTML+PHP+SQL】综合实验:一个简单丑陋的论坛网站(实现用户注册、登录、注销、修改信息、留言等功能)
- [HTTP权威指南]2.HTTP连接管理
- oracle认证视频教程,Oracle认证专家视频教程-OCP全套教程之学习笔记-RMAN-CATALOG
- pandownload复活版下载 | 百度网盘下载速度慢的终极解决方案
- (十进制快速幂+矩阵优化)BZOJ 3240 3240: [Noi2013]矩阵游戏
- 经典黑白搭配 现代简约风格设计美学精神
- 阿里云ECS简易部署教程
- matplotlib学习笔记(六)
- python人脸识别门禁_树莓派人脸识别门禁系统图文教程
- AI正在模仿人类大脑!2021年10篇顶会论文:大脑也在「无监督」学习
- 360全景拍摄教程:鱼眼全景图拍摄方法
- 牛客网-八皇后【搜就完事】