• 26.00 KB
  • 2021-04-22 发布

MATLAB《数学实验》报告14

  • 2页
  • 当前文档由用户上传发布,收益归属用户
  1. 1、本文档由用户上传,淘文库整理发布,可阅读全部内容。
  2. 2、本文档内容版权归属内容提供方,所产生的收益全部归内容提供方所有。如果您对本文有版权争议,请立即联系网站客服。
  3. 3、本文档由用户上传,本站不保证质量和数量令人满意,可能有诸多瑕疵,付费之前,请仔细阅读内容确认后进行付费下载。
  4. 网站客服QQ:403074932
‎  MATLAB《数学实验》报告14‎ ‎  学号 1012011 姓名 成绩 实验内容:实验数据的拟合 一 实验目的 ‎  学会MATLAB软件中利用给定数据进行拟合运算的方法。 二 预备知识 ‎  (1)在MATLAB中,有下面几个命令与拟合相关,他们的含义和调用方法如下:‎ ‎  P=ployfit(x,y,n), C=lsqcurvefit(fun,c0,x,y)‎ ‎  (2) lsqcurvefit用作各种类型曲线的拟合,用最小二乘法寻找符合经验公式的最有曲线。可用非线性函数的你数据拟合。 三 实验内容与要求 (1) 1‎ ‎  ,‎ ‎  已 ‎  知 ‎  数 ‎  据 ‎  x=[1.2,1.4,1.8,2.1,2.4,2.6,3.0,3.3],y=[4.85,5.2,5.6,6.2,6.5,7.0,7.5,8.0],求对x与y进行一次、二次拟合的拟合函数 Matlab命令 结果 1‎ ‎  x=[1.2,1.8,2.1,2.4,2.6,3.0,3.3]; y=[4.85,5.2,5.6,6.2,6.5,7.0,7.5] p1=polyfit(x,y,4); p2=polyfit(x,y,5); p3=polyfit(x,y,6); disp('′aoú'),poly2str(p1,'x') disp('′aoú'),poly2str(p2,'x') disp('áù′aoú'),poly2str(p3,'x') x1=0:.1:3;y1=polyval(p1,x1);y2=polyval(p2,x1); y3=polyval(p3,x1); plot(x,y,'rp',x1,y1,'--',x1,y2,'k-',x1,y3,'*-y') legend('?ao?μ?','′ao?','′ao?','áù′ao?') exam 50拟合点四次拟合五次拟合六次拟合0-50-100-150-200-250-30000.511.522.533.5 (2) 假定某天的气温变化记录如下表,试用最小二乘法找出这一天的气温变化规律 ‎  t/h t/h 0 13 1 14 14 32 2 14 15 31 3 14 16 29 4 14 17 27 5 14 18 25 6 15 19 24 7 16 20 22 8 18 21 20 9 20 22 18 10 22 23 17 11 23 24 16 12 25 T/`C 15 T/`C 31 考虑下列函数,计算误差平方和,并作图比较效果: (1)二次函数 (2)三次函数 (3)四次函数 Matlab命令 结果 2‎ ‎  clear;clc t=0:1:24; T=[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]; %?t′?oˉêy?ao? p2=polyfit(t, T, 2); T2=polyval(p2, t); p2 deltaT2=sum((T2-T).*(T2-T)) %èy′?oˉêy?ao? p3=polyfit(t, T, 3); T3=polyval(p3, t); p3 deltaT3=sum((T3-T).*(T3-T)) %′?oˉêy?ao? p4=polyfit(t, T, 4); T4=polyval(p4, t); p4 deltaT4=sum((T4-T).*(T4-T)) plot(t,T,'--k',t,T2,'*-y',t,T3,'-.rd',t,T4,':go'); legend('1?2a?μ','?t′?oˉêy?ao?','èy′?oˉêy?ao?','′?oˉêy?ao?'); title('1?2a?μoí?÷′aoμμ?±è'); 505202535观测值和各次拟合值的比较观测值二次函数拟合三次函数拟合四次函数拟合30151010152025 3‎