李洪興
(1. 北京師范大學(xué)珠海校區(qū) 應(yīng)用數(shù)學(xué)與交叉科學(xué)研究中心, 廣東 珠海 519085; 2. 大連理工大學(xué) 控制科學(xué)與工程學(xué)院, 遼寧 大連 116024)
為了說(shuō)明問(wèn)題的實(shí)質(zhì),當(dāng)取U=R時(shí),那么這個(gè)函數(shù)可以看作一個(gè)單位脈沖波,見(jiàn)圖1.
圖 1 單位脈沖波
特別地,當(dāng)A=[a,b]?R時(shí),這個(gè)函數(shù)可以視為一個(gè)單位矩形波,見(jiàn)圖2.
圖 2 單位矩形波
(A,B)→A∪B=∪(A,B)=
{u∈U|(u∈A)∨(u∈B)};
(A,B)→A∩B=∩(A,B)=
{u∈U|(u∈A)∧(u∈B)};
A→Ac=c(A)=UA=
{u∈U|u?A}.
χA∪B(u)=χA(u)∨χB(u), ?u∈U;
χA∩B(u)=χA(u)∧χB(u), ?u∈U;
χAc(u)=1-χA(u), ?u∈U,
其中,∨=max,∧=min,即
a∨b=max{a,b}, ?a,b∈R,
a∧b=min{a,b}, ?a,b∈R.
A×B:={(a,b)|a∈A,b∈B}
叫做A與B的直積.特別地,有
U×V={(u,v)|u∈U,v∈V}.
不難驗(yàn)證結(jié)果
χA×B(u,v)=χA(u)∧χB(v)=
χA(u)·χB(v), ?(u,v)∈U×V.
現(xiàn)在引入符號(hào)
Ch(U)={χ|χ:U→[0,1]},
在集合Ch(U)上定義2個(gè)二元代數(shù)運(yùn)算∨、∧和一個(gè)一元代數(shù)運(yùn)算“c”如下:
∨:Ch(U)×Ch(U)→Ch(U),
(χ1,χ2)→χ1∨χ2=∨(χ1,χ2),
(χ1∨χ2)(u)=χ1(u)∨χ2(u), ?u∈U;
∧:Ch(U)×Ch(U)→Ch(U),
(χ1,χ2)→χ1∧χ2=∧(χ1,χ2),
(χ1∧χ2)(u)=χ1(u)∧χ2(u), ?u∈U;
c:Ch(U)→Ch(U),χ→χc=c(χ),
χc(u)=1-χ(u), ?u∈U.
這樣就構(gòu)成一個(gè)代數(shù)系統(tǒng)(Ch(U),∨,∧,c).
不防證明(1)式.事實(shí)上,考慮映射
[(u∈A)∧(u?B)]∨
[(u?A)∧(u∈B)]
為真.無(wú)妨假定(u∈A)∧(u?B)為真,那么有
(χA(u)=1)∧(χB(u)=0),
A={u∈U|χ(u)=1}.
根據(jù)這個(gè)集合A,又可以得到一個(gè)函數(shù)χA∈Ch(U)如下:
χA:U→{0,1},
接下來(lái),證明χA=χ.實(shí)際上,對(duì)于任意的元素u∈U,顯然有
χA(u)=1?u∈A?χ(u)=1.
由此可見(jiàn)χA=χ,即
f(A)=χA=χ,
所以f是個(gè)滿射,從而f是雙射.
f(A∪B)=χA∪B,
f(A∩B)=χA∩B,
f(Ac)=χAc.
只需觀察下面的事實(shí):對(duì)于任意的元素u∈U,有:
[f(A∪B)](u)=χA∪B(u)=
χA(u)∨χB(u)=
[f(A)](u)∨[f(B)](u)=
{[f(A)]∨[f(B)]}(u),
[f(A∩B)](u)=χA∩B(u)=
χA(u)∧χB(u)=
[f(A)](u)∧[f(B)](u)=
{[f(A)]∧[f(B)]}(u),
[f(Ac)](u)=χAc(u)=1-χA(u)=
1-[f(A)](u)=[f(A)]c(u).
可見(jiàn)雙射f確實(shí)保持代數(shù)運(yùn)算,從而f是個(gè)同構(gòu)映射.
因?yàn)?/p>
所以論域U中的一個(gè)集合就是定義在U上的一個(gè)波函數(shù),反之亦然;這就叫做波集二象性.
(2)
事實(shí)上,對(duì)于任意的u∈U有
χA(u)=1?u∈A?
1+0=1.
這種情況足以啟發(fā)產(chǎn)生一種新的思想:在實(shí)際的宏觀物理世界,任何一個(gè)物體可以視為一個(gè)質(zhì)點(diǎn),而該物體一定由大量的微觀粒子所構(gòu)成.如果宏觀質(zhì)點(diǎn)運(yùn)動(dòng)也具有波動(dòng)性,而微觀粒子運(yùn)動(dòng)是具有波粒二象性的,那么宏觀質(zhì)點(diǎn)運(yùn)動(dòng)的波函數(shù)一定是有若干微觀粒子運(yùn)動(dòng)的波函數(shù)構(gòu)成.可喜的是,文獻(xiàn)[1]已經(jīng)證明了宏觀質(zhì)點(diǎn)運(yùn)動(dòng)也具有波動(dòng)性,換言之,宏觀質(zhì)點(diǎn)運(yùn)動(dòng)是具有波質(zhì)二象性的.
從一類具有單輸入單輸出的開(kāi)環(huán)系統(tǒng)談起,見(jiàn)圖3,其中S代表一個(gè)系統(tǒng),x表示在論域
X=[a,b]?R
中取值的輸入變量,y表示在論域
Y=[c,d]?R
中取值的輸出變量.
圖 3 單輸入單輸出的開(kāi)環(huán)系統(tǒng)
如果S是一個(gè)復(fù)雜的不確定系統(tǒng),那么對(duì)這個(gè)系統(tǒng)采用常規(guī)的機(jī)制建模法是困難的,甚至是不可能的.對(duì)于這種情況,可以通過(guò)實(shí)驗(yàn)來(lái)得到一組離散的數(shù)據(jù),由它近似地表示該系統(tǒng)的輸入輸出關(guān)系[2].這組數(shù)據(jù)用符號(hào)IOD來(lái)表達(dá)如下
IOD?{(xi,yi)|i=0,1,…,n}?X×Y,
這里的輸入數(shù)據(jù)和輸出數(shù)據(jù)可以分別表示為:
X0?{xi|i=0,1,…,n},
Y0?{yi|i=0,1,…,n},
a=x0 實(shí)際上,IOD可以看作一個(gè)離散的函數(shù) g0:X0→Y0,xi→g(xi)?yi, i=0,1,…,n. 圖 4 系統(tǒng)的離散響應(yīng) 現(xiàn)在擴(kuò)展函數(shù)g0:X0→Y0的定義如下: g:X→Y, 盡管映射g:X→Y使得系統(tǒng)S對(duì)于每一個(gè)x∈X都有響應(yīng),即 i=0,1,…,n. 如果注意到g(x)=0(?x∈XX0),可見(jiàn)系統(tǒng)S在XX0中的響應(yīng)幾乎無(wú)意義. 于是應(yīng)該思考一個(gè)問(wèn)題:如何利用已知的映射g0:X0→Y0來(lái)構(gòu)造一個(gè)實(shí)用的映射f:X→Y使得系統(tǒng)S對(duì)于X中的每一個(gè)元素都具有實(shí)用價(jià)值的響應(yīng),并且滿足如下的自然條件: f(xi)=g0(xi), ?i∈{0,1,…,n}. 當(dāng)然,眾所周知,插值是解決這個(gè)問(wèn)題的一個(gè)方法;不過(guò)插值方法不具有系統(tǒng)的觀念,也不具有集合運(yùn)算、邏輯和推理的意義[3].下面將采用集合、邏輯和推理的思想來(lái)解決該問(wèn)題. 首先,不難理解,通過(guò)實(shí)驗(yàn)得到的數(shù)據(jù)xi、yi,i=0,1,…,n都具有誤差,于是便有下面的表達(dá)式: xi±δi,yi±εi,δi≥0,εi≥0, i=0,1,…,n, 或者可以寫為下列集合形式(區(qū)間是集合): xi∈[xi-δi,xi+δi], yi∈[yi-εi,yi+εi], i=0,1,…,n. 這意味著數(shù)據(jù)被一步步擴(kuò)展為下面的集合形式[4]: IOD={(xi,yi)|i=0,1,…,n}? {({xi},{yi})|i=0,1,…,n}? {([xi-δi,xi+δi],[yi-εi,yi+εi])| i=0,1,…,n}. 注意每個(gè)誤差區(qū)間的特征函數(shù)形式: χ[xi-δi,xi+δi](x)= i=0,1,…,n. 基于上面討論的波集二象性,每個(gè)單位脈沖波χ{xi}(x)擴(kuò)展為單位矩形波χ[xi-δi,xi+δi](x),見(jiàn)圖5. 圖 5 單位脈沖波到單位矩形波 從單位脈沖波χ{xi}(x)到單位矩形波χ[xi-δi,xi+δi](x),實(shí)際上利用數(shù)據(jù)的誤差獲得了更多的信息.現(xiàn)在計(jì)算來(lái)自誤差區(qū)間[xi-δi,xi+δi]的一種類型的信息量如下 (3) 如圖6所示,系統(tǒng)在位于一些間隔(interspaces)上的輸入量也是沒(méi)有響應(yīng)的. 圖 6 系統(tǒng)在間隔上沒(méi)有響應(yīng) 為了解決這個(gè)問(wèn)題,思想是在保持信息量不變的情況下,考慮使用另外一些形式的波函數(shù)μAi(x),當(dāng)然它們也是某種類型的集合Ai的波函數(shù),見(jiàn)圖7. 所謂保持信息量不變,就是要求滿足條件 圖 7 作為某種集合Ai的波函數(shù)μAi(x) 如果取下面形式的2個(gè)數(shù): 那么就可以得到如下的波函數(shù): 不難驗(yàn)證,上述形式的波函數(shù)已經(jīng)滿足條件: 同時(shí)不難知道這些波函數(shù)適合下面的條件: μAi(X)=[0,1], ?i∈{0,1,…,n}, 這里μAi(X)是函數(shù)μAi關(guān)于論域X的象集,它們不同于χ[xi-δi,xi+δi](X)={0,1}.因此,從波函數(shù)的角度講,是對(duì)集合形式的一種推廣. 基于前面講過(guò)的波集二象性,可以猜測(cè)或想象,波函數(shù)μAi(x)應(yīng)該對(duì)應(yīng)一種新的集合形式,稱之為“Fuzzy集合”. 回顧上一節(jié)的事實(shí) μAi(X)=[0,1], ?i∈{0,1,…,n}, 自然會(huì)給出定義:所謂論域U上的一個(gè)Fuzzy集A是指,對(duì)于任何元素u∈U,存在唯一的一個(gè)實(shí)數(shù)μA(u)∈[0,1]與u對(duì)應(yīng),這里μA(u)叫做元素u關(guān)于A的隸屬度.這意味著得到一個(gè)映射 μA:U→[0,1],u→μA(u), 類似于描述集合的Venn圖,可以給出描述Fuzzy集的泛Venn圖.事實(shí)上,首先用一個(gè)矩形表示論域U.然后,論域U中的元素用具有單位長(zhǎng)度的線段來(lái)表示,而論域U上的一個(gè)Fuzzy集A表示為矩形中的一個(gè)橢圓,見(jiàn)圖8,可見(jiàn): μA(u1)=1,μA(u2)=0, μA(u3)=0.3,μA(u4)=0.6. 圖 8 泛Venn圖 μA(U)={0,1}, 則隸屬函數(shù)便退化為特征函數(shù),這意味著: A={u∈U|μA(u)=1}. 因此,Cantor集合是Fuzzy集合的特例[6].注意到 μA(u0)∈(0,1), ?u0∈U, 這里(0,1)是個(gè)開(kāi)區(qū)間.(4)式所表達(dá)的Fuzzy集 μAi(x),i=0,1,…,n, 就是常用的真Fuzzy集的例子. 回顧代數(shù)系統(tǒng)(Ch(U),∨,∧,c),記 Meb(U)={μ|μ:U→[0,1]}, 定義2個(gè)二元代數(shù)運(yùn)算∨、∧和一個(gè)一元代數(shù)運(yùn)算“c”如下: ∨:Meb(U)×Meb(U)→Meb(U), (μ1,μ2)→μ1∨μ2=∨(μ1,μ2), (μ1∨μ2)(u)=μ1(u)∨μ2(u), ?u∈U; ∧:Meb(U)×Meb(U)→Meb(U), (μ1,μ2)→μ1∧μ2=∧(μ1,μ2), (μ1∧μ2)(u)=μ1(u)∧μ2(u), ?u∈U; c:Meb(U)→Meb(U),μ→μc=c(μ), μc(u)=1-μ(u), ?u∈U, 于是便得到一個(gè)代數(shù)系統(tǒng)(Meb(U),∨,∧,c)[7]. 值得指出的是,從同構(gòu)關(guān)系 A?B?(?u∈U)(μA(u)≥μB(u)); A=B?(A?B)∧(B?A)? (?u∈U)(μA(u)=μB(u)); C=A∪B? (?u∈U)(μC(u)=μA(u)∨μB(u)); C=A∩B? (?u∈U)(μC(u)=μA(u)∧μB(u)); C=Ac?(?u∈U)(μC(u)=1-μA(u)). (6) 其實(shí),(6)式也可以寫為:?(u,v)∈U×V,μA×B(u,v)=μA(u)·μB(v). (7) 正是基于同構(gòu)關(guān)系 現(xiàn)在再探討一下Zadeh的擴(kuò)展原理.設(shè)X和Y為2個(gè)非空論域,考慮映射: f:X→Y,x→y=f(x). Zadeh的擴(kuò)展原理描述的事情是,如何把映射f:X→Y擴(kuò)展為如下的映射: 首先,考慮Zadeh的擴(kuò)展原理的特殊情況,即如何把映射f:X→Y擴(kuò)展為如下的映射: 如圖9所示,熟知,映射f實(shí)際上是一個(gè)關(guān)系f?X×Y.做下面的集變換: A→B=f°(A)=[(A×Y)∩f]Y. (8) 圖 9 Cantor集的集變換 注意,有關(guān)Cantor集合與Fuzzy集合的投影以及集變換的具體內(nèi)容可參考文獻(xiàn)[1]中的第二章.特別地,利用文獻(xiàn)[1]中的Proposition 2.5.1,有下面的命題. 命題 3.1關(guān)于集變換(8)式,有 (9) 證明對(duì)于任意的y∈Y,根據(jù)文獻(xiàn)[1]中的Proposition 2.5.1,有下面的結(jié)果: χB(y)=χf°(A)(y)=χ[(A×Y)∩f]Y(y)= 因此,(9)式為真,參見(jiàn)圖9. 現(xiàn)在轉(zhuǎn)而考慮如何將映射f:X→Y擴(kuò)展為下面的映射 從圖9得到啟發(fā),可以做出如下的Fuzzy集合的集變換: A→B=f*(A)=[(A×Y)∩f]Y. (10) 再利用文獻(xiàn)[1]中的Proposition 2.6.1,有下面的命題. 命題 3.2關(guān)于Fuzzy集變換(10),有 (11) 證明對(duì)于任意的y∈Y,根據(jù)文獻(xiàn)[1]中的Proposition 2.6.1,有下面的表達(dá)式 μB(y)=μf*(A)(y)=μ[(A×Y)∩f]Y(y)= 因此,(10)式為真. 注 3.1表達(dá)式(10)恰好是熟知的Zadeh的擴(kuò)展原理.因?yàn)閆adeh的擴(kuò)展原理是作為原理出現(xiàn)的,故它不需要證明.正如上面的討論,Zadeh的擴(kuò)展原理不是個(gè)原理,而是個(gè)命題,這說(shuō)明在前面給出的擴(kuò)展原理是更為基礎(chǔ)的原理,使用范圍更廣,它把Zadeh的擴(kuò)展原理作為一個(gè)命題證明出來(lái),這是很有意義的事情. (10)式是關(guān)于一元函數(shù)的擴(kuò)展形式,當(dāng)然可以考慮多元函數(shù)的擴(kuò)展形式,以二元函數(shù)為例進(jìn)行討論. 設(shè)X、Y、Z是3個(gè)非空論域,考慮下面的映射: f:X×Y→Z, (x,y)→z=f(x,y). 首先把映射f:X×Y→Z擴(kuò)展為下面的樣子: 根據(jù)文獻(xiàn)[1]中的Definition 2.6.1,得到下面的Fuzzy集變換的形式: R→C=f*(R)=[(R×Z)∩f]Z. (12) 然后,利用文獻(xiàn)[1]中的Proposition 2.6.1,有下面的命題. 命題 3.3關(guān)于Fuzzy集變換(12),有下面的表達(dá)式 (13) 證明對(duì)于任意的z∈Z,文獻(xiàn)[1]中的Proposition 2.6.1,注意事實(shí)f?X×Y×Z,有下面的結(jié)果: μC(z)=μf*(R)(z)=μ[(R×Z)∩f]Z(z)= χf(x,y,f(x,y)))= 便有如下的表達(dá)式 因此,(13)式為真. 再繼續(xù)考慮把映射f:X×Y→Z擴(kuò)展為下面的樣子: (A,B)→C=f*(A,B). (A,B)→C=f*(A,B)= [((A×B)×Z)∩f]Z. (14) 根據(jù)文獻(xiàn)[1]中的Proposition 2.10.2,對(duì)于任意的z∈Z,有下面的結(jié)果 (15) 注 3.2關(guān)于Fuzzy關(guān)系 對(duì)于任何的(x,y)∈X×Y,不難知道下面的事實(shí) μR(x,y)=μA×B(x,y)=μA(x)·μB(y). 因此,(15)式也可以寫為下面的樣子 z∈Z. (16) 還可以考慮另外一個(gè)問(wèn)題.設(shè)X和Y是2個(gè)非空論域,考慮映射: f:X→Y,x→y=f(x). 熟知,Zadeh的擴(kuò)展原理是基于映射f:X→Y獲得下面的映射: 另一方面,基于映射f:X→Y,應(yīng)當(dāng)還考慮如下的逆映射: 為了解作為Fuzzy集合(f-1)*(B)的隸屬函數(shù),先考慮特例:把映射f:X→Y擴(kuò)展為如下的逆映射: 事實(shí)上,首先顯然有 (f-1)°(B)={x∈X|f(x)∈B}, 于是便有下面的表達(dá)式 χ(f-1)°(B)(x)=χ{x∈X|f(x)∈B}(x), ?x∈X. 命題 3.4?x∈X,一定有 χ{x∈X|f(x)∈B}(x)=χB(f(x)). 證明對(duì)于任意的x∈X,不難看出如下的等價(jià)式 χ{t∈X|f(t)∈B}(x)=1? x∈{t∈X|f(t)∈B}? f(x)∈B?χB(f(x))=1. 因此,該命題為真. 根據(jù)命題3.4,利用提出的擴(kuò)展原理,可以得到下面的結(jié)果 μ(f-1)*(B)(x)=μB(f(x)), ?x∈X. (17) 這里 n=1,2,3,… 顯然它們滿足條件 并且 n>m?‖Δ(n)‖<‖Δ(m)‖, ?n,m∈N. 注意到 由此構(gòu)造如下波函數(shù)組(參見(jiàn)圖10): i=1,2,…,n-1, 圖 10 連續(xù)波函數(shù) 從而得到一個(gè)線性子空間序列 利用f構(gòu)造常數(shù) 由此得到基底一個(gè)線性組合 這意味著獲得一個(gè)連續(xù)函數(shù)序列 可以驗(yàn)證 圖 11 分段線性插值函數(shù) 現(xiàn)在仍然使用基函數(shù)組 再利用已知的常數(shù)組 構(gòu)造一組新的常數(shù): ?ε>0, ?N∈N+, ?n,m∈N+,n,m>N? 例 4.1導(dǎo)函數(shù)逼近的仿真.取 f(x)=sin(x)∈C[0,2π], 熟知 f′(x)=cos(x)∈C[0,2π]. 圖12~15是n取不同值時(shí)的近似導(dǎo)函數(shù)(粗線)和真實(shí)導(dǎo)函數(shù)(細(xì)線). 下面討論原函數(shù)的逼近意義.設(shè)函數(shù)g(x)∈R[a,b],考慮[a,b]的等距分割全體: 圖 12 n=20 圖 13 n=40 圖 14 n=50 圖15 n=100 n=1,2,3,… 利用常數(shù) 構(gòu)造線性方程組 這里 是待求的未知量.寫為遞推式: … i=1,2,…,n-1, 不難理解,函數(shù)g(x)的原函數(shù)的一般形式,即原函數(shù)的通解為f(x,C0)+C,其中C為任意常數(shù). 例 5.1原函數(shù)逼近的仿真.取 g(x)=cos(x)∈C[0,2π], 熟知原函數(shù)為 f(x)=sin(x)+C. 取C=0,則 f(x)=sin(x)∈C[0,2π]. 圖16~19是n取不同值時(shí)的近似原函數(shù)(粗線)和真實(shí)原函數(shù)(細(xì)線)的圖像. 圖 16 n=50 圖 17 n=100 圖 18 n=200 圖 19 n=300 3 Fuzzy集合與擴(kuò)展原理
4 導(dǎo)函數(shù)的逼近論意義
5 原函數(shù)的逼近論意義