【正文】
參考文獻(xiàn)[1] 易大義,沈云寶,李有法. 計(jì)算方法(第2版),浙江大學(xué)出版社. .[2] 張琨 高思超 畢靖 編著 MATLAB2010從入門(mén)到精通 電子工業(yè)出版社。 break。Error!x(%d) is out of range!\n39。)。if n~=p fprintf(39。p=length(y0)。程序優(yōu)化由分段線(xiàn)性插值和分段二次插值的原理,x取值在函數(shù)表范圍內(nèi)時(shí),插值結(jié)果有意義,而當(dāng)x取值在函數(shù)表范圍以外,利用分段線(xiàn)性插值公式仍可以進(jìn)行運(yùn)算并得到一個(gè)值,但其結(jié)果不準(zhǔn)確;分段二次插值則無(wú)法找到三個(gè)合適的點(diǎn)以求插值,不予以輸出結(jié)果;若輸入的函數(shù)表x與y的長(zhǎng)度不相等,則無(wú)法插值。 x0=[ ]。 y2=lagrange(x,y,x0)y(1)=Error!x(2) is out of range!y2 = 首先輸入函數(shù)變及待求點(diǎn) x=[ ]。 y=[2 3 4 5]。 x0=[ ]。 y2=lagrange(x,y,x0)Error! Please input again! x=[1 2 3 4]。%輸出插值結(jié)果endendend5 算例分析 測(cè)試示例 x=[1 2 3 4]。y(%d)=%.5f\n39。 end y(i)=s。 for j=1:n if j~=k p=p*(zx0(j))/(x0(k)x0(j))。 break。Error!x(%d) is out of range!\n39。 s=。)。%n,p,m分別為x0,y0,x長(zhǎng)度if n~=p fprintf(39。p=length(y0)。,i,y(i),x0(l1),y0(l1),x0(l),y0(l),x0(l+1),y0(l+1))。fprintf(39。 end end sum=sum+y0(k)*a。 end%輸入正確時(shí),若l還等于零,l=n1 for k=l1:l+1 a=。