|
U - U * | 2 < £ 2 (3.23)
bu erda £ 1 - kuzatuv xatosiga mutanosib ravishda belgilangan taxminiy aniqlik; £ 2 - aniqlik muammo yechimlari, tortishish va maydon U aniqligi bilan mutanosib , asl maydon U va maydon U * o'rtasidagi farq , allaqachon belgilangan tufayli - E qurilish komponentlarini yaqinlashtirish (3.19).
Agar siz toymasin oynada tartibsiz ravishda tarmoqqa tushsangiz - shart qanoatlantiriladi | u | < £ 1 , keyin oynaning markazida U maydonining qiymatini aniqlashda bu nuqtalar e'tiborga olinmaydi. Agar oyna ichidagi barcha nuqtalar uchun yuqoridagi shart bajarilsa, u holda oynaning markazi ostidagi shar mos kelmaydi.
Ko'rinib turibdiki, har bir ursda joylashgan sharlar soni taxminan N , - yo'q S = S 1 , S 2 , S 3 ... nuqtalarning umumiy sonidan har doim kamroq bo'ladi: N Muammoni to'liq hal qilish mezonlari end-busting - lichestva k chuqurlik darajalariga beriladi; (yoki) eritmaning kerakli aniqligiga erishish (3.23); (yoki) takrorlanish jarayonining divergensiyaga moyilligining namoyon bo'lishi F + 1 > F i , F - (uning yakuniy bosqichida).
Manbani taxmin qilish algoritmi quyidagicha ko'rinadi:
U maydonidan doimiy komponentni yo'q qilish .
Simulyatsiya qilingan maydon massivini shakllantirish U : U j - U j , j - 1,2, ..., M
R s boshlang'ich qiymatini tanlash .
Tashqi halqa: sanab darajalari (sirtlari) S (i = 1,2,3 ... k) joylashtirish ann - roksimiruyuschih sohalar (odatda k = 6-8).
Ichki halqa 1: skanerlash hududi R surma oynasi R ' s x R : : tanlash komponenti U , o'rtacha balandlikni kuzatish sirt ta'rifi - N. sohalar uchun bo'linish koordinatalari ( £ n , Z) .
Ichki halqa 2: SLAE yechimi (3.21) Zaydel tomonidan hisob-kitobni yakunlash - shartdan keyin (3.22) yoki oldindan belgilangan takrorlash sonidan keyin.
Ichki halqa 3: Maydonni qayta tiklash U i * dan M tugunlari tartibsiz panjara va chiqarib tashlangan - uni simulyatsiya qilingan maydondan chenie U i + 1 = U i - U *
Fayldagi yozuv samarali taxminiy parametrlarni kiritdi - {^ j , n j , Z j , P j }, j = 1.2, ..., N . ...
Eritmasi smeta (3.23) sifati (mezonlar bajarilishi taqdirda to'liq echimlar hisoblash jarayoni muammoni to'xtaydi).
R S uchun yangi qiymatni tanlash : R S + 1 = R S / 2 va 4-bosqichga o'ting.
Hisobni to'ldirish.
To'g'ri to'rtburchaklar tarmoq tugunlari yoki maydonlardagi tiklanish maydonlari qiymatlari ma'lum fizik va geometrik parametrlarga ega bo'lgan taxminiy tuzilmaning to'g'ridan-to'g'ri muammoli echimlarini o'zgartirish orqali amalga oshiriladi - ekvivalent muhit modeli. F * = F * ( x , y , z , -, Z, P) to'g'ridan-to'g'ri masalani yechish operatorining shakli masalaning gravitatsiyaviy maydoni A g uchun masalaning qo'yilishi bilan aniqlanadi :
• interpolyatsiya bilan F * = F (3.20);
• potentsial V (antiderivativ) ni hisoblashda : F * = 1 / R, bu erda
R = ( K - - x ) 2 + ( n - Y ) 2 + ( Z - z ) 2;
• tortishish maydonining 1-vertikal hosilasini hisoblashda V : F * = [2 ( Z - z) 2 - ( - - x) 2 - ( n -y) 2 ] / R 5 .
Ushbu algoritmni amalga oshirish uchun ANMOD dasturi foydalanish uchun tuzilgan - u vizual tizim ob'ektga yo'naltirilgan dasturlash Delphi 4, Ishga tushirildi - eritish 32-bitli operatsion tizim Windows 95/98/2000 / NT. O'zining qulay foydalanuvchi interfeysi bilan bir qatorda, uni yaratish uchun dasturchidan Windows 32 c uchun Delphi 4 imple-dan foydalangan holda ilovalarni ishlab chiqish uchun minimal kuch talab etiladi - imtiyozlar himoyalangan protsessor rejimini shakllantirish. Bu didga moslab qabul - , printsipial yangi algoritmik hal tezligi va Uve hisoblash oshdi - natijalarini lichenie aniqligini. Xususan, dinamik massivlarda joylashtiriladigan dastlabki ma'lumotlar miqdori bo'yicha deyarli hech qanday oschu - timye cheklovi yo'q , ularning o'lchamlari faqat kompyuter xotirasi hajmiga bog'liq (bir formatli ma'lumotlarni yozishni boshlashda uchun 1 dona kerak bo'ladi) MB).
Boshlang'ich ma'lumotlar vektorlar to'plamidir (x, y, z, U) , ularning har biri od - xlorid nuqtasida asl maydonning fazoviy holatini va amplitudasini tavsiflaydi . Ma'lumotlar matnli fayl ko'rinishida, GIS SURFER formatida (dat - fayl) yoziladi. Gorizontal tekislikda z = 0 yoki krivoli - chiziqli sirtda joylashgan oddiy to'r (grd-fayl) tugunlarida maydon tasvirining hosil bo'lgan qiymatlarini hisoblash jarayonida apriori ko'rsatilgan qiymatlar natijalarni ko'taradi - kommutativ xarakterlanadi. nuqtalar z = z (x, y) . Barcha o'zgartirishlar berilgan koordinatalar bo'yicha amalga oshiriladi - kuchlanish z (3D-versiyada) va ularni hisobga olmasdan (2D-versiyasida).
ANMOD dasturining ishlash protokoli protokol.txt nomli matn fayliga yozilgan (16-tab.).
Muntazam - gidroklorik tarmoqning tugunlarida dastlabki potentsial maydon qiymatlarini tiklashdan tashqari , turli xil hisoblash transformatorlarini, xususan, A g (V z ) uchun tortishish maydonini hisoblashni ta'minlaydi: 1 va 2 vertikal lotin V va V ; gorizontal hosilalar V xz va V yz ; to'liq gorizontal gradient moduli (V x 2 + V 2 ) 05 ; darajadagi sirtning egriligi V A = V - V xx ; tortishish potentsiali V
Bu yerda modellik va oldindan tavsiflovchi ayrim amaliy misollar bor - vakillik algoritmi va dasturi ANMOD.
16-jadval
Do'stlaringiz bilan baham: |
