没有合适的资源?快使用搜索试试~ 我知道了~
利用 Matlab作回归分析.docx
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 197 浏览量
2022-06-23
03:57:29
上传
评论
收藏 146KB DOCX 举报
温馨提示
试读
17页
利用 Matlab作回归分析.docx利用 Matlab作回归分析.docx利用 Matlab作回归分析.docx利用 Matlab作回归分析.docx利用 Matlab作回归分析.docx利用 Matlab作回归分析.docx利用 Matlab作回归分析.docx利用 Matlab作回归分析.docx
资源推荐
资源详情
资源评论
利用 Matlab 作回归分析
一元线性回归模型:
y x ,
N(0, )
2
求得经验回归方程:
ˆ
y
ˆ
ˆ
x
统计量:
总偏差平方和:
n
,其自由度为
n
1
;
SST ( y y)
2
f
i
T
i1
回归平方和:
残差平方和:
n
,其自由度为 ;
SSR
SSE
ˆ
(y y)
f 1
R
2
i
i1
n
,其自由度为 ;
yˆ
f n 2
(y
)
2
i i
E
i1
它们之间有关系:SST=SSR+SSE。
一元回归分析的相关数学理论可以参见《概率论与数理统计教
程》,下面仅以示例说明如何利用 Matlab 作回归分析。
【例 1】为了了解百货商店销售额 x 与流通费率(反映商业活动
的一个质量指标,指每元商品流转额所分摊的流通费用)y 之间
的关系,收集了九个商店的有关数据,见下表 1.试建立流通费
率 y 与销售额 x 的回归方程。
表 1 销售额与流通费率数据
样本点 销售额 x(万元) 流通费率 y
1
2
1.5
4.5
7.0
4.8
3
4
5
6
7
8
9
3.6
3.1
2.7
2.5
2.4
2.3
2.2
10.5
13.5
16.5
19.5
22.5
25.5
【分析】:首先绘制散点图以直观地选择拟合曲线,这项工作可
结合相关专业领域的知识和经验进行,有时可能需要多种尝试。
选定目标函数后进行线性化变换,针对变换后的线性目标函数
进行回归建模与评价,然后还原为非线性回归方程。
【Matlab 数据处理】:
【Step1】:绘制散点图以直观地选择拟合曲线
x=[1.5 4.5 7.5 10.5 13.5 16.5 19.5 22.5 25.5];
y=[7.0 4.8 3.6 3.1 2.7 2.5 2.4 2.3 2.2];
plot(x,y,'-o')
输出图形见图 1。
10
5
0
5 10 15 20 25
图 1 销售额与流通费率数据散点图
根据图 1,初步判断应以幂函数曲线为拟合目标,即选择非
线性回归模型,目标函数为:
y axb (b 0)
其线性化变换公式为:
v ln y , u ln x
线性函数为:
v ln a bu
【Step2】:线性化变换即线性回归建模(若选择为非线性模型)
与模型评价
% 线性化变换
u=log(x)';
v=log(y)';
% 构造资本论观测值矩阵
mu=[ones(length(u),1) u];
alpha=0.05;
% 线性回归计算
[b,bint,r,rint,states]=regress(v,mu,alpha)
输出结果:
b =[ 2.1421; -0.4259]
表示线性回归模型 中:lna=2.1421,b=-0.4259;
v ln a bu
即拟合的线性回归模型为
y 2.1421
0.4259x
;
bint =[ 2.0614 2.2228; -0.4583 -0.3934]
表示拟合系数 lna 和 b 的 100(1-alpha)%的置信区间分别为:
[2.0614 2.2228]和[-0.4583 -0.3934];
r =[ -0.0235 0.0671 -0.0030 -0.0093 -0.0404
-0.0319 -0.0016 0.0168 0.0257]
表示模型拟合残差向量;
rint =[ -0.0700 0.0230
0.0202 0.1140
-0.0873 0.0813
-0.0939 0.0754
-0.1154 0.0347
-0.1095 0.0457
-0.0837 0.0805
-0.0621 0.0958
-0.0493 0.1007]
剩余16页未读,继续阅读
资源评论
春哥111
- 粉丝: 1w+
- 资源: 5万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功