数值逼近课设.docxVIP

南京航空航天大学数值逼近课程设计姓名Ssh学 号保密学 院理学院2012/12/15注：本程序需在matlab2012版本中运行，7.0可能不支持某些语法。1.龙格现象：在[-1,1]区间上构造函数F(x)=1(1+25*x^2)的插值多项式，插值节点是4,8,12个，可以发现随着插值节点的增加，插值函数与原函数相比，在区间边上时震荡较大。Matlab代码如下：syms xf(x)=1/(1+25*x^2)ezplot(f(x),[-1,1])hold onfor k=4:4:12 %插值节点个数从4变到12r(x)=0*x^0p(x)=0*x^0w(x)=1*x^0h=2/ka=[]y=[]fori=1:k+1y(i)=-1+h*(i-1)endfori=1:k+1a(i)=f(y(i))endfori=1:k+1w(x)=(x-y(i))*w(x)endl(x)=diff(w(x))fori=1:k+1r(x)=w(x)/((x-y(i))*l(y(i))) %拉格朗日迭代公式p(x)=p(x)+a(i)*r(x)endezplot(p(x),[-1,1])hold on;end2.remez算法对于上给定的函数的最佳一致逼近多项式的寻求，由于问题的非线性性质，一般求解都比较困难。里米兹于1957年，采用逐次逼近的思想，提出了一个求最佳一致逼近多项式的近似算法，取得了良好的效果。matlab代码如下：syms xf(x)=exp(x);y=[-1,-0.5,0.5,1]; %可改为自己想要的交错点组k=[];ezplot(f(x),-1,1)hold onm=1;a=[1,y(1),y(1)*y(1),1;1,y(2),y(2)*y(2),-1;1,y(3),y(3)*y(3),1;1,y(4),y(4)*y(4),-1];b=[f(y(1));f(y(2));f(y(3));f(y(4))];c=inv(a)*b;n=abs(c(4));p(x)=c(1)+c(2)*x+c(3)*x^2;ezplot(p(x),-1,1)hold ong(x)=f(x)-p(x);h(x)=diff(g(x));ll(x)=diff(h(x));k=[];k(1)=y(1);k(2)=k(1)-h(k(1))/ll(k(1));i=1;while abs(k(i+1)-k(i))10^(-20) %牛顿迭代法求f（x）-p（x）倒数的零点，也就是它的无穷范点。i=i+1;k(i+1)=k(i)-h(k(i))/ll(k(i));endd=abs(g(k(i+1))); %d是g（x）的无穷范数d-n;whileabs(d-n)10^(-10) %或者是for m=1:4,控制迭代次数，此时下面不要“m=m+1“m=m+1;j=1 ; %f(x)-p(x)的极值点来替代交错点组中点while y(j+1)k(i+1)j=j+1;endif g(y(j))*g(k(i+1))0y(j)=k(i+1);elsey(j+1)=k(i+1);enda=[1,y(1),y(1)*y(1),1;1,y(2),y(2)*y(2),-1;1,y(3),y(3)*y(3),1;1,y(4),y(4)*y(4),-1]; %求差值函数时的矩阵b=[f(y(1));f(y(2));f(y(3));f(y(4))];c=inv(a)*b; %求系数和无穷范n=abs(c(4)); %最佳逼近p(x)=c(1)+c(2)*x+c(3)*x^2; %最佳逼近多项式ezplot(p(x),-1,1)hold ong(x)=f(x)-p(x);h(x)=diff(g(x));ll(x)=diff(h(x));k=[];k(1)=y(1);k(2)=k(1)-h(k(1))/ll(k(1));i=1;while abs(k(i+1)-k(i))10^(-20)i=i+1;k(i+1)=k(i)-h(k(i))/ll(k(i));endd=abs(g(k(i+1)));endp(x) n=c(1) %最佳逼近 a=c(2) %ab=c(2) %bc=c(2) %cm %满足条件的迭代次数ans =((5660103570157385*exp(-1))/9007199254740992 + (5899625318394855*exp(1))/9007199254740992 - (707512946269673*exp(1/2))/1125899906842624 - (7374531