主成分分析法编程

丘比特

clear
clc
8 j  t8 ]# W" l, O: C% 读入人口数据（1971－2000年）
* J8 y5 r: y, f7 q# iY=[33815        33981        34004        34165        34212        34327        34344        34458        34498        34476        34483        34488        34513        34497        34511        34520        34507        34509        34521        34513        34515        34517        34519        34519        34521        34521        34523        34525        34525        34527]
5 R* M; n* J9 @% 读入时间变量数据（t＝年份－1970）
T=[1        2        3        4        5        6        7        8        9        10        11        12        13        14        15        16        17        18        19        20        21        22        23        24        25        26        27        28        29        30]
% I# q9 `  E- N* k  ?  D. s% 线性化处理
0 k. T# u0 L7 i' A9 Dfor t = 1:30,
& w3 s, _' L4 m   x(t)=exp(-t);
   y(t)=1/Y(t);
end
% 计算，并输出回归系数B
' n/ v' I5 |) |1 F: |c=zeros(30,1)+1;
8 m) v/ Q0 F7 |, M9 {) B6 iX=[c,x'];
B=inv(X'*X)*X'*y'
8 ^; w# `/ P! |- z- E8 H) |1 D& Wfor i=1:30,
% 计算回归拟合值   
, C% X3 M- x3 R# V) ]: Q3 y    z(i)=B(1,1)+B(2,1)*x(i);
; y: \% ?6 {) X: l% 计算离差
* o5 `3 r7 `$ n+ @) K2 O    s(i)=y(i)-sum(y)/30;
" v9 }+ `7 O4 ^( h% 计算误差   
+ ?. w- s4 P3 v0 K8 _, a, M( I" A  n    w(i)=z(i)-y(i);
7 a4 W+ ]- R( h9 i% s7 X$ e& l$ Oend
6 p9 Y2 s# R0 K% 计算离差平方和S
S=s*s';
# `5 [  ]( W) U" B5 u% 回归误差平方和Q
7 G7 \1 w7 s( w: |Q=w*w';
% 计算回归平方和U
U=S-Q;
! k- G4 {* b! z) s! O/ J% 计算，并输出F检验值
F=28*U/Q
% 计算非线性回归模型的拟合值
for j=1:30,
    Y(j)=1/(B(1,1)+B(2,1)*exp(-j));
6 ^9 X* Q+ T7 i; ^8 N+ @9 dend
% 输出非线性回归模型的拟合曲线（Logisic曲线）
plot(T,Y)

合欢花

这是用的SPSS软件的程序？还是什么呀

臺戶

都值得研究研究 谢谢楼主先

一米阳光的ta

感谢楼主分享。。。

CharlesLei

研究研究，多谢分享

呆呆秀@

谢谢楼主分享!!!