[工学]第二章 Matlab软件的数值计算方法2.ppt

二、函数求极值 1.fmin的命令格式 X=fmin(‘Fun’,x1,x2) 举例35 2.fmins的命令格式 X=fmins(‘Fun’,x0) 举例36 三、函数求零点 命令格式： X=fzero(‘Fun’,x0) 举例37 举例39 表2－6 数值分析的部分MATLAB命令函数 补充：符号计算 ！1、符号对象和符号表达式 在MATLAB中，数值和数值变量用于数值的存储和各种数值计算。而符号常数、符号变量、符号函数、符号操作等则是用来形成符号表达式，严格按照代数、微积分等课程中的规则公式进行运算，并尽可能给出解析表达结果。 1.1 符号对象的生成和使用 数值计算中，数值表达式所用的变量必须事先被赋过值，否则无法计算。Symbolic Math Toolbox 2.1版沿用这种模式，规定：在进行符号运算时， 首先要定义基本的符号对象（可以是常数、变量、表达式），然后利用这些基本符号对象去构成新的表达式，进而从事所需的符号运算。 运算中，凡是由包含符号对象的表达式所生成的衍生对象也都是符号对象。 定义基本符号对象的指令有两个：sym, syms。 f = sym( arg ) f = sym( argn,flagn) argv =sym(‘argv’,flagv) syms(‘argv1’,’argv2’,’argvk’) syms argv1 argv2 argvk 【例1】符号常数形成中的差异 a1=[1/3,pi/7,sqrt(5),pi+sqrt(5)] a2=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) a3=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) a23=a2-a3 a1 = 0.3333 0.4488 2.2361 5.3777 a2 = [ 1/3, pi/7, sqrt(5), 6054707603575008*2^(-50)] a3 = [ 1/3, pi/7, sqrt(5), pi+sqrt(5)] a23 = [ 0, 0, 0, 189209612611719/35184372088832-pi-5^(1/2)] 【例2】把字符表达式转换为符号变量 y=sym(2*sin(x)*cos(x)) y=simple(y) y = 2*sin(x)*cos(x) y = sin(2*x) 说明：符号运算中，如事先没对表达式中的独立符号变量进行定义，软件会自动检查哪些字符是软件的符号函数，哪些是变量符号，且总把在英文字母表中离x最近的字母认为独立符号变量。 【例3】用符号计算验证三角等式。 syms fai1 fai2 y=simple(sin(fai1)*cos(fai2)-cos(fai1)*sin(fai2)) y = sin(fai1-fai2) 【例4】求矩阵的行列式值、逆和特征根 syms a11 a12 a21 a22 A=[a11,a12;a21,a22] DA=det(A),IA=inv(A),EA=eig(A) A = [ a11, a12] [ a21, a22] DA = a11*a22-a12*a21 IA = [ a22/(a11*a22-a12*a21), -a12/(a11*a22-a12*a21)] [ -a21/(a11*a22-a12*a21), a11/(a11*a22-a12*a21)] EA = [ 1/2*a11+1/2*a22+1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)] [ 1/2*a11+1/2*a22-1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)] 1.2 符号计算中的算符和基本函数 由于matlab中采用了重载技术，使得用来构成符号计算表达式的算符和基本函数，无论在形状、名称上，还是在使用方法上，都与数值计算中的算符和基本函数几乎完全相同。 1）基本运算符 + - * \ / 矩阵 .* ./ .\ 元素对元素 ‘ .’ 矩阵 2）关系运算符 符号对象中没有大于、大于等于、小于、小于等于的概念，而只有是否等于的概念。 算符＝＝、～＝分别对算符两边的对象进行相等、不等的比较。当事实为真时，比较结果用1表示，假用0表示。 3）三角函数、双曲函数及它们的反函数 除atan2只用于数值计算外，其余的三角函数（如sin）、双曲函数（如cos