建立M文件夹： %发病人群信息
A=xlsread('c题数据1.xls','sheet1','A2:F58926');X925行-6列 a=size(A);
B=zeros(4,12); for i=1:a(1)
if A(i,4)==2007 if A(i,5)==1
B(1,1)=B(1,1)+1; elseif A(i,5)==2 B(1,2)=B(1,2)+1; elseif A(i,5)==3 B(1,3)=B(1,3)+1; elseif A(i,5)==4 B(1,4)=B(1,4)+1; elseif A(i,5)==5 B(1,5)=B(1,5)+1; elseif A(i,5)==6 B(1,6)=B(1,6)+1; elseif A(i,5)==7 B(1,7)=B(1,7)+1; elseif A(i,5)==8 B(1,8)=B(1,8)+1; elseif A(i,5)==9 B(1,9)=B(1,9)+1; elseif A(i,5)==10 B(1,10)=B(1,10)+1; elseif A(i,5)==11 B(1,11)=B(1,11)+1; elseif A(i,5)==12 B(1,12)=B(1,12)+1;
end 07年各自月份病发者指标 elseif A(i,4)==2008 if A(i,5)==1 B(2,1)=B(2,1)+1; elseif A(i,5)==2 B(2,2)=B(2,2)+1; elseif A(i,5)==3 B(2,3)=B(2,3)+1; elseif A(i,5)==4 B(2,4)=B(2,4)+1; elseif A(i,5)==5 B(2,5)=B(2,5)+1; elseif A(i,5)==6
26
B(2,6)=B(2,6)+1; elseif A(i,5)==7 B(2,7)=B(2,7)+1; elseif A(i,5)==8 B(2,8)=B(2,8)+1; elseif A(i,5)==9 B(2,9)=B(2,9)+1; elseif A(i,5)==10 B(2,10)=B(2,10)+1; elseif A(i,5)==11 B(2,11)=B(2,11)+1; elseif A(i,5)==12 B(2,12)=B(2,12)+1;
end 08年各自月份病发者指标 elseif A(i,4)==2009 if A(i,5)==1 B(3,1)=B(3,1)+1; elseif A(i,5)==2 B(3,2)=B(3,2)+1; elseif A(i,5)==3 B(3,3)=B(3,3)+1; elseif A(i,5)==4 B(3,4)=B(3,4)+1; elseif A(i,5)==5 B(3,5)=B(3,5)+1; elseif A(i,5)==6 B(3,6)=B(3,6)+1; elseif A(i,5)==7 B(3,7)=B(3,7)+1; elseif A(i,5)==8 B(3,8)=B(3,8)+1; elseif A(i,5)==9 B(3,9)=B(3,9)+1; elseif A(i,5)==10 B(3,10)=B(3,10)+1; elseif A(i,5)==11 B(3,11)=B(3,11)+1; elseif A(i,5)==12 B(3,12)=B(3,12)+1;
end 09年各自月份病发者指标 else A(i,4)==2010 if A(i,5)==1 B(4,1)=B(4,1)+1; elseif A(i,5)==2
27
B(4,2)=B(4,2)+1; elseif A(i,5)==3 B(4,3)=B(4,3)+1; elseif A(i,5)==4 B(4,4)=B(4,4)+1; elseif A(i,5)==5 B(4,5)=B(4,5)+1; elseif A(i,5)==6 B(4,6)=B(4,6)+1; elseif A(i,5)==7 B(4,7)=B(4,7)+1; elseif A(i,5)==8 B(4,8)=B(4,8)+1; elseif A(i,5)==9 B(4,9)=B(4,9)+1; elseif A(i,5)==10 B(4,10)=B(4,10)+1; elseif A(i,5)==11 B(4,11)=B(4,11)+1; elseif A(i,5)==12 B(4,12)=B(4,12)+1;
end 10年各自月份病发者指标 end end
C=mean(B)/mean(mean(B)),%季节指数 附录7
主成分分析：
设置步骤如下： （1）打开因子分析.
28
（2）将要分析的变量移入.
（3）选中“初始解”、“相关系数”和“KMO检验和Bartlett球形检验”.
（4）选中“主成分分析方法”、“相关系数矩阵”、“未经过旋转的因子解（因子载荷矩阵）”、“碎石图”、“变量提取数目”.
（5）（如果需要旋转的话）选中“方差最大化旋转方式”、“旋转解”、“因子载荷图”.
点击“OK”就可以输出计算结果.
附录8 程序： x=1:48;
29
y=[1.87 5.01 7.95 11.29 18.49 22.12 26.29 26.61 21.66 16.05 8.93 5.2 0.83 -0.59 7.01 12.1 17.04 21.05 26.81 24.94 22.12 16.82 8.65 2.71 -0.37 5.71 6.17 11.66 16.72 22.5 25.1 25.26 21.72 16.5 7.97 2.54 1.03 3.81 5.44 9.02 17.13 20.76 25.84 26.97 22.67 14.67 8.69 3.08 ];
f=inline('a(1)*cos(a(2)*x+a(3))+a(4)','a','x'); a=lsqcurvefit(f,[100;pi/6;-pi/4;30],x,y) x1=1:0.01:48;
y1=a(1)*cos(a(2)*x1+a(3))+a(4); plot(x,y,'o',x1,y1,'r')
y2=a(1)*cos(a(2)*x+a(3))+a(4); y3=(y2-y)./y;%输出相对误差 y4=abs(y3);%输出绝对相对误差
y5=mean(y4),%输出平均绝对相对误差 y6=max(y4);%输出最大绝对相对误差 t=49:1:60;
y7=a(1)*sin(a(2)*t+a(3))+a(4)
30