A na-eji nyocha nke mmetụta uche chọpụta mmetụta nke nchịkọta ihe ndị nwere onwe ha na mgbanwe dabere n'okpuru ọnọdụ ụfọdụ.
Ọ bụ ụzọ siri ike maka ịchọpụta ka mmepụta ihe nlereanya si emetụta ihe ntinye ihe nlereanya n'ozuzu ya. N'ime ọkwa a, m ga-enye nkọwa ngwa ngwa nke nyocha uche site na iji SALib, ngwugwu nyocha Python sensitivity n'efu.
Uru ọnụọgụgụ nke a maara dị ka ndenye mmetụta uche, na-egosipụtakarị mmetụta ntinye ọ bụla. Enwere ọtụtụ ụdị indices sensitivity:
- Indices-usoro nke mbụ: na-agbakọ ntinye nke otu ihe nlere anya na ọdịiche mmepụta.
- Indices usoro nke abụọ: na-agbakọ ntinye ntinye ihe atụ abụọ maka mgbanwe dị iche iche.
- Ngụkọta usoro ndeksi: na-agụta ntinye ntinye ihe nlereanya na mgbanwe dị iche iche nke mmepụta, na-agụnye ma nsonaazụ mbụ (ntinye na-agbanwe naanị) yana mmekọrịta ọ bụla dị elu.
Kedu ihe bụ SAlib?
SAlib bụ Python dabere Oghere Toolkit maka ime ntule uche. Ọ nwere usoro nrụpụta ọrụ nke gbapụrụ agbawa, nke pụtara na ọ naghị emekọ ihe ozugbo na ụdị mgbakọ na mwepụ ma ọ bụ mgbakọ. Kama, SALib na-ahụ maka ịmepụta ntinye ihe nlereanya (site na otu n'ime ọrụ nlele) na ịgbakọ indices sensitivity (site na otu n'ime ọrụ nyocha) site na mmepụta ihe nlereanya.
Nyocha mmetụta uche SAlib na-ahụkarị nwere usoro anọ:
- Kpebisie ike ntinye ihe nlere anya (parameters) na ihe nlele maka nke ọ bụla.
- Iji mepụta ntinye ihe nlereanya, rụọ ọrụ nlele.
- Nyochaa ihe nlereanya site na iji ntinye ewepụtara wee chekwaa nsonaazụ nlereanya.
- Iji gbakọọ akara nhụsianya, jiri ọrụ nyocha na ntinye.
Sobol, Morris, na FAST bụ ole na ole n'ime ụzọ nyocha uche nke SAlib nyere. Ọtụtụ ihe na-emetụta ụzọ kacha mma maka ngwa enyere, dị ka anyị ga-ahụ ma emechaa. Maka oge a, buru n'uche na naanị ị ga-eji ọrụ abụọ, nlele na nyocha, n'agbanyeghị usoro ọ bụla ị na-eji. Anyị ga-eduzi gị site na ihe atụ bụ isi iji gosi otu esi eji SALib.
Ihe Nlereanya SALib – Sobol' Sensitivity Analysis
N'ihe atụ a, anyị ga-enyocha mmetụta Sobol' nke ọrụ Ishigami, dị ka egosiri n'okpuru. N'ihi nnukwu enweghị usoro ya na enweghị monotonicity, a na-eji ọrụ Ishigami n'ọtụtụ ebe iji nyochaa ejighị n'aka na usoro nyocha uche.
Usoro na-aga dị ka ndị a:
1. Mbubata SAlib
Nzọụkwụ mbụ bụ ịgbakwunye ọba akwụkwọ achọrọ. A na-edobe ihe nlele na nyocha ọrụ SAlib dị iche na modul Python. Ịbubata ihe nlele satịlaịtị na ọrụ nyocha Sobol, dịka ọmụmaatụ, egosiri n'okpuru.
Anyị na-ejikwa ọrụ Ishigami, nke dị ka ọrụ nnwale na SAlib. N'ikpeazụ, anyị na-ebubata NumPy ka SALib na-eji ya echekwa ihe ntinye na ntinye ihe nlereanya na matriks.
2. Ntinye ihe nlereanya
Ekwesịrị ịkọwapụta ntinye ihe nlereanya. Ọrụ Ishigami na-anabata ntinye atọ: x1, x2, na x3. Na SAlib, anyị na-arụ dict nke na-akọwapụta ọnụọgụ ntinye, aha ha, na oke na ntinye ọ bụla, dị ka a hụrụ n'okpuru.
3. Mepụta Samples na Model
A na-emepụta ihe nlele ahụ. Anyị kwesịrị ịmepụta ihe nlele site na iji Saltelli sampler ebe anyị na-eme nyocha uche nke Sobol. N'okwu a, ụkpụrụ param bụ matrix NumPy. Anyị nwere ike ịhụ na matrix bụ 8000 site na 3 site na-agba ọsọ param values.shape. Emebere ihe nlele 8000 na Saltelli sampler. Saltelli sampler na-emepụta ihe atụ, ebe N bụ 1024 (nkeji anyị nyere) na D bụ 3. (ọnụọgụ nke ntinye ihe nlereanya).
Dịka ekwuru na mbụ, SALib anaghị etinye aka na mgbakọ na mwepụ ma ọ bụ nyocha ihe nlere anya. Ọ bụrụ na edere ihe nlereanya a na Python, ị ga-agbanyekarị ntinye nlele ọ bụla wee nyochaa ihe nlereanya:
Enwere ike ichekwa ihe nlele ahụ na faịlụ ederede ma ọ bụrụ na emepụtaghị ihe nlereanya na Python:
Ahịrị ọ bụla dị na param values.txt na-anọchite anya ntinye ụdị. Ekwesịrị ịchekwa mmepụta ihe nlereanya ahụ na faịlụ ọzọ n'ụdị yiri ya, yana otu mmepụta na ahịrị ọ bụla. Mgbe nke ahụ gasịrị, enwere ike ibunye nsonaazụ ya:
N'ihe atụ a, anyị ga-eji ọrụ Ishigami sitere na SAlib. Enwere ike nyochaa ọrụ ule ndị a dị ka ndị a:
4. Mee Analysis
Anyị nwere ike mechaa gbakọọ indices sensitivity ka anyị tinyechara nsonaazụ ihe nlereanya n'ime Python. N'ihe atụ a, anyị ga-eji sobol.analyze gbakọọ akara nke mbụ, nke abụọ na mkpokọta.
Si bụ akwụkwọ ọkọwa okwu Python nwere igodo “S1,” “S2,” “ST,” “S1 conf,” “S2 conf,” na “ST conf.” Igodo _conf na-ejide oge ntụkwasị obi emetụtara, nke a na-atọkarị ya na pasentị 95. Iji wepụta, indices niile, jiri mbipụta paramita isiokwu ka ịkasi obi=Eziokwu. N'aka nke ọzọ, dịka e gosipụtara n'okpuru, anyị nwere ike ibipụta ụkpụrụ nke onye ọ bụla site na Si.
Anyị nwere ike ịhụ na x1 na x2 nwere nghọta nke mbụ, mana x3 apụtaghị na ọ nwere mmetụta ọ bụla nke mbụ.
Ọ bụrụ na ngụkọta usoro ngụkọta dị ukwuu karịa akara nrịbama nke mbụ, mmekọrịta dị elu na-eme n'ezie. Anyị nwere ike ịhụ mmekọrịta ndị a dị elu site n'ile anya n'usoro nke abụọ:
Anyị nwere ike ịhụ na x1 na x3 nwere nnukwu mmekọrịta. Mgbe nke ahụ gasịrị, enwere ike gbanwee nsonaazụ ya ka ọ bụrụ Pandas DataFrame maka ọmụmụ ihe ọzọ.
5. Ịgba nkata
Maka ịdị mma gị, ewepụtara ụlọ ọrụ nhazi ihe. Ọrụ nkata () na-emepụta ihe axis matplotlib maka ngbanwe na-esote.
mmechi
SAlib bụ ngwa nyocha nyocha ọkaibe. Usoro ndị ọzọ dị na SAlib gụnyere Fourier Amplitude Sensitivity Test (FAST), Morris Method, na Delta-Moment Independent Measure. Ọ bụ ezie na ọ bụ ọba akwụkwọ Python, e bu n'obi were ụdị ụdị ọ bụla rụọ ọrụ.
SAlib na-enye interface-ahịrị iwu dị mfe iji mee ihe maka ịmepụta ntinye ihe nlereanya na ịlele mmepụta ihe nlereanya. Lelee akwụkwọ SAlib ịmụtakwu ihe.
Nkume a-aza