sarstestdata1.m

该程序用遗传算法对人工神经网络进行优化

%stn_id	year	month	day	tmax	tmin	hu	cm	prec	空气污染指数
%45007表示香港由于tmax tmin缺失故不纳入测试
%注意：第一列是站台号其中：
%                          57687表示长沙
%                          59287表示广州
dmat=[57687	2003	3	1	102	45	89	94	0	56
57687	2003	3	2	209	82	83	73	5	97
57687	2003	3	3	109	89	92	100	173	76
57687	2003	3	4	46	30	91	100	129	53
57687	2003	3	5	31	7	91	100	8	44
57687	2003	3	6	50	4	81	75	0	50
57687	2003	3	7	116	-13	75	29	0	88
57687	2003	3	8	149	30	69	38	0	95
57687	2003	3	9	144	72	69	75	0	97
57687	2003	3	10	171	34	64	0	31	109
57687	2003	3	11	123	87	79	98	269	90
57687	2003	3	12	95	83	94	100	140	65
57687	2003	3	13	99	84	96	100	104	69
57687	2003	3	14	86	73	96	100	38	67
57687	2003	3	15	76	61	94	100	13	56
57687	2003	3	16	73	59	93	100	70	48
57687	2003	3	17	61	48	91	100	190	48
57687	2003	3	18	56	40	93	100	63	54
57687	2003	3	19	85	47	88	100	36	64
57687	2003	3	20	120	44	87	100	0	94
57687	2003	3	21	175	42	78	100	16	117
57687	2003	3	22	108	70	81	78	18	84
57687	2003	3	23	159	53	78	48	0	80
57687	2003	3	24	219	66	72	48	0	149
57687	2003	3	25	261	119	63	0	0	147
57687	2003	3	26	281	174	63	73	20	89
57687	2003	3	27	191	148	77	95	0	93
57687	2003	3	28	211	146	77	66	5	99
57687	2003	3	29	247	139	76	25	5	94
57687	2003	3	30	313	194	66	45	5	89
57687	2003	3	31	311	224	65	63	0	51
57687	2003	4	1	245	161	85	75	170	61
57687	2003	4	2	118	91	89	100	3	30
57687	2003	4	3	107	92	91	100	135	63
57687	2003	4	4	146	98	85	98	0	46
57687	2003	4	5	187	70	74	35	0	79
57687	2003	4	6	166	122	84	98	32	95
57687	2003	4	7	217	128	80	69	0	84
57687	2003	4	8	247	110	74	58	0	104
57687	2003	4	9	211	159	74	98	138	94
57687	2003	4	10	171	154	96	100	331	61
57687	2003	4	11	162	144	93	100	9	46
57687	2003	4	12	150	135	94	100	278	66
57687	2003	4	13	156	127	86	98	23	64
57687	2003	4	14	222	96	62	34	23	117
57687	2003	4	15	267	127	57	0	23	145
57687	2003	4	16	304	177	60	50	23	106
57687	2003	4	17	316	234	63	63	23	61
57687	2003	4	18	295	241	74	69	23	64
57687	2003	4	19	147	121	84	94	23	40
57687	2003	4	20	137	116	89	100	23	36
57687	2003	4	21	171	122	91	100	1	47
57687	2003	4	22	192	150	91	100	46	80
57687	2003	4	23	154	128	87	100	0	60
57687	2003	4	24	142	127	91	98	252	64
57687	2003	4	25	229	106	77	98	0	73
57687	2003	4	26	277	155	67	6	0	122
57687	2003	4	27	283	176	79	54	61	125
57687	2003	4	28	251	171	89	95	58	96
57687	2003	4	29	213	171	91	100	15	112
57687	2003	4	30	271	156	70	6	0	77
59287	2003	2	1	191	145	80	51	0	102
59287	2003	2	2	146	108	66	98	0	31
59287	2003	2	3	130	104	74	100	0	44
59287	2003	2	4	162	96	61	79	0	31
59287	2003	2	5	186	98	43	31	0	52
59287	2003	2	6	209	78	53	58	0	69
59287	2003	2	7	156	120	87	98	65	93
59287	2003	2	8	229	132	87	83	0	108
59287	2003	2	9	237	149	83	95	0	132
59287	2003	2	10	260	175	78	76	0	80
59287	2003	2	11	263	155	79	76	0	100
59287	2003	2	12	145	105	73	100	1	82
59287	2003	2	13	166	90	74	98	0	49
59287	2003	2	14	202	134	77	98	0	77
59287	2003	2	15	218	168	84	100	2	145
59287	2003	2	16	265	191	85	91	2	99
59287	2003	2	17	185	156	79	98	1	78
59287	2003	2	18	237	171	81	94	0	105
59287	2003	2	19	245	174	86	86	1	150
59287	2003	2	20	203	155	70	95	0	75
59287	2003	2	21	218	164	70	100	1	100
59287	2003	2	22	257	179	80	84	0	108
59287	2003	2	23	263	176	83	69	0	86
59287	2003	2	24	281	178	79	74	0	103
59287	2003	2	25	252	192	74	93	0	100
59287	2003	2	26	258	196	71	93	0	89
59287	2003	2	27	286	173	71	64	0	99
59287	2003	2	28	257	210	73	74	0	82
59287	2003	3	1	279	195	72	53	0	80
59287	2003	3	2	259	201	78	58	0	86
59287	2003	3	3	242	204	81	95	0	81
59287	2003	3	4	231	187	89	98	132	64
59287	2003	3	5	211	173	89	100	6	66
59287	2003	3	6	140	102	77	98	0	67
59287	2003	3	7	165	66	58	54	0	40
59287	2003	3	8	146	95	52	93	24	53
59287	2003	3	9	160	104	76	95	41	77
59287	2003	3	10	187	109	78	93	0	78
59287	2003	3	11	216	135	62	93	0	97
59287	2003	3	12	231	145	69	70	3	97
59287	2003	3	13	265	158	70	64	0	115
59287	2003	3	14	229	163	71	80	0	88
59287	2003	3	15	239	181	78	93	0	85
59287	2003	3	16	267	197	81	74	0	91
59287	2003	3	17	276	209	80	78	8	73
59287	2003	3	18	226	154	82	98	1	61
59287	2003	3	19	147	112	82	100	3	50
59287	2003	3	20	121	100	92	100	83	39
59287	2003	3	21	200	111	78	91	1	44
59287	2003	3	22	220	136	70	93	166	75
59287	2003	3	23	170	150	96	100	125	82
59287	2003	3	24	221	156	89	98	0	75
59287	2003	3	25	268	165	76	69	0	99
59287	2003	3	26	276	162	72	79	0	93
59287	2003	3	27	275	189	78	76	0	88
59287	2003	3	28	251	198	88	89	41	127
59287	2003	3	29	221	189	87	98	41	112
59287	2003	3	30	233	194	90	91	3	113
59287	2003	3	31	261	209	87	98	5	81
59287	2003	4	1	278	226	89	95	0	66
59287	2003	4	2	296	235	87	89	3	67
59287	2003	4	3	292	233	87	93	13	72
59287	2003	4	4	235	190	78	95	0	86
59287	2003	4	5	230	167	61	95	0	36
59287	2003	4	6	200	181	78	98	17	71
59287	2003	4	7	197	171	92	100	71	103
59287	2003	4	8	266	186	87	98	131	66
59287	2003	4	9	216	185	91	100	1	66
59287	2003	4	10	200	175	92	100	26	91
59287	2003	4	11	254	190	92	88	0	109
59287	2003	4	12	268	231	91	90	0	80
59287	2003	4	13	288	238	86	98	193	63
59287	2003	4	14	245	180	71	89	15	53
59287	2003	4	15	236	173	68	74	231	68
59287	2003	4	16	238	196	87	100	1	118
59287	2003	4	17	295	210	78	73	1	123
59287	2003	4	18	305	224	75	90	0	70
59287	2003	4	19	287	231	76	70	0	62
59287	2003	4	20	312	221	74	49	0	69
59287	2003	4	21	323	217	70	64	0	88
59287	2003	4	22	307	230	70	78	0	66
59287	2003	4	23	315	230	73	65	0	88
59287	2003	4	24	329	236	77	76	136	91
59287	2003	4	25	311	213	66	44	0	76
59287	2003	4	26	326	210	63	76	0	71
59287	2003	4	27	310	232	70	83	0	103
59287	2003	4	28	317	225	69	70	0	94
59287	2003	4	29	316	229	68	66	0	71
59287	2003	4	30	292	229	53	83	0	75
59287	2003	5	1	292	192	56	38	0	67
59287	2003	5	2	283	210	69	83	2	84
59287	2003	5	3	271	232	83	98	116	80
59287	2003	5	4	271	223	89	93	116	70
59287	2003	5	5	299	218	84	93	1	74
59287	2003	5	6	320	254	80	90	1	62
59287	2003	5	7	316	266	88	73	278	63
59287	2003	5	8	304	226	85	73	3	60
59287	2003	5	9	291	248	81	98	0	58
59287	2003	5	10	288	245	80	76	1	110
59287	2003	5	11	321	230	68	50	0	106
59287	2003	5	12	312	238	74	83	0	73
59287	2003	5	13	324	244	74	90	0	68
59287	2003	5	14	316	259	82	89	763	68
59287	2003	5	15	322	243	83	93	1	73
59287	2003	5	16	325	263	79	86	0	61
59287	2003	5	17	312	270	81	95	37	61
59287	2003	5	18	262	232	81	98	0	53
59287	2003	5	19	307	231	69	95	0	43
59287	2003	5	20	323	250	68	79	0	77
59287	2003	5	21	332	258	60	40	0	85
59287	2003	5	22	325	241	68	59	0	85
59287	2003	5	23	326	262	78	79	100	72
59287	2003	5	24	334	265	83	90	378	60
59287	2003	5	25	317	242	77	90	0	59
59287	2003	5	26	311	252	76	95	0	70
59287	2003	5	27	315	239	66	79	0	90
59287	2003	5	28	310	253	76	93	0	86
59287	2003	5	29	330	251	76	76	0	99
59287	2003	5	30	339	246	60	54	0	89
59287	2003	5	31	347	241	56	48	0	96
59287	2003	6	1	346	240	56	31	0	99
];
[linetmp columntmp]=size(dmat);
for i=1:linetmp
    if(dmat(i,1)==53463)
        dmat(i,columntmp-1)=dmat(i,columntmp)*6.5;
    elseif(dmat(i,1)==54511)
        dmat(i,columntmp-1)=dmat(i,columntmp)*5.6;
    elseif(dmat(i,1)==57687)
        dmat(i,columntmp-1)=dmat(i,columntmp-1)*4.7;
    else
        dmat(i,columntmp-1)=dmat(i,columntmp-1)*5.0;
    end
end
testdat=dmat(:,5:columntmp);
testdat=unitay(testdat,0);
trainweight
[length column]=size(testdat);
for i=1:length
    output1(i,:)=f(testdat(i,:)*w1'+b1');
    output2(i)=f(output1(i,:)*w2'+b2');
end

%运行后的最终结果为(注意：运行时请先将结果删除)
   0.24414354520175
   0.47391308006893
   0.38043963161321
   0.22823380498376
   0.15782069031407
   0.20251812844892
   0.42201219062549
   0.46613471637063
   0.47732864144177
   0.50918459333065
   0.47751108826034
   0.30842945701358
   0.32038394749228
   0.29259952823254
   0.22914608109694
   0.19818294167284
   0.21781007493287
   0.22162408307438
   0.29333552300873
   0.44573560663716
   0.53676546454100
   0.40523093094783
   0.39509839550290
   0.57600835470532
   0.56646780874209
   0.48274875530578
   0.46570760589437
   0.48202061433172
   0.46683338754221
   0.47952100670596
   0.33802627602948
   0.36072507758346
   0.12963255067057
   0.30793503451010
   0.21848927154546
   0.40043724776798
   0.46427004778462
   0.43073161727194
   0.50452593724764
   0.49334698424308
   0.34975932404023
   0.21386490910266
   0.36074896606104
   0.32185217446149
   0.53308572265301
   0.56589387946092
   0.52267075739385
   0.39284519709565
   0.38868804533408
   0.19408055476668
   0.16425563179023
   0.22272193381683
   0.40764624972135
   0.29460900815007
   0.34770653417603
   0.40216257237844
   0.53560775828836
   0.55012701648717
   0.48495932948638
   0.51273329492145
   0.41261817244814
   0.47944940467293
   0.18064691687252
   0.22596045981897
   0.18739608125673
   0.31601400673411
   0.40363432358111
   0.45636619066354
   0.50515752215195
   0.55480209956176
   0.43293355186292
   0.49593829014502
   0.42310637703714
   0.26409450753219
   0.41101865904999
   0.56771250080771
   0.49239957251255
   0.40742852477269
   0.50484497401724
   0.57380212597850
   0.41176461120553
   0.49767932545172
   0.51409898240047
   0.44827053251465
   0.50626225388629
   0.49944870273598
   0.47345267322609
   0.50046385246752
   0.44476250387752
   0.43902108011520
   0.45012934756301
   0.43409035337199
   0.36848008477711
   0.34726224324645
   0.34489558370106
   0.22875926211052
   0.32459054074685
   0.40663092109685
   0.40799241780832
   0.49593224217345
   0.48796643851328
   0.53081402708698
   0.46071659471635
   0.44966285696992
   0.46934990559073
   0.41132424067002
   0.33681847749673
   0.24827973694849
   0.17643002198613
   0.23932043287397
   0.44250801353269
   0.41534920145146
   0.39041701893278
   0.49484221541769
   0.48569059999290
   0.46496572049015
   0.54893173735025
   0.51898227376811
   0.51631788780766
   0.43457696460132
   0.37396894169380
   0.38746237164090
   0.41031215618370
   0.45261109695146
   0.24436594390342
   0.38676503915839
   0.49291291813161
   0.39297396493741
   0.34578971195306
   0.44989629859016
   0.50945520492362
   0.42678074630402
   0.40393690947541
   0.32314493346744
   0.42599040644699
   0.52991031068166
   0.54536574350303
   0.41712870560159
   0.37011315878434
   0.40461569214927
   0.48022884085401
   0.40406156752435
   0.47676672635623
   0.50470766622899
   0.43729592135283
   0.43490018941154
   0.51636463628417
   0.49615090286550
   0.42523965181832
   0.44892931296389
   0.40427882733123
   0.46397540818995
   0.45603525812794
   0.40877726328829
   0.42125034755900
   0.38471360038743
   0.42077069152424
   0.35473671865590
   0.35605192431410
   0.52142888120499
   0.51916123926926
   0.43075037143946
   0.41745288036501
   0.51472938984774
   0.42738978698204
   0.38282530351655
   0.38547003678682
   0.31845741738041
   0.30180375701686
   0.45241889907266
   0.47413232384987
   0.47253644281947
   0.44290574325639
   0.44032198553926
   0.37265345820347
   0.42098923969386
   0.48894392503180
   0.47428718666914
   0.50828544248732
   0.48909701639417
   0.50652871360874
   0.50922419468583