第4章MATLAB符号计算讲义.ppt

MATLAB基础及应用第2版 主编：于润伟 朱晓慧 第4章 MATLAB符号计算 本章要点 符号计算的基本函数 符号微积分 符号方程求解 4.1 符号函数的计算 一、符号变量和符号矩阵 函数：符号变量名=sym（‘表达式’） 说明：创建一个符号变量。表达式可以是变量、字符、字符串、数学表达式或字符表达式等。 函数：syms 变量名1 变量名2 变量名3 … 说明：一次创建多个符号变量。 例4-1 创建符号变量。 clear a=sym(‘matlab’) b=sym(‘3*x^2+4*x+7’) e=[1 3 5;2 4 6;7 9 11] %创建数值矩阵 m=sym(e) %创建符号矩阵（没加单引号） 例4-2 创建符号变量和符号矩阵。 clear syms A B C // 观察工作窗口 syms a b c d n=[a b c d;b c d a;c d a b;d a c b] 二、常用函数 数值运行结果为数值 符号运算结果为表达式或符号 MATLAB运算中，浮点运算速度最快，符号运算占用时间和内存比较多，但运算结果最精确。在默认情况下，当生成符号变量后，利用符号工具箱中提供的函数进行符号运算。 二、常用函数 1．算术函数 例4-3 计算表达式 与表达式x-1的和、差、积、商河乘方。 2．化简函数 例4-4 对表达式 进行因式分解。 例4-6 对表达式f=x(x(x-8)+6)t,分别将自变量x和t的同类项合并。 clear syms x t f=x*(x*(x-8)+6)*t; %与上例比又一种表示方法 collect(f) %按默认的变量x合并 collect(f,t) %按变量t合并 例4-9 对表达式f=x/y-y/x进行通分。 clear syms x y f=x/y-y/x; [m,n]=numden(f) %如果没有m、n则默认ans为分子 例4-10 对表达式 进行嵌套形式重写。 3．替换函数 函数：[R，SYM]=subexpr（S，‘SYM’） 说明：用变量SYM（字符或字符串）的值代替符号表达式S中重复出现的字符串，R是返回替换后的结果。 函数subs的调用格式： R=subs（S，old，new） 说明：该函数是用新的符号变量new替换原来符号表达式S中的变量old，R是替换后的符号表达式。需要注意的，当变量new是数值形式时，显示的结果虽然是数值，但事实上还是符号变量，可以用vpa函数强制求值。 例4-11 求解并化简三次方程 符号解。 例4-12 求表达式 在x=1时的代数值。 三、可视化符号函数计算器 在MATLAB的命令窗口输入funtool，即可启动可视化符号函数计算器 。 4.2 符号微积分 一、符号极限 例4-13 分别计算表达式 、 、 、 、 。 二、符号求导 函数：diff（s，x，n） 说明：其中s为符号表达式，x为自变量，n为求导的阶数。 例4-14 分别计算表达式 的一阶导数和三阶导数。 clear syms x diff(x^5) diff(x^5,3) %x，n都可以默认，自变量为系统默认，如果有多个变量，选择字母顺序表中最接近x的字母作为自变量。默认n为一阶导数。 三、符号积分  例4-15 分别计算下列表达式的积分。(1) (2) (3) (4) clear syms x y s=(4-3*x^2)^2; int(s) int(x/(x+y),x) int(x/(x+y),y) int(x^2/(x+2),x,1,3) double(ans) %以双精度计算ans 四、积分变换 例4-16 求函数 的傅里叶变换及其逆变换。 4.3 符号方程求解 一、代数方程 例4-19 求解代数方程 例4-20 求解代数方程组 程序： clear syms x y z s1=2*x^2+y^2-3*z-4 s2=y+z-3 s3=x-2*y-3*z [x,y,z]=solve(s1,s2,s3) % 注意[x,y,z]的使用