2013春ME 实验二数据建模实验报告我打野.docVIP

1.编写一个函数实现Lagrangian插值。function f=Lagrangian (x,y,q) m=length(x); n=length(y); s=0; for i=1:n t=ones(1,length(q)); for j=1:n if j~=i t=t.*(q-x(j))/(x(i)-x(j)); end end s=s+t*y(i); end 2.在下面的函数中选择两个，在n个结点上（n不要太大，如5～11）分别用Lagrangian、分段线性、三次样条插值三种插值方法，计算m个插值点的函数值 (m要适中，如50～100)。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n, 再做比较，由此做初步分析。 （1） （2） （3） （4） （1） hold on x=linspace(0,2*pi,5); y=sin(x); p=polyfit(x,y,4); xi=0:0.1:2*pi; yi=polyval(p,xi); plot(xi,yi,’b’) yi=interp(x,y,xi); plot(xi,yi,’—r’) yi=interp(x,y,xi,’spline’); plot(xi,yi,’o’) title(‘n=5三种曲线比较’); legend(‘多项式插值’,’分段线性插值’,’样条插值’); xlabel(‘x轴’);ylabel(‘y轴’) （2） 代码与上题基本相同 x=linspace(-1,1,5); y=sqrt(1-x.*x); (3) x=linspace(-2,2,5); y=(cos(x)).^10; (4) x=linspace(-2,2,10); y=exp(-x.^2); 3. 某天的气温变化如第3题数据表，试用最小二乘法找出这一天的气温变化规律。 第3题数据表（工作表：气温变化） 时刻 t(h) 0 1 2 3 4 5 6 7 8 9 10 11 12 温度 T(℃) 15 14 14 14 14 15 16 18 20 22 23 25 28 时刻 t(h) 13 14 15 16 17 18 19 20 21 22 23 24 　 温度 T(℃) 31 32 31 29 27 25 24 22 20 18 17 16 　 fun2=inline(c(1)*x.^2+c(2)*x+c(3),c,x); x=0:24; y=[15 14 14 14 14 15 16 18 20 22 23 25 28 31 32 31 29 27 25 24 22 20 18 17 16]; c=lsqcurvefit(fun2,[0 0 0],x,y); T=poly2sym(c) T = - (6745190612456979*x^2)/72057594037927936 + (5841826433026835*x)/2251799813685248 + 2368816419385319/281474976710656 4.用电压V=10v的电池给电容器充电，电容器上t时刻的电压为 ，其中是电容器的初始电压，是充电常数。试由第4题数据表中的一组t,V 数据确定 和。 第4题数据表（工作表：电容器） t(s) 0.5 1 2 3 4 5 7 9 V(v) 6.4 6.5 7.3 8.2 8.7 9 9.4 9.6 t=[0.5 1 2 3 4 5 7 9]; v=[6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63]; t1=linspace(0,9,100); p=polyfit(t,v,1); v1=polyval(p,t1); plot(t,v,r+,t1,v1,b-) V=0.4t+10.2 a=ln(10/(10-v0)); b=1/T=0.4 V0= 0.12 T=2.50 5.给定数据见第5题数据表： 第5题数据表（工作表：三次样条） x 0.25 0.30 0.39 0.45 0.53 y 0.5000 0.5477 0.6245 0.6708 0.7280 分别就下列边界条件求三次样条函数并作图。 （1） （2） (1)x=[0.25 0.30 0.39 0.45 0.53]; y=[0.5000 0.5477 0.6245 0.6708 0.7280]; s=spline(x,y); fnplt(s,r) plot(x,y,o,x,y,) s.coefs ans = 0.9224 -0.8943 0.9964 0