模糊聚类分析的应用

1、数学建模论文题目：模糊聚类分析在数学考研真题中的应用 摘 要本文采用模糊聚类分析方法和GM(1,1)灰色预测模型，利用软件matlab求解，预测出出卷者在未来怎样出题以及对考研者的复习指导。 关键词 ：模糊聚类分析 相关系数法 平方法 matlab 时间序列 一、问题的重述在数学建模中，如何用模糊数学中的“模糊聚类分析的方法解决近10年数学考研试题”这一个很模糊的问题？二、模型假设假设本小组从网上下载的考研真题具有真实性。假设从题目中提取的数据是合理的。假设本小组所用的算法在电脑中执行的结果是正确的。三、变量说明函数 - x1极限 - x2连续 - x3一元函数微积分学 - x4向量代数与空间

2、解析几何 - x5多元函数的微积分学 - x6无穷级数 - x7常微分方程 - x8行列式 - x9矩阵 - x10向量 - x11线性方程组 - x12矩阵的特征值和特征向量 - x13二次型 - x14随机事件和概率 - x15随机变量及其概率分布 - x16二维随机变量及其概率分布 - x17随机变量的数字特征 - x18大数定律和中心极限定理 - x19数理统计的基本概念 - x20参数估计 - x214、 模型的准备 首先，本小组对2004-2013年的数学考研试题中的每一道题目进行知识点的标记，然后对所有标记的题目通过知识点进行统计，如下表：表1年份2004200520062007

3、20082009201020112012201324323123424x34211321001x42123258342x51012010101x61432432343x72120101010x81111201111x91010101021x102222342323x111012101210x122111001210x130101210210x141000010111x150101020101x162011303212x172222022232x180110110100x190100000100x200000010000x211012101021其中，表中的数据又分为三

4、类： 高等数学（x1至x8） 线性代数(x9至x14)概率论与数理统计(x15至x21)五、模型的建立与求解 模型一 通过上面模型的准备，下面开始对数进行相应的处理。对上面每一类的数据进行相应的处理得到模糊相似矩阵,下面以高等数学中的知识点为例。步骤如下: (1)提取表1中的x1到x8中的数据，利用相关系数法，构造模糊相似关系矩阵,即。其中利用软件matlab编写程序得到高等数学知识点的模糊相似矩阵：（2） 利用平方法，得到模糊相似矩阵的模糊等价矩阵。利用matlab计算(代码见附录)得到，且。取(见附录)，得到：通过上面的矩阵对上面的x1至x8进行分类得到如下的结果:(x1),（x2,x4,

5、x6）,(x3,),(x5,x7),(x8)下面对所分类的结果进行检验:利用软件matlab的命令把（x2,x4,x6），(x5,x7)的图形分别画出如下： 图1通过图1确定：x5和x7的走向趋势具有相反性，x2与x4具有一致性，均与x6具有相反性。同上面的过程，我们再对线性代数和概率论与数理统计就简单了。线性代数的步骤如下:数据来源于(x9至x14) 模糊相似矩阵模糊等价矩阵取(见附录)，得到:通过上面的矩阵对上面的x9至x14进行分类得到如下的结果:(x9,x14),（x11,x13）,(x10),(x12)下面对所分类的结果进行检验:利用软件matlab的命令把（x9,x14），(x11

6、,x13)的图形分别画出如下： 图2从上面的图2中得出结论: x9与x14具有相反性，x11与x13具有相反性。概率论与数理统计的步骤如下：数据来源于（x15至x21）模糊相似矩阵:模糊等价矩阵：取(见附录)，得到:通过上面的矩阵对上面的x15至x21进行分类得到如下的结果:(x15,x21),（x16),(x17）,(x18),(x19),(x20)下面对所分类的结果进行检验: 图3通过图3得出:x15与x21具有相反性。模型二下面开始通过灰色预测模型（GM(1,1)）来对2014年各知识点出现的频率做预测。记x(0)为原始数列x(0)=(x(0)(k)xk=1,2,，n)=(x(0），x(

7、0），x(0)(n)记x为生成数列x=(x（k)xk=1,2，n)=(x，x，x（n)如果x(0) 与x之间满足下列关系，即称为一次累加生成。5建模步骤a、建模机理b、 把原始数据加工成生成数；c、 对残差（模型计算值与实际值之差）修订后，建立差分微分方程模型；d、 基于关联度收敛的分析；e、 gm模型所得数据须经过逆生成还原后才能用。f、采用“五步建模（系统定性分析、因素分析、初步量化、动态量化、优化）”法，建立一种差分微分方程模型gm(1,1）预测模型。GM（1，1）模型令 x(0)=(x，x，x(n)作一次累加生成， kx(k)= x(m) 消除数据的随机性和波动性m=1有 x=(x，x

8、，x(n)=(x，x+x，x(n-1)+x(n)x可建立白化方程：dx/dt+ax=u 即gm(1,1).该方程的解为： x(k+1)=(x-u/a)e+u/a其中：称为发展灰数；称为内生控制灰数利用matlab软件预测2014年21个变量的值，得到的结果如下表（代码见附录）:各知识点出现的次数预测x1x2x3x4x5x6x7x8x9x10x11x12x13x14x15x16x17x18x19x20x212014331413011301130220001 六、模型结果的分析通过上面模型一和模型二的建立与求解过程，可以得出: x5和x7的走向趋势具有相反性，x2与x4具有一致性，均与x6具有相反

9、性。 x9与x14具有相反性，x11与x13具有相反性。 x15与x21具有相反性。从而可以对考研者在未来一年的考研辅导过程中，进行以下指导: 对x1、x2、x4、x6、x10、x14、x16、x17重点复习，但是在复习x2、x4、x6时，先复习x2和x4,因为x6呈下降趋势，x2、x4与x6具有相反性。在复习x5和x7时重点复习x5,x7可以不复习，因为在预测中x7=0,x5=1;表示复习x5,符合x5和x7具有相反性。在复习x9和x14时重点复习x14，虽然可能会考x9,但可正常复习一下x9即可。在复习x11和x13以及x15和x21的时候，只需复习x13,x11可以简单盖过，符合x11与

10、x13具有相反性;只需复习x21,x15可以简单盖过，符合x15与x21具有相反性。其他的知识点，根据预测的结果得出以下结论：x3 x8 x10均需复习，x18 x19 x20可以简单的复习一下。对于出卷者而言： 在高数方面，重点出x1 x2 x4 x6四个知识点，x7不出，其他的知识点正常出题；在线性代数方面，重点出x10 x14,x11不出，其他的知识点正常出题；在概率论与数理统计方面，重点出x16 x17,再出x21,其他的知识点可以不出。 七、模型的优缺点 1、模糊聚类分析优点：聚类分析模型的优点就是直观，结论形式简明。 缺点：在样本量较大时，要获得聚类结论有一定困难。由于相似系数是根

11、据被试的反映来建立反映被试间内在联系的指标，而实践中有时尽管从被试反映所得出的数据中发现他们之间有紧密的关系，但事物之间却无任何内在联系，此时，如果根据距离或相似系数得出聚类分析的结果，显然是不适当的，但是，聚类分析模型本身却无法识别这类错误。 2、灰色预测模型优点：要求负荷数据少、不考虑分布规律、不考虑变化趋势、运算方便、短期预测精度高、易于检验。缺点：当数据离散程度大，即数据灰度大，预测精度越差。为了解决这一问题，一般提出对历史数据的平滑处理、模型参数修正等方法。 八、参考文献美 MATLAB实用教程(第二版) Holly Moore 著 高会生 刘童娜 议十、附录 高等数学的知识点的代码

12、：1.求模糊相似矩阵的MATLAB程序a=5 5421121; 53210411; 33121321; 43142301; 43340412; 32231300; 33180201; 43031401; 42040411; 34121311;mu=mean(a);sigma=std(a);for i=1:8 for j=1:8r(i,j)=exp(-(mu(j)-mu(i)2/(sigma(i)+sigma(j)2);end endrsave data1 r a2.矩阵合成的MATLAB函数function rtha=hecheng(r);n=length(r);for i=1:n for j

13、=1:nrtha(i,j)=max(min(r(i,:);r(:,j);end end3.求模糊等价矩阵和聚类的程序load data1;r1=hecheng(r)r2=hecheng(r1)r3=hecheng(r2)bh=zeros(8);bh(find(r20.995)=1 线性代数的知识点代码：1.求模糊相似矩阵的MATLAB程序b=1 21211; 021110; 021110; 021110; 030020; 040001; 121100; 032211; 221111; 130001;mu=mean(b);sigma=std(b);for i=1:6 for j=1:6r(i,j

14、)=exp(-(mu(j)-mu(i)2/(sigma(i)+sigma(j)2);end endrsave data2 r b2.矩阵合成的MATLAB函数function rtha=hecheng(r);n=length(r);for i=1:n for j=1:nrtha(i,j)=max(min(r(i,:);r(:,j);end end3.求模糊等价矩阵和聚类的程序load data1;r1=hecheng(r)r2=hecheng(r1)r3=hecheng(r2)bh=zeros(n);bh(find(r20.996)=1 概率论与数理统计知识点代码1.求模糊相似矩阵的MATLA

15、B程序c= 0220001; 1021101; 0121001; 1120001; 0301001; 2021010; 0320000; 1221100; 1130000; 0220001;mu=mean(c);sigma=std(c);for i=1:7 for j=1:7r(i,j)=exp(-(mu(j)-mu(i)2/(sigma(i)+sigma(j)2);end endrsave data3 r c2.矩阵合成的MATLAB函数function rtha=hecheng(r);n=length(r);for i=1:n for j=1:nrtha(i,j)=max(min(r(i,

16、:);r(:,j);end end3.求模糊等价矩阵和聚类的程序load data1;r1=hecheng(r)r2=hecheng(r1)r3=hecheng(r2)bh=zeros(n);bh(find(r20.999)=1高等数学的代码和图形：x=2004:1:2013;x2=4323123424;x4=2123258342;x6=1432432343;x5=1012010101;x7=2120101010;subplot(1,2,1)plot(x,x2,r,x,x4,b,x,x6,k)hold onsubplot(1,2,2)plot(x,x5,b,x,x7,k)线性代数的代码和图形：

17、x=2004:1:2013;x9=1010101021;x14=0101210210;x11=1012101210;x13=1210201210;subplot(1,2,1)plot(x,x9,k,x,x14,b)hold onsubplot(1,2,2)plot(x,x11,b,x,x13,k)概率论与数理统计的代码和图形：x=2004:1:2013;x15=0101020101;x21=1012101021;plot(x,x15,b-v,x,x21,o-k)调用函数代码：function y=gm11(x,n)%x为行向量数据%做一次累加x1=zeros(size(x);for i=1:s

18、ize(x1,2) x1(i)=sum(x(1:i);end%x1的均值数列z1=zeros(size(x);for i=1:size(x1,2)-1 z1(i+1)=0.5*x1(i+1)+0.5*x1(i);endYn=x(2:size(x,2);B=-z1(2:size(z1,2);B(:,2)=1;u=inv(B*B)*B*Yn;a=u(1);b=u(2);%预测x2=zeros(1,n);x2(1)=x(1);for i=1:n-1 x2(1+i)=(x(1)-b/a)*exp(-a*i)+b/a;endx2=0 x2;y=diff(x2);预测结果代码：x1=5 534433443;gm11(x1,11)x2=4 323123424;gm11(x2,11)x3=4 211321001;gm11(x3,11)x4=2 123258342;gm11(x4,11)x5=1 012010101;gm11(x5,11)x6=1 432432343;gm11(x6,11)x7=2 120101010;gm11(x7,11)x8=1 111201111;gm11(x8,11)x9=1 010101021;gm11(x9,11)x10=2 222342323;gm11(x10,11)x11=1 012101