题目

假设向量y=[7613.51  7850.91  8381.86  9142.81 10813.6 8631.43 8124.94 9429.79 10230.81 10163.61 9737.56 8561.06 7781.82 7110.97];

x1=[7666 7704 8148 8571 8679 7704 6471 5870 5289 3815 3335 2927 2758 2591];

x2=[16.22 16.85 17.93 17.28 17.23 17 19 18.22 16.3 13.37 11.62 10.36 9.83 9.25];

是一个表格中三个列向量。

建立的基本模是:

y=Ax1^2+Bx2^2+Cx1+Dx2+Ex1*x2;

利用regress函数求解各个系数的估计值。

分析

regress命令——用于一元及多元线性回归,本质上是最小二乘法。在Matlab中,输入help regress ,会弹出和regress的相关信息:

调用格式:

B = regress(Y,X)
[B,BINT] = regress(Y,X)
[B,BINT,R] = regress(Y,X)
[B,BINT,R,RINT] = regress(Y,X)
B,BINT,R,RINT,STATS] = regress(Y,X)
[...] = regress(Y,X,ALPHA)

参数解释:

B:回归系数,是个向量(“the vector B of regression coefficients in the  linear model Y = X*B”)。
BINT:回归系数的区间估计(“a matrix BINT of 95% confidence intervals for B”)。
R:残差( “a vector R of residuals”)。
RINT:置信区间(“a matrix RINT of intervals that can be used to diagnose outliers”)。
STATS:用于检验回归模型的统计量。有4个数值:判定系数R^2,F统计量观测值,检验的p的值,误差方差的估计。
ALPHA:显著性水平(缺少时为默认值0.05)。

解题

>> y=[7613.51  7850.91  8381.86  9142.81 10813.6 8631.43 8124.94 9429.79 10230.81 10163.61 9737.56 8561.06 7781.82 7110.97]';
>> x1=[7666 7704 8148 8571 8679 7704 6471 5870 5289 3815 3335 2927 2758 2591]';
>> x2=[16.22 16.85 17.93 17.28 17.23 17 19 18.22 16.3 13.37 11.62 10.36 9.83 9.25]';
>> X=[ones(size(y)) x1.^2 x2.^2 x1 x2 x1.*x2];   %size(y)求y的长度,ones全1阵
>> [b,bint,r,rint,stats] = regress(y,X);
>> scatter3(x1,x2,y,'filled')

散列图为:

      

                                

b为对应的参数 b(1)为F(最后那个常数项) ,b(2)为A(第一个参数),b(3)为B,b(4)为C,b(4)为D,b(5)为E。bint为b的95%置信区间。

stats的第三个参数为F检测的P值,p值很小(P<0.001),说明拟合模型有效。(这里的分析技巧,以后有时间的话,会专门总结。)

附加

对散点数据进行拟合

>>  hold on
>>  x1fit = min(x1):100:max(x1);
>>  x2fit = min(x2):1:max(x2);
>> [X1FIT,X2FIT] = meshgrid(x1fit,x2fit);
>> YFIT=b(1)+b(2)*X1FIT.^2+b(3)*X2FIT.^2+b(4)*X1FIT+b(5)*X2FIT+b(6)*X1FIT.*X2FIT;
>>  mesh(X1FIT,X2FIT,YFIT)
>> view(10,10)
>> xlabel('x1')
>>  ylabel('x2')
>> zlabel('y') 

截图为

参考文献:regress命令篇

MATLAB回归分析命令——regress命令相关推荐

  1. MATLAB regress命令

    1 regress命令 用于一元及多元线性回归,本质上是最小二乘法.在Matlab 命令行窗口输入help regress ,会弹出和regress的相关信息,一一整理. 调用格式: B = regr ...

  2. matlab stats里的f值,MATLAB 回归分析regress,nlinfit,stepwise函数

    回归分析 1.多元线性回归 在Matlab统计工具箱中使用命令regress()实现多元线性回归,调用格式为 b=regress(y,x) 或 [b,bint,r,rint,statsl = rege ...

  3. matlab engevalstring,Matlab 執行engEvalString命令介紹 | 學步園

    MATLAB與VC混合編程指導書 在調用 Matlab 引擎之前,首先應在相關文件中加入一行: #include "enging.h" ,該文件包含了引擎 API 函數的說明和所需 ...

  4. Matlab绘图基本控制命令

    图形的控制与表现 (Figure control and representation) MATLAB提供的用于图形控制的函数和命令:    axis:  人工选择坐标轴尺寸.     clf:清图形 ...

  5. 【Matlab】使用Matlab运行Windows命令行命令+实例

    可以使用Matlab的一些命令来帮助程序运行.比如说 ! calc % 打开计算器 ! mspaint % 打开画图 dos calc % 打开计算器 比如一个程序要运行很长时间,而我们又不能一直守在 ...

  6. Matlab最实用画图命令整理(包括Print输出SCI论文高清大图!)

    Matlab最实用画图命令整理 1.关于 gcf 和 gca 2.画多个子图(subplot),以4个子图为例 3.关于figure的一些设置命令整理 4.关于坐标轴显示的设置 5.imagesc函数 ...

  7. matlab中的priniu,MATLAB中的plotyy命令

    MATLAB中的plotyy命令 x = 0:0.01:20; y1 = 200*exp(-0.05*x).*sin(x); y2 = 0.8*exp(-0.5*x).*sin(10*x); [AX, ...

  8. matlab方程近似求根,第七讲MATLAB中求方程的近似根(解)教学目的学习matlab中求根命令.doc...

    第七讲MATLAB中求方程的近似根(解)教学目的学习matlab中求根命令 第七讲 MATLAB中求方程的近似根(解) 教学目的:学习matlab中求根命令,了解代数方程求根求解的四种方法,即图解法. ...

  9. shell 调用matlab,在matlab中调用linux命令 以及后台运行matlab --zz

    由于matlab和linux命令各有所长:matlab比较擅长数值计算,linux命令比较擅长文件操作(matlab也基本上能够完成,但会稍微麻烦一点).所以在很多时候会需要将二者结合起来. 可能更多 ...

最新文章

  1. Mac OS build caffe2 Error:This file was generated by an older version of protoc which is
  2. 《那些年啊,那些事——一个程序员的奋斗史》——59
  3. 在fedora23中安装virtualbox, 然后实现虚拟机irtualbox 或者 vmware 下的xp操作系统
  4. android 相机 全功能,一加7系首个Android 11公测代码暗示了相机应用的诸多功能更新...
  5. html按钮按下效果_CSS+HTMLlt;水滴按钮效果gt;
  6. spike 序列matlab,SPKtool1.0.1 神经信号spike 分类及处理 工具包 matlab 238万源代码下载- www.pudn.com...
  7. 什么叫缺位_一位父亲能够给予孩子的最好礼物是什么?看看这部电影给出的答案...
  8. 在线时序流程图制作工具
  9. mex 1 hello,world
  10. 筛选出c语言成绩前五名的记录,c语言课程论文设计成绩记录簿.doc
  11. 关于标题的多HTML定制
  12. Activity生命周期详解二
  13. 14.性能之巅 洞悉系统、企业与云计算 --- 附录
  14. Java入门到精通全教程
  15. keras深度学习之猫狗分类三(特征提取)
  16. Axure Share ——原型设计工具 Axure ,移动版
  17. python可视化47|最常用5个排序(Ranking)关系图
  18. linux制作U盘启动盘无法使用,完美:使用syslinux制作U盘启动盘
  19. springBoot+itext处理富文本转pdf
  20. 摄像头8mm可以看多远_家用监控摄像头镜头焦距该如何选择_弱电工程

热门文章

  1. 配音秀显示服务器溜去配音啦,配音秀作品恢复方法操作介绍
  2. CSP在线考试环境 | OBS录屏软件下载安装和设置教程
  3. C语言内存分布图----栈空间、堆空间
  4. C语言函数实现ascii码快速查询
  5. Procmon打开后无反应的问题
  6. 自定义封包协议c语言,Socket分包,封包,粘包
  7. 隐私计算头条周刊(12.4-12.10)
  8. html2canvas微信头像没绘制,解决使用canvas生成含有微信头像的邀请海
  9. 2019年全国大学生数学建模E题
  10. 单用户MIMO系统(二):信道信息在发端未知