基本证明就不说啦,基本上就是运用泰勒公式实现对函数的近似运算。
对于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+')%精确的二次导图像
做出图像