(优雅的使用）matlab脚本编程习题5

第5次作业：回归分析

虽然这种题，代码不会差的太多，但是我还是建议借鉴的时候，先理解我的代码，再自己重写一遍，写的时候尽量不要对着我的看。有问题可以评论区问我~

1.一金属块加热到300°F，在接下来7小时内的温度如下表：
编程完成如下功能：

①画出time(时间)-temp(温度)散点图（scatter）；

②使用三次多项式拟合、在同一坐标上用plot画出time(时间)-temp(温度)曲线（带legend）；

③求9小时的温度

2022.9.30更新，期末考完才发现这里一直弄错了，要三次拟合，我给当成三项了，应该再加个time^3的项

time=[0,0.5,1.0,1.5,2.0,2.5,3.0,3.5,4.0,5.0,6.0,6.1,7.0];
temp=[300,281,261,244,228,214,202,191,181,164,151,149,141];
scatter(time,temp);
Time=[ones(13,1) time' (time.^2)'];
[b,bint,r,rint,stats]=regress(temp',Time);
y=b(1)+b(2)*time+b(3)*time.^2;
plot(time,y,'b-')
hold on
plot(time,temp,'r.')
legend('y=b1+b2*time+b3*time^2','temp')
y1=b(1)+b(2)*9+b(3)*81;

fprintf('9小时的温度预计为：%.2f',y1)

2.“综合打分” 是 去年 体育老师根据这15名同学的体重、肺活量、50m短跑、1分钟仰卧起坐、跳远成绩、1000米成绩、1分钟跳绳、引体向上、坐位体前屈等等数据综合评价打出的分数。今年 因为学校器材有限，体育老师只测了这15名同学的三项指标：跳远成绩、1000米成绩、1分钟跳绳。现在体育老师想要知道这三项指标能不能 线性合成 成最后的 今年 的综合分数?

原题来自：这篇文章
因为是多项式拟合，画图就注意把每项画出来

x1=[180,201,205,208,213,217,218,222,226,230,233,238,240,2242,253];
    x2=[280,240,226,224,220,217,225,221,211,213,199,198,195,186,183];
    x3=[153,170,162,160,162,165,170,168,169,179,172,172,175,181,176];
Y=[60,75,70,70,75,75,85,80,80,85,90,90,90,95,95]';
X=[ones(15,1),x1',x2',x3'];
[b,bint,r,rint,stats]=regress(Y,X);
y=[60,75,70,70,75,75,85,80,80,85,90,90,90,95,95];
y2=b(1)+b(2).*x1+b(3).*x2+b(4).*x3;
plot(x1,'b.')
hold on
plot(x2,'b*')
plot(x3,'k.')
plot(y,'g.')
plot(y2,'r-')
axis([0 15 0 300])

legend('跳远成绩','1000米成绩','1分钟跳绳个数','综合成绩','拟合成绩')