[工学]MATLAB符号计算.pptVIP

作业： 1.化简 2. syms a b x X Y k=sym(’3’); z=sym(c*sqrt(delta)+y*sin(theta)); f=a*z*X+(b*x^2+k)*Y; 试写出下面命令的结果。 findsym(f) findsym(f,1) findsym(f,2) findsym(f,3) 3.设A是一个符号矩阵（定义如下），试指出findsym(A,1) 的输出结果，并由此得出结论： findsym 确定自由变量时，是对整个矩阵进行的，还是对各个矩阵元素分别进行？ syms a b t u v x y; A=[a+b*x,sin(t)+u; x*exp(-t),log(y)+v]; 4.(1)设 求复合函数 f(g(z)); (2)设 求复合函数 f(g(z))。 5.替换函数subs的使用。 (1)试指出下面程序中的f1、f2、f3、f4、f5的值。 syms a x; f=a*sin(x)+5; f1=subs(f,sin(x),sym(y)) f2=subs(f,{a,x},{2,sym(pi/3)}) f3=subs(f,{a,x},{2,pi/3}) f4=subs(subs(f,a,2),x,0:pi/6:pi) f5=subs(f,{a,x},{0:6,0:pi/6:pi}) syms x y; f=2*x+y; subs(f,{x,y},{x+y,x-y}) (2)试指出下面替换的结果，并说明理由。 4.4 常见积分变换函数 常见的积分变换有傅立叶变换、拉普拉斯变换和Z变换。 1．傅立叶(Fourier)变换 在MATLAB中，进行傅立叶变换的函数是： fourier(f):默认得到关于w的函数 fourier(f,v):指定函数为v fourier(f,u,v)：f有多个变量 2．傅立叶(Fourier)反变换 f = ifourier(F) f = ifourier(F,u) f = ifourier(F,v,u) 例题： 2．拉普拉斯(Laplace)变换 在MATLAB中，进行拉普拉斯变换的函数是： laplace(fx,x,y)： 求函数f(x)的拉普拉斯像函数F(y)， 默认为F(s)。 ilaplace(Fw,s,y)：求拉普拉斯像函数F(s)的原函数f(y)， 默认为f(x) 例9-5 计算y=x3a的拉普拉斯变换及 逆变换。 3．Z变换 当函数f(x)呈现为一个离散的数列f(n)时，对数列f(n)进行z变换的MATLAB函数是： ztrans(fn,n,z)：求fn的Z变换像函数F(z)。 iztrans(Fz,z,n)：求Fz的z变换原函数f(n)。 例 求数列 的Z变换。 4.6 符号方程求解 4.6.1 符号代数方程求解 在MATLAB中，求解用符号表达式表示的代数方程可由函数solve实现，其调用格式为： solve(s)： 求解符号表达式s的代数方程， 其中求解变量为默认变量。 solve(s,v)： 求解符号表达式s的代数方程， 其中的求解变量为v。 solve(s1,s2,…,sn,v1,v2,…,vn)： 求解符号表达式s1,s2,…,sn组成的代数 方程组，求解变量分别v1,v2,…,vn。 syms a b c x f=sym(a*x*x+b*x+c=0) solve(f) ans = [ 1/2/a*(-b+(b^2-4*c*a)^(1/2))] [ 1/2/a*(-b-(b^2-4*c*a)^(1/2))] solve(1+x=sin(x)) ans = -1.9345632107520242675632614537689 4.6.2 符号常微分方程求解 在MATLAB中，用大写字母D表示导数。例如，Dy表示y，D2y表示y，Dy(0)=5表示y(0)=5。 D3y+D2y+Dy-x+5=0表示微分方程y+y+y-x+5=0。符号常微分方程求解可以通过函数dsolve来实现，其调用格式为： dsolve(e,c,v) 该函数求解常微分方程e在初值条件c下的特解。参数v描述方程中的自变量，省略时按缺省原则处理，若没有给出初值条