數(shù)值逼近課程設(shè)計(jì)-其他專業(yè)-全文預(yù)覽
【正文】 } u=p[m1]/g[m1]。 g[k1]=(xxx0)/(x[k1]x[mi2])。i=n1。j++){ k=m。 g[i]=u。 } while(1==1){ u=。 for(k=0。 double x[21],g[21],d,t,u,s,xx,x0,h,yy。 參考文獻(xiàn)：數(shù)值分析算法描述與習(xí)題解答（徐士良著） 二 、 程序代碼 (C語言 實(shí)現(xiàn) )： include include double remezf(double x) { double y。n=839。原函數(shù) 39。b.39。k39。 y6=lagrange(x3,y3,x)。 y5=lagrange(x2,y2,x)。 y4=lagrange(x1,y1,x)。 end y=y+y1*y0(i)/prod(x0(i)x0(ij))。y=zeros(size(x))。 程序代碼 (matlab實(shí)現(xiàn) )： function y=lagrange(x0,y0,x) ii=1:length(x0)。 for j=1:length(ij),y1=y1.*(xx0(ij(j)))。 x=1::1。 x=1::1。 x=1::1。 plot(x,y,39。,x,y5,39。) legend(39。,39。) 三 、 運(yùn)行結(jié)果（運(yùn)行過程）： 輸入 runge(4,8,12)，可得插值圖像： 作業(yè) 二 Remez 算法 求函數(shù) f(x)= 在 [1,1]上的二次多項(xiàng)式逼近。 } void remez(double a,double b,double p[],int n,double eps) { int i,j,k,m。 d=+35。 x[k]=(b+a+(ba)*t)/。i++){ p[i]=remezf(x[i])。j=n1。
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號(hào)-1