实验3Matlab符号运算及求函数极值汇编.pdfVIP

学号： 班级： 姓名： 实验 3 Matlab 符号运算及求函数极值 一、实验目的和要求 掌握用 Matlab 软件进行符号运算以及求函数的极值。 二、实验环境 Windows系列操作系统， Matlab 软件。 三、实验内容 1. 用 MATLAB进行符号运算； 2. 编程求函数的极值。 四、实验步骤 3. 开启软件平台—— Matlab ，开启 Matlab 编辑窗口； 4. 根据求解步骤编写 M文件； 5. 保存文件并运行； 6. 观察运行结果 ( 数值或图形 ) ； 7. 根据观察到的结果和体会写出实验报告。 五、示例 1．计算一元函数的极值 2 3x 4x 4 例 1 求 y 2 的极值 x x 1 解 首先建立函数关系： syms x y=(3*x^2+4*x+4)/( x^2+x+1); 然后求函数的驻点： dy=diff(y); xz=solve(dy) xz= [0] [-2] 知道函数有两个驻点 x =0 和 x =-2 ， 1 2 接下来我们通过考察函数的图形，则它的极值情况和许多其它特性是一目了然的。而借助 MATLAB 的作图功能，我们很容易做到这一点。 例 2 画出上例中函数的图形 学号： 班级： 姓名： 解 syms x y=(3*x^2+4*x+4)/( x^2+x+1); 得到如下图形 ezplot(y) 2．计算二元函数的极值 MATLAB 中主要用 diff 求函数的偏导数 ,用 jacobian 求 Jacobian 矩阵。 diff(f,x,n) 求函数 f 关于自变量 x 的 n 阶导数。 jacobian(f,x) 求向量函数 f 关于自变量 x(x 也为向量 )的 jacobian 矩阵。 可以用 help diff, help jacobian 查阅有关这些命令的详细信息 4 2 例 1 求函数 z x 8xy 2 y 3 的极值点和极值 . 首先用 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 z 3 z 即 4x 8 y , 8x 4 y 再求解方程，求得各驻点的坐标。一般方程组的符号解用 x y solve 命令，当方程组不存在符号解时， solve 将给出数值解。求解方程的 MA TLAB 代码为： 学号： 班级： 姓名： 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; x=-5:0.2:5; y=-5:0.2:5; [X,Y]=meshgrid(x,y); Z=X.^4-8