数学建模社区-数学中国

标题: 急求帮助！！！ [打印本页]

作者: won19911129 时间: 2012-5-14 12:13
标题: 急求帮助！！！
恳求高手看一下这个程序错在哪里了？数据拟合求多元回归方程的。
X1=[2018.6 1354.58 2347.46 2884.62 871.51 669 1063.6 846.78 1047 86.265 870.04 810.91 625 432.4 1404.8 1270.08 1035.8 217.79 810.6 330];
X2=[427.2 180 270 275.97 160.84 88.5 70 160 205 156.89 183.25 120 133 61.6 257.9 212 190 21.94 100 35];
X3=[16410.54 11917.3 6340.5 82402.95 1159 13237 53100 9983 8488 7446.2 16596 6597 21473 8034 12390 7434.4 1953 350 12948 14000];
X4=[8205.27 3900.1 2346.2 27467.35 399.7 6612.5 18030.1 3331.5 2998.3 2482.1 5196.3 3196.4 7521.5 3624.5 4430.6 1618.1 750.3 119.2 4316.2 2800.3];
Y=[114025 89.375 92.75 70.2967 90.875 58.58333 92.5 104.9583 89.95833 102.875 93.875 102.2917 61.95833 66 87.5 59.54167 58.16667 124.625 88.41667 129.8333];
n=20;m=4;
X=[ones(n,1),x1',x2',x3,x4'];
[b,bint,r,rint,s]=regress(Y’,X,0.05);
B,bint,r,rint,s,

作者: zjqylcy 时间: 2012-5-14 18:26
好歹也要把错误提示贴出来吧

作者: won19911129 时间: 2012-5-14 18:35
标题: RE: 急求帮助！！！

zjqylcy 发表于 2012-5-14 18:26
好歹也要把错误提示贴出来吧

出现了这种情况，麻烦看看呗，嘿嘿。??? Undefined function or variable 'x1'.（? ? ?未定义的函数或变量的x1”）

作者: 0xvirtual 时间: 2012-5-15 07:39
我当时校赛很头疼，同求高手指点

作者: liuxingwuye 时间: 2012-5-15 08:20
路过、、、、

作者: won19911129 时间: 2012-5-15 09:14

zjqylcy 发表于 2012-5-14 18:26
好歹也要把错误提示贴出来吧

出现了这种情况，麻烦看看呗，嘿嘿。??? Undefined function or variable 'x1'.（? ? ?未定义的函数或变量的x1”）

作者: zjqylcy 时间: 2012-5-16 16:24
matlab区分大小写的把x1改成X1

作者: liwenhui 时间: 2012-5-17 12:56
试试这个：

x1=[2018.6 1354.58 2347.46 2884.62 871.51 669 1063.6 846.78 1047 86.265 870.04 810.91 625 432.4 1404.8 1270.08 1035.8 217.79 810.6 330];
x2=[427.2 180 270 275.97 160.84 88.5 70 160 205 156.89 183.25 120 133 61.6 257.9 212 190 21.94 100 35];
x3=[16410.54 11917.3 6340.5 82402.95 1159 13237 53100 9983 8488 7446.2 16596 6597 21473 8034 12390 7434.4 1953 350 12948 14000];
x4=[8205.27 3900.1 2346.2 27467.35 399.7 6612.5 18030.1 3331.5 2998.3 2482.1 5196.3 3196.4 7521.5 3624.5 4430.6 1618.1 750.3 119.2 4316.2 2800.3];
Y=[114025 89.375 92.75 70.2967 90.875 58.58333 92.5 104.9583 89.95833 102.875 93.875 102.2917 61.95833 66 87.5 59.54167 58.16667 124.625 88.41667 129.8333];
n=20;m=0.05;
X=[ones(n,1),x1',x2',x3',x4'];
[b,bint,r,rint,s]=regress(Y',X,m)
rcoplot(r,rint);

运行结果：

b =

  1.0e+004 *

-1.8715
-0.0018
0.0227
-0.0003
0.0011

bint =

  1.0e+004 *

-3.6451 -0.0979
-0.0044 0.0008
0.0058 0.0396
-0.0007 -0.0000
0.0001 0.0020

r =

  1.0e+004 *

3.9949
0.0639
-0.4277
-0.7673
-0.2531
-1.4972
0.8595
-0.4295
-1.2417
-1.6630
-0.6819
-0.5679
-0.8065
0.1067
-2.0157
0.0950
-0.7306
1.7648
0.8153
3.3818

rint =

  1.0e+004 *

3.9925 3.9973
-3.5827 3.7105
-3.0295 2.1741
-2.6854 1.1508
-3.9239 3.4178
-4.3178 1.3234
-2.1716 3.8906
-4.1414 3.2825
-4.8625 2.3790
-4.5270 1.2011
-4.2461 2.8823
-4.1287 2.9929
-4.3757 2.7626
-3.3590 3.5725
-5.4481 1.4166
-3.3651 3.5552
-4.3713 2.9102
-1.5273 5.0568
-2.8350 4.4657
0.9033 5.8604

s =

  1.0e+008 *

0.0000 0.0000 0.0000 3.0362

残差图：