利用matlab做回归分析

合金强度y与其中含碳量x有密切关系，如下表

xx	0.10	0.11	0.12	0.13	0.14	0.15	0.16	0.17	0.18	0.20	0.21	0.23
yy	42.0	41.5	45.0	45.5	45.0	47.5	49.0	55.0	50.0	55.0	55.5	60.5

根据此表建立y(x)y(x)。并对结果作可信度进行检验、判断xx对yy影响是否显著、检查数据中有无异常点、由xx的取值对yy作出预测。

解：在x−yx-y平面上画散点图，直观地知道y与x大致为线性关系。用命令polyfit(x,y,1)可得y=140.6194x+27.0269y=140.6194x+27.0269。

x=0.1:0.01:0.18;x=[x,0.2,0.21,0.23]';
y=[42,41.5,45,45.5,45,47.5,49,55,50,55,55.5,60.5]';
polyfit(x,y,1)
plot(x,y,'r*')

作回归分析用命令[b,bint,r,rint,ststs]=regress(y,x,alpha)

可用help查阅此命令的具体用法

残差及置信区间可以用

rcoplot(r,rint)
x=0.1:0.01:0.18;x=[x,0.2,0.21,0.23]';
y=[42,41.5,45,45.5,45,47.5,49,55,50,55,55.5,60.5]';
X=[ones(12,1),x];
[b,bint,r,rint,stats]=regress(y,X,0.05);
b,bint,stats,rcoplot(r,rint)

结果含义为
β0=27.0269,β1=140.6194β_0=27.0269,β_1=140.6194
β0β_0的置信区间是[22.3226,31.7313][22.3226,31.7313]
β1β_1的置信区间是[111.7842,169.4546] [111.7842,169.4546]
R2=0.9219,F=118.0670,p<10−4.R_2=0.9219,F=118.0670, p

RR是衡量y与x的相关程度的指标,称为相关系数.RR越大,xx与yy关系越密切.通常RR大于0.90.9才认为相关关系成立.

FF是一统计指标,pp是与FF对应的概率,当p<0.05p时,回归模型成立.

此例中p=0<10−4<0.05p=0 ,所以,所得回归模型成立。

观察所得残差分布图,看到第8个数据的残差置信区间不含零点,此点视为异常点,剔除后重新计算。

此时键入：

X(8,:)=[];y(8)=[];
[b,bint,r,rint,stats]=regress(y,X);
b,bint,stats,rcoplot(r,rint)X(8,:)=[];y(8)=[];[b,bint,r,rint,stats]=regress(y,X);b,bint,stats,rcoplot(r,rint)

得:
b=27.0992 137.8085b =27.0992\ 137.8085
bint=23.8563 30.342 117.8534 157.7636 bint = 23.8563 \ 30.342 \ 117.8534 \ 157.7636
stats=0.9644244.05710.0000 stats = 0.9644 \quad244.0571\quad0.0000

可以看到：置信区间缩小；R2R^2、FF变大，所以应该采用修改后的结果。建立的回归预测方程为：y=27.0992+137.8085xy=27.0992+137.8085x

利用matlab做回归分析相关推荐

简单使用matlab做回归分析、拟合分析
简单使用matlab做回归分析.拟合分析前言一元线非线性回归例子多元线性回归例子多元线非线性回归前言这里是简单做一元非线性,及多元线性. 一元:就是只有一个X.比如:y=x^2+x+1 ...
如何利用matlab做BP神经网络分析（包括利用matlab神经网络工具箱）
如何利用matlab做BP神经网络分析(包括利用matlab神经网络工具箱) 转载:https://blog.csdn.net/xgxyxs/article/details/53265318 最近一段 ...
利用matlab做dsp实验,华工DSP实验一利用Matlab进行频谱分析
实验一利用Matlab进行频谱分析一.对无限长时域离散信号进行频谱分析 1．代码: function test1(N) N = input('Type in the length = '); k = ...
我的spark学习之路（三）：利用spark做回归分析
spark的机器学习库(MLlib)下有简单的回归分析方法,今天只说最简单的线性回归,spark提供有两个回归分析库(mllib和ml),我学习的时候在网上也查了不少资料,有一个奇怪的现象是网上关于s ...
MATLAB做回归分析
根据以下帖子整理: https://blog.csdn.net/heres_/article/details/78375651 http://blog.sina.com.cn/s/blog_4a1ae ...
利用Matlab做二元huffman,信息论与编码课程作业_huffman编码的matlab_实现
信息论与编码课程作业 --霍夫曼编码求信源熵和存储前后的信息量的变化一:设计目的: 1.学习离散信源平均信息量的计算方法. 2.理解和掌握huffman 编码的基本原理,实现对信源符号的huffma ...
利用Matlab做二维图形和三维图形
二维图形 1.曲线图 Matlab作图是通过描点.连线来实现的,故在画一个曲线图形之前,必须先取得该图形上的一系列的点的坐标(即横坐标和纵坐标),然后将该点集的坐标传给Matlab函数画图. 命令为: ...
利用matlab做一个不定二重积分的计算
计算如下不定二重积分实验程序:syms x y f=(x*exp(-x*y))//进行二次简单的不定积分,分别进行计算,计算二次就好了c=int(f,'x')int(c,'y')实验结果: y=ds ...
利用matlab程序自动生成feko近场云图及matlab程序实现自动化feko运行
Feko是一款非常实用的电磁仿真软件,目前许多研究的验证及数据的制作都需要借助Feko来实现.但是对于刚入门的同学来说,feko的数据获取比较麻烦,需要在每一次运行结束后手动调节参数再执行下一个运行程 ...

利用matlab做回归分析

利用matlab做回归分析相关推荐

最新文章

热门文章