MATLAB求导及线性方程组求解.docVIP

实验四 物理与电子工程系 电气工程及其自动化 2008126132—吴林 1. 的极限存在吗？画出其图像。 解：A.syms x limit(((x-2)/abs(x-2)),x,2) ans = NaN B.ezplot(((x-2)/abs(x-2)),[-10,15]) 2. 求方程dx/dt = -2x + 8 的解并在同一个图上绘制常量C1=1,2,...,5 时的图象。 解：dsolve(Dx=-2*x+8) ans = 4+exp(-2*t)*C1 ezplot(Dx=-2*x+8) 3．绘制下列函数的图像 （1）， 解：f = sym(sin(x) + x^2); ezplot(f,[0,2*pi]); （2）， 解：f = sym(x^3 + 2*x^2 + 1); ezplot(f,[-2 2]); 4．计算下列各式 （1），求 解：y = sym(sin(2*x)*log(x)); diff(y) ans = 2*cos(2*x)*log(x)+sin(2*x)/x （2），求，， 解： f = x*y*log(x+y); fx = diff(f,x) fx = y*log(x+y)+x*y/(x+y) fy = diff(f,y) fy = x*log(x+y)+x*y/(x+y) f2xy = diff(fx,y) f2xy = log(x+y)+y/(x+y)+x/(x+y)-x*y/(x+y)^2 （3）， 解：A.syms t y = log(1+t); int(y) ans = log(1+t)*(1+t)-1-t B.int(y,0,27) ans = 56*log(2)+28*log(7)-27 5．计算下列各式 （1） 解：symsum(sym((3/n)^n),1,inf) （2）在 0 附近的Taylor 展开 解：taylor(sym(sin(x))) 6．绘制函数在，上的表面图。 解：syms x y z = 1/(2*pi)*exp(-(x^2+y^2)); ezsurf(x,y,z,[-3,3,-3,3]); 7．求解线性方程组 解：[x,y] = solve(sym(2*x+3*y=1),sym(3*x+2*y=-1)) 8．对符号表达式，进行如下变换 （1）关于的拉普拉斯变换 （2）关于的傅立叶变换 （3）分别关于和的z变换 （1） 解：syms x y z = x*exp(-(x^2+y^2)); syms u v fourier(z,x,u) ans = -1/2*i*pi^(1/2)*u*exp(-y^2-1/4*u^2) （2） 解： laplace(z,y,v) ans = 1/2*x*exp(-x^2)*pi^(1/2)*exp(1/4*v^2)*erfc(1/2*v) （3） 解：A.ztrans(z,x,u) ans = -u*diff(ztrans(exp(-x^2-y^2),x,u),u) B.ztrans(z,y,v) ans = x*ztrans(exp(-x^2-y^2),y,v)