Kusanthula kwachidziwitso kumagwiritsidwa ntchito kuti adziwe momwe kusonkhanitsa kwazinthu zodziyimira pawokha kumadalira pazikhalidwe zina.
Ndi njira yamphamvu yodziwira momwe zotulutsa zachitsanzo zimakhudzidwira ndi zolowa zamtundu uliwonse. Mu positi iyi, ndikuwonetsa mwachidule kusanthula kwachangu pogwiritsa ntchito SALib, phukusi laulere la Python sensitivity kusanthula.
Nambala yomwe imadziwika kuti sensitivity index, nthawi zambiri imayimira kukhudzika kwa chilichonse. Pali mitundu ingapo ya sensitivity indices:
- Ma index a dongosolo loyamba: amawerengera zopereka zachitsanzo chimodzi pakusintha kotulutsa.
- Ma indices a dongosolo lachiwiri: amawerengera zopereka zamitundu iwiri kuti zitheke.
- Mlozera wa dongosolo lonse: zimawerengera zomwe zathandizira zachitsanzo pakusiyana kotulutsa, kuphatikiza zonse zoyambira (zolowetsazo zimasinthasintha zokha) ndi kulumikizana kulikonse kwapamwamba.
Kodi SALib ndi chiyani?
SALIB ndi Python-based gwero lotseguka zida zoyesera zowunikira. Ili ndi kayendedwe ka ntchito kosiyana, zomwe zikutanthauza kuti sizimalumikizana mwachindunji ndi masamu kapena masamu. M'malo mwake, SALib imayang'anira kupanga zolowetsa zachitsanzo (kudzera m'chitsanzo chimodzi) ndikuyika ma index a sensitivity (kudzera m'modzi mwa ntchito zowunikira) kuchokera pazotsatira zachitsanzo.
Kuwunika kodziwika bwino kwa SALib kumakhala ndi njira zinayi:
- Dziwani zolowetsa zachitsanzo (magawo) ndi mtundu wachitsanzo cha chilichonse.
- Kuti mupange zolowetsa zachitsanzo, gwiritsani ntchito chitsanzo.
- Yesani chitsanzocho pogwiritsa ntchito zolowetsa zomwe zapangidwa ndikusunga zotsatira zachitsanzo.
- Kuti muwerenge ma sensitivity indices, gwiritsani ntchito kusanthula pazotulutsa.
Sobol, Morris, ndi FAST ndi ochepa chabe mwa njira zowunikira zomwe zimaperekedwa ndi SALib. Zinthu zambiri zimakhudza njira yomwe ili yabwino kwambiri pakugwiritsa ntchito, monga tiwonera mtsogolo. Pakadali pano, kumbukirani kuti muyenera kugwiritsa ntchito ntchito ziwiri zokha, zitsanzo ndi kusanthula, mosasamala kanthu za njira yomwe mumagwiritsa ntchito. Tikuwongolerani chitsanzo choyambirira chofotokozera momwe mungagwiritsire ntchito SALib.
Chitsanzo cha SALib - Sobol' Sensitivity Analysis
Mu chitsanzo ichi, tiwona mphamvu ya Sobol ya ntchito ya Ishigami, monga momwe tawonetsera pansipa. Chifukwa cha kusagwirizana kwake kwakukulu komanso kusakhalapo kwa monotonicity, ntchito ya Ishigami imagwiritsidwa ntchito kwambiri poyesa njira zowunika za kusatsimikizika komanso kukhudzidwa.
Masitepe amapita motere:
1. Kuitanitsa SALib
Gawo loyamba ndikuwonjezera malaibulale ofunikira. Zitsanzo ndi kusanthula ntchito za SALib zimasungidwa mosiyana m'ma module a Python. Kulowetsa chitsanzo cha satellite ndi ntchito zowunikira za Sobol, mwachitsanzo, zikuwonetsedwa pansipa.
Timagwiritsanso ntchito ntchito ya Ishigami, yomwe imapezeka ngati ntchito yoyesera ku SALib. Pomaliza, timatumiza NumPy monga SALib imagwiritsa ntchito kusunga zolowetsa ndi zotuluka mu matrix.
2. Kulowetsa kwachitsanzo
Kenako zolowetsa zachitsanzo ziyenera kufotokozedwa. Ntchito ya Ishigami imavomereza zolowetsa zitatu: x1, x2, ndi x3. Ku SALib, timapanga lamulo lomwe limafotokoza kuchuluka kwa zolowetsa, mayina awo, ndi malire pazolowera zilizonse, monga tawonera pansipa.
3. Pangani Zitsanzo ndi Chitsanzo
Zitsanzozo zimapangidwa. Tiyenera kupanga zitsanzo pogwiritsa ntchito Saltelli sampler popeza tikuchita kafukufuku wa Sobol sensitivity. Pankhaniyi, ma param ndi matrix a NumPy. Titha kuwona kuti matrix ndi 8000 ndi 3 poyendetsa param values.shape. Zitsanzo 8000 zidapangidwa ndi Saltelli sampler. Saltelli sampler imapanga zitsanzo, kumene N ndi 1024 (parameter yomwe tinapereka) ndi D ndi 3. (chiwerengero cha zolowetsa chitsanzo).
Monga tanena kale, SALib sichita nawo masamu kapena ma computational model. Ngati chitsanzocho chalembedwa mu Python, nthawi zambiri mumadutsa pazolowera zilizonse ndikuyesa chitsanzo:
Zitsanzo zitha kusungidwa ku fayilo yamawu ngati mtunduwo sunapangidwe mu Python:
Mzere uliwonse wa param values.txt umayimira chitsanzo chimodzi. Zotulutsa zachitsanzo ziyenera kusungidwa ku fayilo ina mwanjira yofananira, ndikutulutsa kumodzi pamzere uliwonse. Pambuyo pake, zotsatirazo zikhoza kuikidwa ndi:
Muchitsanzo ichi, tigwiritsa ntchito Ishigami kuchokera ku SALib. Ntchito zoyesererazi zitha kuyesedwa motere:
4. Pangani Analysis
Titha kuwerengera ma index a sensitivity titatha kutsitsa zotsatira zachitsanzo ku Python. Muchitsanzo ichi, tidzagwiritsa ntchito sobol.analyze kuti tiwerenge zizindikiro zoyamba, zachiwiri, ndi zonse.
Si ndi dikishonale ya Python yokhala ndi makiyi "S1," "S2," "ST," "S1 conf," "S2 conf," ndi "ST conf." Makiyi a _conf amakhala ndi nthawi yodalirika yolumikizana, yomwe nthawi zambiri imayikidwa 95 peresenti. Kuti mutulutse, ma indices onse, gwiritsani ntchito mawu osakira parameter print to console=Zowona. Kapenanso, monga tawonera pansipa, titha kusindikiza zikhalidwe za Si.
Titha kuwona kuti x1 ndi x2 zili ndi chidwi choyitanitsa koyamba, koma x3 sikuwoneka kuti ili ndi zotsatira zoyambira.
Ngati zizindikiro za dongosolo lonse ndi zazikulu kwambiri kuposa zizindikiro zoyamba, kuyanjana kwapamwamba kumachitikadi. Titha kuwona kuyanjana kwapamwambaku poyang'ana zizindikiro zachiwiri:
Titha kuwona kuti x1 ndi x3 zimalumikizana kwambiri. Pambuyo pake, zotsatira zake zitha kusinthidwa kukhala Pandas DataFrame kuti mupitirize kuphunzira.
5. Kukonza chiwembu
Kuti mukhale omasuka, zida zoyambira zopangira ma chart zimaperekedwa. Chiwembu () chimapanga zinthu za matplotlib axis kuti zisinthe.
Kutsiliza
SALib ndi chida chowunikira kwambiri cha sensitivity. Njira zina mu SALib zikuphatikiza Mayeso a Fourier Amplitude Sensitivity Test (FAST), Morris Method, ndi Delta-Moment Independent Measure. Ngakhale ili laibulale ya Python, idapangidwa kuti izigwira ntchito ndi mitundu yamtundu uliwonse.
SALib imapereka mawonekedwe osavuta kugwiritsa ntchito pamzere wamalamulo popanga zolowetsa zachitsanzo ndikuwunika zotuluka. Onani Zolemba za SALib kudziwa zambiri.
Siyani Mumakonda