matlab计算微积分

多项式微积分

多项式求值

polyval()

a=[9,-5,3,7];//系数
x=-2:0.01:5;
f=polyval(a,x);
plot(x,f,'LineWidth',2);

求导

polyder()

a=[9,-5,3,7];//系数
a1=polyder(a);//b是导函数的系数
polyval(a1,x0)//返回x0点的导数值

不定积分

polyint()

a=[9,-5,3,7];
a1=polyint(a,3);//返回原函数的系数

数值微积分

求导

diff(x)返回x向量相邻项作差形成的向量

在某一点的导数

x0=pi/2;
detaX=0.01;
x=[x0,x0+detaX];
y=[sin(x0),sin(x0+detaX)];
y1=diff(y)./diff(x);//导数值

导函数

x=0:0.01:2*pi;
y=sin(x);
y1=diff(y)./diff(x);
plot(x,y1);//画出导函数

高阶导

y2=diff(y1)./diff(x(1:end-1));//y1向量的长度比x少1
y3=...
......

定积分

integral()

f=@(x) x.^2+2*x+1;
integral(f,0,2);//从0到2对f作定积分

多元函数定积分

integral2()
integral3()

f=@(x,y) y.*sin(x)+x.*cos(y);
integral2(f,pi,2*pi,0,pi);//x:pi~2pi,y:0~pi

f=@(x,y,z) ...;
integral3(f,0,pi,0,1,-1,1);

符号微积分

符号微分

diff(f)
diff(f,x)

符号积分

int(f)
int(f,x)
int(f,x,a,b)