рд╕рд░рд▓ рдкреНрд░рд╡рд╛рд╣ рдореЙрдбрд▓рд┐рдВрдЧ

рд╕рд┐рд░реНрдл рджреЛ рд╢рдмреНрджреЛрдВ рдореЗрдВ, рдЬреИрд╕реА рдЪреАрдЬреЗрдВ



рдФрд░ рднреА рдмрд╣реБрдд рдХреБрдЫред рдпрд╣ рд╕рдм рдПрдХ "рд╕рд░рд▓ рдзрд╛рд░рд╛" рдХрд╣рд╛ рдЬрд╛рддрд╛ рд╣реИред

рдЬрд╛рд╣рд┐рд░ рд╣реИ, рд╕реНрдЯреЗрд╢рди рдкрд░ рдХреЙрд▓ рдХреА рд╕рдВрдЦреНрдпрд╛, рд░рд╛рдЬрдорд╛рд░реНрдЧ рдкрд░ рдХрд╛рд░реЗрдВ рдкреНрд░рддреНрдпреЗрдХ рдмрд┐рдВрджреБ рдкрд░ рд╕рдордп рдкрд░ рдпрд╛рджреГрдЪреНрдЫрд┐рдХ рд╣реЛрддреА рд╣реИрдВред рдпрд╣ рджреЗрдЦрддреЗ рд╣реБрдП рдХрд┐ рд╕реВрдЪреА рд╕реЗ рд╕рднреА рддреАрди рдЪреАрдЬреЗрдВ рдкреВрд░реЗ рд╕рдорд╛рдЬ рдХреЗ рд▓рд┐рдП рдФрд░ рдЙрд╕рдХреЗ рдЗрд╖реНрдЯрддрдо рдЕрд╕реНрддрд┐рддреНрд╡ рдХреЗ рд▓рд┐рдП рдХрд╛рдлреА рдорд╣рддреНрд╡рдкреВрд░реНрдг рд╣реИрдВ, рдпрд╣ рд╕реАрдЦрдирд╛ рдЖрд╡рд╢реНрдпрдХ рд╣реИ рдХрд┐ рдЗрд╕ рдкреНрд░рд╡рд╛рд╣ рдХреЛ рдХреИрд╕реЗ рдореЙрдбрд▓ рдХрд┐рдпрд╛ рдЬрд╛рдПред

рд╡рд┐рд╢реЗрд╖ рд░реВрдк рд╕реЗ рдкреНрд░рд╛рд╕рдВрдЧрд┐рдХ рд╕реВрдЪреА рдкрд░ рдЕрдВрддрд┐рдо рдЖрдЗрдЯрдо рд╣реИ, рдХреНрдпреЛрдВрдХрд┐ рдЬреЛ рдЙрдкрдпреЛрдЧрдХрд░реНрддрд╛ рдХрдВрдкреНрдпреВрдЯрд░ рдХреА рд╕рд░рд▓ рдореВрд▓ рдмрд╛рддреЗрдВ рдирд╣реАрдВ рдЬрд╛рдирддреЗ рд╣реИрдВ рд╡реЗ рдЙрдиреНрд╣реЗрдВ рдЕрд╡рд╛рд╕реНрддрд╡рд┐рдХ рдкреНрд░рд╛рдкреНрдд рдХрд░рдиреЗ рдХрд╛ рдкреНрд░рдмрдВрдзрди рдХрд░рддреЗ рд╣реИрдВ, рдФрд░ рдкреНрд░рд╡рд╛рд╣ рдорд╛рдкрджрдВрдбреЛрдВ рдХреЛ рдЬрд╛рдирддреЗ рд╣реБрдП, рдЖрдк рдЙрдирдХреА рд╕рдВрдЦреНрдпрд╛ рдХреА рднрд╡рд┐рд╖реНрдпрд╡рд╛рдгреА рдХрд░ рд╕рдХрддреЗ рд╣реИрдВред

рд╕рдмрд╕реЗ рд╕рд░рд▓ рдкреНрд░рд╡рд╛рд╣ рдХреЗ рдЧрдгрд┐рддреАрдп рдореЙрдбрд▓ рдХреЛ рдЗрд╕рдореЗрдВ рд╕реЗ рдХрдИ рдЧреБрдгреЛрдВ рдХреА рдЖрд╡рд╢реНрдпрдХрддрд╛ рд╣реЛрддреА рд╣реИред

  1. Stationarityред рдЗрд╕рдХрд╛ рдорддрд▓рдм рдпрд╣ рд╣реИ рдХрд┐ рдПрдХ рдирд┐рд╢реНрдЪрд┐рдд рд╕рдордп рдкрд░ рдХрд┐рд╕реА рд╡рд┐рд╢реЗрд╖ рд╕рдВрдЦреНрдпрд╛ рдореЗрдВ рд╣реЛрдиреЗ рд╡рд╛рд▓реА рдШрдЯрдирд╛рдУрдВ рдХреА рд╕рдВрднрд╛рд╡рдирд╛ рдХреЗрд╡рд▓ рдЕрдиреБрднрд╛рдЧ рдХреА рд▓рдВрдмрд╛рдИ рдкрд░ рдирд┐рд░реНрднрд░ рдХрд░рддреА рд╣реИ рдФрд░ рдпрд╣ рдирд┐рд░реНрднрд░ рдирд╣реАрдВ рдХрд░рддреА рд╣реИ рдХрд┐ рдпрд╣ рдЦрдВрдб рдХрд╣рд╛рдБ рд╕реНрдерд┐рдд рд╣реИред
  2. рдЕрднрд╛рд╡ рдХреЗ рдмрд╛рджред рдЗрд╕рдХрд╛ рдЕрд░реНрде рдирд┐рдореНрди рд╣реИ: рдпрджрд┐ рджреЛ рд╕рдордп рдЕрдВрддрд░рд╛рд▓ рдПрдХ рджреВрд╕рд░реЗ рдХреЛ рдирд╣реАрдВ рдХрд╛рдЯрддреЗ рд╣реИрдВ, рддреЛ рдЙрдирдореЗрдВ рд╕реЗ рдПрдХ рдкрд░ рдЧрд┐рд░рдиреЗ рд╡рд╛рд▓реА рдШрдЯрдирд╛рдУрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ рдХрд┐рд╕реА рдЕрдиреНрдп рд╕рд╛рдЗрдЯ рдкрд░ рдЖрдиреЗ рд╡рд╛рд▓реА рдШрдЯрдирд╛рдУрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ рдкрд░ рдирд┐рд░реНрднрд░ рдирд╣реАрдВ рдХрд░рддреА рд╣реИред
  3. рд╕рд░реНрд╡рд╕рд╛рдзрд╛рд░рдгред рдПрдХ рдкреНрд░рд╡рд╛рд╣ рд╕рд╛рдорд╛рдиреНрдп рд╣реИ рдпрджрд┐ рджреЛ рдпрд╛ рджреЛ рд╕реЗ рдЕрдзрд┐рдХ рдШрдЯрдирд╛рдУрдВ рдХреЗ рдПрдХ рдкрд░реНрдпрд╛рдкреНрдд рдЫреЛрдЯреЗ рд╕рдордп рдЕрдВрддрд░рд╛рд▓ рдореЗрдВ рдЧрд┐рд░рдиреЗ рдХреА рд╕рдВрднрд╛рд╡рдирд╛ рдирдЧрдгреНрдп рд╣реИред


рдзреНрдпрд╛рди рджреЗрдВ рдХрд┐ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рдШрдЯрдирд╛рдУрдВ рдХреЛ рд╕рдордп рдЕрдВрддрд░рд╛рд▓ рдореЗрдВ рдХреИрд╕реЗ рд╡рд┐рднрд╛рдЬрд┐рдд рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдЗрд╕ рдкрд░ рдирд┐рд░реНрднрд░ рдХрд░рддреЗ рд╣реБрдП, рдпрд╣ рдирд┐рд░реНрднрд░ рдХрд░ рд╕рдХрддрд╛ рд╣реИ рдХрд┐ рдпрд╣ рдзрд╛рд░рд╛ рд╕рдмрд╕реЗ рд╕рд░рд▓ рд╣реИ рдпрд╛ рдирд╣реАрдВред

рдЙрджрд╛рд╣рд░рдг рдХреЗ рд▓рд┐рдП, рдкреНрд░рддрд┐ рджрд┐рди рд░рд╛рдЬрдорд╛рд░реНрдЧ рдкрд░ рдХрд╛рд░реЛрдВ рдХрд╛ рдкреНрд░рд╡рд╛рд╣ рд╕рдмрд╕реЗ рд╕рд░рд▓ рдирд╣реАрдВ рд╣реИ, рдХреНрдпреЛрдВрдХрд┐ рдпрд╣ рд╕реНрдерд┐рд░ рдирд╣реАрдВ рд╣реИ: рджрд┐рди рдХреЗ рджреМрд░рд╛рди рд╕реБрдмрд╣ рдФрд░ рд╢рд╛рдо рдХреА рддреБрд▓рдирд╛ рдореЗрдВ рдХрдо рдХрд╛рд░реЗрдВ рд╣реИрдВ, рдЬрдмрдХрд┐ рднреАрдбрд╝ рдХреЗ рдШрдВрдЯреЗ рдореЗрдВ рд╕рдорд╛рди рдкреНрд░рд╡рд╛рд╣ рдкрд╣рд▓реЗ рд╕реЗ рд╣реА рд╕рдмрд╕реЗ рд╕рд░рд▓ рдорд╛рдирд╛ рдЬрд╛ рд╕рдХрддрд╛ рд╣реИред

рдпрд╣рд╛рдВ рдореБрдЦреНрдп рдмрд╛рдд рдпрд╣ рд╣реИ рдХрд┐ рд╣рдорд╛рд░реЗ рд▓рд┐рдП рдмреНрдпрд╛рдЬ рдХреА рдорд╛рддреНрд░рд╛ рдЕрд╕рддрдд рдкреЙрдЗрд╕рди рд╡рд┐рддрд░рдг рдХрд╛ рдкрд╛рд▓рди рдХрд░рддреА рд╣реИред

рдПрд▓реНрдЧреЛрд░рд┐рдереНрдо


рд╣рдореЗрдВ рдкреИрд░рд╛рдореАрдЯрд░ рдХрд╛ рдорд╛рди рджрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ - рддрдХрдиреАрдХреА рд╕рд╣рд╛рдпрддрд╛ рдХреЗ рд▓рд┐рдП рдХреЙрд▓ рдХреА рдФрд╕рдд рд╕рдВрдЦреНрдпрд╛, рдЙрджрд╛рд╣рд░рдг рдХреЗ рд▓рд┐рдП, рдФрд░ рд╣рдореЗрдВ рдПрдХ рдпрд╛рджреГрдЪреНрдЫрд┐рдХ рдЪрд░ X рдЙрддреНрдкрдиреНрди рдХрд░рдирд╛ рд╣реЛрдЧрд╛ рдЬреЛ рд╡рд┐рддрд░рдг рдХрд╛рдиреВрди рдХрд╛ рдкрд╛рд▓рди рдХрд░рддрд╛ рд╣реИ Po (a)

рд╕реА # рднрд╛рд╖рд╛ рдХреЗ рд░реИрдВрдбрдо рд╡рд░реНрдЧ рдореЗрдВ, рдЬрд┐рд╕реЗ рдореИрдВрдиреЗ рдПрдХ рдЙрджрд╛рд╣рд░рдг рдХреЗ рд░реВрдк рдореЗрдВ рдЪреБрдирд╛, рд╡рд╣рд╛рдБ рдирдореВрдирд╛ () рд╡рд┐рдзрд┐ рд╣реИ, рдЬреЛ 0 рд╕реЗ 1 рддрдХ рдПрдХ рдпрд╛рджреГрдЪреНрдЫрд┐рдХ рд╕рдВрдЦреНрдпрд╛ рдЙрддреНрдкрдиреНрди рдХрд░рддреА рд╣реИ, рд╕рдорд╛рди рд░реВрдк рд╕реЗ рдЗрд╕ рдЦрдВрдб рдкрд░ рд╡рд┐рддрд░рд┐рдд рдХреА рдЬрд╛рддреА рд╣реИред

рддреЛ рд╣рдо рдПрдХ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛ рд╕реЗ рдкреВрд░реНрдгрд╛рдВрдХ рдХреИрд╕реЗ рдкреНрд░рд╛рдкреНрдд рдХрд░реЗрдВ? рдЧрдгрд┐рддреАрдп рдЧрдгрдирд╛рдУрдВ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рдирд╛ рдЬреЛ рд╕рдВрднрд╛рд╡реНрдпрддрд╛ рд╕рд┐рджреНрдзрд╛рдВрдд рдореЗрдВ рд╡рд┐рд╢реЗрд╖рдЬреНрдЮреЛрдВ рдХреЗ рджреГрд╖реНрдЯрд┐рдХреЛрдг рд╕реЗ рдмрд╣реБрдд рд╕рд░рд▓ рд╣реИрдВ, рдирд┐рдореНрдирд▓рд┐рдЦрд┐рдд рдХрдерди рд╕рд╛рдмрд┐рдд рд╣реЛ рд╕рдХрддрд╛ рд╣реИред

рдорд╛рди рд▓реАрдЬрд┐рдП рдХрд┐ рд╣рдорд╛рд░реЗ рдкрд╛рд╕ рдЕрдВрддрд░рд╛рд▓ рдкрд░ рд╕рдорд╛рди рд░реВрдк рд╕реЗ рд╡рд┐рддрд░рд┐рдд рдПрдХ рдирд┐рд╢реНрдЪрд┐рдд рдпрд╛рджреГрдЪреНрдЫрд┐рдХ рд╕рдВрдЦреНрдпрд╛ рд╣реИ [0; 1]ред рдлрд┐рд░ рд╕рдмрд╕реЗ рдЫреЛрдЯреА рд╕рдВрдЦреНрдпрд╛ X, рдЬреИрд╕реЗ рдХрд┐ рд╣рдорд╛рд░реЗ рдпрд╛рджреГрдЪреНрдЫрд┐рдХ рд╕рдВрдЦреНрдпрд╛рдУрдВ рдХрд╛ рдЙрддреНрдкрд╛рдж X, exp (-a) рд╕реЗ рдХрдбрд╝рд╛рдИ рд╕реЗ рдХрдо рд╣реИ, рдПрдХ рдкреИрд░рд╛рдореАрдЯрд░ рдХреЗ рд╕рд╛рде Poisson рд╡рд┐рддрд░рдг рдХрд╛ рдкрд╛рд▓рди рдХрд░рддрд╛ рд╣реИ

рдкрд░рд┐рдгрд╛рдорд╕реНрд╡рд░реВрдк, рд╣рдореЗрдВ рдирд┐рдореНрдирд▓рд┐рдЦрд┐рдд рдПрд▓реНрдЧреЛрд░рд┐рдереНрдо рдорд┐рд▓рддреЗ рд╣реИрдВ:

class Poisson: Random { public static uint Generate(double a){ uint X = 0; double Prod = 1; double U = base.Sample(); Prod *= U; while (Prod >= Math.exp(-a)) { X++; U = base.Next(); Prod *= u; } return X; } } } 


рдЗрд╕ рдПрд▓реНрдЧреЛрд░рд┐рдереНрдо рдХреА рдХрдореНрдкреНрдпреВрдЯреЗрд╢рдирд▓ рдЬрдЯрд┐рд▓рддрд╛ рд╣реЗ (рдП) рд╣реИ, рд╣рд╛рд▓рд╛рдВрдХрд┐ рдпрд╣рд╛рдВ рд╣рдореЗрдВ рдмрд╛рд░-рдмрд╛рд░ рдПрдХ рдпрд╛рджреГрдЪреНрдЫрд┐рдХ рдЪрд░ рдЙрддреНрдкрдиреНрди рдХрд░рдирд╛ рд╣реЛрдЧрд╛, рдЬрдм, рд╡рд╛рд╕реНрддрд╡ рдореЗрдВ, рдЖрдк рд╕рд┐рд░реНрдл рдПрдХ рдХрд░ рд╕рдХрддреЗ рд╣реИрдВ:

  public static uint Generate(double a){ uint X = 0; double Prod = Math.Exp(-a); double Sum = Prod; double U = base.Sample(); while (U > Sum) { X++; Prod *= a/Convert.ToDouble(X); Sum += Prod; } return X; } 


рдЗрд╕ рдПрд▓реНрдЧреЛрд░рд┐рдереНрдо рдХреА рдЬрдЯрд┐рд▓рддрд╛ рднреА рдУ (рдП) рд╣реИред

Source: https://habr.com/ru/post/In204324/


All Articles