3_0、最小二乘法及录像机计数器的用途问题.ppt

最小二乘拟合 补充知识 最小二乘法 最小二乘拟合的MATLAB实现 在MATLAB中，用polyfit函数来求得最小二乘拟合多项式的系数，再用polyval函数按所得的多项式计算所给出的点上的函数近似值。 polyfit函数的调用格式为： [P,S]=polyfit(X,Y,m) 函数根据采样点X和采样点函数值Y，产生一个m次多项式（最小二乘拟合中m=1）P及其在采样点的误差向量S。其中X,Y是两个等长的向量，P是一个长度为m+1的向量，P的元素为多项式系数。 polyval函数的功能是按多项式的系数计算x点多项式的值。 例子 函数polyfit是指用一个多项式函数来对已知数据进行拟合，我们以下列数据为例介绍这个函数的用法： x=0:0.1:1； y=[ -0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2 ] 为了使用polyfit，首先必须指定我们希望以多少阶多项式对以上数据进行拟合，我们选择二阶多项式进行拟合。 P= polyfit (x, y, 2) P= -9.8108 20.1293 -0.0317 函数返回的是一个多项式系数的行向量，写成多项式形式为： 为了比较拟合结果，我们绘制两者的图形： plot(X,Y,bo); hold on; plot(X,Z,‘r*-’); %给出拟合曲线 legend(‘原始数据’,‘拟合数据’); %图形标注 大作业 问题：配件厂为装配线生产若干种部件，轮换生产不同的部件时因更换设备要付生产准备费（与生产数量无关），同一部件的产量大于需求时因积压资金、占用仓库要付存贮费。该厂生产能力非常大，即所需数量可在很短时间内产出（即生产能力远大于需求）。今已知某一部件的日需求量100件，生产准备费5000元，存贮费每日每件1元。 （1）、如果不允许出现缺货，试安排该产品的生产计划：即多少天生产一次（称为生产周期），每次产量多少，可使总费用最小。 （2）、如果允许缺货，缺货时因失去销售机会而遭受损失，每天每件产品缺货损失费c3 ，但缺货数量需在下次生产（订货）时补足。试安排该产品的生产计划。 （3）、假设该配件厂已进入稳定发展状态，现有该厂某年前十个月的销售收入，试预测该厂剩下两个月的销售收入。 * 当问题的机理非常不清楚难以直接利用其他知识来建模时，一个较为自然的想法是直接利用数据进行曲线拟合，找出变量之间的近似依赖关系即函数关系。 设经实际测量已得到n组数据（xi , yi），i=1,…, n。将数据画在平面直角坐标系中，见图。如果建模者判断这n个点很象是分布在某条直线附近，令该直线方程为y=ax+b，进而利用数据来求参数a和b。由于该直线只是数据近似满足的关系式，故 yi-(axi+b)=0一般不成立，但我们希望 最小 此式对a和b的偏导数均 为0，解相应方程组，求得： y=ax+b y O (xi ,yi) x 其中 和 分别为xi和yi的平均值 如果建模者判断变量间的关系并非线性关系而是其他类型的函数，则可作 变量替换使之转化为线性关系或用类似方法拟合。 问 题 在一次使用中录像带已经转过大半，计数器读数为 4450，问剩下的一段还能否录下1小时的节目？ 要求 不仅回答问题，而且建立计数器读数与 录像带转过时间的关系。 思考 计数器读数是均匀增长的吗？ 建模实例： 录像机计数器的用途 经试验，一盘标明180分钟的录像带从头走到尾，时间用了184分，计数器读数从0000变到6061。 录像机计数器的工作原理 主动轮 压轮 0000 左轮盘 右轮盘 磁头 计数器 录像带 录像带运动方向 录像带运动 右轮盘半径增大 右轮转速不是常数 录像带运动速度是常数 计数器读数增长变慢 问题分析 观察 计数器读数增长越来越慢！ 模型假设 录像带的运动速度是常数 v ； 计数器读数 n与右轮转数 m成正比，记 m=kn； 录像带厚度（加两圈间空隙）为常数 w； 空右轮盘半径记作 r ； 时间 t=0 时读数 n=0 . 建模目的 建立时间t与读数n之间的关系 （设v,k,w ,r为已知参数） 模型建立 建立t与n的函数关系有多种方法 1. 右轮盘转第 i 圈的半径为r+wi, m圈的总长度 等于录像带在时间t内移动的长度vt, 所以 2. 考察右轮盘面积的 变化，等于录像带厚度 乘以转过的长度，即 3. 考察t到t+dt录像带在 右轮盘缠绕的长度，有 模型建立 思 考 3