[2019/5/31]学习MATLAB第四节，求函数的导数和定积分的几个方法

一.问题
1.绘制函数$f(x)=(20\ast x^3-7\ast x^2+5\ast x+10)\ast (4x^2+12\ast x-3)$;及$f^{\prime }(x)$的图形

a = [20,-7,5,10];
b = [4,12,-3];
x = -2:0.01:1;
c = conv(a,b);  %提取两个函数的乘积后的系数值
f = polyval(c,x);%求出函数f(x)
hold on;
plot(x,f,'b--');
xlabel('x');
ylabel('f(x)');
f1 = polyval(polyder(c),x);%polyder(c)表示提取函数f(x)求导后的系数
plot(x,f1,'r');
hold off;
legend('f(x)','f^\prime(x)');%显示f(x)和f‘(x)

结果：
2.求f(x)的积分,利用polyint(p,n)函数
如$f(x)=5\ast x^4-2\ast x^2+1$，求$\int f(x)d(x)$

p=[5 0 -2 0 1];
a = polyint(p, 3);%polyint(p,n)表示提取函数积分后的系数值，n是常数项
polyval(a,7);%计算 ∫f(7)d(x)

3.利用diff()函数(计算向量内两元素的差值)，使用导数定义求$f^{\prime }(x)$，
定义：$f^{\prime }(x)={\lim }_{h\to 0}\frac{f(x0+h)-f(x0)}{h}$
例：$f(x)=sin(x),求f^{\prime }(x0),其中x0=\pi /2，h=0.1$

x0 = pi/2; h = 0.1;
x = [x0,x0+h];
y = [sin(x0),sin(x0+h)]; 
m = diff(y)./diff(x)

结果：

m =

   -0.0500

注：在线LaTex公式编辑器在CSDN博客上插入数学公式，http://private.codecogs.com/latex/eqneditor.php

二.问题
已知$f(x)=e^{-x}\ast sin(x^2/2),求当h=0.1、0.01、0.001时，f^{\prime }(x)的图形$

for i = 1:4
    h = power(10,-i);
    x = 0:h:2*pi;
    y = exp(-x).*sin(x.^2/2);
    m = diff(y)./diff(x);
    hold on;
    %plot(x,y);
    plot(x(1:end-1),m);%x(1:end-1)表示x(1)、x（2）...导函数比原函数少一个元素
    hold off;
end;

结果：
怎么二次求导？

m = diff(y)./diff(x);	%求f'(x)
m2 = diff(m)./diff(x(1:end-1));	%求f''(x)

三.求定积分的几个方法(求得是近似值)
①.中值规则
${\int }_{x0}^{x3}f(x)dx\approx hf0+hf1+hf2=h{\sum }_{i=0}^{n-1}fi$

其中$f0=f(\frac{x0+x1}{2})$、$f1=f(\frac{x1+x2}{2})$、$f2=f(\frac{x2+x3}{2})$

例：求$A={\int }_0^{2}4x^{3}dx=(2{)}^4-(0{)}^4=16$

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

结果：

s = 
	15.9950

注解：
x(1:end-1)表示 x0,x1,x2…
x(2:end)表示 x1,x2,x3…
(x(1:end-1)+x(2:end))./2表示$(\frac{x0+x1}{2})、(\frac{x1+x2}{2})、(\frac{x2+x3}{2})...$
midpoint 表示求$f(\frac{x0+x1}{2})$

②.求梯形面积规则

${\int }_{x0}^{x3}f(x)dx\approx h(\frac{f0+f1}{2}+\frac{f1+f2}{2}+\frac{f2+f3}{2})$

例：求$A={\int }_0^{2}4x^{3}dx=(2{)}^4-(0{)}^4=16$

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

结果：

s =

   16.0100

③.Simpson（辛普森） 法则

${\int }_{x0}^{x3}f(x)dx\approx \frac{h}{3}(f0+4f1+f2)$

例：求$A={\int }_0^{2}4x^{3}dx=(2{)}^4-(0{)}^4=16$

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

结果：

s =
	16

推广：

④integral()函数取定积分的准确值

求：${\int }_0^2\frac{1}{x^3-2x-5}dx$

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

结果：

ans =

   -0.4605

注：integral(y,a,b)中，y表示原函数，a、b分别表示积分上下限；
integral2(y,a,b,c,d);二次积分；a、b分别表示内层积分上下限，c、d表示外层积分上下限；