| 查看: 6583 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求助 用matlab拟合三元线性回归方程及分析
|
|||
|
用matlab拟合三元线性回归方程，我自己弄了一个如下 x1=[23.00 23.00 23.00 23.00 23.00 28.00 28.00 28.00 28.00 28.00]'; x2=[0.30 0.30 0.30 0.30 0.30 0.10 0.10 0.10 0.10 0.10]'; x3=[0.00 0.69 1.09 1.38 1.61 0.00 0.69 1.09 1.38 1.61]'; y=[-0.91 -0.35 -0.05 0.15 0.29 -0.58 -0.04 0.23 0.41 0.54]; x=[ones(10,1) x1 x2 x3]; >> [b,bint,r,rint,stats]=regress(y,x); ??? Error using ==> or Matrix dimensions must agree. Error in ==> regress at 70 wasnan = (isnan(y) | any(isnan(X),2)); >> b,bint,stats,rcoplot(r,rint) b = 0 -0.0115 -1.7878 0.6751 bint = 0 0 -0.0168 -0.0062 -2.2782 -1.2973 0.5701 0.7800 stats = 0.9747 134.6812 0.0000 0.0059 问题： 1.我编写的程序对吗？ 2.输入程序中，总是会出现warning的那段文字，不知道怎么修改？ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB用于三元回归参数
已经有3人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- 求助SAS 多元非线性回归拟合编程
已经有3人回复
- 求助EXCEL 用LINEST函数求解三元线性回归方程
已经有3人回复
- 菜鸟求助：关于Origin多元线性回归拟合y=x1*（a*x2+b*x3+c)
已经有10人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
- 【求助】matlab曲面拟合表达式
已经有10人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
【答案】应助回帖
★ ★
臭水沟: 金币+2, 谢谢交流~~ 2012-07-11 12:12:08
臭水沟: 金币+2, 谢谢交流~~ 2012-07-11 12:12:08
|
Most likely you have too many predictor variables, and too few observations. As an alaogy, imagine trying to fit a cubic polynomial regression with only two observations. REGRESS will make a choise about which coefficients to set to zero, but it can't possibly know what you really want. 变量太多，你可以试一下，分别使用x(:,1:2), x(:,1:3), x(:,2:3)算一下，使用3列时总会出现一个b=0。 |
4楼2012-07-11 10:38:00
2楼2012-07-11 09:50:45
|
请问y写成列向量后又出现问题了，这是怎么回事啊？？ >> x1=[23.00 23.00 23.00 23.00 23.00 28.00 28.00 28.00 28.00 28.00]'; x2=[0.30 0.30 0.30 0.30 0.30 0.10 0.10 0.10 0.10 0.10]'; x3=[0.00 0.69 1.09 1.38 1.61 0.00 0.69 1.09 1.38 1.61]'; >> y=[-0.91 -0.35 -0.05 0.15 0.29 -0.58 -0.04 0.23 0.41 0.54]; >> x=[ones(10,1) x1 x2 x3]; >> [b,bint,r,rint,stats]=regress(y,x); Warning: X is rank deficient to within machine precision. > In regress at 82 >> b,bint,stats,rcoplot(r,rint) |
3楼2012-07-11 10:17:45
5楼2012-07-11 11:29:51
