數(shù)值計(jì)算課程設(shè)計(jì)報告(插值法)-資料下載頁
2025-08-05 06:42本頁面
【正文】 function f=Hermite(x,y,y_1,x0) syms t。 f=。 if(length(x)==length(y)) if(length(y)==length(y_1)) n=length(x)。 else disp(39。y和y的導(dǎo)數(shù)的維數(shù)不相等！39。)。 return。 end end else disp(39。x和y的維數(shù)不相等！)。 return。 end for i=1:n h=。 a=。 for j=1:n if(j~=i) h=h*(tx(j))^2/((x(i)x(j))^2)。 a=a+1/(x(i)x(j))。 end end f=f+h*((x(i)t)*(2*a*y(i)y(i)y_1(i))+y(i))。 if(i==n) if(nargin==4) f=subs(f,39。t39。,x0)。 else f=vpa(f,6)。 end end End3 拉格朗日插值的MATLAB實(shí)現(xiàn) 在MATLAB編輯窗口中輸入 x=[2,0,4,5]。 y=[5,1,3,1]。 f=Language(x,y,1) f= 運(yùn)行程序時調(diào)用的Language函數(shù)如下： function f=Language(x,y,x0) syms t。 if(length(x)==length(y)) n=length(x)。 else disp(39。x和y的維數(shù)不相等！39。)。 return。 end f=。 for(i=1:n) l=y(i)。 for(j=1:i1) if(j~=i) l=l*(tx(j))/((x(i)x(j))。 end。 for(j=i+1:n) l=l*(tx(j))/((x(i)x(j))。 end f=f+1。 simplify(f)。 if(i==n) if(nargin==3) f=subs(f,39。t39。,x0)。 else f=vpa(f,6)。 end end End4 牛頓插值法的MATLAB實(shí)現(xiàn) 在MATLAB編輯窗口中輸入 x=[1 4]。 y=[1 16]。 f=Newton(x,y,) f=4運(yùn)行程序時調(diào)用的Newton函數(shù)如下： function f=Newton(x,y,x0) syms t。 f=。 if(length(x)==length(y)) n=length(x)。 c(1:n)=。 else disp(39。x和y的維數(shù)不相等！39。)。 return。 end f=y(1)。 y1=0。 l=1。 for(i=1:n1) for(j=i+1:n) y1(j)=(y(j)y(i))/(x(j)x(i))。 end c(i)=y1(i+1)。 l=l*(tx(i))。 f=f+cc(i)*l。 simplify(f)。 y=y1。 if(i==n1) if(nargin==3) f=subs(f,39。t39。,x0)。 else f=collect(f)。 f=vpa(f,6)。 end end end指導(dǎo)教師評語：指導(dǎo)教師(簽字)： 年 月 日各組員課程設(shè)計(jì)成績姓名學(xué)號成績高育坤1309064043王冬妮1309064044韓建1309064046李婧1309064047
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1