matlab学习（3）

%% 单变量函数求极限:声明符号变量、描述函数、直接求解
syms x;
f = sin(x)/x;
limit(f,x,0)

%% 复杂一点得极限求解
syms x a b;
f = x * (1 + a/x)^x*sin(b/x);
L = limit(f,x,inf)

%% 单边极限求解
syms x;
limit( (exp(x^3)-1) / (1-cos(sqrt(x - sin(x)))), x, 0, ‘right’)

%% 极限及拓展知识
x = -0.1:0.001:0.1;
y = (exp(x.^3)-1)./(1-cos(sqrt(x- sin(x))));
plot(x, y, ‘-’, [0], [12], ‘o’); % 在点（0，12）的地方标注原点

%%
syms x;
limit( (exp(x.^3)-1)./(1-cos(sqrt(x- sin(x)))), x, 0)

%% 正切函数的单边极限
syms t;
f = tan(t);
L1 = limit(f, t, pi/2, ‘left’)
L2 = limit(f, t, pi/2, ‘right’)

%% 序列的极限 factorial(n) n的阶乘
syms n;
f = n^(2/3)*sin(factorial(n))/(n+1);
F = limit(f,n, inf)

%%
syms x n;
f = natan(1/(n(x^{2+1)+x))*tan(pi/4+x/2/n)}n;
limit(f, n, inf)

%% 区间函数的极限运算
syms x n real;
f = x^n;
L1 = limit(f, x, inf)
L2 = limit(f, n, inf)

% 分段区间 piecewise
% L1 = piecewise(n == 0, 1, 0 < n, Inf, n < 0, 0)
% L2 = piecewise(x == 1, 1, 1 < x, Inf, x < -1, limit(x^n, n, Inf), x in Dom::Interval(-1, 0) | x in Dom::Interval(0, 1), 0)

%% 底层limit函数调用
% a, b > 0;
syms a b positive, syms x;
f = asin(8x^2)+bcos(2x-2);
% f 函数的极限
% symengine 符号运算引擎
% limit 使用 limit 函数来进行求解，是MuPAD下面调用的
% ‘Intervals’ 求取区间极限
% ‘x=infinity’ x趋向于无穷大
L = feval(symengine, ‘limit’, f, ‘x=infinity’, ‘Intervals’)

%% 分段函数绘图 调用 piecewise()
f = piecewise(‘abs(x) > 1.1’, …
‘1.1*sin(x)’, ‘abs(x)<=1.1’, ‘x’);
syms x;
x0 = -3:0.01:3;
f1 = subs(f, x, x0); % 将f函数中的x变量全部使用x0进行替换
plot(x0, f1)

%% 多变量函数的极限 累极限（是有次序的）
syms x a; syms y positive;
f = exp(-1/(y^2 + x^2))*sin(x)2/x^2 …
*(1+1/y²⁾(x+a^2*y2);
L1 = limit(limit(f,x,1/sqrt(y)),y, inf)
% L2 = limit(limit(f,x,x0),y, y0)

%% 函数的导数和高阶导数
% diff(fun, x); diff(fun, x, n) 自变量为唯一符号变量，可以省去x

syms x;
f = sin(x)/(x^2 + 4*x + 3);
f1 = diff(f);
f4 = diff(f, x, 4);

tic, diff(f, x, 50);toc

% 绘制原函数和其一阶函数
fplot(f, [0,5]);
hold on;
fplot(f1, [0, 5]);

%% 复合泛函求导 泛函:一个函数是另外一个函数的函数
syms t f(t);
G = simplify(diff(t^2*sin(t)*f, t, 3))

%% 矩阵函数的求导
syms x;
H = [ 4sin(5x), exp(-4x^2);
3x^2+4x+1, sqrt(4x^2+2) ];
H1 = diff(H, x, 3)

%% 参数方程的求导
syms t;

y = sin(t)/(t+1)^3;
x = cos(t)/(t+1)^3;
f = paradiff(y, x, t, 3);

[n, d] = numden(f); % 单独提取分子和分母

% 对分子分母进行分别化简
F = simplify(n)/simplify(d)