scshuju2.m

最小二乘法和主成分回归法对同一数据进行分析

function[X1 Y1 X2 Y2]=scshuju2(n,e)

%产生100个（x，y）数据
%先通过B来构造y，再反过来验证B是否与设定值一致
B=[1 2 3 4 5 6 7 8 9 10]'; %B=[a0 a1 a2 a3 a4 a5 a6 a7 a8 a9]
%e=randn(1,n)';%e是白噪声1、2、5、10、20、50、100、200、500、1000、2000、5000
%e=randn(1,n); e=e/std(e); e=e-mean(e); t=0; s=sqrt(15); e=t+s*e;e=e';
x1=rand(n,1);x3=chi2rnd(5,n,1);x2=0.4*x1+1.2*x3;%根据前三个x的关系建立,x1为rand随机分布，x3为自由度为5的 χ 2 随机分布
x4=raylrnd(5,n,6);%产生rayleigh分布的随机数,从x4到x9
X0=[x1 x2 x3 x4 ones(n,1)];
X1=X0([1:n/2],:);%训练集
X2=X0([(n/2+1):n],:);%校正集
Y0=X0*B+e;
Y1=Y0([1:n/2],:);
Y2=Y0([(n/2+1):n],:);
t1=1:50;