1、附件(fjin)：程序(chngx)1%回歸方程1系數(xsh)求解%programx1=6859 7702 8472 9422 10493;x2=2366 2475 2622 2936 3255;X=x1,x2y=3055 3385 3765 4245 4776;X=ones(5,1) x1 x2;b=regress(y,X) %回歸函數程序2計算回歸方程1的解x1=6859 7702 8472 9422 10493;x2=2366 2475 2622 2936 3255;g=-532.3585 0.3648 0.4557;for i=1:5 k(i)=g(1)+x1(i)*g(2)+x2(

2、i)*g(3);endk程序3%回歸方程2系數求解%programclc clearsh=52378 1417.4*108 275.06481*108 2.11154*108 5000 35791 331.9*108 89.55621*108 0.52307*108 4600 14436 318.2*108 150.89366*108 0.47096*108 3500 9312 395.2*108 166.91618*108 3.42308*108 4500 10085 33.8*108 8.47082*108 0.5555*108 3100 7701 312.6*108 153.68677*

3、108 1.52673*108 3400;x1=sh(:,1)./100; %對數據作了相應處理 x2=sh(:,2)./109;x3=sh(:,3)./108;x4=sh(:,4)./106;X=x1,x2,x3,x4;y=sh(:,5);X=ones(size(x1),x1,x2,x3,x4; format long eng %使用長浮點b,bint,r,rint,stxts=regress(y,X) %對數據作回歸Q=r*r;sit=Q/(6-2) 程序(chngx)4%計算(j sun)回歸方程2的解clearc=52378 1417.4*108 275.06481*108 2.111

4、54*108; 35791 331.9*108 89.55621*108 0.52307*108; 14436 318.2*108 150.89366*108 0.47096*108; 9312 395.2*108 166.91618*108 3.42308*108; 10085 33.8*108 8.47082*108 0.5555*108; 7701 312.6*108 153.68677*108 1.52673*108;x1=c(:,1)./100;x2=c(:,2)./109;x3=c(:,3)./108;x4=c(:,4)./106;g=2119.66 7.15 23.84 -5.2

5、2 5.10;for i=1:6 k(i)=g(1)+x1(i)*g(2)+x2(1)*g(3)+x3(i)*g(4)+x4(i)*g(5);endk程序(chngx)5%回歸方程3系數求解clcclearsh1=1659.18 1617.88 498.878804757.73 895.841293.82414.038814938.34 341.381118.30372.478945243.60 102.741010.35397.608915452.75 58.60994.08467.319575885.21;x1=sh1(:,1);x2=sh1(:,2);x3=sh1(:,3);x4=sh1

6、(:,4);X=x1,x2,x3,x4;y=sh1(:,5);X=ones(size(x1),x1,x2,x3,x4;format long engx,xint,r,rint,stats=regress(y,X)Q=r*r;sit=Q/(5-2)程序6%計算回歸方程3的解clearc= 1659.18 1617.88 498.87880 %對目標影響的矩陣 895.841293.82414.03881 341.381118.30372.47894 102.741010.35397.60891 58.60994.08467.31957;x1=c(:,1);x2=c(:,2);x3=c(:,3);

7、x4=c(:,4);g=637.59 -1.14 1.27 3.29 2.63;for i=1:5 %求解(qi ji)方程 k(i)=g(1)+x1(i)*g(2)+x2(1)*g(3)+x3(i)*g(4)+x4(i)*g(5);endk程序(chngx)7%作圖程序(chngx)a=1990200 1510 13 686 29 1996 2000 4839 41 1926 251 1997 2000 5160 39 2090 247 1998 3500 5425 65 2162 251 1999 3200 5854 55 2210 265 2002 4218 7703 56 2476 2

8、07 2003 4424 8472 52 2622 161 2004 4785 9422 51 2936 161;x=a(:,1);y1=a(:,2);y2=a(:,3);y3=a(:,4);y4=a(:,5);y5=a(:,6);plot(x,y1,r,x,y2,b,x,y4,g)grid onxlabel(年)ylabel(元)title()legend(學費,城鎮,鄉村);gtext(學費);gtext(城鎮);gtext(鄉村);程序8%作圖程序a=1990200 1510 13 686 29 1996 2000 4839 41 1926 251 19972000 5160 39 20

9、90 247 19983500 5425 65 2162 251 19993200 5854 55 2210 265 20024218 7703 56 2476 207 20034424 8472 52 2622 161 20044785 9422 51 2936 161;x=a(:,1);y3=a(:,4);y5=a(:,6);plot(x,y5,r,x,y3,b)grid onxlabel(年)ylabel(%)title()legend(學費(xufi)/農村,學費/城鎮);gtext(學費(xufi)/農村);gtext(學費/城鎮);程序(chngx)9%作圖程序clcclearnu

10、m1= xlsread(biao.xls, A7:H20);nian=num1(:,1);country=num1(:,3);sh=num1(:,5);jk=num1(:,6);plot(nian,country,nian,sh,nian,jk)legend(國家投入,社會辦學,捐款);gtext(國家投入);gtext(社會辦學);gtext(捐款);grid on程序10%計算加權平均數a=6859 7702 8472 9422 10493;b=2366 2475 2622 2936 3255;c=0.38 0.39 0.41 0.42 0.43;d=0.62 0.61 0.59 0.58

11、 0.57;sum1=0;for i=1:5 sum(i)=a(i)*c(i)+b(i)*d(i); sum1=sum(i)+sum1endsum ; av=sum1/5; f=sum*3*0.25表1 “十五”期間(qjin)城鄉居民可支配收入情況表 年分項目2001年2002年2003年2004年2005年城鎮居民人均收入（元）685977028472942210493農村居民人均收入（元）23662475262229363255城鎮人口比（%）3839414243農村人口比（%）6261595857人均可支配收入加權平均數（元）40744514502156606367“十五”期間人均可支

12、配收入（元）5125 表4 單位(dnwi)（元） 地區(人均GDP)(當地財政收入)(國家財政性教育經費)(社會捐資和集資辦學經費)(學費)上海523781417.397275.064812.111545000天津35791331.850789.556210.523074200-5000黑龍江14436318.2056150.893660.470963000-4000湖南9312395.2651166.916183.423084000-5000青海1008533.82228.470820.55552800-3300云南7701312.6490153.686771.526732800-4000

13、以Y為被解釋(jish)變量即學費，以 ，為解釋變量建立多元回歸模型如下：假設Y為隨機變量(su j bin lin)，且服從正態分布，列出下列回歸方程：因為(yn wi)我們選擇了上海、天津、黑龍江、湖南、青海、云南六個地區(dq)。所以i取16。用最小二乘法對模型求解:， ， 表7 單位（元） 院校助學貸款獎學金勤工儉學住宿費(學費)元低等收入家庭1659.181617.88498.878804757.73中低等收入家庭895.841293.82414.038814938.34中等收入家庭341.381118.30372.478945243.60較高等收入家庭102.741010.3539

14、7.608915452.75高等收入家庭58.60994.08467.319575885.21同樣，以Y為被解釋變量，以助學貸款，獎學金，勤工儉學，住宿費，建立多元回歸模型：我們先假設Y為隨機變量，且服從正態分布，列出下列方程。因為我們考慮五類學生家庭情況，所以i取之為1、2、3、4、5。， ， 模型的求解(qi ji)方法同上代入數據(shj)用Matlab軟件(run jin)編程求解（詳見附件：程序5）得出： 代入回歸方程：y=637.59-1.14+1.27+3.29+2.63將變量代入得出以下結果 (詳見程序附件：程序6)：表8 單位（元）低等收入家庭中低等收入家庭中等收入家庭中高等

15、收入家庭高等收入家庭學費4756.55350.25879.86226.66679.9附錄一 模糊綜合分析程序disp(請輸入判斷矩陣A(n階);A=input(A=);n,n=size(A);x=ones(n,100);y=ones(n,100);m=zeros(1,100);m(1)=max(x(:,1);y(:,1)=x(:,1);x(:,2)=A*y(:,1);m(2)=max(x(:,2);y(:,2)=x(:,2)/m(2);p=0.0001;i=2;k=abs(m(2)-m(1);while kp i=i+1; x(:,i)=A*y(:,i-1); m(i)=max(x(:,i); y(:,i)=x(:,i)/m(i); k=abs(m(i)-m(i-1);enda=sum(y(:,i);w=y(:,i)/a;t=m(i);disp(權向量);disp(w);disp(最大特征值);disp(t); %以下是一致性檢驗CI=(t-n)/(n-1);RI=0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 1.59;CR=CI/RI(n);if CR0.10 disp(此矩陣的一