matlab中int积分不准_MATLAB在高等数学中的应用（极限、导数、积分）

(一)求极限

matlab中求极限的命令为

limit(expr, x, a)
limit(expr, a)
limit(expr)
limit(expr, x, a, 'left')
limit(expr, x, a, 'right')

其中，limit(expr, x, a)表示求符号表达式expr关于符号变量x趋近于a时的极限；limit(expr)表示求默认变量趋近于0时的极限。

关于left和right则表示左极限和右极限。

clc,clear;
syms x
b = limit(sqrt(1+x^2)-1)/(1-cos(x))

b =
0

clc,clear;
syms x a
b = limit((1+a/x)^x,x, inf)

b =
exp(a)

(二)求导数

matlab的求导数命令为

diff(expr)
diff(expr, v)
diff(expr, sym('v'))
diff(expr, n)
diff(expr, v, n)
diff(expr, n, v)

其中：diff(expr)表示求表达式expr，关于默认变量的1阶导数；diff(expr, v, n)和diff(expr, n, v)都表示求表达式expr关于符号变量v的n阶导数。

clc, clear;
syms y(x)
y = log((x+2)/(1-x));
d3y = diff(y,3)
d3y = simplify(d3y)%对符号函数进行化简
pretty(d3y)%分数线居中显示

d3y =
(2*(1/(x - 1) - (x + 2)/(x - 1)^2)*(x - 1))/(x + 2)^3 - (2*(2/(x - 1)^2 - (2*(x + 2))/(x - 1)^3))/(x + 2) - (2*(1/(x - 1) - (x + 2)/(x - 1)^2))/(x + 2)^2 + (2*(2/(x - 1)^2 - (2*(x + 2))/(x - 1)^3)*(x - 1))/(x + 2)^2 + ((6/(x - 1)^3 - (6*(x + 2))/(x - 1)^4)*(x - 1))/(x + 2)
d3y =
-(18*(x^2 + x + 1))/(x^2 + x - 2)^3
2
(x  + x + 1) 18
- ---------------
2         3
(x  + x - 2)

(三)求极值

求函数f(x) = x^3+6*x^2+8*x-1的极值点，并画出函数的图形。

解：对这类问题一般的求解方式为对f(x)求导，然后令导函数等于零，解方程则可求得函数f(x)的极值点。

clc,clear;
syms x
y = x^3+6*x^2+8*x-1;
dy = diff(y);
dy_zero = solve(dy)
dy_zero_num = double(dy_zero)%变成数值类型
ezplot(y)%符号函数画图

dy_zero =
- (2*3^(1/2))/3 - 2
(2*3^(1/2))/3 - 2
dy_zero_num =
-3.154700538379252
-0.845299461620748

(四)求积分

1.求不定积分

Matlab求符号函数不定积分的命令为

int(expr)

int(expr, v)

例：求不定积分 

clc,clear;

syms x

I = int(1/(1+sqrt(1-x^2)))

pretty(I)

结果为：

I =

(x*asin(x) + (1 - x^2)^(1/2) - 1)/x

2

x asin(x) + sqrt(1 - x ) - 1

----------------------------

x

2.求定积分

1)求定积分的符号解

Matlab求fuha符号函数的定积分命令为

int(expr, a, b)

int(expr, v, a, b)

例：求定积分

clc,clear;

syms x

I = int(cos(x)*cos(2*x), -pi/2, pi/2)

结果为：

I =

2/3