离散元课件二
《离散数学》课件-第四章 二元关系
R2= R • R={<1,1>，<2,2>，<1,3>，<2,4>， <3,5>}
R3=R2 • R={<1,2>，<2,1>，<1,4>，<2,3>， <2,5>}
R4= R3 • R={<1,1>，<2,2>，<1,5>，<2,4>，
从关系图来看关系的n次幂
R：
1
2
3
4
5
R2：
1
2
3
4
5
R2就是从R的关系图中的任何一个结点x出发，长 为2的路径，如果路径的终点是y，则在R2 的关系 图中有一条从x到y的有向边。其他以次类推：
R3：
1
2
3
4
5
R4：
1
2
3
4
5
定理 设|A|=n，R A×A，则必有i，j∈N， 0≤i＜j≤2n2，使得Ri=Rj。
＝R5，R7＝R6•R＝R5，…，Rn＝R5 (n＞5) 故Rn{R0,R1,R2,R3,R4,R5}。
S0=IA,S1＝S,
S2＝S•S＝{<a,c>,<b,d>,<c,e>,<d,f>}， S3＝S•S•S＝S2•S＝{<a,d>,<b,e>,<c,f>}， S4＝S3•S＝{<a,e>,<b,f>}， S5＝S4•S＝{<a,f>}， S6＝S5•S＝Φ， S7＝Φ， …， 故，Sn{S0,S1,S2,S3,S4,S5,S6}
离散完整ppt课件5.2-3共23页文档
定义 非空集合 S 和 S 上 k 个一元或二元运算 f1, f2, … , fk 组成的系统称为一个代数系统, 简称代 数，记做 V=<S, f1, f2, … , fk>.
S 称为代数系统的载体， S 和运算叫做代数系 统的成分. 有的代数系统定义指定了S中的特殊 元素，称为代数常数, 例如二元运算的单位元. 有时也将代数常数作为系统的成分.
6
积代数
定义 设 V1=<S1,o>和 V2=<S2,>是代数系统，其中 o 和 是二元运算. V1 与 V2 的 积代数 是V=<S1S2,∙>,
<x1,y1>, <x2,y2>S1S2 , <x1,y1> ∙ <x2,y2>=<x1ox2, y1y2>
例3 V1=<Z,+>, V2=<M2(R), ∙ >, 积代数< ZM2(R),o> <z1,M1>, <z2,M2>ZM2(R) , <z1,M1> o <z2,M2> = <z1+z2, M1∙M2>
单同态、满同态、同构 自同态
同态映射的性质
9
同态映射的定义
定义 设 V1=<S1,∘>和 V2=<S2,>是代数系统，其中 ∘ 和 是二元运算. f: S1S2, 且x,yS1, f (x∘y) = f(x) f( y), 则称 f 为V1到 V2 的同态映射，简称同态.
10
更广泛的同态映射定义
f (x∘y)=f(x)f(y), f (x∙y)=f(x)◊f(y), f (∆ x)=∇f(x) 则称 f 为V1到 V2 的同态映射，简称同态.
离散完整ppt课件3.1-3共41页
证明 X=Y
命题演算法 等式代入法 反证法 运算法
以上的 X, Y 代表集合公式
18
命题演算法证 XY
任取 x ， xX … xY
例3 证明AB P(A)P(B) 任取x xP(A) xA xB xP(B) 任取x xA {x}A {x}P(A) {x}P(B) {x}B xB
13
例1
F:一年级大学生的集合
S：二年级大学生的集合
R：计算机系学生的集合
M：数学系学生的集合
T：选修离散数学的学生的集合
L：爱好文学学生的集合
P：爱好体育运动学生的集合
所有计算机系二年级学生都选修离散数学
数学系一年级的学生都没有选修离散数学
数学系学生或爱好文学或爱好体育运动 只有一、二年级的学生才爱好体育运动 除去数学和计算机系二年级学生外都不 选修离散数学3.2 集合的基本运算
集合基本运算的定义
文氏图（John Venn） 例题 集合运算的算律 集合包含或恒等式的证明
10
集合基本运算的定义
并 交 相对补 对称差
绝对补
AB = { x | xA xB } AB = { x | xA xB } AB = { x | xA xB } AB = (AB)(BA)
由已知包含式通过运算产生新的包含式 XY XZYZ, XZYZ
= (AB)(AB) A = EA
11
文氏图表示
12
关于运算的说明
运算顺序： 和幂集优先，其他由括号确定 并和交运算可以推广到有穷个集合上，即
A1A2…An= {x | xA1xA2…xAn} A1A2…An= {x | xA1xA2…xAn} 某些重要结果 ABA AB AB=（后面证明） AB= AB=A
离散数学第2章 关系(祝清顺版)
离散数学
关系矩阵的实例
例9 设A={3, 5, 6, 9}, A上的二元关系
R={<x, y｜x＞y},
试求出关系矩阵。
[解] 关系的集合表示为:
R＝{9, 3, 9, 5, 9, 6, 6, 3, 6, 5, 5, 3}.
关系矩阵为: 0 1 MR= 1 1
关系的三种表示方法： 集合表达式 关系矩阵
关系图
关系矩阵和关系图可以表示有限集合上的关系。
离散数学
第二章 二元关系
2007年8月20日
关系矩阵
设给定集合A={a1,a2,…,an}，B={b1,b2,…,bm}，R为从A到B
的一个二元关系，构造一个n×m矩阵。用集合A的元素标注矩
阵的行，用集合B的元素标注矩阵的列，对于aiA和bjB，令
n2 n2
个. 不
每一个子集代表一个A上的二元关系，所以A上有 2 同的二元关系。 |Ai|=mi，则A1×A2×…×An上有 2 二元关系。
离散数学 第二章 二元关系 2007年8月20日
m1m2…mn
个不同的
常用的关系
定义 对任意集合A，定义 (1) 空关系 (2) 全域关系 EA={<x, y>|x∈A且y∈A}=A×A (3) 恒等关系 IA={<x, x>|x∈A} (4) 小于或等于关系：LA={<x, y>|x, y∈A且x≤y}, 其中 AR。 (5) 整除关系：DA={<x, y>|x, y∈A且x整除y}, 其中 AZ* , Z*是非零整数集 (6) 包含关系：R＝{<x, y>|x, y∈A且xy}, 其中A是集 合族。
离散数学课件第六章(第2讲)
《定理》：设*是S上的二元运算,且x S,对任一m,n
I+有（1）xmxn=xm+n
（2）(xm)n=xmn
证明： （1） xmxn= (xm x) x… x = (xm+1 x) x… x
n
n-1
=….= xm+n
（2）(xm)n= xm … xm= xm+m xm … xm=…=xmn
n
例：设M= {0º,60º,120º,240º,300º,180º}表示平面上几何图形 顺时针旋转的六种位置，定义一个二元运算*，对M中任一 元素a,b有a*b=图形旋转(a+b)的角度，并规定当旋转到 360º时即为0º，试验证<M ,*>是一个群。
* 0º 60º 120º 180º 240º 300º 0º 0º 60º 120º 180º 240º 300º 60º 60º 120º 180º 240º 300º 0º 120º 120º 180º 240º 300º 0º 60º 180º 180º 240º 300º 0º 60º 120º 240º 240º 300º 0º 60º 120º 180º 300º 300º 0º 60º 120º 180º 240º
例： <I ,max>，其中max(x1,x2)取二者之大值；<I ,min>， 其中min(x1,x2)取二者之小值，均不为独异点（不存在幺 元）。<N ,max>则为独异点，其中 e =0
《定义》：设< S ,* >是一半群，TS，且*在T上是封闭的， 那么< T ,* >也是半群，称< T ,* >是< S ,* >的子半群。
离散元课件
转动方程：
转动方程可以表示为
dωi Ii Ti dt j 1
式中，I i 与ωi 分别为颗粒 i 的转动惯量与角速度，对于 球形颗粒 I i为
ki
2 I i mi Ri2 5
二 基本原理-球形颗粒元离散元法
接触模型
综述：
关于接触力的计算模型已有大量的研究成果，目前仍旧是 一个活跃的研究领域，特别是对于切向力的计算方法。
二 基本原理
离散化模型
图1 颗粒元与块体元示意图
二 基本原理-球形颗粒元离散元法
俞 缙
bugyu0717@
华侨大学岩土工程研究所
二 基本原理
根据离散化模型中所采用的单元种类分别介绍离 散元法的基本原理：
颗粒元
• 二维圆盘单元 • 三维圆球单元
块体元
• 多边形单元 • 多面体单元
对于理想散体颗粒（无粘连）：采用 Hertz 理论描述法向作用，而采 用Mindlin与Deresiewicz理论描述切向作用； 对于存在粘连的散体颗粒：法向接触力根据在 Hertz 理论基础上考虑 粘连力的JKR（Johnson-Kendall-Roberts）理论确定，切向接触力增量 则根据把 Savkoor 和 Briggs 理论与 Mindlin 和 Deresiewicz理论相结合形 成的理论确定。
1988 年 Cundall 所在的 ITASCA 咨询公司推出针对三维块体元的 3DEC程序。至此，离散元的理论体系基本形成。
一 历史由来及研究现状
早期的离散单元法
Cundall称之为“Distinct Element Method”，随着该方法的推广， 有的学者称其为“ Discrete Element Method” ，缩写形式均为 DEM。 最初，离散元的研究对象主要是岩石等非连续介质的力学行为， 它的基本思想是把不连续体分离为刚性元素的集合，使各个刚性 元素满足运动方程，用时步迭代的方法求解各刚性元素的运动方 程，继而求得不连续体的整体运动形态。离散元法允许单元间的 相对运动，不一定要满足位移连续和变形谐调条件，计算速度快， 所需存储空间小，尤其适合求解大位移和非线性问题。
离散完整ppt课件2.1-2共25页
▪ 个体词 ▪ 谓词 ▪ 量词 ▪ 一阶逻辑中命题符号化
1
基本概念——个体词、谓词、量词
个体词（个体）: 所研究对象中可以独立存在的具 体或抽象的客体
个体常项：具体的事物，用a, b, c表示 个体变项：抽象的事物，用x, y, z表示 个体域: 个体变项的取值范围
有限个体域，如{a, b, c}, {1, 2} 无限个体域，如N, Z, R, … 全总个体域: 宇宙间一切事物组成
12
原子公式
定义 设R(x1, x2, …, xn)是任意的n元谓词，t1,t2,…, tn 是任意的n个项，则称R(t1, t2, …, tn)是原子公式. 原子公式是由项组成的n元谓词. 例如，F(x,y), F(f(x1,x2),g(x3,x4))等均为原子公式
13
合式公式
定义 合式公式（简称公式）定义如下： (1) 原子公式是合式公式. (2) 若A是合式公式，则 (A)也是合式公式 (3) 若A, B是合式公式，则(AB), (AB), (AB), (AB)也是合式公式 (4) 若A是合式公式，则xA, xA也是合式公式 (5) 只有有限次地应用(1)~(4)形成的符号串是合 式公式.
15
公式的解释与分类
给定公式 A=x(F(x)G(x)) 成真解释: 个体域N, F(x): x>2, G(x): x>1
代入得A=x(x>2x>1) 真命题 成假解释: 个体域N, F(x): x>1, G(x): x>2
(2) x (F(x)G(x))
这是两个基本公式, 注意这两个基本公式的使用.
7
一阶逻辑中命题符号化(续)
例3 在一阶逻辑中将下面命题符号化
离散数学第2章ppt课件
C
n
A k A 1A 2 A n
k 1
二、集合的并 (Union)
3、性质
1)幂等律 A∪A =A
2)零律
A∪U =U
3)同一律 A∪ =A
4)交换律 A∪B =B∪A
5)结合律 A∪(B∪C) =(A∪B)∪C
二、集合的并 (Union)
3、性质
， 6)
若A⊆B，C⊆D，则A∪C
是集合，没有元素
有1个元素的集合
2) ∈{}， {}
五、特殊集合
1、空集
定理 空集是任一集合A的子集，即 ⊆A。
下列命题是否为真。
1)√⊆;
2) ∈ ; 3) ⊆{}; 4) ∈{} 。
√
√
五、特殊集合
1、空集
推理 空集是唯一的。(绝对唯一)
证明： 设1，2是两个空集, 则1 2，且2 1，
证明唯一性 一般采用反
1、符号表示法
通常用大写字母A, B, C, …代表集合; 用小写字母a, b, c, …代表元素。
1)如果a是集合A的一个元素, 则记为 a∈A, 读做“a属于A”,或 “a在集合A中”。
2)如果a不是集合A的一个元素, 则记为 a∈A, 读做“a不属于A”,或 “a不在集合A中”。
注：任一元素, 对某一集合而言, 或属于该集合, 或不属于该集合, 二者必居其一, 且只居其一。
1) 若b∈A，则b是不给自己刮脸的人， 而由题意，b只给集合A中的人刮脸。 ∴b 要给b 刮脸， 即b ∈ A。
理发师问题
在一个很僻静的孤岛上，住着一些人家，岛上只 有一位理发师，该理发师专给那些并且只给那些自己 不刮脸的人刮脸。那么，谁给这位理发师刮脸？
方世昌离散数学第三版教材课件第3章二元关系(可编辑)
方世昌离散数学第三版教材课件第3章二元关系31 基本概念 32 关系的合成 33 关系上的闭包运算 34 次序关系 35 等价关系和划分 31 基本概念311 关系关系的数学概念是建立在日常生活中关系的概念之上的让我们先看两个例子例31-1 设 A abcd 是某乒乓球队的男队员集合 B efg 是女队员集合如果A和B元素之间有混双配对关系的是a 和gd和e我们可表达为 R 〈ag〉〈de〉这里R表示具有混双配对关系的序偶集合所有可能具有混双配对关系的序偶集合是 A×B 〈xy〉x∈A∧y∈B 〈ae〉〈af〉〈ag〉〈be〉〈bf〉〈bg〉〈ce〉〈cf〉〈cg〉〈de〉〈df〉〈dg〉例31-2 设学生集合A1 abcd 选修课集合A2 日语法语成绩等级集合A3 甲乙丙如果四人的选修内容及成绩如下 a 日乙 b 法甲c 日丙 d 法乙我们可表达为S 〈a日乙〉〈b法甲〉〈c日丙〉〈d法乙〉这里S表示学生和选修课及成绩间的关系而可能出现的全部情况为 A1×A2×A3 〈xyz〉x∈A1∧y∈A2∧z∈A3 〈a日甲〉〈a 日乙〉〈a日丙〉〈a法甲〉〈a法乙〉〈a法丙〉〈b日甲〉〈b日乙〉〈b日丙〉〈b 法甲〉〈b法乙〉〈b法丙〉〈c日甲〉〈c日乙〉〈c日丙〉〈c法甲〉〈c法乙〉〈d法丙〉定〈c法丙〉〈d日甲〉〈d日乙〉〈d日丙〉〈d法甲〉〈d法乙〉义31―1 1 A×B的子集叫做A到B的一个二元关系2 A1×A2××An n≥1 的子集叫做A1×A2××An上的一个n元关系3 从定义可看出关系是一个集合所有定义集合的方法都可用来定义关系例31-1和例31-2是列举法的例子一个谓词Px1x2xn 可以定义一个n元关系R R 〈x1x2xn〉P x1x2xn 例如实数R上的二元关系＞可定义如下＞〈xy〉x∈R∧y∈R∧x＞y 反之一个n元关系也可定义一个谓词当n 1时R 〈x〉P x 称为一元关系它是一重组集合表示论述域上具有性质P的元素集合其意义与R xP x 相同仅记法不同而已例如设P x 表示x是质数论述域是N则质数集合可表示为〈x〉｜P x 或x｜P x 关系也可归纳地定义自然数上的小于关系可定义如下1 基础〈01〉∈＜2 归纳如果〈xy〉∈＜那么i 〈xy1〉∈＜ ii 〈x1y1〉∈＜ 3 极小性对一切xy∈Nx＜y当且仅当〈xy〉是由有限次应用条款 1 和 2 构成定义31―2 设R是的子集如果R 则称R为空关系如果则称R为全域关系现在定义关系相等的概念在关系相等的概念中不仅需要n重组集合相等还需其叉积扩集也相同定义31―3设R1是上的n元关系R2是上的m元关系那么R1 R2当且仅当n m且对一切i1≤i≤nAi Bi并且R1和R2是相等的有序n重组集合 312 二元关系最重要的关系是二元关系本章主要讨论二元关系今后术语关系都指二元关系若非二元关系将用三元或n元一类术语指出二元关系有自己专用的记法和若干新术语设 A x1x2x7 B y1y2y6 R〈x3y1〉〈x3y2〉〈x4y4〉〈x6y2〉 A到B的二元关系R可如图31―1那样形象地表示〈x3y1〉∈R也可写成x3Ry1称为中缀记法读做x3和y1有关系R中缀记法常用来表示诸如＜＞等关系例如〈35〉∈＜通常写作3＜5 A叫做关系R的前域B叫做关系R的陪域 D R x｜y 〈xy〉∈R 叫做关系R的定义域R R y｜x 〈xy〉∈R 叫做关系R的值域关系是序偶的集合对它可进行集合运算运算结果定义一个新关系设R和S是给定集合上的两个二元关系则R∪SR∩SR-S 等可分别定义如下x R∪S y xRy∨xSy x R∩S y xRy∧xSy x R-S y xRy∧xy x y xRy 例31-3平面上的几何图形是平面R2的子集也是一种关系设参看图31―2 R1 〈xy〉｜〈xy〉∈R2∧x2y2≤9 R2 〈xy〉｜〈xy〉∈R2∧1≤x≤3 ∧0≤y≤3 R3 〈xy〉｜〈xy〉∈R2∧x2y2≥4 则 R1∪R2 〈xy〉｜〈xy〉∈R2∧ x2y2 ≤9∨ 1≤x≤3∧0≤y≤3 R1∩R3〈xy〉｜〈xy〉∈R2∧ x2y2 ≤9∧x2y2≥4 R1-R3 〈xy〉｜〈xy〉∈R2∧ x2y2≤9∧ L x2y2≥4 〈xy〉｜〈xy〉∈R2∧ x2y2≥4 313 关系矩阵和关系图表达有限集合到有限集合的二元关系时矩阵是一有力工具定义31―4 给定集合A a1a2am 和B b1b2bn 及一个A到B的二元关系R 使例31-4 设A a1a2 B b1b2b3 R 〈a1b1〉〈a2b1〉〈a1b3〉〈a2b2〉则其关系矩阵为例31-5 设A 1234 A上的二元关系R 〈xy〉｜x＞y 试求出关系矩阵解R 〈41〉〈42〉〈43〉〈31〉〈32〉〈21〉例31-6 设 A 12345 R 〈12〉〈22〉〈32〉〈34〉〈43〉其图示如图31―3所示图中结点5叫做孤立点利用关系R的图示也可写出关系R 314 关系的特性在研究各种二元关系中关系的某些特性扮演着重要角色我们将定义这些特性并给出它的图示和矩阵的特点定义31―5 设R是A上的二元关系 1如果对A中每一xxRx那么R是自反的即 A上的关系R是自反的x x∈A→xRx A 123 R1 〈11〉〈22〉〈33〉〈12〉是自反的其关系图和关系矩阵的特点如图31―4所示 2 如果对A中每一xxRx那么R是反自反的即 A上的关系R是反自反的 x x∈A→xRx 例如 A 123 R2 〈21〉〈13〉〈32〉是反自反的其关系图和关系矩阵的特点如图31―5所示有些关系既不是自反的又不是反自反的如图31―6 例如R3 〈11〉〈12〉〈32〉〈23〉〈33〉 3 如果对每一xy∈AxRy蕴含着yRx那么R是对称的即A上的关系R是对称的x y x∈A∧y∈A∧xRy→yRx 例如A 123 R4 〈12〉〈21〉〈13〉〈31〉〈11〉是对称的其关系图和关系矩阵的特点如图31―7所示 4 如果对每一xy∈AxRyyRx蕴含着x y那么R是反对称的即A上的关系R是反对称的x y x∈A∧y∈A∧xRy∧yRx→x y 例如A 123 R5 〈12〉〈23〉是反对称的其关系图和关系矩阵的特点如图31―8所示 5 如果对每一xyz∈AxRyyRz蕴含着xRz那么R是传递的即A上的关系R是传递的x y z x∈A∧y∈A∧z ∈A∧xRy∧yRz→xRz 例如A 1234R5 〈41〉〈43〉〈42〉〈32〉〈31〉〈21〉是传递的其关系图和关系矩阵如图31―10所示例31-7 1 任何集合上的相等关系是自反的对称的反对称的和传递的但不是反自反的 2 整数集合I上关系≤是自反的反对称的可传递的但不是反自反的和对称的关系＜是反自反的反对称的可传递的但不是自反的和对称的 3 设 ab 试考察上的下列关系 i 关系与有同样长度是自反的对称的可传递的但不是反自反的和反对称的 ii xRy当且仅当x是y的真词头这里R是反自反的反对称的可传递的但不是自反的和对称的 iii xRy当且仅当x的某真词头是y的一个真词尾这里R既不是自反的又不是反自反的因为aaRaa但abRab既不是对称的也不是反对称的并且不是传递的 4 非空集合上的空关系是反自反的对称的反对称的和传递的但不是自反的空集合上的空关系则是自反的反自反的对称的反对称的和可传递的 5 基数大于1的集合上的全域关系是自反的对称的和传递的但不是反自反的和反对称的例如图31―11所示的关系 321 关系的合成前边已经指出关系是序偶的集合因此可以进行集合运算本节介绍一种对关系来说更为重要的运算合成运算假设R1是A到B的关系R2是B到C的关系参看图32-1合成关系R1R2是一个A到C的关系如果在关系图上从a∈A到c∈C有一长度路径中弧的条数为2的路径其第一条弧属于R1其第二条弧属于R2那么〈ac〉∈R1R2合成关系R1R2就是由〈ac〉这样的序偶组成的集合其第一条弧属于R1其第二条弧属于R2那么〈ac〉∈R1R2合成关系R1R2就是由〈ac〉这样的序偶组成的集合定义32―1 设R1是从A到B的关系R2是从B到C的关系从A到C的合成关系记为R1R2定义为 R1R2 〈ac〉｜a∈A∧c∈C∧b〔b∈B∧〈ab〉∈R1∧〈bc〉∈R2〕例32-11 如果R1是关系是的兄弟R2是关系是的父亲那么R1R2是关系是的叔伯R2R2是关系是的祖父 2 给定集合A 1234 B 234 C 123 设R是A到B的关系S是B到C的关系 R 〈xy〉｜xy 6 〈24〉〈33〉〈42〉S 〈yz〉｜y-z 1 〈21〉〈32〉〈43〉则R·S 〈23〉〈32〉〈41〉如图32―2所示 3 设A 12345 R和S都是A上二元关系如果 R 〈12〉〈34〉〈22〉 S 〈42〉〈25〉〈31〉〈13〉则R·S 〈15〉〈32〉〈25〉 S·R 〈42〉〈32〉〈14〉 R·S ·R 〈32〉 R· S·R 〈32〉 R·R 〈12〉〈22〉 S·S〈45〉〈33〉〈11〉 4 设R是A到B的二元关系IAIB分别是A和B上的相等关系则IA·R R·IB R 5 如果关系R的值域与关系S的定义域的交集是空集则合成关系R·S是空关系下边介绍合成关系的性质定理32―1 设R1是从A到B的关系R2和R3是从B到C的关系R4是从C到D的关系那么 1 R1 R2∪R3 R1R2∪R1R3 2 R1 R2∩R3 R1R2∩R1R3 3 R2∪R3 R4 R2R4∪R3R4 4 R2∩R3 R4 R2R4∩R3R41 2 3 部分的证明留作练习我们仅证明 2 部分证先证明公式因为〈ac〉∈R1 R2∩R3 b〔〈ab〉∈R1∧〈bc〉∈R2∩R3 〕b〔〈ab〉∈R1∧〈bc〉∈R2∧〈bc〉∈R3 〕b〔〈ab〉∈R1∧〈bc〉∈R2 ∧〈ab〉∈R1∧〈bc〉∈R3 〕b〔〈ab〉∈R1∧〈bc〉∈R2〕∧b〔〈ab〉∈R1∧〈bc〉∈R3〕〈ac〉∈R1R2∧〈ac〉∈R1R3 〈ac〉∈R1R2∩R1R3 即〈ac〉∈R1 R2∩R3 〈ac〉∈R1R2∩R1R3 所以R1 R2∩R3 R1R2∩R1R3 再证包含可能是真包含举反例证明如果 A a B b1b2b3 C c A到B的关系R1〈ab1〉〈ab2〉 B到C的关系R2 〈b1c〉〈b3c〉 B到C的关系R3〈b2c〉〈b3c〉那么R1 R2∩R3 R1R2∩R1R3 〈ac〉此时R1 R2∩R3 ≠R1R2∩R1R3证毕定理32―2 设R1R2和R3分别是从A到BB到C和C到D的关系那么 R1R2 R3 R1 R2R3 证先证 R1R2R3 R1 R2R3 设〈ad〉∈ R1R2 R3那么对某c∈C〈ac〉∈R1R2和〈cd〉∈R3因为〈ac〉∈R1R2存在b∈B使〈ab〉∈R1和〈bc〉∈R2因为〈bc〉∈R2和〈cd〉∈R3得〈bd〉∈R2R3所以〈ad〉∈R1 R2R3 这样就证明了 R1R2 R3 R1 R2R3 R1 R2R3 R1R2 R3的证明是类似的留给读者自证上述证明也可用等价序列表达 322 关系R的幂当R是A上的一个关系时R可与自身合成任意次而形成A上的一个新关系在这种情况下RR常表示为R2RRR表示为R3等等我们能归纳地定义这一符号如下定义32―2设R是集合A上的二元关系n∈N那么R的n次幂记为Rn定义如下 1R0是A上的相等关系R0 〈xx〉｜x∈A 2 Rn1 Rn·R 定理32―3 设R是A上的二元关系并设m和n是N的元素那么 1Rm·Rn Rmn 2 Rm n Rmn 可用归纳法证明请读者自证定理32―4 设｜A｜ nR是集合A上的一个关系那么存在i和j使Ri Rj而0≤i＜j≤证A上的每一二元关系是A×A的子集因为｜A×A｜ n2｜ρ A×A ｜因此A上有个不同关系所以R的不同的幂不会超过个但序列R0R1 有项因此R的这些幂中至少有两个是相等的证毕定理32―5 设R是集合A上的一个二元关系若存在i和ji＜j使Ri Rj记d j-i那么 1 对所有k≥0Rik Rjk 2 对所有km≥0Rimdk Rik 3 记S R0R1R2Rj-1 那么R的每一次幂是S的元素即对n∈NRn∈S 证 1 和 2 部分用归纳法证明留作练习3 对于 c 设n∈N如果n＜j那么根据S的定义Rn∈S假设n≥j那么我们能将n表示为imdk这里k＜d根据 b 部分得Rn Rik因为ik＜j这证明了Rn∈S定理中的ij在实用时宜取最小的非负整数以保证S中无重复元素例32-2 设 A abcd R 〈ab〉〈cb〉〈bc〉〈cd〉其关系图如图32―3所示则R0 〈aa〉〈bb〉〈cc〉〈dd〉 R2 〈ac〉〈bb〉〈bd〉〈cc〉 R3〈ab〉〈ad〉〈bc〉〈cb〉〈cd〉 R4 〈ac〉〈bb〉〈bd〉〈cc〉它们的关系图如图32―4所示由于R4 R2根据定理32―5 c 对所有n∈NRn∈ R0R1R2R3 可见不必再算了事实上易证R5 R3R6 R4 R2用归纳法可得R2n1 R3和R2n R2这里n≥1 323 合成关系的矩阵表达定理32―6 设X x1x2xm Y y1y2yn Z z1z2zp R是X到Y的关系MR 〔aij〕是m×n矩阵S是Y到Z的关系MS 〔bij〕是n×p矩阵则MR·S 〔cij〕 MR·MS这里证因为如果存在某k使aik和bki都等于1则cij 1但aik和bkj都等于1意味着xiRyk和ykSzj所以xi R·S zj可见如此求得的MR·S确实表达了R·S的关系因此上述等式是正确的如果不仅存在一个k使aik和bki都是1此时cij仍为1只是从xi到zj不止一条长度为2的路径但等式仍然正确上段的论证已隐含了不止一个k的情况本定理说明合成关系矩阵可用关系矩阵布尔矩阵的乘法表达例32-3设X 12 Y abc Z αβ R 〈1a〉〈1b〉〈2c〉 S 〈aβ〉〈bβ〉则定理32―7 关系矩阵的乘法是可结合的证利用关系合成的可结合性证明 MR·MS ·MT MR·S·MT M R·S ·T MR· S·T MR·MS·T MR· MS·MT 不仅合成关系可用关系矩阵表达而且关系的集合运算也可用关系矩阵表达设R和S是X到Y上的二元关系MR 〔aij〕MS 〔bij〕cij是运算后所得新关系之关系矩阵的元素则 MR∩S MR∧MS cij aij∧bij MR∪S MR∨MS cij aij∨bij cij aij MR-S MR∧ cij aij∧ bij 331 逆关系在讨论闭包运算时要用到逆关系的概念因此我们先介绍逆关系定义33―1设R是从A到B的二元关系关系R的逆或叫R的逆关系记为是一从B到A的二元关系定义如下例33-11 I上的关系2 集合族上的关系的逆是关系3 空关系的逆是空关系4 B×A即A×B的全域关系的逆等于B×A的全域关系定理33―1设R是从A到B的关系而S是从B到C 的关系则定理33―2 设RR1和R2都是从A到B的二元关系那么下列各式成立 332 关系的闭包运算关系的闭包运算是关系上的一元运算它把给出的关系R扩充成一新关系R′使R′具有一定的性质且所进行的扩充又是最节约的定义33―2设R是A上的二元关系R的自反对称传递闭包是关系R′使 i R′是自反的对称的传递的ii R′R iii 对任何自反的对称的传递的关系R〃如果R〃R那么R〃R′ R的自反对称和传递闭包分别记为r R s R和t R 由定义可以看出R的自反对称传递闭包是含有R并且具有自反对称传递性质的最小关系如果R已经是自反的对称的传递的那么具有该性质并含有R的最小关系就是R自身下一定理说明这一点定理33―4设R是集合A上的二元关系那么 a R是自反的当且仅当r R R b R是对称的当且仅当s R R c R是传递的当且仅当t R R 证 a 如果R是自反的那么R具有定义33―2对R′所要求的性质因此r R R反之如果r R R那么根据定义33―2的性质 i R是自反的b 和 c 的证明是类似的略构造R的自反对称和传递闭包的方法就是给R补充必要的序偶使它具有所希望的特性下面我们用关系图来说明如何实现这一点定理33―5 设R是集合A上的二元关系那么r R R ∪E 这里E是A上相等关系在本节中均如此证设R′ R∪E显然R′是自反的且R′R余下只需证明最小性现假设R〃是A上的自反关系且R〃R因R〃是自反的所以R〃E又R〃R所以R〃R∪E R′这样定义33―2都满足所以R′ r R 证毕设G是集合A上二元关系R的关系图我们把G的所有弧都画成有来有往即如果有从a到b的弧那么也有从b到a的弧就得到了R的对称闭包的有向图下一定理体现了这一想法定理33―7 设R 是集合A上的二元关系那么例33-2 a 整数集合I 上的关系＜的自反闭包是≤对称闭包是关系≠传递闭包是关系＜自身b 整数集合I上的关系≤的自反闭包是自身对称闭包是全域关系传递闭包是自身 c E的自反闭包对称闭包和传递闭包都是 E d ≠的自反闭包是全域关系对称闭包是≠≠的传递闭包是全域关系e 空关系的自反闭包是相等关系对称闭包和传递闭包是自身 f 设R是I上的关系xRy当且仅当y x1那么t R 是关系＜定理33―8设R是集合A上的二元关系这里A有n个元素那么证设〈xy〉∈t R 于是必存在最小的正整数k使〈xy〉∈Rk现证明k≤n若不然存在A的元素序列x a0a1a2ak-1ak y使xRa1a1Ra2ak-1Ry因k＞na0a1ak中必有相同者不妨设ai aj0≤i＜j≤k于是xRa1a1Ra2ai-1RaiajRaj1ak-1Ry 成立即〈xy〉∈Rs这里s k- j-i但这与k是最小的假设矛盾于是k≤n又〈xy〉是任意的故定理得证例33-3 设A abcd R如图33―1 a 所示则t R R∪R2∪R3∪R4如图33―1 b 所示本例即是32-2 定理33―9 1 如果R是自反的那么s R 和t R 都是自反的 2 如果R是对称的那么r R 和t R 都是对称的 3 如果R是传递的那么r R 是传递的定理33―10 设R是集合A上的二元关系那么 1 rs R sr R 2 rt R tr R 3 ts R st R 2 注意到ER RE R 和对一切n∈NEn E可得 34 次序关系 341 偏序集合定义34―1 如果集合A上的二元关系R是自反的反对称的和传递的那么称R为A上的偏序称序偶〈AR〉为偏序集合如果R是偏序〈AR〉常记为〈A ≤〉≤是偏序符号由于≤难以书写通常写作≤读做小于或等于因为小于或等于也是一种偏序故不会产生混乱R是偏序时aRb就记成a≤b 如果R是集合A上的偏序则 R 也是A上的偏序如果用≤表示R 可用≥表示R〈A≤〉和〈A ≥〉都是偏序集合并互为对偶例34-1 1 〈I≤〉是偏序集合这里≤表示整数中的小于或等于关系 2 〈ρ A 〉是偏序集合这里是集合间的包含关系 3 A 2468 D代表整除关系M代表整倍数关系则 D 〈22〉〈44〉〈66〉〈88〉〈24〉〈26〉〈28〉〈48〉 M 〈22〉〈44〉〈66〉〈88〉〈42〉〈62〉〈82〉〈84〉〈AD〉〈AM〉都是偏序集合且互为对偶例2 a P 1234 〈P≤〉的哈斯图为图34―2 b A 236122436 〈A整除〉的哈斯图为图34―3 c A 1212 〈A整除〉的哈斯图为图34―4 定义34―2 设〈A≤〉是一偏序集合B是A的子集 a 元素b∈B是B的最大元素如果对每一元素x∈Bx≤b b 元素b∈B是B的最小元素如果对每一元素x∈Bb≤x 例3考虑在偏序整除下整数1到6的集合其哈斯图为图34―5 a 如果B 1236 那么1是B的最小元素6是B的最大元素 b 如果B 23 因为2和3互相不能整除那么B没有最小元素和最大元素 c 如果B 4 那么4是B的最大元素也是B的最小元素定理34―1 设〈A≤〉是一偏序集合且B A如果B有最大最小元素那么它是唯一的证假设a和b都是B的最大元素那么a≤b和b≤a从≤的反对称性得到a b当a和b都是B的最小元素时证明是类似的定义34―3设〈A≤〉是一偏序集合B是A的子集 a如果b∈B且B中不存在元素x使b≠x且b≤x那么元素b∈B叫做B的极大元素b 如果b∈B且B中不存在元素x使b≠x且x≤b那么元素b∈B叫做B的极小元素定义34―4设〈A≤〉是一偏序集合B是A的子集a 如果对每一b∈Bb≤a那么元素a∈A叫做B的上界如果对每一b∈Ba≤b那么元素a∈A叫做B的下界 b 如果a是一上界并且对每一B的上界a′有a≤a′那么元素a∈A叫做B的最小上界记为lub如果a是一下界并且对每一B的下界a′有a′≤a那么元素a∈A叫做B的最大下界记为glb 例34-4 a 考虑偏序集合〈〈11〉〈10〉〈01〉〈00〉≤〉这里≤按〈 ab〉≤〈cd〉a≤c∧b≤d 规定其哈斯图如图34―6 如果B 〈10〉那么〈10〉是B的最小和最大元素也是B的极大和极小元素B的上界是〈10〉和〈11〉〈10〉是最小上界B的下界是〈00〉和〈10〉〈10〉是最大下界 b 考虑偏序集合〈I≤〉设B 2i｜i∈N那么B既没有最大元素和极大元素也没有上界和最小上界B的最小元素和极小元素是0B的下界集合是 i｜i∈I∧i≤0 0是最大下界 c 考虑在偏序集合〈 256101530 整除〉其哈斯图如图34―7设B是全集合 256101530 那么2和5都是B的极小元素但B没有最小元素集合B没有下界所以没有最大下界元素30是B的最大元素极大元素上界最小上界定理34―2 如果〈A≤〉是非空有限的偏序集合则A的极小大元素常存在最大下界和最小上界也可能存在或不存在但如果它们存在则是唯一的定理34―3 设〈A≤〉是偏序集合且B A 如果B的最小上界最大下界存在那么是唯一的下述定理描述了存在于诸特异元素之间的某些关系定理34―4 设〈A≤〉是偏序集合B是A的子集 a 如果b是B的最大元素那么b是B的极大元素 b 如果b是B的最大元素那么b是B的lub c 如果b是B的一个上界且b∈B那么b是B的最大元素证明可由最大元素极大元素和lub的定义直接得出故略去另外读者不难给出表达最小元素极小元素和glb间关系的定理 342 拟序集合定义34―5如果集合A上的二元关系R是传递的和反自反的那么R叫做A上的拟序〈AR〉称为拟序集合常借用符号＜表示拟序拟序是反对称的虽然定义中没有明确指出但容易证明这一点因为如果xRy和yRx由R的传递性得xRx但这与R的反自反性矛盾所以xRy∧yRx常假于是xRy∧yRx→x y常真即R是反对称的例34-5 a 实数集合中的＜是拟序关系 b 集合族中的真包含是拟序关系拟序集合和偏序集合是紧密相关的唯一区别是相等关系E下述定理将说明这一点定理34―5在集合A上 a 如果R是一拟序那么rR R∪E是偏序 b 如果R是一偏序那么R-E是一拟序 343线序集合和良序集合如果≤是一偏序或a≤b或b≤a我们说a和b是可比较的偏序集合中的元素不一定都可比较所以叫偏序下面介绍的都是可比较的情况定义34―6在偏序集合〈A≤〉中如果每一ab∈A或者a≤b或者b≤a那么≤叫做A上的线序或全序这时的序偶〈A≤〉叫做线序集合或链例34-6 a P a ab abc 〈P〉是线序集合其哈斯图如图34―8所示 b 〈I≤〉是线序集合其哈斯图不完全如图34―9所示 c 设S是区间套的集合〔0a ｜a∈R 则〈S〉是线序集合 d 〈 1236 整除〉不是线序集合如果A是多于一个元素的集合那么〈ρ A 〉不是线序集合定义34―7如果A上的二元关系R是一线序且A的每一非空子集都有一最小元素那么R叫做A上的良序序偶〈AR〉叫做良序集合定理34―6〈N≤〉是良序集合证我们必须证明N的每一非空子集S在关系≤之下都有一最小元素因为S非空所以在S中可以取一个数n显然S中所有不大于n的数形成非空集T S如果T有最小数那么这最小数就是S中的最小数但从0到n只有n1个自然数于是T中所含的数最多是n1个所以T有最小数因此定理成立例34-8 a 每一有限线序集合是良序的 b 线序集合〈I≤〉不是良序集合因为I的某些子集诸如I自身不包含最小元素 c 关系≤是实数R的线序但不是良序例如子集A 01〕无最小元素如果A中的a是最小元素那么也在A中而≤a且不相等这与假设a是线序关系≤下A的最小元素矛盾2 应用N上的良序定义出Nn上的良序例如n 2时N2上的次序关系可如下定义〈ab〉〈cd〉a＜c∨ a c∧b d 〈N2〉是良序集合关系严格小于可如下定义〈ab〉＜〈cd〉〈ab〉≤〈cd〉∧〈ab〉≠〈cd〉类似地应用I上的线序能定义出线序集合〈In≤〉 3 应用字母表∑上的线序可定义出∑上的通常叫词典序的线序定义34―8 设∑是一有限字母表指定了字母表序线序如果xy∈∑ a x是y的词头或 b x zu和y zv这里z∈∑是x和y的最长公共词头且在字母表序中u的第一个字符前于v的第一个字符那么x≤y≤叫做词典序4 由于〈N〉和有限线序集合都是良序集合可应用它们定义出∑上的一个良序通常叫标准序定义34―9设∑是一有限字母表指定了字母表序‖x‖表示x∈∑的长度如果xy∈∑ a ‖x‖＜‖y‖或b ‖x‖‖y‖且在∑的词典序中x前于y那么x≤y ≤叫做标准序不论在词典序和标准序下∑的每一元素都有直接后继者设∑ abc 且a≤b≤cx∈∑在标准序下 xa和xb的直接后继者分别是xb和xc xc的直接后继者是ya这里y是x的直接后继者在词典序下x的直接后继者是xa在标准序下 xb和xc的直接前趋分别是xa和xb xa的直接前趋是yc这里y是x的直接前趋在词典序下 xa的直接前趋是x非a结尾的串都无直接前趋例如babaab但有无限个前趋 345 数学归纳法的推广前章我们把数学归纳法第一第二原理看作是自然数域上的一个推理规则本小节我们把它推广到一般的良序集合对任一个自然数n我们先取0如果n≠0取0的后继者1如果n≠1再取1的后继者2如此进行下去最终会得出n 给定一个良序集合如果对它的任一元素x我们先取该集合的最小元素m0如果x≠m0取m0的后继者m1如果x≠m1再取m1的后继者m2如此以往最终会得出x那么就称这样的良序集合是像自然数的例 8 1 设∑ ab 良序集合〈∑标准序〉是像自然数的因为定长的串的个数有限给定任一个串x在x之前的串的个数有限所以从∧开始反复取后继者终可得出x 2 良序集合〈N×N≤〉不像自然数这里≤按上一小节规定因为有许多元素没有直接前趋例如〈50〉就是这样因而有无限个元素前于〈50〉所以从〈00〉开始反复地取后继者不可能取得〈50〉像自然数的良序集合可以应用数学归纳法第一原理因为第一原理是建立在后继运算上而这种良序集合的每一元素都可通过重复地取后继者得到设m0是该良序集合〈S≤〉的最小元素S x 是元素x的后继者则推理规则如下对不像自然数的良序集合不能应用数学归纳法第一原理因为这种良序集合的有些元素不能由后继运算得到但对它可应用数学归纳法第二原理第二原理是建立在良序集合上的适用于一切良序集合设〈S≤〉是良序集合＜表示≤-E 即x＜y表示x≤y且x≠y 则推理规则如下下面证明良序集合上这个推理规则是有效的假设我们能证明前提例34-10〈Q≤〉是线序集合现说明在此线序集合中第二原理不是有效推理规则设谓词P x 表示x小于或等于5 i 当x≤5时 y〔y＜x→P y 〕是真P x 也真所以是真综合 i 和 ii 得在论述域Q上 x 〔 y y＜x→P y →P x 〕是真但结论 x P x 是假这说明第二原理不能应用于线序集合〈Q≤〉 35 等价关系和划分 351 等价关系二元关系的另一重要类型是等价关系其定义如下定义35―1 如果集合A上的二元关系R是自反的对称的和传递的那么称R是等价关系设R是A上的等价关系abc是A的任意元素如果aRb 即〈ab〉∈R 通常我们记作a～b读做a等价于b 定义35―2 设k是一正整数而ab∈I如果对某整数ma-b m·k那么a和b是模k等价写成a≡b modk 整数k叫做等价的模数定理35―1模k等价是任何集合A I上的等价关系证如果A 例35-1 c 已指出它是等价关系如果A≠则 i 自反的因为对任一aa-a 0·k得出a≡a modk ii 对称的因为a≡b mod k 时存在某m∈I使a-b m·k于是b-a-m·k 因此 b≡a mod k iii 传递的设a≡b mod k 和b≡c mod k 那么存在m1m2∈I 使a-b m1k和b-c m2·k 将两等式两边相加得a-c m1m2 ·k所以a≡c mod k 例1 a 同学集合A abcdefgA中的关系R是住在同一房间这是等价关系因为 i 任一个人和自己同住一间具有自反性 ii 若甲和乙同住一间则乙和甲也同住一间具有对称性 iii 若甲和乙同住一间乙和丙同住一间则甲和丙也同住一间具有传递性现假设a和b同住一间def同住一间c住一间则 R 〈aa〉〈ab〉〈ba〉〈bb〉〈cc〉〈dd〉〈ee〉〈ff〉〈de〉〈ed〉〈ef〉〈fe〉〈df〉〈fd〉其有向图如图35―1所示 b 数中的相等关系集合中的相等关系命题演算中的关系等都是等价关系 c 空集合中的二元关系R是等价关系因为i x x∈→xRx ii x y〔x∈∧y∈∧xRy→yRx〕iii x y z〔x∈∧y∈∧z∈∧xRy∧yRz→xRz〕都无义地真所以R是等价关系集合A上的全域关系R A×A是等价关系模数等价是整数域或其子集上的等价关系并且是等价关系中极为重要的一类定理 35-1 模k等价是任何集合A I上的等价关系证如果A 例35-1 3 已指出它是等价关系如果A≠则 i 自反的因为对任一aa-a 0·k得出a≡a mod k ii 对称的因为a≡b modk 时存在某m∈I使a-b m·k于是b-a -m·k因此b≡amodk iii 传递的设a≡b modk 和b≡c modk 那么存在m1m2∈I使a-b m1k和b-c m2·k将两等式两边相加得a-c m1m2 ·k所以a≡c modk 例35-2 a 若R是I上模4等价关系则〔0〕4 -8-4048 〔1〕4 -7-3159 〔2〕4 -6-22610 〔3〕4 -5-13711 b 若R是I上模2等价关系则〔0〕2 -4-2024 〔1〕2 -3-1135 每一集合中的数相互等价 c 时钟是按模12方式记数的设备13点钟和1点钟有相同的记数定义35―3 设R是集合A上等价关系对每一a∈Aa关于R的等价类是集合 x｜xRa 记为〔a〕R简记为〔a〕称a为等价类〔a〕的表示元素如果等价类个数有限则R的不同等价类的个数叫做R的秩否则秩是无限的对每一a∈A等价类〔a〕R非空因为a∈〔a〕R 例 35-3 1 如图35―2设A abcdef R 〈aa〉〈bb〉〈cc〉〈ab〉〈ba〉〈ac〉〈ca〉〈bc〉〈cb〉〈dd〉〈ee〉〈de〉〈ed〉〈ff〉则等价关系R的等价类如下〔a〕〔b〕〔c〕 abc 〔d〕〔e〕 de 〔f〕 f等价关系R的秩是3 2 I上模4等价的等价类是〔0〕4〔1〕4〔2〕4〔3〕4 参看例2 a I上模2等价的等价类是〔0〕2 〔1〕2 参看例2 b3 集合A上相等关系的秩等于A的元素个数定理35―2 设R是非空集合A上的等价关系aRb 当且仅当〔a〕〔b〕证充分性因为a∈〔a〕〔b〕即a∈〔b〕所以aRb 定理35―3设R是集合A上的等价关系则对所有ab∈A或者〔a〕〔b〕或者〔a〕∩〔b〕证如果A 断言无义地真现设A≠若〔a〕∩〔b〕≠则存在某元素c∈〔a〕和c∈〔b〕根据定理35―2得〔a〕〔c〕〔b〕又因〔a〕和〔b〕都非空〔a〕∩〔b〕和〔a〕〔b〕不能兼得因而定理得证定义35― 4 给定非空集合A和非空集合族πA1A2Am 如果那么称集合族π是A的覆盖定理35―4设R是集合A上的等价关系则证先证定理35―5设R1和R2是集合A上的等价关系那么R1 R2当且仅当〔a〕R1｜a∈A 〔a〕R2｜a∈A 证必要性因为R1 R2所以对任意a∈A有〔a〕R1 x｜xR1a x｜xR2a。
岩石力学-离散元法
数据结构一般是静态的，所谓静态数据 结构，就是一个数据结构存在期间，其结构 的特性不变。
动态数据结构是使用线性链表（一组任 意的存贮单元）存放数据元素的一种结构， 它可以动态地增加或减少结点。
3、参数的选择和本构模型 ①阻尼
岩块在运动中一般不发生弹跳，这是由于运动 时其动能转化成热能而耗散掉的缘故。
②力和位移的计算循环
离散元法可采用动态松弛法求Biblioteka ，其基本运动 方程为
mu(t) cu(t) ku(t) f (t) （5－129）
力－位移的关系
力F
运动方程 u F m
位移 u
力边界条件
计算循环
位移边界条件
③分格检索
分格检索实际上是块体接触点的识别问 题，这是离散元法所特有的一个问题。
A、粘性阻尼 其自由振动微分方程为
mu c u ku 0
a、欠阻尼振动系统 b、临界兹尼振动系统 c、过阻尼振动系统
（5－143）
第五节 离散元法
离散元法是一种显式求解的数值方法。 1、离散元法的基本方程 ①物理方程－力和位移的关系
假定块体之间的法向力 Fn 正比于它们之间位 移的“叠合”Un ，即
Fn KnUn （5－124） 式中 Kn 为法向刚度系数 ② 运动方程－牛顿第二定律
2、离散元法的计算实施 ①动态松弛法
动态松弛法是把非线性静力学问题转化为动力 学问题求解的一种数值方法。
离散数学 二元关系 PPT课件
▪ 常见的几种特殊的二元关系
▪≤ ≥ ＜ ＞ ＝ ▪| ▪ 集合之间的关系 ： ＝ ≠
2020/7/15
20
计算机科学学院 刘芳
7.2.2 二元关系的表示
1.集合表示法
2.关系矩阵(matrix of relation)
▪ 设A＝{a1,a2,…,am} ，B＝{b1,b2,…,bn}，R是A到B的一个二
所以， （A∩B）×C＝（A×C）∩（B×C）成立。
2020/7/15
11
计算机科学学院 刘芳
7.1.3 有序 n 元组和 n 阶笛卡尔积
定义：
▪ n个元素x1，x2，…，xn组成的有序序列,记做：
＜x1,x2,…,xn＞
▪ 称为n重组（n元序偶、n元组）。
约定:
▪ ＜x1，x2，…, xn-1, xn＞= ＜<x1，x2，… ,xn-1 >,xn＞
1 1 1 0 1 0 0 0 0 0 0 0 0 1 0 0
2020/7/15
28
计算机科学学院 刘芳
7.2.2 二元关系的表示
关系的表示方法
▪ 关系R的集合表达式 ▪ 关系矩阵MR ▪ 关系图GR
三者均可以唯一相互确定。
2020/7/15
29
计算机科学学院 刘芳
7.3 关系的运算
7.3.1 关系的定义域、值域 和 域
例:
▪ (1) R＝{<x,y> | x,yN, x＋y＜3}
＝{<0,0>, <0,1>, <0,2>, <1,0>, <1,1>, <2,0>}
▪ (2) C＝{<x,y> | x,yR, x2＋y2＝1}
离散数学讲义ppt课件
课程概况
教材：
《离散数学（第三版）》，耿素云等编著 清华大学出版社，2004年3月
参考书：
（1） 《离散数学(第二版)》及其配套参考书《离散 数学题解》作者：屈婉玲，耿素云，张立昂 清华大学出版社
（2） 《离散数学》焦占亚主编 电子工业出版社 2005年1月
2
课程概况
选修课/必修课：选修 周学时：3（学时） 上课周：1－16周 总学时：48（学时）
3
课程内容及学时安排
第一篇 数理逻辑（14学时）
第一章 命题逻辑（8） 第二章 谓词逻辑（6）
第二篇 集合论（12学时）
第三章 集合（4） 第四章 二元关系与函数（8）
第四篇 图论（14学时）
第七章 图论（8） 第八章 一些特殊图（4） 第九章 树 （2）
4
课程考核
第四篇 代数系统（8学时）
第5、6章 图论（8）
所以，伊勒克持拉既知道并且又不知道这个人是她的 哥哥。
20
NO.3 M：著名的理发师悖论是伯特纳德·罗素提出的。一个理发 师的招牌上写着： 告示：城里所有不自己刮脸的男人都由我给他们刮脸，我 也只给这些人刮脸。 M：谁给这位理发师刮脸呢？ M：如果他自己刮脸，那他就属于自己刮脸的那类人。但 是，他的招牌说明他不给这类人刮脸，因此他不能自己来 刮。 M：如果另外一个人来给他刮脸，那他就是不自己刮脸的 人。但是，他的招牌说他要给所有这类人刮脸。因此其他 任何人也不能给他刮脸。看来，没有任何人能给这位理发 师刮脸了！
P
Q
PQ
P
0
0
0
1
0
1
0
1
1
0
0
0
1
1
1
离散数学完整版课件全套ppt教学教程最全整套电子讲义幻灯片(最新)
1.1 命题及联结词
运算符“析取” 与汉语的“或”几乎一致但有 区别：哪些老师讲离散数学？有人回答如下：
(16)“讲离散数学的老师是杨老师或吴老师”， 分解为
“讲离散数学的老师是杨老师”或 “讲离散数学的老师是吴老师”， 这两个原子命题有可能都是对的， 这种“或”称为“可同时为真的或”， 或简称为“可兼或”。 这种“或”表示可表 示为“析取”
1.1 命题及联结词
定义1.4条件:当p是1 ,q是0时,pq为0,即10 为0,其他情况为1。
逻辑运算符“如果…那么”， 如老妈说：“如果期终考了年级前10 名，那么奖励1000元”。 p:期终考了年级前10名 q:奖励1000元 则上面的语句表示为pq。 先考虑值为0即假的情况： 当p为1即“期终考了年级前10”， 且q为0即“没有奖励1000元” 这时老妈的话是假话空话，
这个例题有点不正点！ “郎才当且仅当女貌”，
可以表示为“郎才女貌”
1.2命题公式
对错明确的陈述语句称为命题，其真值t/f 0/1 C运算：加+、减-、乘x、除/、余数%， 命题逻辑：合、析、否定、条件、双条件(版) C语言中用变量x表示某些数，如x*x+x+10是表达式，
命题逻辑中用变量p,q,r表示任意命题，由命题常元与 此类变量所构成表达式，称为“命题公式”。
无论p/q取何值，这两个公式的值，与前面各例 不同，此表是将运算结果写在联结词的下方！
1.3 等值式
定义1.3.1等值： 对于合法的命题公式A、B， 若无论其中的命题变元取何值，A 、B值总相等， 称为两个公式等值，记为AB (边播边板)
目的:
1．掌握离散数学五大核心内容（集合论、数 理逻辑、代数结构、图论、组合数学）的基本概 念、基本理论、基本方法，训练提高学生的概括 抽象能力、逻辑思维能力、归纳构造能力，培养 学生严谨、完整、规范的科学态度和学习思维习 惯。
离散元ppt课件
一 历史由来及研究现状
离散单元法 的研究现状
离散元理论研究的发展
近30年来，离散元法的应用领域在不断扩大，它自身的内涵也发 生了变化，以致于目前很难对离散元法给出一个严格的定义。下 面，我们从离散元法的离散模型特点及便于甄别与其它数值计算 方法的关系的角度给予离散元法一个比较宽松的定义。
一 历史由来及研究现状
一 历史由来及研究现状
产生背景
散粒岩土材料在自然界中普遍存在
从本质上讲，岩土材料都是由离散的、尺寸不一、形状各异的颗 粒或块体组成的，例如，土就是松散颗粒的堆积物，同样，天然 岩体也是由被结构面切割而成的大小不一、形态各异的岩石块体 所组成。散粒岩土材料的力学特性有着重要的工程应用，如泥砂 的沉淀，土堤、土（岩）坡、铁路道渣等的稳定性研究，散粒岩 土材料的力学特性研究是岩土力学中最基本的、也是最重要的问 题之一。
离散单元法 的研究现状
离散元理论研究的发展
• 数值方法通常将实际具有无限自由度的介质近似为具有有限 自由度的离散体(或网络)的计算模型(有限离散模型)进行计算。 有限离散模型具有三个要素：单元(或网络)、节点和节点间 的关联。
一 历史由来及研究现状
离散单元法 的研究现状
离散元理论研究的发展
• 离散元单元的形状有形形色色，但它只有一个基本节点(取单 元的形心点)，是一种物理元(physicalelement)．这种单元与有 限元法、边界元法等数值方法采用的由一组基本节点联成的 单元(一般称为网络元，mesh element)相比有明显的不同。
一 历史由来及研究现状
早期的离散单元法
Cundall称之为“Distinct Element Method”，随着该方法的推广， 有的学者称其为“Discrete Element Method”，缩写形式均为 DEM。 最初，离散元的研究对象主要是岩石等非连续介质的力学行为， 它的基本思想是把不连续体分离为刚性元素的集合，使各个刚性 元素满足运动方程，用时步迭代的方法求解各刚性元素的运动方 程，继而求得不连续体的整体运动形态。离散元法允许单元间的 相对运动，不一定要满足位移连续和变形谐调条件，计算速度快， 所需存储空间小，尤其适合求解大位移和非线性问题。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
二基本原理球形颗粒元离散元法
接触模型
两个处于接触颗粒单位法向和切向向量：
单位法向向量 n Ri / Ri
单位切向向量 t Vij Vij • n n Vij Vij • n n
单位切向量之所以通过两个颗粒的相对速度来计算， 是因为接触力与粘性阻尼力的方向与相对速度的方 向相同。
二基本原理球形颗粒元离散元法
二 基本原理-球形颗粒元离散元法
运动描述
处于一个理想散体中的任意一个颗粒,具有6个自由 度,3个平动自由度与三个转动自由度,可通过Newton 第二定律分别描述。
二基本原理球形颗粒元离散元法
运动描述
平动方程：
mi
dVi dt
ki
Fc,ij Fd ,ij
j 1
mi g
式中，mi与 Vi 分别为颗粒 i 的质量和速度。t 为时间，
,
Fct ,ij s Fcn,ij Fct ,ij s Fcn,ij
3
Fct,ij
δt 为颗粒 i 与 j 间的累积切向位移矢量
δt N δt N1 Vij n n t
二 基本原理-球形颗粒元离散元法
接触模型
切向接触力模型 —阻尼力 ：
Fdt,ij ct Vij n n
式中，ct 为切向粘性接触阻尼系数。
二 基本原理-球形颗粒元离散元法
计算模型总结
运动方程
mi g 为颗粒的重力，Fc,ij 与 Fd ,ij 分别为颗粒 i 与 j 的
接触力与粘性接触阻尼力，ki 为所有与颗粒接触的颗
粒总数。
二基本原理球形颗粒元离散元法
运动描述
接触力的分解：
颗粒 i与 j间的接触力可分解为法向与切向接触力，即
Fc,ij Fcn,ij Fct,ij
同理，粘性接触阻尼力也可分解为法向与切向分量形 式，即
式中，
k
为法向弹簧刚度。
n
Hale Waihona Puke 二 基本原理-球形颗粒元离散元法
接触模型
法向接触力计算模型 —法向粘性接触阻尼力：
Fdn,ij cn Vij • n n
式中，cn为法向粘性接触阻
尼系数。
二 基本原理-球形颗粒元离散元法
接触模型
切向接触力计算模型 —综述：
处于接触中的两个颗粒的切向作用，从本质上讲，是 一种摩擦行为，按照摩擦机理，摩擦力包括：滑动摩 擦、滚动摩擦与静摩擦，其中滑动摩擦与静摩擦属于 切向摩擦力；滚动摩擦是由于法向接触应力的不均匀 分布产生的。介绍两个切向接触力模型： Coulomb准则 Mindlin与Deresiewicz切向接触力模型
Fd ,ij Fdn,ij Fdt,ij
二基本原理球形颗粒元离散元法
运动描述
接触力产生的力矩：
颗粒间的接触力作用在两个颗粒的接触点上，而不是作 用在颗粒的中心，所以这些接触力（除法向接触力Fcn,ij 外）将会对颗粒产生力矩 Ti ，
Ti Ri Fct,ij Fdt,ij
式中，Ri为从颗粒i 的质心指向接触点的矢量，其幅值为
对于理想散体颗粒（无粘连）：采用Hertz理论描述法向作用，而采 用Mindlin与Deresiewicz理论描述切向作用； 对于存在粘连的散体颗粒：法向接触力根据在Hertz理论基础上考虑 粘连力的JKR（Johnson-Kendall-Roberts）理论确定，切向接触力增量 则根据把Savkoor和Briggs理论与Mindlin和Deresiewicz理论相结合形 成的理论确定。
二 基本原理-球形颗粒元离散元法
岩土工程研究所
刘军
二 基本原理
根据离散化模型中所采用的单元种类分别介绍离 散元法的基本原理：
颗粒元
• 二维圆盘单元 • 三维圆球单元
块体元
• 多边形单元 • 多面体单元
二 基本原理-球形颗粒元离散元法
基本假设
假定速度和加速度在每个时间步长内为常量 ； 选取的时间步长应该足够小以至于在单个时间步长 内扰动的传播不会超过当前与之相邻的粒子 。
R（i 颗粒的半径）。
二基本原理球形颗粒元离散元法
运动描述
转动方程：
转动方程可以表示为
Ii
dωi dt
ki
Ti
j 1
式中，Ii与ωi 分别为颗粒 i 的转动惯量与角速度，对于
球形颗粒 Ii为
Ii
2 5
mi Ri2
二 基本原理-球形颗粒元离散元法
接触模型
综述：
关于接触力的计算模型已有大量的研究成果，目前仍旧是 一个活跃的研究领域，特别是对于切向力的计算方法。
二 基本原理-球形颗粒元离散元法
接触模型
切向接触力计算模型 —Coulomb准则：
在离散元模拟中，一般用Coulomb准则这种简单的形 式描述，静摩擦的详细刻画需要涉及切向位移甚至可能 要考虑时间依赖效应。
Fct ,ij
Fct ,ij ,
s Fcn,ij
,
Fct ,ij s Fcn,ij Fct ,ij s Fcn,ij
式中，μs 为静摩擦系数，切向摩擦力的方向为与相对滑
动的趋势相反。
二 基本原理-球形颗粒元离散元法
接触模型
切向接触力模型 —Mindlin与Deresiewicz模型 ：
式中
3
Fct,ij
μs Fcn,ij δt
1
1
min
δt , δt,max δt,max
2
δt
2v
δt,max μs 2 1 v δn
接触模型
两个处于接触颗粒接触点的相对速度：
Vij V j Vi ω j R j ωi Ri
法向相对速度为
Vn,ij Vij • n n
切向相对速度为
Vt,ij Vij Vij • n n
Vt,ij Vij n n
或者写为
二 基本原理-球形颗粒元离散元法
接触模型
法向接触力计算模型 —Hertz模型：
式中
Fcn,ij
4 3
E*
R*
3/ n
2n
E*
E 2(1 v2 )
R* 1 1 Ri R j
δn为颗粒i与j接触时的侵入深度
n Ri R j R j Ri
二基本原理球形颗粒元离散元法
接触模型
法向接触力计算模型—
Cundall模型：
Fcn,ij kn nn
mi
dVi dt
ki
Fc,ij Fd ,ij
j 1
mi g
Ii
dωi dt
ki
Ti
j 1
接触力的计算
• 法向接触力
Fcn,ij
4 3
E*
R*
3/ n
2n
Fcn,ij kn nn
Fdn,ij cn Vij • n n
• 切向接触力
Fct ,ij
Fct ,ij , s Fcn,ij