线性拟合matlab程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

clc
clear all
x=[2002 2003 2004 2005 2006 2007 2008 2009 2010]';
y=[1573.52 1813.38 2129.33 2462.22 2644.45 2812.12 2938.74 3037 3217.06]';
X=[ones(length(x),1),x];                %构造自变量观察值矩阵
[b,bint,r,rint,stats]=regress(y,X);     %线性回归建模与评价
b,stats
%作图
figure;
plot(x,y,'-+');                          %数据散点图
legend('第二产业能源消费量');
figure(2)
rcoplot(r,rint);                        %作残差分析图
hold on;
figure(3)
x0=[x' 2011 2012 2013 2014 2015]';
z=b(1)+b(2)*x0;
plot(x,y,'*',x0,z,'-r^');
legend('第二产业能源消费量预测图','第二产业能源消费量回归线图');             %预测及作回归线图
%检验
SSE=sum((y-(b(1)+b(2)*x)).^2)           %残差平方和
STD=sqrt(SSE/(length(x)-2))             %标准误差
DELTA=2*STD                             %0.05显著水平下的边际误差
ci=[z-DELTA,z+DELTA]                    %0.95置信区间

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：MATLAB程序 MATLAB matla atlab Lab 能源 建模 matlab legend figure

不错。

不错，借鉴借鉴

如果的滞后的X，那么如何构造观察矩阵呢，如有X（—2），X（-4）呢！

我也想问与楼上相同的问题！！！