Matlab实现数值微分

基本证明就不说啦，基本上就是运用泰勒公式实现对函数的近似运算。

对于f(x)=-u‘’(x),u=sinx这个例子，代码如下

a(1)=2;
a(2)=-1;
for i=3:99
    a(i)=0;
end
for i=1:97
    for j=1:99
        C(i,j)=0;
    end
end
for i=1:97
     C(i,i)=-1;
     C(i,i+1)=2;
     C(i,i+2)=-1;
end
for i=1:97
    c(i)=0;
end
c(98)=-1;
c(99)=2;
A=[a;C;c];
x=0:0.01:1;
u=sin(x);
f=-diff(u,2);
b(1)=sin(0)+1/100^2*f(1);
b(99)=sin(1)+1/100^2*f(99);
for i=2:98
    b(i)=1/100^2*f(i);
end
U=inv(A)*b';
plot(U,'r*')%近似的二次导图像
hold on
plot(sin(x),'b+')%精确的二次导图像

做出图像