MATLAB学习系列——微积分

符号变量

1.定义符号变量
%定义单个或多个符号变量：
syms  x  y  z  t. 
%定义单个符号变量也可以用：
sym  x或：x='x'.
2.定义符号表达式函数
syms x;
f=sin(x)+2*x-8;
x=1;
eval(f)
3.符号方程定义
fs='x^2+3*x+2';
x=1;
eval(fs)

符号表达式运算

syms x;
fx=x^2+exp(5*x)+3;
gx=sin(x)-4*x^2;
z1=fx+gx       %加法
z2=fx*gx       %乘法
compose(fx,gx)   %嵌套函数f(g(x))

%反函数
%求y=sin（x-1)的反函数
syms x;
y=sin(x-1);
g=finverse(y)

符号与数转换

digits(6)     %定义有效位数为6位
sym(pi,'d')   %用十进制表示pi

求解代数方程

重点介绍第二种方法(solve())和第三种零点法

%1.利用求解多项式的方法
%求解5x^2-2x-4=0
a=[5 -2 -4];
x=roots(a)

%2.调用solve()
%求解ax^2+bx+c=0
syms a b c x;
F=a*x^2+b*x+c;
x=solve(F)

%3.利用求零点的方法
%求解2*x^3-6x-2=0的根
% 调用格式：x=fzero(f,x0)
% x=fzero(f,[a,b]) 
% x返回使函数值为零的点；
% f为符号表达式；x0为初值点；[a,b]是一个初始搜索区间，
% 应在区间端点处函数值异号，否则会有错误信息. 
g='exp(x)-4*x^2';
x=fzero(g,5)

函数极限

其中 inf是无穷的意思。
limit(f,x,a) ---------------$\underset{x\to a}{\lim }f(x)$
limit(f,a)------------------默认变量x活唯一符号变量
limit(f,x,a,‘right’)-------右极限
limit(f,x,a,‘light’)-------左极限

syms x;
f=sin(x)/x;
limit(f,inf)
limit(1/x,x,0,'left')

函数微分

求解f=ax^3+x^2-bx-c

syms x a b c;
f=a*x^3+x^2-b*x-c;
diff(f)       %求一阶微分
diff(f,a)     %关于变量a对于f的微分
diff(f,2)     %求2阶导数
diff(f,a,2)

Taylor级数展开

syms x;
y=sin(x);
yt=taylor(y,x,0,'order',6)
       %y对于x在x=0处展开6项

一元极值

1.函数fminbnd(‘f ’,a,b)是对f(x)在(a,b)上搜索极小值点.
2.若求f(x)的极大值点，需对－f(x)取极小值点. 所得即为f(x)的极大值点.
3.x返回极值点；fval返回极值；exitflag若返回1表示搜索成功. 若返回－1则未成功搜索到极小值点；若返回0表示已经达到函数评价或迭代的最大次数. output返回包含优化信息的结构输出

实例：

syms x;
f='x^3-x^2-x+1';
[x1,favl]=fminbnd( f,-2,2);

函数积分

syms x c;
f=sin(x)*x*c;
int(f)    %求不定积分
int(f,c)  %求对c的不定积分
int(f,0,2*pi) %求在[0,2*pi]上定积分

广义积分实例

syms x;
f=exp(-x^2);
a=int(f,-inf,inf)