matlab經(jīng)典算法程序---回歸分析-文庫吧
2025-09-17 16:24 本頁面
【正文】 (i =1,2,… ,n), 如何利用MATLAB軟件實現(xiàn)以上的統(tǒng)計計算？ MATLAB軟件實現(xiàn) 使用命令 regress實現(xiàn)一元線性回歸模型的計算 b = regress (Y, X) 或 [b, bint, r, rint, stats] = regress(Y, X, alpha) ??????????????????????nn yyxx???11Y,11X回歸系數(shù) a， b以及它們的置信區(qū)間 殘差向量 e=YY及它們的置信區(qū)間 相關(guān)系數(shù) R2， F統(tǒng)計量和與 χ0對應的概率 p。 殘差及其置信區(qū)間可以用 rcoplot(r,rint)畫圖。 默認值是 引例求解 輸入： () x=[1097 1284 1502 1394 1303 1555 1917 2051 2111 2286 2311 2021 2435 2625 2948 3155 3372]。 y=[698 872 988 807 738 1025 1316 1539 1561 1765 1762 1960 1902 2021 2446 2736 2825]。 X=[ones(size(x39。)),x39。],pause [c,cint,r,rint,stats]=regress(y39。,X,),pause rcoplot(r,rint) 輸出： c = (參數(shù) a) (參數(shù) b) cint = ( a的置信區(qū)間 ) ( b的置信區(qū)間 ) r = [ ]’(殘差向量 ) rint=（略）（參見殘差分析圖） stats = (R2) ( F ) ( P{χ0} ) xbay ??? ??5 10 154002000200400600Residual Case Order PlotResidualsCase Number第 12個數(shù)據(jù)點異常，可刪出 預測 x1(1)=3372。() for i=1:5 x1(i+1)=*x1(i)。%未來五年國民收入以 %的 速度遞增 y1(i+1)=+*x1(i+1)。%鋼材的預 測值 end x1, y1 結(jié)果 x1 = y1 = Z 如果從數(shù)據(jù)的散點圖上發(fā)現(xiàn) y與 x沒有直線關(guān)系，又如何計算？ 例如，試分析年齡與運動（旋轉(zhuǎn)定向）能力 年齡 17 19 21 23 25 27 29 第 1人 第 2人 15 20 25 301820222426283032假設(shè)模型 ??? ????),0(~。23221???NaxaxaY 一元多項式回歸在 matlab 軟件中用命令polyfit實現(xiàn)。如前面的例子，具體計算如下： 輸入 ： () x1=17:2:29。x=[x1,x1]。 y=[ ]。 [p,S]=polyfit(x,y,2)。p 注意： x， y向量的維數(shù)要一致。 S是一個數(shù)據(jù)結(jié)構(gòu)，用于其它函數(shù)的計算。 計算 y的擬合值： 輸入： [Y,delta]=polyconf(p,x,S)。Y 結(jié)果： Y= 擬合效果圖： 15 20 25 301520253035用 polytool(x,y,2)還可以得到一個交互式畫面。 Export Parameters Parameters CI Prediction Prediction CI Residuals All ??? ????),0(~。23221???NaxaxaY在工作空間中，輸入 yhat，回車，得到預測值。 實驗內(nèi)容 確定企業(yè)年設(shè)備能力與年勞動生產(chǎn)率的關(guān)系 某市電子工業(yè)公司有 14個所屬企業(yè) ， 各企業(yè)的年設(shè)備能力與年勞動生產(chǎn)率統(tǒng)計數(shù)據(jù)如下表 。試分析企業(yè)年設(shè)備能力與年勞動生產(chǎn)率的關(guān)系 。若該公司計劃新建一個設(shè)備能力為 /人的企業(yè) ， 估計勞動生產(chǎn)率將為多少 ？ 企業(yè) 設(shè)備能力
教學課件相關(guān)推薦
鄂ICP備17016276號-1