求y=a+b1x1+b2x2+b3x3

求y=a+b1x1+b2x2+b3x3

我大概扫描了一下你的资料，
其实说简单一点，就是给你一些4维空间的点：
(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*x1+c*y1+d*z1不一定会等于y1，这只是y1的近似，下面就是要让这个近似最接近！
也就是
方差=[(a+b*x1+c*y1+d*z1)-p1]^2+[(a+b*x2+c*y2+d*z2)-p2]^2+……+[(a+b*xn+c*yn+d*zn)-pn]^2要最小！！
由数学分析的知识，你要把其中的a，b，c，d看成未知数，然后求最小值，那么关于a，b，c，d的函数对各个未知量的偏导数必须等于0

(这里再注释一下，对a求偏导数的时候，只把a看作未知数，其它的数看作常数，就像你对二次函数求导一样简单）

这样你会得到4个方程(分别是对a，b，c，d求偏导数后，再令偏导数＝0得到的)那么4个方程解4个未知数，利用软件来解就可以了！
下面是matlab解答：
>> M=[1 000001 1.9400 0.2200 12.41 1.00
2 000156 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 0.49 %'川投长钢'
9 000693 2.1400 0.2000 9.36 1.55 %'聚友网络'
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 1.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.1375 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.71 %'耀皮玻璃'
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]; %'南京化纤'
>> p=M(:,6);X=[ones(50,1),M(:,3:5)];
>> B=regress(p,X);
>> a=B(1)

a =

-0.6518

>> b=B(2)

b =

0.8351

>> c=B(3)

c =

1.4092

>> d=B(4)

d =

0.0343

>>
具体数据可能有点误差，不过和你的那个答案很接近，呵呵，我只能做这么多了

最后祝你工作顺利

<<数理统计>>，汪荣鑫编写，西安交通大学出版社出版，P198页到P208页。
上面有具体公式。
这种东西，说起来也挺麻烦，最好自己看书