求y=a+b1x1+b2x2+b3x3。
其实简单来说,就是给你四维空间里的一些点:
(x1,y1,z1,p1),(x2,y2,z2,p2),……,(xn,yn,zn,pn)
让你找到回归方程P = A+B * X+C * Y+D * Z。
下面的过程就是求a,b,c,d,让这个回归方程最大限度的满足以上n点!
而最大的意义就是当你把这些点代入方程时,使p和pi的方差最小!
我先解释一下:
p = A+B * x 1+C * y1+D * z 1可能不等于y1,只是y 1的一个近似值。这里是让这个近似值最接近!
即
方差=[(a+b * x 1+c * y 1+d * z 1)-p 1]2+[(a+b * x2+c * y2+d * z2)-p2]2。!
从数学分析的知识来说,你要把A,B,C,D看成未知数,然后求最小值,那么函数关于A,B,C,D的偏导数一定等于0。
(这里注意,当你对A求偏导数时,你只把A作为一个未知数,其他数都是常数,就像你对一个二次函数求导一样简单。)
这样你会得到四个方程(对A,B,C,D取偏导数,然后使偏导数= 0)。然后,你可以用软件求解四个未知数!
以下是matlab解:
& gt& gtm =[1 000001 1.9400 0.2200 12.41.00
2000156 5.2600 0.2900 12.90 4.45% '安苏'
3 000501 1.9500 0.0100 4.91 0.84% '鄂乌商A '
4 000507 2.7800 0.0600 5.46 0.98% '广东华福'
5 000532 1.4900 0.1400 8.35 1.29% '广东华电A '
6 000543 3.1700 0.1500 7.92 2.60% '万能电力'
7 000550 2.3900 0.3300 11.23 1.43% '江铃汽车'
8 000569 1.1980 0.0250 4.71.49% '川投长钢'
9 000693 2.1400 0.2000 9.36 1.55% ' MySpace '
10 000782 2.3500 0.1100 6.13 1.72% '美达股份'
11 000803 1.5100 0.0300 13.76 1.54% '美亚股份'
12 000948 3.5000 0.2900 12.15 2.30% '南天资讯'
13 600149 1.7700 0.0950 7.81.1.75% '邢台轧辊'
14 600175 4.3100 0.0400 8.90 4.61% '包华工业'
15 600179 2.4200 0.0200 5.44 0.80% '黑化股'
16 600180 3.2900 0.3300 10.93 3.33% '九发股份'
17 600237 2.5200 0.3100 9.94 2.38% '同丰电子'
18 600281 2.9030 0.0620 6.51.16% '泰华股份'
19 600286 3.1100 0.2090 9.70 2.38% '国光瓷业'
20 600292 4.7160 0.4200 13.57 4.72% '九龙电力'
21 600305 2.9968 0.2155 12.53 2.48% '恒顺醋业'
22 600306 2.4600 0.1000 10.47 1.40% '商业城'
23 600311 3.4990 0.2430 13.89 2.80% '扬农化工'
24 600486 3.4990 0.2430 13.89 2.80% '扬农化工'
25 600530 4.4900 0.2600 14.07 4.60% '交通大学'
26 600599 2.2100 0.1600 13.92 2.21% '浏阳花炮'
27 600614 1.8800 0.2330 16.34 1.77% '胶带股'
28 600616 1.8500 0.1700 12.12 2.00% '首食'
29 600618 2.4220 0.0055 11.92 1.78% '氯碱化工'
30 600619 2.6000 0.0440 10.02 2.34% '李海股份'
31 600620 1.7000 0.3020 11.16 1.53% '天辰股份'
32 600621 1.7000 0.1 375 6.92 2.00% '上海金陵'
33 600630 3.9560 0.1190 5.58 4.00% '龙头股'
34 600632 3.4600 0.2400 6.53 2.35% '华联百货'
35 600636 2.5100 0.3086 8.77 2.34% '三爱富人'
36 600637 4.2944 0.1404 12.82 3.12% '广播信息'
37 600638 3.8819 0.0965 7.85 3.60% '新黄浦'
38 600655 3.6920 0.2030 7.77 3.69% '豫园商城'
39 600662 3.2500 0.3600 8.59 3.50% '强生控股'
40 600665 1.7480 0.0120 6.05 1.10% '胡昶特钢'
41 600673 1.4900 0.0200 7.86 0.96% '梁成股份'
42 600683 1.4300 0.1400 6.71.1.08% '宁波华联'
43 600688 2.0270 0.1260 3.90 1.22% '上海石化'
44 600788 4.7700 0.1000 6.42 2.80% '达尔曼'
45 600814 1.6400 0.0090 6.00 1.10% '杭州白洁'
46 600819 3.4400 0.2270 14.05 2.7438+0% '耀皮玻璃'
47 600825 2.0190 0.4340 12.09 2.05% '华联超市'
48 600827 3.3640 0.2500 12.18 2.80% '友谊股份'
49 600874 1.3100 0.2200 7.20 1.42% '创业与环保'
50 600889 3.0800 0.2000 9.18 1.92];% '南京化纤'
& gt& gtp=M(:,6);X=[ones(50,1),M(:,3:5)];
& gt& gtb =回归(p,X);
& gt& gta=B(1)
a =
-0.6518
& gt& gtb=B(2)
b =
0.8351
& gt& gtc=B(3)
c =
1.4092
& gt& gtd=B(4)
d =
0.0343
& gt& gt
具体数据可能会有误差,但是和你的答案很接近,呵呵,我只能做这么多了。
最后,祝你工作顺利。