《物元可拓源碼》由會員分享,可在線閱讀,更多相關(guān)《物元可拓源碼(4頁珍藏版)》請?jiān)谘b配圖網(wǎng)上搜索。
1、精選優(yōu)質(zhì)文檔-----傾情為你奉上
物元可拓法于80年代由我國蔡文教授創(chuàng)立,目前已廣泛應(yīng)用于新產(chǎn)品構(gòu)思與設(shè)計(jì)、優(yōu)化決策、控制、識別與評價(jià)等各個領(lǐng)域,無論在理論還是在實(shí)踐上都發(fā)揮了越來越重要的作用。
物元是描述事物的名稱、特征及量值3個基本元素的簡稱,在形式上可記為M=(N,c,v)=(N,c,c(N))。其中M、N、c、v分別是Matter、Name,?Character,?Value的縮寫??赏丶鲜怯藐P(guān)聯(lián)度將模糊集合的[0,1]閉合區(qū)間連續(xù)取值拓廣到(-∞,+∞)實(shí)數(shù)軸,以表達(dá)物元的量值為實(shí)軸上的一點(diǎn)時(shí)符合要求的程度。物元分析是研究物元及其變化并用以解決矛盾問題的規(guī)律和方法,可拓學(xué)是用
2、形式化的工具,從定性和定量兩個角度去研究解決矛盾問題的規(guī)律和方法。物元可拓法結(jié)合二者,是將辨證邏輯和形式邏輯相結(jié)合的可拓邏輯,豐富了事物的內(nèi)涵,客觀地反映了物質(zhì)世界的真實(shí)狀態(tài)。
本次選用評價(jià)因子污染貢獻(xiàn)率方法來確定權(quán)系數(shù)。主要計(jì)算程序:
Dim?sRow?As?Integer,?sCol?As?Integer????'起始的行與列
Dim?i?As?Integer,?j?As?Integer??????????'循環(huán)變量
Dim?Xj?As?Double????????????????????????'定義實(shí)測值
Dim?Aij?As?Double,?Bij?As?Double??????
3、??'定義標(biāo)準(zhǔn)域區(qū)間
Dim?Apj?As?Double,?Bpj?As?Double????????'定義節(jié)域變量
Dim?YZS?As?Integer??????????????????????'定義評價(jià)因子個數(shù)
Dim?DJS?As?Integer??????????????????????'定義評價(jià)等級數(shù)
'得到起始行列值
sRow?=?InputBox("請輸入監(jiān)測數(shù)據(jù)第一個數(shù)的行號!",?"輸入行號",?0)
sCol?=?InputBox("請輸入監(jiān)測數(shù)據(jù)第一個數(shù)的列號!",?"輸入列號",?0)
YZS?=?InputBox("請輸入評價(jià)因子個數(shù)!",?"輸入因子個數(shù)",?0)
4、
DJS?=?InputBox("請輸入評價(jià)等級個數(shù)!",?"輸入評價(jià)等級數(shù)",?0)
'插入標(biāo)記列文字
With?Sheets("sheet1")
??For?i?=?1?To?DJS
??????Cells(sRow?+?DJS?+?2?+?i,?sCol?-?1).Value?=?"關(guān)聯(lián)函數(shù)k_等級"?&?i
??Next?i
??Cells(sRow?+?2?*?DJS?+?3,?sCol?-?1).Value?=?"X/S"
??Cells(sRow?+?2?*?DJS?+?4,?sCol?-?1).Value?=?"歸一化權(quán)重"
??For?i?=?1?To?DJS
??????C
5、ells(sRow?+?2?*?DJS?+?4?+?i,?sCol?-?1).Value?=?"關(guān)聯(lián)度K_等級"?&?i
??Next?i
??Cells(sRow?+?3?*?DJS?+?5,?sCol?-?1).Value?=?"可拓指數(shù)"
??
??'按列循環(huán)計(jì)算
??For?j?=?sCol?To?sCol?+?YZS?-?1
????'賦初值
????Xj?=?Cells(sRow,?j).Value????????????'實(shí)測值
????Apj?=?Cells(sRow?+?1,?j).Value???????'可拓域最小值
????Bpj?=?Cells(sRow?+?DJS?
6、+?2,?j).Value?'可拓域最大值
????
????For?i?=?1?To?DJS
??????'對aij,bij賦值
??????Aij?=?Cells(sRow?+?i,?j).Value
??????Bij?=?Cells(sRow?+?i?+?1,?j).Value
??????
??????'按條件選擇公式計(jì)算關(guān)聯(lián)度
??????If?Xj?>?Aij?And?Xj?
7、JS?+?2,?j).Value?=?-(Aij?-?Xj)?/?(Bij?-?Aij)
????????Else
??????????Cells(sRow?+?i?+?DJS?+?2,?j).Value?=?-(Xj?-?Bij)?/?(Bij?-?Aij)
????????End?If
????????
??????Else??????????'xj<>Xij?點(diǎn)x位于本標(biāo)準(zhǔn)之外
??????
????????If?Xj?
8、????????????Cells(sRow?+?i?+?DJS?+?2,?j).Value?=?(Aij?-?Xj)?/?(Apj?-?Aij)
??????????Else
????????????Cells(sRow?+?i?+?DJS?+?2,?j).Value?=?(Aij?-?Xj)?/?(2?*?Xj?-?Bpj?-?Aij)
??????????End?If
????????
????????ElseIf?Xj?>?Bij?Then?'x位于標(biāo)準(zhǔn)的右邊,此時(shí)有x>(ai+bi)/2
????????
??????????If?Xj?<=?(Apj?+?Bpj)?/?2?Then
9、
????????????Cells(sRow?+?i?+?DJS?+?2,?j).Value?=?(Xj?-?Bij)?/?(Apj?+?Bij?-?2?*?Xj)
??????????Else
????????????Cells(sRow?+?i?+?DJS?+?2,?j).Value?=?(Xj?-?Bij)?/?(Bij?-?Bpj)
??????????End?If
??????????
????????End?If
??????End?If
????Next?i
??Next?j
??
????'計(jì)算X/S
????For?j?=?sCol?To?sCol?+?YZS?-?1
?
10、???????Dim?a?As?Double
????????a?=?0
????????For?i?=?1?To?DJS?+?2
????????????a?=?a?+?Cells(sRow?+?i,?j)
????????Next?i
????????Cells(sRow?+?2?*?DJS?+?3,?j).Value?=?Cells(sRow,?j).Value?*?(DJS?+?2)?/?a
????Next?j
????'計(jì)算權(quán)重
????'計(jì)算x/s的總和
????a?=?0
????For?i?=?sCol?To?sCol?+?YZS?-?1
????????a?=?a?+?Ce
11、lls(sRow?+?2?*?DJS?+?3,?i)
????Next?i
????
????For?j?=?sCol?To?sCol?+?YZS?-?1
????????Cells(sRow?+?2?*?DJS?+?4,?j).Value?=?Cells(sRow?+?2?*?DJS?+?3,?j).Value?/?a
????Next?j
????
????'計(jì)算關(guān)聯(lián)度
????Cells(sRow?+?2?*?DJS?+?4,?sCol?+?YZS)?=?"綜合關(guān)聯(lián)度"
????For?i?=?1?To?DJS
????????For?j?=?sCol?To?sCol?+?YZS?-?1
12、
????????????Cells(sRow?+?2?*?DJS?+?4?+?i,?j).Value?=?Cells(sRow?+?DJS?+?2?+?i,?j).Value?*?Cells(sRow?+?2?*?DJS?+?4,?j).Value
????????Next?j
????????Dim?k?As?Integer
????????a?=?0
????????For?k?=?sCol?To?sCol?+?YZS?-?1
???????????a?=?a?+?Cells(sRow?+?2?*?DJS?+?4?+?i,?k)?'綜合關(guān)聯(lián)度累加
????????Next?k
?????
13、???Cells(sRow?+?2?*?DJS?+?4?+?i,?sCol?+?YZS).Value?=?a
????Next?i
????'計(jì)算可拓指數(shù)
????'找最小與最大關(guān)聯(lián)度
????Dim?Kmax,?Kmin?As?Double
????Kmax?=?Cells(sRow?+?2?*?DJS?+?4?+?1,?sCol?+?YZS).Value
????Kmin?=?Kmax
????For?i?=?2?To?DJS
??????If?Kmax?
14、=?Cells(sRow?+?2?*?DJS?+?4?+?i,?sCol?+?YZS).Value
??????End?If
??????If?Kmin?>?Cells(sRow?+?2?*?DJS?+?4?+?i,?sCol?+?YZS).Value?Then
????????Kmin?=?Cells(sRow?+?2?*?DJS?+?4?+?i,?sCol?+?YZS).Value
??????End?If
????Next?i
????
????Dim?KXP()?As?Double
????ReDim?KXP(DJS)?As?Double
????For?i?=?1?To?DJS
??
15、????KXP(i)?=?(Cells(sRow?+?2?*?DJS?+?4?+?i,?sCol?+?YZS).Value?-?Kmin)?/?(Kmax?-?Kmin)
????Next?i
????Dim?FZ,?FM?As?Double
????For?i?=?1?To?DJS
??????FZ?=?FZ?+?i?*?KXP(i)
??????FM?=?FM?+?KXP(i)
????Next?i
????Cells(sRow?+?3?*?DJS?+?5,?sCol).Value?=?FZ?/?FM
End?With
我做地下水水質(zhì)評價(jià)時(shí)用到的公式是這樣的:
矩跟關(guān)聯(lián)度的公式是
16、通常用到的那個,即:
p(xj,xij)=|xj-0.5(aij+bij)|-0.5(bij-aij)
關(guān)聯(lián)函數(shù)計(jì)算公式:
當(dāng)xj屬于xij時(shí),ki(xj)=-p(xj,xij)/|xij|
當(dāng)xj不屬于xij時(shí),ki(xj)=p(xj,xij)/[p(xj,xrj)-p(xj,xij)]
綜合關(guān)聯(lián)度計(jì)算公式:
Kj(p)=a1*kj(x1)+a2*kj(x2)+...+an*kj(xn)
權(quán)重用的污染因子貢獻(xiàn)率法:
ai=(xi/Si)/(x1/S1+x2/S2+...xn/Sn)
可拓指數(shù)用公式:
KXP=(Kj(X)-Min(Kj(X)))/(Max(Kj(X))-Min(Kj(X)))
由于涉及的內(nèi)容或者思路一同,可能選用的關(guān)聯(lián)度計(jì)算公式,權(quán)重計(jì)算公式不同,酌情參考。
參考源程序里最初的部分,后面的變量說明有詳細(xì)解釋,定義要評價(jià)的參數(shù)個數(shù),評價(jià)的等級等,相應(yīng)的在excel表格里的參數(shù)個數(shù)與評價(jià)等級的行數(shù)要與之對應(yīng)。
PS:上面的代碼里因子個數(shù)是通過跳出的對話框手動輸入進(jìn)去的。
專心---專注---專業(yè)