matlab数学实验

1、matlab 数学实验管理 数学实验 实验报告班级姓名实验 1：MATLAB的数值运算【实验目的】（ 1）掌握 MATLAB 变量的使用（ 2）掌握 MATLAB 数组的创建，（ 3）掌握 MA TLAB 数组和矩阵的运算。（ 4）熟悉 MATLAB 多项式的运用【实验原理 】矩阵运算和数组运算在MA TLAB中属于两种不同类型的运算，数组的运算是从数组元素出发， 针对每个元素进行运算， 矩阵的运算是从矩阵的整体出发， 依照线性代数的运算规则进行。【实验步骤 】（ 1）使用冒号生成法和定数线性采样法生成一维数组。（ 2）使用 MA TLAB 提供的库函数 reshape，将一维数组转换为二维和

2、三维数组。（ 3）使用逐个元素输入法生成给定变量，并对变量进行指定的算术运算、关系运算、逻辑运算。（4）使用 MA TLAB绘制指定函数的曲线图，将所有输入的指令保存为M 文件。【实验内容 】（ 1）在 0,2*pi 上产生 50 个等距采样数据的一维数组，用两种不同的指令实现。0:(2*pi-0)/(50-1):2*pi或linspace(0,2*pi,50)（ 2）将一维数组A=1:18 ，转换为2×9 数组和 2× 3× 3 数组。reshape(A,2,9)ans =Columns 1 through 7135246789111012131415 1716

3、 18reshape(A,2,3,3)ans(:,:,1) =135246ans(:,:,2) =791181012ans(:,:,3) =131517141618matlab 数学实验（ 3）A=0 2 3 4 ;1 3 5 0,B=1 0 5 3;1 5 0 5，计算数组 A、 B 乘积，计算 A&B,A|B,A,A=B,A>B 。A.*Bans=00151211500A&Bans =00111100A|Bans =11111111Aans =10000001A=Bans =00001000A>=Bans =01011010tt（ 4）绘制 y= 0.5 e3-

4、t*t*sin(t),t=0,pi并标注峰值和峰值时间,添加标题 y= 0.5 e3 -t*t*sint ，将所有输入的指令保存为M 文件。a=0.5b=1/3t=0:0.001:piy=a*exp(b*t)-t.*t.*sin(t)y_max,t_max=max(y)t_text='t=',num2str(t(t_max)y_text='y=',num2str(y_max)max_text=char('maximum',t_text,y_text)tit='y=a*exp(',num2str(b),'t)-t*t*sin

5、(t)'hold onplot(t,y,'y.')plot(t(t_max),y_max,'r')text(t(t_max)+0.3,y_max+0.1,max_text)title(tit),xlabel('t'),ylabel('y'),hold offmatlab 数学实验【实验心得与总结】通过这次试验让我了解常用简单函数的功能，学会利用函数解决一些；数值计算和符号计算的实际问题 ; 利用 Matlab 的 help 命令查询一些函数的功能。利用MA TLAB可以让繁琐的计算问题变得更加简单化，如矩阵运算等。实验 2

6、：MATLAB绘图【实验步目的 】利用 MTALAB画墨西哥帽子，及参数方程的图像【实验原理 】(1)二维绘图命令：plot(x,y) 函数(2)三维绘图命令中三维曲线：plot3(x,y,z),(3)利用 mesh 函数画三维的网格表面的。【实验内容 】（含参考程序、实验结果及结果分析等）x2 cos(t)画出函数图形yt 30 t 10 。zt方程：x 2cos(t)yt 30t10zt【参考程序】matlab 数学实验>> t=0:0.1:4*pi;>> plot3(2*cos(t),t.3,t)【实验结果 】22sin xy7.5,7.5z f (x, y)2,

7、 x y2画出曲面x的图像。方程：sin x2y27.5,7.5, y 7.5,7.5z f (x, y), xx2y2【参考程序】x = -7.5:0.5:7.5;y = x;xx, yy = meshgrid(x, y);R = sqrt(xx.2 + yy.2) + eps;z = sin(R)./R;surf(xx, yy, z)【实验结果 】matlab 数学实验10.50-0.510510050-5-5-10-10【实验心得与总结】Matlab 的常见错误： Inner matrix dimensions must agree1、因为在 Matlab 的输入变量是矩阵，参与运算的矩

8、阵维数必须对应，矩阵相应元素的运算必须全部加dot（点），例 2中方程如果这样输入：x=2*(cos(t)+t*sin(t)，就会出现该错误.2、 mesh函数是用来画三维的网格表面的。三维空间中的一个点是用(x,y,z) 来表示的，mesh就是把这些点之间用网格连接起来。实验 3：MATLAB微积分问题的计算【实验目的】 利用 MTALAB求解二重积分、勒展开式及级数求和。【实验原理 】1 利用 int(int(f,x,a,b),y,c,d)db函数求二重积分计算累次积分f ( x, y)dxdyca2 利 用 泰 勒 函 数 taylor（ f,n,x,a) 来 求f(x,y) 的 n-1

9、阶泰勒展开式f ( x)n 1 f (k ) (a) ?( x a)k ；k 0k!n23.利用函数 symsum(f,k,n1,n2) 来求级数的和函数f (k)kn1【实验内容 】（含参考程序、实验结果及结果分析等）21x1xydydx求02 x。【参考程序 】>> syms x y>> z=x*y;>> f=int(int(z,y,2*x,x2+1),x,0,1)【 实验结果 】f =1/12matlab 数学实验将 f （x） =ln x 展开为幂为（ x-2) 的 5 阶泰勒展开式。【 参考程序 】>> syms x n;>>

10、; f=(-1)n*x(n+1)/(n+1);>> symsum(f,n,1,inf)【 实验结果 】ans =log(1+x)-x级数求和( 1) n xn 1, x ( 1,1) 。n 1n 1【参考程序 】>> syms x n;>> f=(-1)n*x(n+1)/(n+1);>> symsum(f,n,1,inf)【 实验结果 】ans = log(1+x)-x【实验心得与总结】1、在实验过程中，要是一句程序结束后加了分号，则说明，不要求执行程序时输出执行结果；2、在 matlab 中是区别大小写的，如果 N 写成 n 会出现 Undef

11、ined function or variable 'n'.Undefined function or variable 'n'. 的错误提示 .实验 4： MATLAB优化计算【实验目的 】掌握应用matlab 求解无约束最优化问题的方法【实验原理与方法 】1：标准形式：minn f (X )x R其中 f : RnR为n元函数2无约束优化问题的基本算法一最速下降法（共轭梯度法）算法步骤： 给定初始点X 0E n ，允许误差0, 令 k=0；计算fX k；检验是否满足收敛性的判别准则：fX k，若满足，则停止迭代，得点X *X k ，否则进行；令 S kfX

12、k，从 X k 出发，沿Sk 进行一维搜索，即求k 使得：min f X kSkf X kk Sk ；0令 X k1X kk Sk ， k=k+1 返回 .matlab 数学实验最速下降法是一种最基本的算法，它在最优化方法中占有重要地位. 最速下降法的优点是工作量小，存储变量较少, 初始点要求不高；缺点是收敛慢，最速下降法适用于寻优过程的前期迭代或作为间插步骤，当接近极值点时，宜选用别种收敛快的算法. 牛顿法算法步骤：(1)选定初始点 X 0E n ，给定允许误差0 ，令 k=0；(2)求f X k,2f X k1, 检验：若f X k, 则停止迭代 ,X *X k . 否则 , 转向 (3)

13、；(3)令Sk 2 fX k1f X k （牛顿方向）；(4)X k 1X kSk ,kk1, 转回 (2).如果 f 是对称正定矩阵A 的二次函数，则用牛顿法经过一次迭代就可达到最优点 , 如不是二次函数，则牛顿法不能一步达到极值点，但由于这种函数在极值点附近和二次函数很近似, 因此牛顿法的收敛速度还是很快的 .牛顿法的收敛速度虽然较快，但要求 Hessian 矩阵要可逆， 要计算二阶导数和逆矩阵， 就加大了计算机计算量和存储量.【实验内容 】1. 求 f = 2 e x sin x 在0<x<8中的最小值与最大值主程序为 wliti1.m:f='2*exp(-x).*s

14、in(x)'fplot(f,0,8);%作图语句xmin,ymin=fminbnd (f, 0,8)f1='-2*exp(-x).*sin(x)'xmax,ymax=fminbnd (f1, 0,8)运行结果：xmin = 3.9270xmax = 0.7854ymin = -0.0279ymax = 0.64482. 对边长为 3 米的正方形铁板，在四个角剪去相等的正方形以制成方形无盖水槽，问如何剪法使水槽的容积最大？先编写 M文件 fun0.m 如下 : function f=fun0(x)2) x解：设剪去的正方形的边长为x，则水槽的容积为：(3 2xf=-(3-

15、2*x).2*x;min y=- (3 2x2 )x ， 0<x<1.5主程序为建立无约束优化模型为：wliti2.m:x,fval=fminbnd('fun0',0,1.5);xmax=xfmax=-fval运算结果为 : xmax = 0.5000,fmax =2.0000.即剪掉的正方形的边长为0.5 米时水槽的容积最大 , 最大容积为2立方米.实验 5： MATLAB图论问题计算【实验目的 】了解用 Matlab 软件求解图论模型及层次分析模型的方法。【实验内容与原理】内容： 1.某城市要建立一个消防站， 为该市所属的七个区服务， 如图所示 问应设在那个区，

16、才能使它至最远区的路径最短。matlab 数学实验2.某矿区有七个矿点，如图所示已知各矿点每天的产矿量q(v j ) （标在图的各顶点上） 现要从这七个矿点选一个来建造矿厂问应选在哪个矿点， 才能使各矿点所产的矿运到选矿厂所在地的总运力（千吨公里）最小原理：利用层次分析法和图论方法模型的一般概念， 理解建立层次分析法和图论方法模型的一般方法，初步学会建立层次分析法和图论方法模型以解决实际问题。【操作方法与步骤】步骤： 1.(1)用 Floyd 算法求出距离矩阵D= ( dij )(2)计算在各点 vi 设立服设施的最大服务距离S(vi )S(vi ) max dij i1,2,1 jS(vk

17、)min S(vi )(3) 求出顶点 vk ，使1i建立M文件a=0 3 inf inf inf inf inf;3 0 2 inf 18 2.5 inf;inf 2 0 6 2 inf inf;inf inf 6 0 3 inf inf;inf 18 2 3 0 4 inf;inf 2.5 inf inf 4 0 1.5;inf inf inf inf inf 1.5 0;D,R=floyd(a)点击运行则 vk 就是要求的建立消防站的地点此点称为图的中心点2. （ 1）求距离阵 D= ( dij ) （ 2） 计算各顶点作为选矿厂的总运力m(vi )m(vi )q(v j ) dijj

18、1i 1,2,matlab 数学实验m(vk )min m(vi )，（ 3）求 vk 使1 i（ 4）建立 M 文件a=0 3 inf inf inf inf inf;3 0 2 inf inf 4 inf;inf 2 0 6 2 inf inf;inf inf 6 0 1 inf inf;inf inf 2 1 0 4 inf;inf 4 inf inf 4 0 1.5;inf inf inf inf inf 1.5 0;D,R=floyd(a)q=3,2,7,1,6,1,4;m=0;for i=1:7for j=1:7m=m+q(i)*D(i,j);endmm=0;end点击运行（ 5）

19、则 vk 就是选矿厂应设之矿点此点称为图G 的重心或中位点【实验结果与分析】实验结果与分析：1.S(v1)=10, S(v2)=7, S(v3)=6, S(v4)=8.5, S(v5)=7, S(v6)=7, S(v7)=8.5S(v3)=6, 故应将消防站设在v3 处。2.0358778.53025445.55203267.5D8530156.57421045.57465401.58.55.57.56.55.51.50由上述公式可得：m(v1)=38.5*3=115.5, m(v2)=23.5*2=47, m(v3)=23.5*7=164.5, m(v4)=28.5, m(v5)=23.5*

20、6=141, m(v6)=27.5, m(v7)=35*4=140再求其中的最小值，m(v6)=27.5 ，则 v6 就是选矿厂应设之矿点实验 6：MATLAB计算机模拟 计算matlab 数学实验【实验目的 】学会用数学软件 matlab和蒙特卡洛方法估计积分值， 并于其中应用概率论中的概率密度等知识点。【实验问题 】估计积分值，并对误差进行估计。【实验要求 】针对要估计的积分选择适当的概率分布设计蒙特卡洛估计算法；利用计算机产生所选分布的随机样本的估计积分值；通过计算平均误差对估计结果进行评价。【实验过程分析】（ x为运行结果平均值， y为样本方差）估计程序如下：clc;clear;m=1

21、0;n=10000;d=0;e=0;for i=1:md=0;a=rand(1,n);for j=1:nb=a(j)+2;c=b2;e(i)=d+c/n;d=e(i);endfprintf('e=%.8fn',e(i)endp=sum(e)/m;for j=1:m;s(j)=(e(j)-p)2;endq=sum(s);fprintf('x=%.8fny=%.8fn',p,q);结果为：e=6.34879520e=6.34068140e=6.35081124e=6.31353632e=6.35586630e=6.33058791e=6.32419121e=6.33

22、707454matlab 数学实验e=6.30357011e=6.35063255x=6.33557468y=0.003042432.估计程序变动如下：b=a(j)*3;c=b*sin(b)*3;结果为：e=3.12211717e=3.11373037e=3.07484948e=3.08660758e=3.10052243e=3.10475698e=3.13762746e=3.16481618e=3.11552000e=3.09615989x=3.11167076y=0.006692233.估计程序变动如下：b=exp(-a(j)2/2);c=b/n*(2*pi)0.5;e(i)=d+c/2;

23、结果为：e=0.88617655e=0.88538972e=0.88635209e=0.88575809e=0.88653705e=0.88606366e=0.88634011e=0.88613926e=0.88573325e=0.88644154x=0.88609313y=0.000001324.估计程序变动如下：b=exp(a(j)2);c=b/n;e(i)=d+c;结果为：e=1.46211146e=1.46154792matlab 数学实验e=1.46327379e=1.46256348e=1.46318297e=1.46235828e=1.46241378e=1.46316145e=

24、1.46203052e=1.46280489x=1.46254485y=0.000003205.估计程序变动如下：b=a(i)*4;c=1/(1+b2)0.5);e(i)=d+c*4/n;结果为：e=1.98511173e=1.02167881e=1.26713031e=0.98770837e=1.14216662e=1.75642022e=1.97055988e=1.96227794e=1.83229787e=1.06190231x=1.49872541y=1.90228258【实验结果与分析】通过对实验所得平均值与真实值的比较， 可以看出实验结果与真实值相比非常接近， 而且样本方差很小，从

25、而说明概率分布的选取比较适当，计算机实验很准确。实验 7：MATLAB与马尔科夫预测模型【实验目的 】基于 matlab编程应用马尔可夫预测模型【实验原理 】马尔可夫通过实践认为：世界上无论是社会领域还是自然领域，有一类事物的变化过程只与事物的近期状态有关，与事物的过去状态无关，这类事物的性质称为无后效性。例如，事物， 从初始状态 (0)起，变动一次后为 (1) ，变动 n 次后为 (n)，则 (n) 仅与 (n-1)有关，与 n-1 以后的各次变动无关。马尔可夫链：如果n 个连续变动的事物，在变动的过程中，其中任一次变动的结果都具有无后效性， 那么，这 n 个连续变动事物的集合， 就叫做马尔

26、可夫链， 这类事物的演变过程就叫做马尔可夫过程。【实验内容 】1.农业收成变化预测考虑某地区农业收成变化的三个状态，即“丰收”、“平收”和“欠收” 。记E1 为“丰收”matlab 数学实验状态， E2 为“平收”状态，E3 为“欠收”状态。下表给出了该地区1965 2004年期间农业收成的状态变化情况。试计算该地区农业收成变化的状态转移概率矩阵，并进行预测。使用 matlab 实现如下：P=0.2000 0.4667 0.3333;0.5385 0.1538 0.3077;0.3636 0.4545 0.1818;% 读入状态转移概率矩阵x=0,1,0;% 读入初始状态概率向量 (2004

27、年的农业收成状态 )fori=1:11% 预测今后11 年 (2005 2015)的农业收成状态y= x*Piend运行结果如下：y =0.53850.15380.3077y =0.30240.41480.2827y =0.38670.33340.2798y =0.35860.35890.2823y =0.36770.35090.2813y =0.36480.35340.2817y =0.36570.35260.2815y =0.36540.35290.2816y =0.36550.35280.2815y =0.36540.35280.2815y =0.36540.35280.28152.市场

28、占有率预测某厂对某产品的市场占有率和销售情况进行了调查：一月份共销售了 50 万件，其中普通、一级、特级品分别为 35、10、 5 万件。二月份中，一月份买普通品的顾客 25%的顾客转买一级品， 8%的顾客转买特级品；一月份买一级品的顾客 10%转买特级品， 3%转买普通品；一月份买特级品的顾客 2%买普通品， 15% 转买一级品。请预测以后月份各个等级产品的市matlab 数学实验场占有率。035, 10, 5 0.7,0.2,0.1由所给的资料可知50 50 50使用 matlab 实现如下：P=0.67 0.25 0.08;0.03 0.87 0.1;0.02 0.15 0.83;%读入

29、状态转移矩阵x=0.7 0.2 0.1;% 读入初始状态概率向量（一月份各等级产品的市场占有率）fori=1:11% 预测今年剩余11 个月各产品等级的市场占有率y= x*Piend运行结果如下：y =0.47700.36400.1590y =0.33370.45980.2065y =0.24150.51440.2441y =0.18210.54450.2734y =0.14380.56030.2959y =0.11910.56780.3131y =0.10310.57070.3262y =0.09270.57120.3361y =0.08600.57050.3435y =0.08160.56

30、940.3490y =0.07870.56810.3532结论：顾客对普通品的需求有减少的趋势，对一级品和特级品的需求有增加的趋势，因此，可以调整相应等级产品的产量。实验 8：基于 MATLAB的灰色预测模型【实验目的 】实验目的：掌握灰色预测模型及其应用基本内容：灰色预测模型的提出，建模以及实现代码。【实验内容 】灰色系统 ( Grey System) 理论是我国著名学者邓聚龙教授20世纪 80年代初创立的一种兼备软硬科学特性的新理论。该理论将信息完全明确的系统定义为白色系统，将信息完全不明确的系统定义为黑色系统，将信息部分明确、 部分不明确的系统定义为灰色系统。由于客观世界中，诸如工程技术

31、、社会、经济、农业、环境、军事等许多领域，大量存在着信息不完全的情况。要么系统因素或参数不完全明确，因素关系不完全清楚；要么系统结构不完全知道，系统的作用原理不完全明了等，从而使得客观实际问题需要用灰色系统理论来解决。灰色预测是应用灰色模型GM(1， 1) 对灰色系统进行分析、建模、求解、预测的过程。由于灰色建模理论应用数据生成手段，弱化了系统的随机性，使紊乱的原始序列呈现某种规律，规律不明显的变得较为明显， 建模后还能进行残差辨识，即使较少的历史数据， 任意随机分布，也能得到较高的预测精度。因此，灰色预测在社会经济、管理决策、农业规划、气象生态等各个部门和行业都得到了广泛的应用。一、 GM（

32、1，1）模型建立matlab 数学实验设有 k个原始非负样本序列统理论采用了独特的数据预处理方式，对序列为揭示系统的客观规律，灰色系 进行一阶累加生成，即 AGO生成，由此得生存数列：GM(1,1) 模型的原始形式为：为的紧邻值生存序列其中GM(1,1) 模型的基本形式：若 =为参数列，且,则GM(1,1) 模型的最小二乘估计参数列满足关于的白化方程也叫影子方程为：定理： 白化方程的解也称时间响应函数为：GM(1,1) 模型的时间响应序列为matlab 数学实验实际预测值为二、模型检验为确保所建立的GM(1,1) 模型有较高的预测精度，还需要进行以下检验（1）求出及之残差 e(k) 、相对误差

33、和平均相对误差：,，.（2）求出原始数据平均值x:.三、残差修正模型记，当然我们知道并不一定全为负或者全为正，这时我们令, 同时令, 则是一个非负序列，我们可以用方法来建立它的GM模型，求解，得出其预测值, 而后还原残差预测值, 最后用修正原来的预测值，得到修正后的预测值.四、实现代码clear allX0=input('请输入序列矩阵');%输入数据请用如例所示形式：48.757.17 68.76 92.15或者43823050,44649620,45793750 46613 48526270 49713050，该向量为原始向量X0n=length(X0);for i=2:n

34、 %开始进行建模可行性分析Q(i)=X0(i-1)/X0(i);endQ(1)=;ma=max(Q);mi=min(Q);if ma>exp(2/(n+1)disp('序列无法进行灰色预测');returnmatlab 数学实验elseif mi<exp(-2/(n+1)disp('序列无法进行灰色预测');returnelsedisp('序列可以进行灰色预测');endclear Q ma mi %检验结束X1=cumsum(X0);%累加生成算子向量X1Z1=ones(n-1,2);for i=1:(n-1)Z1(i,1)=-(X

35、1(i)+X1(i+1)/2;Z1(i,2)=1;%均值生成算子 Z1endZ1T=Z1'% 均值生成算子矩阵Z1的转置 Z1Tfor j=1:n-1Y(j)=X0(j+1);endYT=Y'A=inv(Z1T*Z1)*Z1T*YT;%最小二乘估计计算参数a、 ua=A(1);%Z1 参数 au=A(2);% 系统给定参数 ut=u/a;t_test=input('请输入需要预测个数：');i=1:t_test+n;X1S(i+1)=(X0(1)-t).*exp(-a.*i)+t;% 计算时间响应序列，得出估计累加向量 X1S X1S(1)=X0(1);X0S(

36、1)=X0(1);for j=n+t_test:-1:2X0S(j)=X1S(j)-X1S(j-1);% 计算 X1S的逆累加向量 X0S,还原 X0得到估计值 endfor i=1:nQ(i)=X0S(i)-X0(i);%求残差E(i)=abs(Q(i)/X0(i);%求相对误差endAVG=sum(E)/(n-1);% 求平均相对误差av=input('请输入允许的平均相对误差if AVG>=av;% 如果平均相对误差大于');% 输入如 0.1 ，或 0.05 等形式，不要用av%，则进入残差GM模型5%这类形式clear cn Q1 CZ1 CZ1T CY CYT

37、 CA ca cu ct Q1Scn=length(Q);Q1=cumsum(Q);% 累加生成算子向量Q1CZ1=ones(cn-1,2);for i=1:(n-1)CZ1(i,1)=-(Q1(i)+Q1(i+1)/2;CZ1(i,2)=1;%均值生成算子 CZ1matlab 数学实验endCZ1T=CZ1'%均值生成算子矩阵CZ1的转置 CZ1Tfor j=1:cn-1CY(j)=Q(j+1);endCYT=Y'CA=inv(CZ1T*CZ1)*CZ1T*CYT;%最小二乘估计计算参数ca、 cuca=CA(1);%CZ1参数 acu=CA(2);% 系统给定参数cuct=

38、cu/ca;i=1:t_test+cn;Q1S(i+1)=(Q(1)-ct).*exp(-ca.*i)+ct;%计算时间响应序列，得出估计累加向量Q1SX1S=X1S+Q1S;%将残差拟合值加入，提高精度for j=n+t_test:-1:2X0S(j)=X1S(j)-X1S(j-1);%计算 X1S的逆累加向量endclear avfor i=1:cnQ(i)=X0S(i)-X0(i);%求残差E(i)=abs(Q(i)/X0(i);%求相对误差endAVG=sum(E)/(n-1);% 求平均相对误差endx=1:n;xs=2:n+t_test;yn=X0S(2:n+t_test);plo

39、t(x,X0,'r',xs,yn,'*-b');%作图X0Sdisp('disp('百分平均相对误差为：',num2str(AVG*100),'%');拟合值为：',num2str(X0S(1:n+t_test);disp(A);【实验目总结 】1、 灰色预测模型的应用范围有限，在使用时不能超过其能预测的范围2、 灰色预测模型的模型建立，以及检验非常重要，可以减少预测带来的误差实验 9：基于 MATLAB的模糊聚类分析【实验目的】1. 加强对 MA TLAB 软件使用的能力。2. 加强对模糊数学学习的理解。3. 基

40、于 MATLAB 对数据挖掘中使用模糊数学进行聚类分析学会使用 MATLAB 中的模糊工具箱。【实验方案】1. 实验原理首先，我们要了解一般聚类分析分为三个步骤： （1）数据标准化，（2）标定，（ 3）聚类。然后我们在来了解下什么是模糊数学：matlab 数学实验（ 1）模糊数学是1956 年，美国加利福尼亚大学控制论专家扎德（L.A.Zadeh ）教授提出的。是一门研究和处理模糊性现象的数学方法。（ 2）对于有限论域 U x1 , x2 K xn 构造映射 A(x) ：U 0,1 ，确定 U 上的模糊子集 A ，映射 A(x) 称为 A 的隶属函数，它表示x 对 A 的隶属程度。（ 3）模糊

41、集的基本运算：相等：A=BA(x) = B( x)；包含：A BA(x) B(x)；并： A B 的隶属函数为(AB)( x)=A(x)B(x)=max A(x)， B(x） ；交： A B 的隶属函数为(AB)( x)=A(x)B(x)=min A(x)，B(x） ；余： Ac 的隶属函数为Ac(x) = 1- A(x).2. 基于 MATLAB 利用模糊数学实现聚类分析。（ 1） 数据标准化。设论域 X = x1, x2, ,xn 为被分类对象 ,每个对象又由 m 个指标表示其性状 : xi = xi 1, xi2, xim, i = 1, 2, n。,根据模糊矩阵的要求，将数据压缩到0,

42、1 区间上。通常需要作如下几种变换：平移? 标准差变换xikxikxk (i1,2,.,n; k1,2,., m)sk1 n1n平移 ? 极差变换xkxk )2其中xik , skn i(xikxikn i 1|1 i n1xikmin xikinmin xik | 1 in对数变换max xik |1xiklg xik (i1,2, n; k1,2, m)取对数以缩小变量间的数量级。（ 2） 标定相似性度量，又称标定，就是根据实际情况，按一定准则或某一种方法，给论域 X 中的元素两两之间都赋以 0,1 内的一个数，称为相似系数。它的大小表征两个元素彼此接近或相似的程度。用 rij 表示元素

43、xi 与 xj 的相似系数，其中：xi = xi1, xi 2,xim,i = 1, 2,n,xj = xj1, xj 2,xjm,j = 1, 2,n,rij 0,1数量积法1,ijm夹角余弦法其中 M max(xik x jk )rij1 mmxikx jk , iji jk 1相关系数法Mk 1xik x jkrijk1mmm2| xik2x j |最大最小法xikxxjki | | x jkrijmk 1k 1k1mmx jk。( xik)2( x jkxj )2( xikxi)rijk1k 1k1m一般海明距离法x jk)( xikk1matlab 数学实验m一d般(x欧i,式x距j)离| xik