【導(dǎo)讀】本文對(duì)針對(duì)某地區(qū)私家車保有量增長(zhǎng)及調(diào)控問(wèn)題進(jìn)行了探討和研究。對(duì)于問(wèn)題一，要求分析影響汽車保有量的因素并對(duì)汽車保有量進(jìn)行預(yù)測(cè)。析，最后預(yù)測(cè)出該地區(qū)2020年的私人汽車保有量參考值為282萬(wàn)。加，年存款利率增加可使該地區(qū)私人汽車保有量減少。我們對(duì)各項(xiàng)廢氣排放進(jìn)行量化，建立各項(xiàng)廢棄排放指標(biāo)與公交車和私。大限度廢氣排放量建立關(guān)系，然后根據(jù)模型以及該地區(qū)的實(shí)際情況進(jìn)行調(diào)控。保有量作了預(yù)測(cè)，建立了較為科學(xué)嚴(yán)謹(jǐn)?shù)挠?jì)量經(jīng)濟(jì)模型及主成分分析預(yù)測(cè)模型，而事實(shí)表明，隨著中國(guó)人均GDP的穩(wěn)健增長(zhǎng)，近年來(lái)，我國(guó)的家用。汽車銷量以兩位數(shù)的增速急劇擴(kuò)大。少，與經(jīng)濟(jì)發(fā)展程度、居民收入以及道路建設(shè)等有著密切的聯(lián)系。私人汽車保有量與一個(gè)國(guó)家或地區(qū)的社會(huì)經(jīng)濟(jì)發(fā)展的有關(guān)數(shù)據(jù)有著密切關(guān)。中數(shù)據(jù)除因素城市交通干線噪音均值以外，其它的數(shù)據(jù)呈高度相關(guān)性。總體有逐年遞增的趨勢(shì)。然后考慮人均國(guó)內(nèi)生產(chǎn)總值、全社會(huì)固定資產(chǎn)投資總額等因素，對(duì)汽

　　

【正文】 0 07 6 99 9 6 7 4 4 9 5 55 相關(guān)系數(shù)矩陣： x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x1 0 5 6 6 7 0 9 5 2 5 x2 5 0 6 6 1 9 7 3 9 2 x3 6 6 0 0 0 4 4 4 4 7 x4 6 6 0 0 1 8 3 9 6 4 x5 7 1 0 1 0 3 7 3 1 5 x6 0 9 4 8 3 0 2 9 6 5 28 x7 9 7 4 3 7 2 0 5 7 8 x8 5 3 4 9 3 9 5 0 6 6 x9 2 9 4 6 1 6 7 6 0 2 x10 5 2 7 4 5 5 8 6 2 0 附表 y = + 1x 2R = () () DW = ； y = + 2x 2R = () () DW = ； y = + 3x 2R = () () DW = ； y = + 4x 2R = () () DW = y = + 5x 2R = 0. 929 () () DW = ； y = + 7x 2R = () () DW = ； y = + 8x 2R = () () DW = ； y = + 9x 2R = () () DW = ； y = + 10x 2R = () () DW = ； y = + 11x 2R = () () DW = 。 附錄 function std=cwstd(vector) vector=[27000 2658 15948 737 16316 29 30619 325 2763 17130 789 18600 33282 423 2801 17866 894 19886 33689 2887 18961 1015 20249 41020 2920 19688 1198 21626 43344 3495 24465 1361 23544 46030 3495 24814 1710 24941 13 53887 4885 36149 2100 25936 59271 5376 43008 2314 26596 29 64507 6091 51946 2500 28494 70597 7305 65745 2614 29628 79221 1345 8188 70346 2897 30063 113 ]。 cwsum=sum(vector,1)。 %對(duì)列求和 [a,b]=size(vector)。 %矩陣大小 ,a為行數(shù) ,b為列數(shù) for i=1:a for j=1:b std(i,j)= vector(i,j)/cwsum(j)。 end end xlswrite(39。F:\subei\39。,std) 附錄 模型程序 clc for x=2020:2020 t0=1995。 s1=log10()+log10()+log10()+log10()。 s2=log10()+log10()+log10(13)+log10()。 s3=log10(29)+log10()+log10()+log10(113)。 m=4。 b=((s3s2)/(s2s1))^(1/m)。 a1=((s2s1)*(b1)/(b*(b^m1)^2))。 a=10^a1。 k1=(1/m*(s1a1*b*(b^m1)/(b1)))。 k=10^k1。 30 y=k*a^(b^(xt0)) end 附錄 97 年 07 年保有量觀察值與預(yù)測(cè) 對(duì)比的 matlab 程序 clc clear x=[1 2 3 4 5 6 7 8 9 10 11 12]。 y=[ 13 29 113]。 hold on plot(x,y,39。*39。) xx=0:11 y1=*.^(.^x)。 hold on plot(xx,y1) title(39。97年 07年保有量觀察值與預(yù)測(cè)值對(duì)比圖 39。)。 xlabel(39。年份 39。)。 ylabel(39。私人汽車保有量（萬(wàn)輛） 39。)。 legend(39。觀察值 39。,39。預(yù)測(cè)值 39。)。 grid on 附錄 clc clear for x1=1:10 x0=input(39。nianfen:39。) 。 x=x01995。 y1=+*x+*x^2+*x^3 y2=*x+*x^*x^3 y3=+*x+*x^*x^3 y5=+**x^2+*x^3 y9=+*x+*x^*x^3 y10=*x+*x^*x^3 y11=+**x^2+*x^3 zi =+*y1+**y3+***y10+*y11 end 附錄 function result=cwfac(vector)。 vector=[ 31 1 1 1 1 1 1 1 1 1 1]。 fprintf(39。相關(guān)系數(shù)矩 陣 :\n39。) std=CORRCOEF(vector) %計(jì)算相關(guān)系數(shù)矩陣 fprintf(39。特征向量 (vec)及特征值 (val)： \n39。) [vec,val]=eig(std) %求特征值 (val)及特征向量 (vec) newval=diag(val) 。 [y,i]=sort(newval) 。 %對(duì)特征根進(jìn)行排序， y為排序結(jié)果， i為索引 fprintf(39。特征根排序： \n39。) for z=1:length(y) newy(z)=y(length(y)+1z)。 end fprintf(39。%g\n39。,newy) rate=y/sum(y)。 fprintf(39。\n貢獻(xiàn)率： \n39。) newrate=newy/sum(newy) sumrate=0。 newi=[]。 for k=length(y):1:1 sumrate=sumrate+rate(k)。 newi(length(y)+1k)=i(k)。 32 if sumrate break。 end end %記下累積貢獻(xiàn)率大 85%的特征值的序號(hào)放入 newi中 fprintf(39。主成分?jǐn)?shù)： %g\n\n39。,length(newi))。 fprintf(39。主成分載荷： \n39。) for p=1:length(newi) for q=1:length(y) result(q,p)=sqrt(newval(newi(p)))*vec(q,newi(p))。 end end %計(jì)算載荷 disp(result) xlswrite(39。F:\subei\39。,vec) 附 錄 matlab 程序 x=[6 6 6 6 8 13 13]。 y=[20249 21626 23544 24941 25936 26596 28494 19886 29628 30063 16316 18600]。 plot(x,y,39。o39。) axis([6,14,10000,35000])。 hold on y1=polyfit(x,y,2) y2=polyval(y1,x) plot(x,y2) hold off xlabel(39。存款準(zhǔn)備金率 39。)。 ylabel(39。人均可支配收入 39。)。 legend(39。實(shí)際值 39。,39。二次擬合 39。) 附錄 matlab 程序 x0=[33689 41020 43344 46030 53887 59271 64507 33282 70597 19223 27000 30619]。 x=x039。 y=[20249 21626 23544 24941 25936 26596 28494 19886 29628 30063 16316 18600]。 plot(x,y,39。o39。) 33 hold on y1=polyfit(x,y,1)。 y2=polyval(y1,x) plot(x,y2) hold off xlabel(39。人均國(guó)內(nèi)生產(chǎn)總值 39。)。 ylabel(39。人均可支配收入 39。)。 legend(39。實(shí)際值 39。,39。一次擬合 39。)。 附錄 與居民存儲(chǔ)款余額關(guān)系的 matlab 程序 x=[ ]。 y=[ ]。 plot(x,y) xlabel(39。一年存款利率 39。)。 ylabel(39。居民存儲(chǔ)款余額 39。)。 附錄 matlab 程序 x=[33689 41020 43344 46030 53887 59271 64507 33282 70597 19223 27000 30619]。 y=[ ]。 plot(x,y) xlabel(39。人均國(guó)內(nèi)生產(chǎn)總值 39。)。 ylabel(39。居民存儲(chǔ)款余額 39。)。 附錄 2回歸模型一的 matlab 程序 X=load(39。F:\subei\39。)。 x1=X(:,1)。%人均國(guó)內(nèi)生產(chǎn)總值 q=X(:,2)。%存款準(zhǔn)備金率 y=X(:,3)。%居民儲(chǔ)蓄款余額 subplot(1,2,1) plot(x1,y) subplot(1,2,2) plot(q,y) T=[ones(12,1) x1.^3 p.^2]。 [b,bint,r,rint,stats]=regress(y,T)。%回歸分析 b,stats 附錄 2回歸模型二的 matlab 程序 X=load(39。F:\subei\39。)。 x1=X(:,1)。%人均國(guó)內(nèi)生產(chǎn)總值 q=X(:,2)。%一年存款利率 y=X(:,3)。%居民儲(chǔ)蓄款余額 plot(x1,y)%作圖分析 plot(q,y) T=[ones(12,1) x1 1./q]。 [b,bint,r,rint,stats]=regress(y,T)。%回歸分析 34 b,stats 附錄 Y=[ 1 ]。%因變量 X=[