Lecture5_符号运算基础.ppt

利用diff也可以对有参数方程确定的函数、隐含数求导，求多元函数的偏导数。 [例] 设 ，求u的一阶偏导数。 方法1： syms x y z;f=sqrt(x^2+y^2+z^2); du_dx=diff(f,x), du_dy=diff(f,y), du_dz=diff(f,z), 方法2： syms x y z;f=sqrt(x^2+y^2+z^2); jacobian(f,[x,y,z]) [例]设 ，求 syms x y; z=x^6-3*y^4+2*x^2*y^2; D2z_d2x=diff(z,x,2) D2z_d2y=diff(z,y,2) D2z_dxdy=diff(diff(z,y),x) [例] 设 ，求 syms a t; dx_dt=diff(a*(t-sin(t)),t); dy_dt=diff(a*(1-cos(t)),t); dy_dx=dy_dt/dx_dt [例]设 ，求 分析： 对 两端微分可得： 程序： syms x y; F=log(x)+exp(-y/x)-exp(1); dy_dx=-diff(F,x)/diff(F,y) [练习] 设 求 3.4.3 符号级数求和－symsum symsum(S,v,a,b) 求符号表达式S当v从a到b时的级数和。当v省略时，由findsym(S,1)确定。 [例] 求 syms n x; symsum([1/2^n, 1/n^2, x^n/(n*2^n)],n,1,inf) [练习]求 3.4.4 符号积分－int 函数int既可求定积分，也可求不定积分。格式： int(S,v,a,b) 求符号表达式S关于变量v 在区间[a,b]上的定积分 int(S,v) 求符号表达式S关于变量v 的不定积分 当v省略时，自动由findsym(S,1) 确定。 [例]设 ,求 syms x; f=cos(x-pi/6)*sin(x+pi/6); S=int(f,x,0,2*pi) double(S) [例]求 syms t x; F=int(1/log(t),t,0,x) [例]计算 syms x; int(abs(x-1),0,2) [例]求 syms a b x; f=[a*x,b*x^2;1/x,sin(x)]; int(f) [例]求二次积分 syms x y; F= int(int(x*y,y,2*x,x^2+1),x,0,1) 3.5 积分变换 3.5.1 符号卷积 [例] 已知系统的刺激响应 ，求输入为 时的输出。 说明：由于 ，故输出为 syms T t tao; ut=exp(-t);ht=exp(-t/T)/T; uh_tao=subs(ut,t,tao)*subs(ht,t,t-tao); %被积函数 yt=int(uh_tao,tao,0,t); yt=simple(yt) 3.5.2 傅里叶变换及其反变换—fourier和i fourier fourier的调用格式： F=fourier(f,t,w) 返回函数f(t)的傅里叶变换F(w) ifourier的调用格式类似： f=ifourier(F,w,t) 返回函数F(w)的逆傅里叶变换f(t) 当参数t缺省时，f的默认自变量是x;当 w缺省时，F默认的自变量是w. 3.5.3 拉普拉斯变换及其反变换—laplace和i laplace 格式： F= laplace(f,t,s) 返回函数f(t)的拉普拉斯变换F(s) f=ilaplace(F,s,t) 返回函数F(s)的逆拉普拉斯变换f(t) 当参数t缺省时，f的默认自变量是t;当 s缺省时，F默认的自变量是s. 3.5.4 Z 变换及其反变换—ztrans和iztrans 格式： F= ztrans(f,n,z) 返回序列f(n)的Z变换F(z) f=iztrans(F,z,n) 返回函数F(z)的逆Z变换f(n) 当参数n缺省由时，f的默认自变量是n;当 z缺省时，F默认的自变量是z.