MATLAB知识：极限、导数、微分方程的数值解

MATLAB知识：极限、导数、微分方程的数值解
极限
MATLAB提供了求极限函数limit()，函数调用格式为：y = limit(fun,x,x0)。
其中，y为返回的函数极限值；
fun为要求解的函数；
x为函数自变量；
x0位函数自变量的取值，x趋近于x0。
MATLAB求解极限问题有专门的函数，编程如下：
clc,clear,close all
syms x a
I1=limit(‘(sin(x)-sin(3x))/sin(x)‘,x,0)
I2=limit(’(tan(x)-tan(a))/(x-a)‘,x,a)
I3=limit(’(3x-5)/(x^3*sin(1/x2))’,x,inf)
运行程序输出结果如下：
I1 =
-2

I2 =
tan(a)^2 + 1

I3 =
3
导数
diff是求微分最常用的函数，其输入参数既可以是函数表达式，也可以是符号矩阵。
MATLAB常用的格式是：diff(f, x, n)；
其中，f关于x求n阶导数。
MATLAB求解导数问题，编程如下：
clc,clear,close all
syms x y
f=sym('exp(-2x)cos(3x^(1/2))')
diff(f,x)
g=sym(‘g(x,y)’) %建立抽象函数。
f=sym(‘f(x,y,g(x,y))’) %建立复合抽象函数。
diff(f,x)