matlab 符号运算 数值,MATLAB中的微积分运算（数值&符号）

显然这个函数是单词differential(微分)的简写，用于计算微分。实际上准确来说计算的是差商。

如果输入一个长度为n的一维向量，则该函数将会返回长度为n-1的向量，向量的值是原向量相邻元素的差，于是可以计算一阶导数的有限差分近似。

(1)符号微分

1.常用的微分函数

函数:diff(f)     求表达式f对默认自变量的一次微分值

diff(f,x)  求表达式f对自变量x的一次积分值

diff(f,n)  求表达式f对默认自变量的n次微分值

diff(f,t,n)求表达式f对自变量t的n次微分值

>> x=1:10

x =

1 2 3 4 5 6 7 8 9 10

>> diff(x)

ans =

1 1 1 1 1 1 1 1 1

例1:求矩阵中各元素的导数

求矩阵[1/(1+a)  (b+x)/cos(x)

1/(x*y)   exp(x^2)]

对x的微分,可以输入以下命令

A = sym('[1/(1+a),(b+x)/cos(x);1,exp(x^2)]');

B = diff(A,'x')

可得到如下结果:

例2:求偏导数

求

的偏导数。

syms x y;

f = x*exp(y)/y^2;

fdx = diff(f,x)

fdy = diff(f,y)

可得到如下结果：

例3:求复合函数的导数

求

的导数

sym('x');

y = 'x*f(x^2)'

y1 = diff(y,'x')

得到结果如下：

例4:求参数方程的导数

对参数方程

求导

syms a b t

f1 = a*cos(t);

f2 = b*sin(t);

A = diff(f2)/diff(f1) %此处代入了参数方程的求导公式

B = diff(f1)*diff(f2,2)-diff(