第3章 MATLAB符号运算317510155.pptVIP

第3章 MATLAB符号运算 3.1 符号运算基础 3.1.2 建立符号对象 【例3.1】符号变量（常量）与数值变量（常量）的区别 3.1.3 符号表示式的运算 2. 函数运算 2）符号表达式系数的提取coeffs(s[,x])函数：该函数返回多项式中按指定变量升幂顺序排列的系数，若没有指定变量，则返回所有项的常系数，且按离字符“X”近原则确定主变量。 3）通分 4）反函数 5）复合函数 6）书写格式美化 7）符号表达式变量的确定 3.3 微积分 【例3.23】 2. 导数 3. 积分 【例3.26】已知f(x)=ax2+bx+c，求f(x)的积分 3.4 方程求解 （2） 非线性方程的数值求解 fzero 的另外一种调用方式： 【例3.33】 2. 常微分方程 clear syms a x limit(1/x,x,0) ans= NaN limit(1/x,x,0,left) ans= -Inf limit(1/x,x,0,right) ans= Inf limit(((x+a)/(x-a))^x,inf) ans= exp(2*a) 求函数的导数（即微分）最基本的方法莫过于极限方法： 【例3.24】 clear syms t x limit((cos(x+t)-cos(x))/t,t,0) ans= -sin(x) MATLAB 提供了专门求导数的diff函数。 格式1：diff(f) 功能：表达式f对默认自变量的一阶导数； 格式2： diff(f, t) 功能：求表达式f对自变量t的一阶导数； 格式3： diff(f,n) 功能：求表达式f对默认自变量的n阶导数； 格式4： diff(f,t,n) 功能：求表达式f对自变量t的n阶导数。 【例3.25】 clear syms a b c x f=sym(a*x^2+b*x+c) f1=diff(f) %对默认自变量求导数 f1= 2*a*x+b diff(f,2) %求X二阶导数，等价于diff(f1) ans= 2*a 格式1：int(f) 功能：求表达式f对默认自变量的积分值； 格式2： int(f, t) 功能：求表达式f对自变量t的不定积分值； 格式3： int(f, a, b) 功能：求表达式f对默认自变量的定积分值，积分区间为[a,b]； 格式4： int(f, t, a, b) 功能：求表达式f对自变量t的定积分值，积分区间为[a,b] clear syms a b c x f=sym(a*x^2+b*x+c) int(f) %表达式f(x)的不定积分，自变量为x ans= 1/3*a*x^3+1/2*b*x^2+c*x int(f,x,0,2) %表达式f在(0,2)的定积分，自变量为x ans= 8/3*a+2*b+2*c int(f,a) %表达式f(x)的不定积分，自变量为a ans= 1/2*a^2*x^2+b*x*a+c*a int(int(f,a),x) ans= 1/6*a^2*x^3+1/2*b*a*x^2+c*a*x 【例3.27】求积分 【例3.28】求积分 clear int(sqrt(x)/(1+x)^2,1,inf) ans= 1/4*pi+1/2 clear int(int(x*exp(-x*y), x),y) ans= 1/y*exp(-x*y) 1. 代数方程 代数方程的求解由函数solve实现： 格式1： solve(f) 功能：求方程式f关于默认自变量的解 格式2： s=solve(f,v) 功能：求方程关于指定自变量的解； 说明：1）f 可以是用字符串表示的方程，或符号表达式； 2）若 f 中不含等号，则表示解方程 f=0。 【例3.29】解方程 x3-3x+1=0 syms x; f=x^3-3*x+1; s=solve(f,x) s=solve(x^3-3*x+1,x) %s=solve(x^3-3*x+1=0,x) 注意：roots(p)是求多项式的所有零点，p 是多项式系数向量。 【例3.30】求下列代数方程的根 1） ax2+bx+c=0 2）1+x=sin(x) syms a b c x f=sym(a*x*x+b*x+c=0) solve(f) ans = [ 1/2/a*(-b+(b^2-4*c*a)^(