《数值计算方法》实验报告 1 温度分布的曲线拟合 学号:XX 姓名:XXX 1
实验描述 美国洛杉矶郊区11 月8 日的温度(华氏温度)如表1 所示
采用24 小时制
表 1 温度数据 时间,p
温度 时间,a
温度 1 66 1 58 2 66 2 58 3 65 3 58 4 64 4 58 5 63 5 57 6 63 6 57 7 62 7 57 8 61 8 58 9 60 9 60 10 60 10 64 11 59 11 67 午夜 58 正午 68 要求:1
线性的最小二乘拟合 2
曲线的最小二乘抛物线拟合; 3
三次样条插值拟合 4
T7 的三角多项式拟合 5
有4 个控制点的贝塞尔曲线拟合 2
实验内容 一、线性最小二乘拟合 定理5
1(最小二乘拟合曲线)设1{(,)}Nkkkx y 有N 个点,其中横坐标1{ }Nk kx 是确定的
《数值计算方法》实验报告 2 最小二乘拟合曲线 yAx B (1) 的系数是下列线性方程组的解,这些方程称为正规方程: 211111NNNkkkkkkkNNkkkkxAxBx yxANBy (2) 核心代码为: %求 方 程 组 am=b的 根 m=a\b; x 1=1:0
1:24; y 1=m(1)*x 1+m(2); %绘 图 , 其 中(x ,y )为已知点, 用红色的 星号表示, y 1为拟合曲线 plot(x ,y ,'*r',x 1,y 1) grid on legend('已知点','最小二乘拟合') 主要算法为: (1)
输入x,y; (2)
求正规方程的系数21Nkkx,1Nkkx,1Nkky,1Nkkkx y (3)
解正规方程组am=b