2MATLAB符号计算介绍.ppt

MATLAB 符号运算 Matlab 符号运算介绍 Matlab符号运算是通过符号数学工具箱（Symbolic Math Toolbox）来实现的。 符号表达式的建立 相关函数 findsym: 查找符号表达式中的符号变量 相关函数 subs：符号替换 例： syms x y a f=2*x+y; x=3,y=4; subs(f) subs(f,x,’a’) 符号矩阵 使用sym函数直接生成 符号矩阵的基本运算 符号矩阵的基本运算与数值矩阵的基本运算相类似。 六大常见符号运算 因式分解、展开、合并、简化及通分等 展开函数： expand 多项式展开 合并同类项： collect collect(f,v): 按指定变量 v 的次数合并系数； collect(f): 合并 f 中的默认自变量的各项系数。 简化函数： simple 和 simplify simple(f): 对 f 尝试多种不同的算法简化， 返回其中最短的简化形式； [R,HOW]=simple(f): R为f的最短简化形式， HOW中记录的为简化过程中使用的主要方法。 simplify(f): 简化函数 例：化简 分式通分： numden [N,D]=numden(f): N为通分后的分子，D为通分后的分母 六大常见符号运算 因式分解、展开、合并、简化及通分等 例：求极限 syms h n x L=limit((log(x+h)-log(x))/h,h,0) M=limit((1-x/n)^n,n,inf) 六大常见符号运算 因式分解、展开、合并、简化及通分等 例：设 y=sin(ax),求 syms a x y=sin(a*x) A=diff(y,x) B=diff(y,a) C=diff(y,x,2) D=diff(y,a,2) 六大常见符号运算 因式分解、展开、合并、简化及通分等 例：求积分 syms x f=(x^2+1)/(x^2-2*x+2)^2; I=int(f) g=cos(x)/(sin(x)+cos(x)); J=int(g,x,0,pi/2) h=exp(-x^2); K=int(h,x,0,inf) 六大常见符号运算 因式分解、展开、合并、简化及通分等 例：求级数 ，以及其前10项的部分和。 syms n S=symsum(1/n^2,n,1,inf) S10=symsum(1/n^2,n,1,10) 六大常见符号运算 因式分解、展开、合并、简化及通分等 1. 傅立叶(Fourier)变换 在MATLAB中，进行傅立叶变换的函数是： fourier(fx,x,t) 求函数f(x)的傅立叶像函数F(t)； ifourier(Fw,t,x) 求傅立叶像函数F(t)的原函数f(x)。 例1 求函数的傅立叶变换及其逆变换。 命令如下： syms x t; y=abs(x); Ft=fourier(y,x,t) %求y的傅立叶变换 fx=ifourier(Ft,t,x) %求Ft的傅立叶逆变换 2. 拉普拉斯(Laplace)变换 在MATLAB中，进行拉普拉斯变换的函数是： laplace(fx,x,t) 求函数f(x)的拉普拉斯像函数F(t)； ilaplace(Fw,t,x) 求拉普拉斯像函数F(t)的原函数f(x)。 例2 计算y=x2的拉普拉斯变换及其逆变换. 命令如下： x=sym(x);y=x^2; Ft=laplace(y,x,t) %对函数y进行拉普拉斯变换 fx=ilaplace(Ft,t,x) %对函数Ft进行拉普拉斯逆变换 * * MATLAB 符号运算 符号对象的建立：sym 和 syms 例： syms x y z x=sym(‘x’); y=sym(‘y’); z=sym(‘z’); == f2=sym(’sin(x)+cos(x)’) syms x f1=sin(x)+cos(x) f3=’sin(x)+cos(x)’ 用这种方法创建的符号表达式对空格很敏感，不要在字符间随意添加空格！ 推荐！ 符号变量与符号常量 a=sym(‘a’); b=sym(‘5’); MATLAB 符号运算 findsym(f)按字母顺序列出符号表达式 f 中的所有自由变量findsym(f,N)列出 f 中距离 x 最近的 N 个自由变量（i，j 除外） 默认自变量＝findsym(f,1) a=sym(a);x=sym(x);k=sym(3); f=k*x+a; findsym(f) 例： an