数据拟合试题及答案.docVIP

... ... 1．用给定的多项式，如 y=x3-6x2+5x-3 ，产生一组数据 (xi,yi ，i=1,2, ? ,n再),在 yi 上添加随 机干扰 (可用 rand 产生 (0,1)均匀分布随机数 ,或用 rands 产生 N(0,1)分布随机数 )，然后用 xi 和添加了随机干扰的 yi 作的 3 次多项式拟合， 与原系数比较。 如果作 2 或 4 次多项式拟合， 结果如何？ 用 电压V=10 伏 的 电 池 给 电 容 器 充 电 ， 电 容 器 上 t时刻 的 电压为 v(t ) V (V V )exp( t / ) ，其中 V0 是电容器的初始电压， 是充电常数。试由下面 0 一组 t，V 数据确定 V ， 。 0 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 解：程序如下： x=1:0.5:10; y=x.^3-6*x.^2+5*x-3; y0=y+rand; f1=polyfit(x,y0,1)% 输出多项式系数 y1=polyval(f1,x);% 计算各x点的拟合值 plot(x,y,+,x,y1) grid on title( 一次拟合曲线); figure(2); f2=polyfit(x,y0,2)%2 次多项式拟合 y2=polyval(f2,x); plot(x,y,+,x,y2); grid on title( 二次拟合曲线 ); figure(3); f4=polyfit(x,y0,4)%4 次多项式拟合 y3=polyval(f4,x); plot(x,y,+,x,y3) grid on title( 四次拟合曲线 ); figure(4); f6=polyfit(x,y0,6)%6 次多项式拟合 y4=polyval(f6,x); plot(x,y,+,x,y4) grid on title( 六次拟合曲线 ); 运行结果如下：依次为各个拟合曲线的系数（按降幂排列） f1 =43.2000 -149.0663 f2 = 10.5000 -72.3000 89.8087 f4 =0.0000 1.0000 -6.0000 5.0000 -2.5913 f6 = 0.0000 -0.0000 0.0000 1.0000 -6.0000 5.0000 -2.4199 运行后，比较拟合后多项式和原式的系数， 发现四次多项式系数与原系数比较接 近，四次多项式的四次项系数很小。 作图后， 发现一次和二次多项式的图形与原 函数的差别比较大， 属于欠拟合的情况， 而四次多项式和六次多项式符合得比较 好。作图如下： 解 ： 据 题 意 分 析 如 下 : 电 容 器 充 电 的 数 学 模 型 已 经 建 立 。 v(t ) V (V V )exp( t / ) (已知 V=10) 可见，v(t)与τ成指数变化关系，所以在 0 通过 曲线拟合的时候，使用指数曲线 y＝a x。( 非线性拟合 ) 。首先进行变量 1ea 2 代换在程序中用 v1 代替 v(t) ，t0 代替τ,v 2 是拟合后的曲线方程： 对 v(t ) 10 (10 v )exp( t / ) 变形后取对数，有ln(10 v(t)) ln(10 v0) ( t / ) 0 令y=ln(10- v(t ) ) ,f 1=ln(10- 0 ) ,f2= -1/t0, 则 v0=10-exp(f(2)) ， t0= -1/ f(1) 。 编写程序如下： t=[0.5 1 2 3 4 5 7 9]; v1=[6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63]; y=log(10- v1); f=polyfit(t,y,1) t0=-1/f(1) v0=10-exp(f(2)) v2=10-(10-v0)*exp(-t/t0); plot(t,v1,rx,t,v2,k:) grid on xlabel(时间 t(s)),ylabel(充电电压(V)); title(电容器充电电压与时间 t 的曲线); 程序运行输出结果如下： f =-0.2835 1.4766 t0 = 3.5269 v0 =5.6221 即电容器的初始电压为 v0 =5.6221 ，τ=3.5629。