erepca.asv

使用PCA算法的故障诊断MATLAB仿真程序

%原始数据集data0[],159×6.
data0=[0.05	61	65	1.15	2.5		7.7625;
0.3	60	62.5	1.1	2.4		6.8625;
0.2	59	62.5	0.9	1.7		9.5625;
0.3	63	66	0.9	1.6		7.65;
0.2	62.5	65	0.8	 2		7.65;
0.3	60	64	1	1.55		8.4375;
0.1	59	64	1.1	2		8.6625;
0.05	62.5	65	0.9	 2.1		6.4125;
0.1	58.5	63	1	2.1		8.4375;
0.15	62	63	1.05	2.1		7.4250;
0.05	58.5	63	1.2	 2.1		7.9875;
0.05	62.5	63.5	1.05	2.1		7.200;
0.15	59	63.5	0.9	1.6		8.2125;
0.05	61	65	1	1.75		6.8625;
0.35	60	63	1.2	2.65		7.875;
0.25	60	64	1	1.85		7.0875;
0.1	 58.5	62.5	0.85	1.9		8.4375;
0.15	59	63.5	0.9	1.85		7.2000;
0.1	61	64	0.9	1.85		7.7625;
0.15	56.5	63	0.95	1.8		8.325;
0.2	62	63	1	2.2		10.0125;
%12月份
0.35	58	62.5	1.05	1.9		7.3125;
0.4	58	62.5	1.05	2.2		7.5375;
0.05	59	65	1.1	 2		8.2075;
0.35	60	65	1	1.75		8.1;
0.25	55.5	63.5	1.2	 2.5		9.7875;
0.05	59.5	63	1.2	 2.1		8.2125;
0.45	59.8	63	0.8	 1.8		7.875;
0.8	 60	64.5	0.8	1.65		7.425;
0.1	61	64.8	0.9	1.7		8.55;
0.3	57.5	62.5	1.05	2.35		9.5625;
0.2	60.5	62.5	1	2.75		7.0875;
0.2	58	63	1.1	2.15		7.2;
0.9	58	62	0.9	2		7.3125;
0.3	60.5	65	1.1	2.15		6.975;
0.15	58.5	62	1	2.1		7.0875;
0.2	60	63	1.15	2.5		7.2;
0.25	60	64	1.05	1.65		6.4125;
0.15	57	64	0.95	2.65		7.2;
0.1	58.5	64	1.2	2.4		7.2;
0.15	60.5	64	0.9	2.1		7.425;
0.15	60	63	1	2.6		6.975;
0.25	60.5	64	1.1	2.1	7.2;
0.1	60.5	64	1.1	1.9		6.8625;
0.15	60.5	61	0.85	1.6		5.0625;
0.1	60	62.5	1	2.15		5.2875;
0.2	55	62.5	0.95	2.1	9.45;
0.95	57.5	61.5	1.1	1.8		6.6375;
0.5	58	61	0.8	1.9		6.4125;
0.3	57.5	62.5	0.75	1.7		6.525;
0.15	59.5	62.5	1	2		6.75;
%2003年1月份
0.35	59	64.5	0.95	1.95		7.9875;
0.1	60.5	65	1.05	1.8		7.9875;
0.35	60	64.5	0.85	1.8		8.4375;
0.1	58.5	64	1.05	2.05		8.55;
0.15	61	64	0.9	1.9		7.875;
0.05	59	65	0.9	1.7		7.9875;
0.25	59	64	1.2	2.1		8.4375;
0.2	61	63.5	0.85	1.5		7.2;
0.25	59.5	63.5	0.7	1.55		7.5375;
0.2	58	62	1.05	1.7		7.3125;
0.15	60	65	0.95	1.9		8.4375;
0.05	61.5	65	0.75	1.7		7.425;
0.35	59	64	0.9	1.6		8.1;
0.1	59	63	1	1.9		7.9875;
0.15	61	65	1.2	1.9 	6.525;
0.25	60.5	66.5	1.05	1.5		8.2125;
0.1	62	65	1	1.8		6.975;
0.15	59	62	1	2.15		7.875;
0.15	60	65	1	2		7.425;
0.05	61.5	63.5	0.7	1.9		7.5375;
0.05	60	66	1	2.45		8.775;
0.05	56.5	65	1.05	2.15	8.1;
0.05	60.5	65	1.1	1.95		7.0875;
1	60	61.5	1.3	2.65		7.3125;
1.95	57	61	1.1	2.35		8.1;
0.75	59	62.5	1.1	2.2		8.6625;
0.1	59.5	66.5	1.8	3.2		7.65;
1.15	59	61.5	1.55	2.5		6.6375;
1.25	56	61	1.2	2.1		7.5375;
0.5	59	64	1.2	2.5		7.2;
1.45	60	62.5	1.6	3.4		4.3875;
0.05	58.5	63	1.5	3.3		9.45;
0.15	59	62.5	1	2.2		6.6375;
0.05	57.5	62.5	1.8	3.6		7.9875;
0.2	57.5	62.5	1.9	3.6		8.4375;
1.4	56	63.5	1.3	2.1		5.9625;
0.05	56	65	2.05	3.8		7.3125;
0.1	55	63	1.9	3.7		7.7825;
0.5	61	65	1.35	2.8		6.8625;
0.65	57	63.5	1.5	2.3		7.65;
0.15	55.5	61.5	1.8	3.5		7.5375;
0.25	56	63	1.95	3.1		8.325;
0.1	58.5	65	1.5	3.15		7.0875;
0.65	58.5	64	1.4	2.8		6.975;
0.05	59.5	66	1.5	3.1		9.45;
0.35	60.5	66.5	0.9	 2.3 8.4375;
%12  月份 
0.25	58.5	65	1.2	2.4		7.2;
0.8	61	63.5	0.95	1.95		6.8625;
0.95	55.5	65.5	1.6	2.9		7.875;
0.5	59	63	1.1	2.2		6.4125;
0.6	53.5	64	1.6	3.1		7.5375;
0.25	56	61.5	1.95	3.3		7.3125;
0.2	57.5	60	1	2.1		6.8625;
0.45	58	64	1.4	3.6		7.9875;
0.5	57.5	63.5	1.1	2.15		7.3125;
0.4	56	63	1.6	3.3		7.0875;
0.25	55	64.5	1.6	3.3		7.525;
0.1	55	63.5	1.1	2.1		8.775;
0.4	57.5	63.5	1.25	2.4		8.4375;
0.65	53.5	62.5	1.2	2.6		8.4375;
0.6	58	63	1.6	3.7		7.425;
0.3	56	61.5	1.7	3.6		8.325;
0.3	60	65	1.25	2.8		6.4125;
0.5	57	64	1.3	2.45		7.425;
0.25	57.5	65	1.45	3		7.425;
0.45	57.5	63	1.7	3.2		7.425;
0.3	59	62	1.6	3.55		6.975;
0.05	56	65.5	1.2	1.9		6.975;
0.1	58.5	65	1.4	2.8		6.6375;
0.2	56	64.5	1.3	2.8		7.425;
0.4	55	62	1.5	3.6		8.1;
0.15	57	64	1.3	3.4		7.7625;
0.1	60	65.5	1.3	2.3		5.7375;
0.05	60	65	1.6	3.2		6.075;
0.2	58	68	1.15	2.6		6.55;
0.15	60	67	1.15	2.15		6.8625;
0.1	57	64.5	1.4	2.6		7.7625;
0.5	57.5	64	1.05	2.3		7.2;
1.1	58	63	1.05	2.2		7.2;
0.05	59.5	66	1.4	2.65		6.6375;
%2003 -1
0.15	59.5	64.5	1.3	2.15		6.4125;
0.25	59	65.5	1.4	2.4		5.625;
0.35	60	64	1.1	2.7		8.4375;
0.1	58.5	64	1.1	2.5		7.875;
0.35	61	66	1.45	2.7		7.3125;
0.1	60.5	65.5	1.2	2.3		6.075;
0.15	60	65.5	1.5	2.65		8.4375;
0.15	60.5	65	1.3	3		6.975;
0.3	59	65.5	1.2	2.2		8.6625;
0.45	58.5	64	1.05	2.35		7.9875;
0.15	58.5	65.5	1.2	2.05		7.9875;
0.2	58	64	1.35	2.35		6.975;
0.5	58	63.5	1.1	2.2		7.65;
0.1	58	64	1.1	2.05		7.425;
1.5	57.5	62.5	1.45	3.25		7.0875;
0.3	57	64	1.8	3.1		8.775;
0.1	55	64	1.45	2.9		8.1;
0.1	54	60	1.4	2.7		7.7625;
0.4	55	59	1.4	2.9		6.75;
%12月份
0.2	55	62	1.55	3		8.6625;
0.3	57	61.5	1.5	3.1		7.875;
0.75	54.5	63	1.85	3.2		8.4375;
0.5	58.5	61.5	1.55	2.9		7.3125;
0.4	57.5	61.5	1.35	3.1		8.55;
0.5	56.5	62	1.25	1.9		7.7625;
%2003-1
0.1	59	65.5	1.5	2.8		8.4375;
0.25	57	66	1.6	2.6		7.9875;
0.2	55.5	64.5	1.9	3.4		7.7625
];
%  
data01=cov(data0);
data011=cov(data0);

[pc,latent,explained]=pcacov(data0);
[pc1,latent1,explained1]=pcacov(data01);
[pc11,latent11,explained11]=pcacov(data011);