数学建模社区-数学中国

标题: 主成分分析法编程 [打印本页]

作者: 丘比特 时间: 2012-9-6 21:02
标题: 主成分分析法编程
clear
clc
% 读入人口数据（1971－2000年）
Y=[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]
% 读入时间变量数据（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]
% 线性化处理
for t = 1:30,
x(t)=exp(-t);
y(t)=1/Y(t);
end
% 计算，并输出回归系数B
c=zeros(30,1)+1;
X=[c,x'];
B=inv(X'*X)*X'*y'
for i=1:30,
% 计算回归拟合值
z(i)=B(1,1)+B(2,1)*x(i);
% 计算离差
s(i)=y(i)-sum(y)/30;
% 计算误差
w(i)=z(i)-y(i);
end
% 计算离差平方和S
S=s*s';
% 回归误差平方和Q
Q=w*w';
% 计算回归平方和U
U=S-Q;
% 计算，并输出F检验值
F=28*U/Q
% 计算非线性回归模型的拟合值
for j=1:30,
Y(j)=1/(B(1,1)+B(2,1)*exp(-j));
end
% 输出非线性回归模型的拟合曲线（Logisic曲线）
plot(T,Y)

作者: 合欢花 时间: 2013-8-6 10:45
这是用的SPSS软件的程序？还是什么呀

作者: 臺戶 时间: 2014-8-27 14:35
都值得研究研究谢谢楼主先

作者: 一米阳光的ta 时间: 2014-9-8 16:49
感谢楼主分享。。。

作者: CharlesLei 时间: 2014-9-9 21:55
研究研究，多谢分享

作者: 呆呆秀@ 时间: 2017-1-20 16:02
谢谢楼主分享!!!

作者: 呆呆秀@ 时间: 2017-1-20 20:51
有用啊。。。。。

作者: 呆呆秀@ 时间: 2017-1-20 20:51
谢谢楼主分享!!!

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)