第五讲MATLAB的符号计算功能.pptVIP

fplot(function,[a,b]) %画出函数f(x)在[a,b]上的图形 例 画出函数在[0，8]上的图形 f=2*exp(-x)*sin(x); %定义函数 fplot(f,[0,8]) 五、符号函数画图 fzero(fun,[a,b],[ ],options) 在求零点之前先要确定一个根的隔离区间，也就是说，[a,b]区间内有且仅有函数的一个零点，这里采用的是对分法求零点。 fzero(fun,x0) %用牛顿法求零点。 例 求函数的零点 fx=x^3-x+1; a=fzero(fx,-1) 六、求符号函数的零点 七、求微分方程的解 dsolve(‘Dy=1+y^2’) %求通解 ans = tan(t+C1) dsolve(‘Dy=1+y^2’,‘y(0)=1’)%求特解 ans = tan(t+1/4*pi) 求二阶微分方程的解 y=dsolve(D2y-2*Dy-3*y=0) y = C1*exp(3*t)+C2*exp(-t) y=dsolve(D2y-2*Dy-3*y=0,y(0)=0,y(1)=1) y = 1/(exp(3)-exp(-1))*exp(3*t)-1/(exp(3)-exp(-1))*exp(-t) y=simple(y) %化简表达式 y = (exp(3*t)-exp(-t))/(exp(3)-exp(-1)) ezplot(y,[-6,2]) 画出二阶微分方程的解的图形 求微分方程组的解 [f,g]=dsolve(Df=3*f+4*g,Dg=-4*f+3*g) f = exp(3*t)*(cos(4*t)*C1+sin(4*t)*C2) g = -exp(3*t)*(sin(4*t)*C1-cos(4*t)*C2) 练习 1、求下列符号函数的零点 （1） （2） 建议先画出函数的图形，观察根的范围。 2、求微分方程的符号解 （1） （2） （3） 希望提出指导与建议 MATLAB第五讲 MATLAB的符号计算功能 一、符号表达式和符号矩阵的创建 1、 符号变量的赋值 2、符号常量 3、符号矩阵的创建 4、符号方程的创建 5、将数字矩阵转化为符号矩阵 6、符号表达式的升幂 7、符号表达式的合并 8、变量代换 二、符号函数的微积分 1、求导数 2、求积分 3、 求泰勒级数taylortool(f) 三、求符号方程的精确解和近似解 1、求解单个符号方程 2 、对代数方程组求解 四、符号矩阵的基本运算 1、加、减、乘 2、求逆运算 3、求符号矩阵的除法 4、求矩阵的特征值和特征向量 5、求符号矩阵的行列式 6、求符号矩阵的约当矩阵 7、求符号矩阵的奇异值 五、符号函数画图 六、求符号函数的零点 七、求微分方程的解 一、入门 Symbolic工具包 MATLAB有一个符号计算工具包叫作 Symbolic Math Toolbox 其中有60多个专用函数。包括微积分、线性代数、方程求解、多项式的简约与展开、特殊数学函数等。 sym的指令在C:\MATLAB6p1\toolbox\symbolic，如果在搜索路境中没有设定则要添加。 符号表达式和符号矩阵的创建 符号表达式是数字、函数、变量的MATLAB字符串，或字符串数组。 符号运算是指使用已知的规则和给定的符号恆等式求解符号方程。 生成符号表达式用引号或sym函数 例 M=[a,b;c,d] ??? Undefined function or variable a. M=sym(‘[a,b;c,d]’) M=sym([a,b;c,d]) M = [ a, b] [ c, d] 1、符号变量的赋值 f1=sin(x)^2’；f2=exp(-x^2/2)’；f3=1/(1+x^2)’； f1 = sin(x)^2 f2 = exp(-x^2/2) f3 = 1/(1+x^2) 符号常量 没有变量的符号表达式叫作符号常量 a2=3 a2 = 3 a2+1 ans = 52 a2是一个符号常量，它是用ASCII码来存储的，‘3’的ASCII码是51因此a2+1得到的是52而不是4 （1）sym命令 M1=sym([sin(x),cos(x);-cos(x),sin(x)]) M1 = [ sin(x), cos(x)] [ -cos(x), sin(x)] （2）直接输入法