matlab ppt 第三讲.pptVIP

第三讲 MATLAB的符号计算 一 符号对象的定义 二 符号计算中的基本运算和函数 三 符号数值任意精度控制和运算 四 符号表达式中自由变量的确定 五 符号对象的操作与转换 六 符号函数绘图 一 符号对象的定义 sym syms 三 符号数值任意精度控制和运算 符号运算与数值运算的区别主要有以下几点： 数值型运算每一次运算都会有一定的截断误差，重复的多次数值运算就可能会造成很大的累积误差。 符号运算不会出现截断误差，符号运算可以得出完全的封闭解或任意精度的数值解。 符号运算的时间较长，而数值型运算速度快。 四 符号表达式中自由变量的确定 MATLAB将基于以下原则选择一个自由变量： 小写字母i和j不能作为自由变量。 符号表达式中如果有多个字符变量，则按照以下顺序选择自由变量：首先选择x作为自由变量；如果没有x，则选择在字母顺序中最接近x的字符变量；如果与x相同距离，则在x后面的优先。 大写字母比所有的小写字母都靠后。 五 符号对象的操作与转换 1 符号表达式的操作 2 置换操作 3 复合函数和反函数 4 微积分 5 符号代数方程/微分方程求解 例: 求 解： syms a b x; f=[a*x, b*x^2; 1/x, sin(x)]; int(f) 5) 符号积分变换 Fw=fourier(ft, t, w)---求时域函数ft的傅立叶变换Fw ft=ifourier(Fw, w, t)---求频域函数Fw的傅立叶反变换ft Fs=laplace(ft, t, s)---求时域函数ft的Laplace变换Fs ft=ilaplace(Fs, s, t)---求频域函数Fs的Laplace反变换ft FZ=ztrans(fn, n, z)---求时域系列fn的Z变换FZ fn=iztrans(FZ, z, n)---求频域系列FZ的Z变换fn 5 符号代数方程求解 这里所讲的一般代数方程包括线性，非线性和超越方程等，求解指令为solve. S=solve(‘eq1’, ‘eq2’,…, ‘eqn’, ‘v1’, ‘v2’, …,’vn’)---方程组指定变量的解(推荐格式) S=solve(exp1, exp2,…, expn, v1, v2,…,vn)---方程组指定变量的解(可用格式) 说明： ● ‘eq1’, ‘eq2’,…, ‘eqn’或是字符串表达的方程，或是字符串表达式. ‘v1’, ‘v2’, …,’vn’是字符串表达的求解变量名． ● exp1, exp2,…, expn只能是符号表达式． v1, v2,…,vn是求解的符号变量． ● 如eq1, eq2,…,eqn是不含等号的表达式，则指令是对eq1=0, eq2=0, …, eqn=0求解． ● S是一个结构数组，如要显示求解结果，必须采用S.V1, S.V2,…,S.Vn的援引 方式. 例 求(x+2)x=2的解 解： syms x; S=solve(‘(x+2)^x=2’, ’x’) 符号微分方程求解 S=dsolve(‘eqn1’, ‘eqn2’, …)---符号微分方程组求解指令格式． 说明： ● 输入宗量包括三部分内容：微分方程，初始条件，指定独立变量．其中微分方程是必不可少的．其余可以没有． ● 对于独立变量，若要指定独立变量，总是由全部输入宗量的最后一个指令，默认的独立变量是t ● 字母D代表微分算子，即d/dt,D2代表d^2/dt^2, dny代表对y(t)的n阶微分． ● 初始条件可以写成格式：y(a)=b, Dy(c)=d,在解中可以出现任意常数符C1,C2等． ● 输出S为结构对象． 例 求dx/dt=y, dy/dt=-x的解 S=dsolve(‘Dx=y, Dy=-x’); disp([blanks(12), ‘x’, blanks(21), ‘y’], disp([S.x, S.y]) x y [ cos(t)*C1+sin(t)*C2, -sin(t)*C1+cos(t)*C2] 六 符号函数画图 ezplot(f)---在默认区间(-2pi, 2pi)绘制f=f(x)的函数图． ezplot(f, [a, b])---在区间(a, b)绘制f=f(x)的函数图． fplot(fun, lims)---绘制fun函数在区间lims=[xmin,xmax]的图形 ezplot(x^3-10*x, [0,10]) 例 * * 符号表达式: 非数字的数学表达式。由符号常量,符号变量,符号函数,符号操作等构成