MATLAB教学课件第3章MATLAB符号计算3.ppt

Dsolve 参数说明 (1)Eqn 微分方程输入方法： 例： 命令：Y=dsolve(‘D2y-2*Dy-3*y=0’,’x’) (2)condition 初始条件 y(a)=b;Dy(a)=c;D2y(a)=d; 最多可接受12个输入参量 例3.25 求微分方程的通解。 命令如下： y=dsolve(Dy-(x^2+y^2)/x^2/2,x) %解(1)。方程的右端为0时可以不写 y=dsolve(Dy*x^2+2*x*y-exp(x),x) %解(2) y=dsolve(Dy-x/y/sqrt(1-x^2),x) %解(3) 例3.26 求微分方程的特解。 命令如下： y=dsolve(Dy=2*x*y^2,y(0)=1,x) %解(1) y=dsolve(Dy-x^2/(1+y^2),y(2)=1,x) %解(2) 3.6 积分变换 1. 傅立叶(Fourier)变换 在MATLAB中，进行傅立叶变换的函数是： fourier(fx,x,t) 求函数f(x)的傅立叶像函数F(t)。 ifourier(Fw,t,x) 求傅立叶像函数F(t)的原函数f(x)。 2. 拉普拉斯(Laplace)变换 在MATLAB中，进行拉普拉斯变换的函数是： laplace(fx,x,t) 求函数f(x)的拉普拉斯像函数F(t)。 ilaplace(Fw,t,x) 求拉普拉斯像函数F(t)的原函数f(x)。 3. Z变换 对数列f(n)进行z变换的MATLAB函数是： ztrans(fn,n,z) 求fn的Z变换像函数F(z) iztrans(Fz,z,n) 求Fz的z变换原函数f(n) 例3.29求函数的傅立叶变换及其逆变换。 命令如下： syms x t; y=abs(x); Fw=fourier(y,x,w) %求y的傅立叶变换 fx=ifourier(Fw,w,x) %求Ft的傅立叶逆变换 命令如下： f = exp(-x^2) Ft=fourier(f) fx=ifourier(Ft) 例3.30计算y=x2的拉普拉斯变换及其逆变换. 命令如下： x=sym(x);y=x^2; Ft=laplace(y,x,s) %对函数y进行拉普拉斯变换 fx=ilaplace(Ft,s,x) %对函数Ft进行拉普拉斯逆变换 例 求数列 fn=e-n的Z变换及其逆变换。 命令如下： syms n z fn=exp(-n); Fz=ztrans(fn,n,z) %求fn的Z变换 f=iztrans(Fz,z,n) %求Fz的逆Z变换 第3章 MATLAB符号计算 2010-9-1 本章内容 3.1 符号计算基础 3.2 微分运算 3.3 积分运算 3.4 泰勒级数 3.5 求解方程 3.6 积分变换 3.1 符号计算基础 3.1.1 符号对象 符号表达式是代表数字、函数、算子和变量的MATLAB字符串，或者字符串组； 符号运算是使用已知的规则和给定符号恒等式求解符号方程的实践； 符号变量和非符号变量区别：非符号变量参与计算之间必须赋值；而符号变量无须赋值。 1 建立符号变量和符号常数 (1) sym函数 S=sym(A) sym函数用来建立单个符号变量，例如，a=sym(‘a’)建立符号变量a，此后，用户可以在表达式中使用变量a进行各种运算。 定义符号常量 如 pi=sym(pi)；a=sym(3/4) 例3.1考察符号变量和数值变量的差别。 在 MATLAB命令窗口，输入命令： a=sym(a);b=sym(b);c=sym(c);d=sym(d); %定义4个符号变量 w=10;x=5;y=-8;z=11; %定义4个数值变量 A=[a,b;c,d] %建立符号矩阵A B=[w,x;y,z] %建立数值矩阵B det(A) %计算符号矩阵A的行列式 det(B) %计算数值矩阵B的行列式 例3.2 比较符号常数与数值变量在代数运算时的差别。 在 MATLAB命令窗口，输入命令： pi1=sym(‘pi’);k1=sym(‘8’);k2=sym(‘2’);k3=sym(‘3’); % 定义符号常数 pi2=pi;r1=8;r2=2;r3=3; % 定义数值变量 sin(pi1/3) % 计算符号表达式值 sin(pi2/