拟合作业拟合作业.doc

新乡学院数学与信息科学系 实验报告 实验项目名称 拟合实验 所属课程名称 数学实验 实 验 类 型 综合性实验 实 验 日 期 2013-6-7 班 级 数学与应用数学二班 学 号 11111012037 姓 名 宋甜甜 成 绩 一、实验概述： 【实验目的】 熟练的掌握用Matlab解决拟合问题，用图像来清晰地比较出拟合的好处； 了解插值与拟合的关系； 掌握用拟合来解决最优解问题。 【实验原理】 1．已知一组（二维）数据，即平面上 n个点（xi,yi) i=1,…n, 寻求一个函数（曲线）y=f(x), 使 f(x) 在某种准则下与所有数据点最为接近，即曲线拟合得最好。 曲线拟合问题最常用的解法——线性最小二乘法的基本思路 作多项式f(x)=a1xm+ …+amx+am+1拟合,可利用已有程序: a=polyfit(x,y,m) 其中（输出拟合多项式系数a=[a1, …,am , am+1] (数组)）；输入同长度的数组x，y；拟合多项式次数） 多项式在x处的值y可用以下命令计算： y=polyval（a，x） 【实验环境】 Matlab7.1 Window XP 二、实验内容： 【实验方案】 对多项式进行拟合运算，比较次数不同的拟合结果的差别 利用拟合解决实际问题。 【实验过程】（实验步骤、记录、数据、分析） 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次多项式拟合，结果如何？ x=1:0.05:10; y=x.^3-6*x.^2+5*x-3; y0=y+rand; f0=polyfit(x,y0,1) y1=polyval(f0,x); plot(x,y,+,x,y1) title(一次拟合曲线); f0 = 42.0345 -149.6663 f1=polyfit(x,y0,3) y1=polyval(f1,x); plot(x,y,+,x,y1) title(三次拟合曲线); f1 = 1.0000 -6.0000 5.0000 -2.5140 f2=polyfit(x,y0,2) y2=polyval(f2,x); plot(x,y,+,x,y2); title(二次拟合曲线)； f3=polyfit(x,y0,4) y3=polyval(f3,x); plot(x,y,+,x,y3) title(四次拟合曲线) f3 = 0.0000 1.0000 -6.0000 5.0000 -2.5140 三次拟合比较符合，通过观察二次与四次多项式拟合图像知多项式次数越高 拟合越符合。 2. 用电压V=10伏的电池给电容器充电，电容器上t时刻的电压为v(t)=v-(v-v)e, 其中V是电容器的初始电压, 是充电常数。试由下面一组t，V数据确定V， 。 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 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,bx,t,v2,k:) xlabel(时间t(s)), ylabel(充电电压(V)); title(电容器充电电压与时间t的曲线); f = -0.2835 1.4766 t0 = 3.5269 v0 =