第二节 MATLAB的数学运算.ppt

2.6 高等数学中的问题 1. 求极限 设 ，求当x→0, x→1+, x→1-,x →∞时函数的极限。 syms x f=(exp(3*x)-1)/x; limit(f,x,0) %求x→0的极限 ans = 3 limit(f,x,1,‘right’) %求x→1+的极限 ans = exp(3) - 1 limit(f,x,1,left) %求x→1-的极限 ans = exp(3) - 1 limit(f,x,inf) %求x→∞的极限 ans = Inf 2. 求导数 设 ，求 。 syms x n; y=x^(n)*sqrt(1-2*x); diff(y) %求 ans = n*x^(n - 1)*(1 - 2*x)^(1/2) - x^n/(1 - 2*x)^(1/2) %注意： diff(y) 与diff(y,x)和diff(y,x,1)得到的结果完全一样 diff(diff(y)) %求 ans = n*x^(n - 2)*(1 - 2*x)^(1/2)*(n - 1) - (2*n*x^(n - 1))/(1 - 2*x)^(1/2) - x^n/(1 - 2*x)^(3/2) 注意：diff(diff(y)) 与diff(y,x,2)得到的结果完全一样。 另外，如n也为变量，我们可以求得偏导数 如 diff(y,n,1) %求 ans = x^n*log(x)*(1 - 2*x)^(1/2) diff(y,n,2) %求 ans = x^n*log(x)^2*(1 - 2*x)^(1/2) diff(diff(y,x,1),n,1) %求 ans = x^(n - 1)*(1 - 2*x)^(1/2) - (x^n*log(x))/(1 - 2*x)^(1/2) + n*x^(n - 1)*log(x)*(1 - 2*x)^(1/2) 3. 求积分 syms x; int((sin(x)+cos(x))/sqrt(sin(x)-cos(x)),x) ans = -(2*cos(x) - 2*sin(x))/(sin(x) - cos(x))^(1/2) %求不定积分 %求定积分 syms x; f=(exp(x)-1)^4*exp(x); int(f,x,0,1) ans = (exp(1) - 1)^5/5 %求二重积分 syms x y; f=x^2*exp(-y^2) int(int(f,y,x,1),x,0,1) %先y后x ans = 1/6 - 1/(3*exp(1)) int(int(f,x,0,y),y,0,1) %先x后y ans = 1/6 - 1/(3*exp(1)) syms x y z; f=x*y*z; int(int(int(f,z,0,1-x-y),y,0,1-x),x,0,1) ans = 1/720 %求三重积分 ，其中Ω是由平面 x+y+z=1与三个坐标面围成的四面体。 详细分解可有： A=[10 2 5;1 8 6.3;4 5.6 55] A = 10.0000 2.0000 5.0000 1.0000 8.0000 6.3000 4.0000 5.6000 55.0000 rref(A) ans = 1 0 0 0 1 0 0 0 1 可以看出，此处矩阵A的秩为3。 2、矩阵的秩 MATLAB中利用rank()函数求得矩阵的秩。调用格式为： k=rank(A) 3、矩阵的范数 在MATLAB中利用norm()函数求矩阵或者矢量的范数。调用格式为 n=norm(A) %只计算矩阵的2范数，与下面p=2的效果相同； n=norm(A,p) %根据p的不同，用来计算矩阵的不同范数。 P=1为1范数，即矩阵A列矢量和的最大值； P=2为2范数，即为矩阵A的最大奇异值； P=-inf,为负无穷范数，即行矢量和的最小值； P=inf,为正