matlab --微分和积分

1.微分：

>> a=[9,-5,3,7];
>> x=-2:0.1:5;
>> f=polyval(a,x);
>> plot(x,f,'LineWidth',2);

a=fx=9x^3
-5x^2+ 3*x+7;
polyval 求多项式的值
plot

>> p=[5 0 -2 0 1];
>> polyder(p)

polyder求导数

conv([20 -7 5 10],[4 12 -3])

向量卷积：
向量乘法；
2.积分

>> p=[5 0 -2 0 1];
>> polyint(p,3)

ans =

    1.0000         0   -0.6667         0    1.0000    3.0000

polyint(a,c) 求函数的积分
c为常数项

poyval(polyint(p,3),7)

积分且x=7的值
3.diff()

x=[ 1 2 5 2 1];
diff(x);
ans =
	1 3 -3 -1

>> x=[1 2];y=[5 7];
>> slope=diff(y)./diff(x)

slope =

     2

求【1，5】 【2，7】之间的导数值

>> x=0:0.1:2*pi;
>> t=0:0.1:2*pi;
>> y=sin(t);
>> re=diff(y)./diff(x)

sinx求导
4.二次和三次微分

x=-2:0.005:2;
y=x.^3;
m=diff(y)./diff(x);
m2=diff(m)./diff(x(1:end-1));
plot(x,y,x(1:end-1),m,x(1:end-2),m2);
xlabel('x','Fontsize',18);
ylabel('y','Fontsize',18);

二次微分
5 积分

h=0.05;
x=0:h:2;
midpoint=(x(1:end-1)+x(2:end))./2;
y=4*midpoint.^3;
s=sum(h*y)

积分：4x^3在0-2的积分 midpoint法

y=@(x) 1./(x.^3-2*x-5);
integral(y,0,2)

1/x^3-2*x-5 在0-2的积分
@用到function

• Double and Triple Integrals
• f=@(x,y) y.*sin(x)+x.*cos(y); intergral12(f,pi,2*pi,0,pi)
  二重积分

三重积分

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