數(shù)值逼近課程設(shè)計(jì)-展示頁
2025-06-14 22:50本頁面
【正文】 。n=439。原函數(shù) 39。m:39。b.39。r39。k39。 y=1./(1+25*x.^2)。 y6=lagrange(x3,y3,x)。 y3=1./(1+25*x3.^2)。 y5=lagrange(x2,y2,x)。 y2=1./(1+25*x2.^2)。 y4=lagrange(x1,y1,x)。 y1=1./(1+25*x1.^2)。 end y=y+y1*y0(i)/prod(x0(i)x0(ij))。y1=1。y=zeros(size(x))。數(shù) 值 逼 近 課 程 設(shè) 計(jì) 報(bào) 告 080710205 邱明敏 作業(yè) 一 多項(xiàng)式插值的 Runge 現(xiàn)象 對(duì)于 Runge函數(shù) f(x)= ，在 [1,1]上作等距節(jié)點(diǎn)插值，分別取 n=4,n=8,n=12,編出程序，畫出此插值的圖像。 程序代碼 (matlab實(shí)現(xiàn) )： function y=lagrange(x0,y0,x) ii=1:length(x0)。 for i=ii ij=find(ii~=i)。 for j=1:length(ij),y1=y1.*(xx0(ij(j)))。 end function runge(m1,m2,m3) x1=1+2*[0:m1]/m1。 x=1::1。 x2=1+2*[0:m2]/m2。 x=1::1。 x3=1+2*[0:m3]/m3。 x=1::1。 x=1::1。 plot(x,y,39。,x,y4,39。,x,y5,39。,x,y6,39。) legend(39。,39。,39。,39。) 三 、 運(yùn)行結(jié)果（運(yùn)行過程）： 輸入 runge(4,8,12)，可得插值圖像： 作業(yè) 二 Remez 算法 求函數(shù) f(x)= 在 [1,1]上的二次多項(xiàng)式逼近。 y=exp(x)。 } void remez(double a,double b
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1