{"id":373319,"date":"2024-05-21T05:33:25","date_gmt":"2024-05-21T05:33:25","guid":{"rendered":"http:\/\/savepearlharbor.com\/?p=373319"},"modified":"-0001-11-30T00:00:00","modified_gmt":"-0001-11-29T21:00:00","slug":"","status":"publish","type":"post","link":"https:\/\/savepearlharbor.com\/?p=373319","title":{"rendered":"<span>\u041c\u0435\u0442\u0440\u0438\u043a\u0430 \u0432 \u043c\u0430\u0448\u0438\u043d\u043d\u043e\u043c \u043e\u0431\u0443\u0447\u0435\u043d\u0438\u0438 \u0441\u043b\u043e\u0436\u043d\u044b\u0445 \u0441\u0438\u0441\u0442\u0435\u043c, \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c \u0438 \u043f\u0440\u043e\u0433\u0440\u0430\u043c\u043c\u043d\u044b\u0439 \u043a\u043e\u0434<\/span>"},"content":{"rendered":"<div><!--[--><!--]--><\/div>\n<div id=\"post-content-body\">\n<div>\n<div class=\"article-formatted-body article-formatted-body article-formatted-body_version-2\">\n<div xmlns=\"http:\/\/www.w3.org\/1999\/xhtml\">\n<p>\u0414\u043b\u044f \u0440\u0435\u0448\u0435\u043d\u0438\u044f \u0437\u0430\u0434\u0430\u0447 \u043c\u0430\u0448\u0438\u043d\u043d\u043e\u0433\u043e \u043e\u0431\u0443\u0447\u0435\u043d\u0438\u044f \u043f\u0440\u0435\u0434\u043b\u0430\u0433\u0430\u0435\u0442\u0441\u044f \u043c\u0435\u0442\u0440\u0438\u043a\u0430, \u043e\u0441\u043d\u043e\u0432\u0430\u043d\u043d\u0430\u044f \u043d\u0430 \u0444\u043e\u0440\u043c\u0443\u043b\u0435 \u043e\u0442\u043d\u043e\u0448\u0435\u043d\u0438\u044f \u0441\u0438\u0433\u043d\u0430\u043b\u0430 \u043a \u0448\u0443\u043c\u0443, SNR:<\/p>\n<figure class=\"full-width\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/habrastorage.org\/r\/w1560\/getpro\/habr\/upload_files\/364\/f4e\/3a4\/364f4e3a41723f85cbf507bd23be5ce4.png\" width=\"749\" height=\"371\" data-src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/364\/f4e\/3a4\/364f4e3a41723f85cbf507bd23be5ce4.png\"\/><\/figure>\n<p>\u0414\u043b\u044f \u043f\u0440\u0438\u043b\u043e\u0436\u0435\u043d\u0438\u0439 \u0432 \u0441\u0444\u0435\u0440\u0435 \u0438\u043d\u0444\u043e\u0440\u043c\u0430\u0446\u0438\u043e\u043d\u043d\u044b\u0445 \u0442\u0435\u0445\u043d\u043e\u043b\u043e\u0433\u0438\u0439, \u0438\u0437\u043c\u0435\u0440\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0439 \u0442\u0435\u0445\u043d\u0438\u043a\u0435, \u0433\u0434\u0435 \u0441\u0443\u0449\u0435\u0441\u0442\u0432\u0443\u0435\u0442 \u043f\u0440\u043e\u0431\u043b\u0435\u043c\u0430 \u0432\u044b\u0431\u043e\u0440\u0430 \u043e\u0431\u044a\u0435\u043a\u0442\u0438\u0432\u043d\u043e\u0433\u043e \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u0435\u043d\u043d\u043e\u0433\u043e \u043a\u0440\u0438\u0442\u0435\u0440\u0438\u044f \u043f\u043e\u043c\u0435\u0445\u043e\u0443\u0441\u0442\u043e\u0439\u0447\u0438\u0432\u043e\u0441\u0442\u0438 \u0443\u0441\u0442\u0440\u043e\u0439\u0441\u0442\u0432 \u043f\u0440\u0438\u043c\u0435\u043d\u044f\u0435\u0442\u0441\u044f \u043e\u0442\u043d\u043e\u0448\u0435\u043d\u0438\u0435 \u0441\u0438\u0433\u043d\u0430\u043b\u0430 \u043a \u0448\u0443\u043c\u0443, SNR. <\/p>\n<p>\u041f\u0440\u0435\u0434\u043b\u043e\u0436\u0435\u043d\u043d\u0430\u044f \u0444\u043e\u0440\u043c\u0443\u043b\u0430 SNR \u0438\u043d\u0432\u0430\u0440\u0438\u0430\u043d\u0442\u043d\u0430 \u043e\u0442\u043d\u043e\u0441\u0438\u0442\u0435\u043b\u044c\u043d\u043e \u043b\u044e\u0431\u044b\u0445 \u043b\u0438\u043d\u0435\u0439\u043d\u044b\u0445 \u043f\u0440\u0435\u043e\u0431\u0440\u0430\u0437\u043e\u0432\u0430\u043d\u0438\u0439 \u043d\u0430\u0431\u043e\u0440\u043e\u0432 <strong>a(i)<\/strong> (\u0438\u043d\u0432\u0430\u0440\u0438\u0430\u043d\u0442\u043d\u043e\u0441\u0442\u044c \u043f\u043e \u043c\u0430\u0441\u0448\u0442\u0430\u0431\u0443 \u0438 \u0441\u0434\u0432\u0438\u0433\u0443) \u0438 \u043d\u0435 \u0437\u0430\u0432\u0438\u0441\u0438\u0442 \u043e\u0442 \u0447\u0438\u0441\u043b\u0430 \u0433\u0440\u0430\u043d\u0443\u043b\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u044f <strong>M<\/strong> \u0434\u043b\u044f \u043c\u043d\u043e\u0436\u0435\u0441\u0442\u0432\u0430 \u0442\u043e\u0447\u0435\u043a \u043d\u043e\u0440\u043c\u0430\u043b\u044c\u043d\u043e\u0433\u043e \u0433\u0430\u0443\u0441\u0441\u043e\u0432\u0430 \u0440\u0430\u0441\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u044f. \u041a\u0432\u0430\u0434\u0440\u0430\u0442\u0438\u0447\u043d\u044b\u0435 \u0444\u043e\u0440\u043c\u044b S \u0438 N \u043f\u043e\u043b\u043e\u0436\u0438\u0442\u0435\u043b\u044c\u043d\u043e \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u044b \u0434\u043b\u044f \u043f\u0440\u043e\u0438\u0437\u0432\u043e\u043b\u044c\u043d\u044b\u0445 \u043d\u0430\u0431\u043e\u0440\u043e\u0432 <strong>a(i).<\/strong><\/p>\n<p>\u0412\u044b\u0432\u043e\u0434 \u0444\u043e\u0440\u043c\u0443\u043b\u044b \u043e\u0441\u0443\u0449\u0435\u0441\u0442\u0432\u043b\u044f\u0435\u0442\u0441\u044f \u0438\u0437 \u0444\u0440\u0430\u043a\u0442\u0430\u043b\u0430 \u00ab\u043f\u044b\u043b\u044c \u041a\u0430\u043d\u0442\u043e\u0440\u0430\u00bb. \u0418\u0437 \u0432\u044b\u0432\u043e\u0434\u0430 \u0444\u043e\u0440\u043c\u0443\u043b\u044b \u0438 \u043e\u043f\u044b\u0442\u0430 \u043f\u0440\u0438\u043c\u0435\u043d\u0435\u043d\u0438\u044f, \u0443\u043c\u0435\u0441\u0442\u043d\u043e \u043f\u0440\u0438\u043b\u043e\u0436\u0435\u043d\u0438\u0435 \u0432 \u0441\u043b\u043e\u0436\u043d\u043e\u0439 \u0441\u0438\u0441\u0442\u0435\u043c\u0435, \u043a\u043e\u0433\u0434\u0430 \u0443\u043f\u043e\u0440\u044f\u0434\u043e\u0447\u0435\u043d\u043d\u043e\u0441\u0442\u044c \u043f\u0440\u0438\u0447\u0438\u043d\u043d\u043e-\u0441\u043b\u0435\u0434\u0441\u0442\u0432\u0435\u043d\u043d\u044b\u0445 \u0441\u0432\u044f\u0437\u0435\u0439 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u0432\u043d\u0443\u0442\u0440\u0438 \u0441\u043b\u043e\u0436\u043d\u043e\u0439 \u0441\u0438\u0441\u0442\u0435\u043c\u044b \u043d\u0435 \u0437\u043d\u0430\u0447\u0438\u043c\u0430, \u0430 \u0437\u0430\u043c\u0435\u0442\u043d\u044b \u043d\u0435\u043b\u0438\u043d\u0435\u0439\u043d\u044b\u0435 \u0442\u0435\u043d\u0434\u0435\u043d\u0446\u0438\u0438 \u043a \u0441\u0430\u043c\u043e\u043e\u0440\u0433\u0430\u043d\u0438\u0437\u0430\u0446\u0438\u0438 \u0434\u0430\u043d\u043d\u044b\u0445. <\/p>\n<p>\u041e\u0442\u043d\u043e\u0448\u0435\u043d\u0438\u0435 \u0441\u0438\u0433\u043d\u0430\u043b\u0430 \u043a \u0448\u0443\u043c\u0443 \u043e\u0431\u0440\u0430\u0442\u043d\u043e \u043f\u0440\u043e\u043f\u043e\u0440\u0446\u0438\u043e\u043d\u0430\u043b\u044c\u043d\u043e \u0447\u0438\u0441\u043b\u0443 \u0441\u0442\u0435\u043f\u0435\u043d\u0435\u0439 \u0441\u0432\u043e\u0431\u043e\u0434\u044b \u0441\u043b\u043e\u0436\u043d\u043e\u0439 \u0441\u0438\u0441\u0442\u0435\u043c\u044b. \u0425\u0430\u0440\u0430\u043a\u0442\u0435\u0440\u043d\u044b \u0434\u043b\u044f \u0441\u0438\u0441\u0442\u0435\u043c\u044b \u043c\u0430\u043a\u0441\u0438\u043c\u0430\u043b\u044c\u043d\u044b\u0435 \u0438 \u043c\u0438\u043d\u0438\u043c\u0430\u043b\u044c\u043d\u044b\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f SNR, \u043f\u043e\u0441\u043b\u0435\u0434\u043e\u0432\u0430\u0442\u0435\u043b\u044c\u043d\u043e\u0435 \u0438\u0445 \u0447\u0435\u0440\u0435\u0434\u043e\u0432\u0430\u043d\u0438\u0435. <\/p>\n<p>\u0414\u043b\u044f \u043f\u0440\u0438\u043c\u0435\u043d\u0435\u043d\u0438\u044f \u043c\u0435\u0442\u0440\u0438\u043a\u0438 \u0432 \u043c\u0430\u0448\u0438\u043d\u043d\u043e\u043c \u043e\u0431\u0443\u0447\u0435\u043d\u0438\u0438 \u043a\u0440\u0438\u0442\u0438\u0447\u0435\u0441\u043a\u0438 \u0432\u0430\u0436\u043d\u0430 \u043e\u043f\u0442\u0438\u043c\u0438\u0437\u0430\u0446\u0438\u044f \u043f\u043e \u0432\u0440\u0435\u043c\u0435\u043d\u0438 \u0440\u0430\u0441\u0447\u0451\u0442\u0430. \u0418\u0441\u0445\u043e\u0434\u044f \u0438\u0437 \u043e\u043f\u044b\u0442\u0430, \u043f\u0440\u0435\u0434\u043b\u0430\u0433\u0430\u0435\u0442\u0441\u044f \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u0439 \u043f\u0440\u043e\u0433\u0440\u0430\u043c\u043c\u043d\u044b\u0439 \u043f\u0440\u043e\u0434\u0443\u043a\u0442. <\/p>\n<pre><code class=\"python\"># SNR, \u0430\u0432\u0442\u043e\u0440 \u0410\u043b\u0435\u043a\u0441\u0430\u043d\u0434\u0440 \u041c\u0430\u043a\u0430\u0440\u0438\u043a(\u0441) import numpy as np from numba import njit  @njit(fastmath=True) def calcChisl(m):     n=m.size-1     s=0     for i in range (2, n - 1):         s = s + m[i] * (2 * m[i] - m[i + 2] - m[i - 2])         #'s = s + massivOfElements(i] * (massivOfElements2(i] - massivOfElements(i + 2] - massivOfElements(i - 2])              #'\u0434\u043b\u044f \u0437\u0430\u043c\u043a\u043d\u0443\u0442\u043e\u0433\u043e \u043a\u043e\u043d\u0442\u0443\u0440\u0430     s = s + m[0] * (2 * m[0] - m[2] - m[n - 1])     s = s + m[1] * (2 * m[1] - m[3] - m[n])      s = s + m[n - 1] * (2 * m[n - 1] - m[0] - m[n - 3])     s = s + m[n] * (2 * m[n] - m[1] - m[n - 2])     return s      @njit(fastmath=True) def calcZnam(m):     n=m.size-1     s=0     for i in range (2, n - 1):         s = s + m[i] * (6 * m[i] - 4 * m[i + 1] + m[i + 2] + m[i - 2] - 4 * m[i - 1])         # 's = s + massivOfElements(i] * (massivOfElements6(i] - massivOfElements4(i + 1] + massivOfElements(i + 2] + massivOfElements(i - 2] - massivOfElements4(i - 1])              #'\u0434\u043b\u044f \u0437\u0430\u043c\u043a\u043d\u0443\u0442\u043e\u0433\u043e \u043a\u043e\u043d\u0442\u0443\u0440\u0430     s = s + m[0] * (6 * m[0] - 4 * m[1] + m[2] + m[n - 1] - 4 * m[n])     s = s + m[1] * (6 * m[1] - 4 * m[2] + m[3] + m[n] - 4 * m[0])      s = s + m[n - 1] * (6 * m[n - 1] - 4 * m[n] + m[0] + m[n - 3] - 4 * m[n - 2])     s = s + m[n] * (6 * m[n] - 4 * m[0] + m[1] + m[n - 2] - 4 * m[n - 1])      #        '\u0434\u043b\u044f \u043d\u0435\u0437\u0430\u043c\u043a\u043d\u0443\u0442\u043e\u0433\u043e \u043a\u043e\u043d\u0442\u0443\u0440\u0430     #        's = s + m(1] * (3 * m(1] - 4 * m(2] + m(3])     #        's = s + m(2] * (7 * m(2] - 4 * m(1] + -4 * m(3] + m(4])      #        's = s + m(n - 1] * (7 * m(n - 1] - 4 * m(n - 2] - 4 * m(n] + m(n - 3])     #        's = s + m(n] * (3 * m(n] - 4 * m(n - 1] + m(n - 2])     return s  def calcSNR(arr):     #print(arr[0:10])     z=calcZnam(arr)     #if z!=0:     a=calcChisl(arr)\/z     #else:     #    a=0     if a&lt;0:print ('znam=',znam)     if a&lt;0:print ('Chisl=',znam)     return a      def snrVectorSlow(arr,deslength):     vSNR=np.zeros(arr.size-deslength)      for i in range(arr.size-deslength):         b=arr[i:i+deslength]         # vSNR = np.append(vSNR, SNR.calcSNR(b])         vSNR[i]= calcSNR(b)     return vSNR      @njit(fastmath=True) def snrChislVector(m):     chisl=np.zeros(m.size)     n=m.size-1     for i in range (2, n - 1):         chisl[i] = m[i] * (2 * m[i] - m[i + 2] - m[i - 2])         #'s = s + massivOfElements(i] * (massivOfElements2(i] - massivOfElements(i + 2] - massivOfElements(i - 2])              #'\u0434\u043b\u044f \u0437\u0430\u043c\u043a\u043d\u0443\u0442\u043e\u0433\u043e \u043a\u043e\u043d\u0442\u0443\u0440\u0430     chisl[0] = m[0] * (2 * m[0] - m[2] - m[n - 1])     chisl[1] = m[1] * (2 * m[1] - m[3] - m[n])      chisl[n-1] = m[n - 1] * (2 * m[n - 1] - m[0] - m[n - 3])     chisl[n] = m[n] * (2 * m[n] - m[1] - m[n - 2])     return chisl      @njit(fastmath=True) def snrZnamVector(m):          n=m.size-1          k=m          tempV=np.zeros(m.size)     for i in range (2, n - 1):         tempV[i] = k[i] * (6 * k[i] - 4 * k[i + 1] + k[i + 2] + k[i - 2] - 4 * k[i - 1])          tempV[0] = k[0] * (6 * k[0] - 4 * k[1] + k[2] + k[n - 1] - 4 * k[n])     tempV[1] = k[1] * (6 * k[1] - 4 * k[2] + k[3] + k[n] - 4 * k[0])      tempV[n - 1] = k[n - 1] * (6 * k[n - 1] - 4 * k[n] + k[0] + k[n - 3] - 4 * k[n - 2])     tempV[n] = k[n] * (6 * k[n] - 4 * k[0] + k[1] + k[n - 2] - 4 * k[n - 1])          znam=tempV #*9     return znam  @njit(fastmath=True) def snrVector(arr,deslength):     vSNR=np.zeros(arr.size-deslength)     b=arr[0:deslength]     chisl=snrChislVector(b)     znam=snrZnamVector(b)     vSNR[0]=np.sum(chisl)\/np.sum(znam)          n=deslength-1     #d=arr[0:deslength]     for k in range(1,arr.size-deslength):         if (k%1000==0):print(k)         d=arr[k:k+deslength]         #chisl=np.delete(chisl,0)                  chisl=chisl[1:chisl.size]  #\u0442\u0430\u043a \u0431\u044b\u0441\u0442\u0440\u0435\u0435         chisl = np.append(chisl, 0) # \u043f\u0440\u043e\u0441\u0442\u043e \u0432 \u043a\u043e\u043d\u0435\u0446 \u0434\u043e\u0431\u0430\u0432\u043b\u044f\u0435\u043c \u044d\u043b\u0435\u043c\u0435\u043d\u0442                  chisl[1] = d[1] * (2 * d[1] - d[3] - d[n])         i = n - 2         chisl[i] = d[i] * (2 * d[i] - d[i + 2] - d[i - 2])         chisl[n - 1] = d[n - 1] * (2 * d[n - 1] - d[0] - d[n - 3])         chisl[n] = d[n] * (2 * d[n] - d[1] - d[n - 2])                           znam=znam[1:znam.size]  #\u0442\u0430\u043a \u0431\u044b\u0441\u0442\u0440\u0435\u0435         znam = np.append(znam, 0) # \u043f\u0440\u043e\u0441\u0442\u043e \u0432 \u043a\u043e\u043d\u0435\u0446 \u0434\u043e\u0431\u0430\u0432\u043b\u044f\u0435\u043c \u044d\u043b\u0435\u043c\u0435\u043d\u0442                  znam[0] = d[0] * (6 * d[0] - 4 * d[1] + d[2] + d[n - 1] - 4 * d[n])         znam[1] = d[1] * (6 * d[1] - 4 * d[2] + d[3] + d[n] - 4 * d[0])          znam[i] = d[i] * (6 * d[i] - 4 * d[i + 1] + d[i + 2] + d[i - 2] - 4 * d[i - 1])         znam[n - 1] = d[n - 1] * (6 * d[n - 1] - 4 * d[n] + d[0] + d[n - 3] - 4 * d[n - 2])         znam[n] = d[n] * (6 * d[n] - 4 * d[0] + d[1] + d[n - 2] - 4 * d[n - 1])                  vSNR[k]=np.sum(chisl)\/np.sum(znam)                            return vSNR          @njit(fastmath=True) def snrVectorFastest(arr,deslength):     vSNR=np.zeros(arr.size-deslength)     d=arr[0:deslength]     chisl=np.sum(snrChislVector(d))     znam=np.sum(snrZnamVector(d))          kkk=chisl\/znam          vSNR[0]=kkk          n=deslength-1          v=np.zeros(6)          # chisl[i] = m[i] * (2 * m[i] - m[i + 2] - m[i - 2])          for k in range(1,arr.size-deslength):         # if (k%1000==0):print(k)                  v[0] = d[0] * (2 * d[0] - d[2] - d[n-1])         v[1] = d[1] * (2 * d[1] - d[3] - d[n])         v[2] = d[2] * (2 * d[2] - d[4] - d[0])   #\u0435\u0433\u043e \u043d\u0435 \u043e\u0431\u044f\u0437\u0430\u0442\u0435\u043b\u044c\u043d\u043e...c\u0447\u0438\u0442\u0430\u0442\u044c \u043f\u043e\u0442\u043e\u043c                  i = n - 2         # v[3] = d[i] * (2 * d[i] - d[i + 2] - d[i - 2])  #\u0435\u0433\u043e \u0443\u0434\u0430\u043b\u044f\u0442\u044c \u043d\u0435 \u043d\u0430\u0434\u043e...!!! \u043d\u043e \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u043d\u0430\u0434\u043e         v[4] = d[n - 1] * (2 * d[n - 1] - d[0] - d[n - 3])         v[5] = d[n] * (2 * d[n] - d[1] - d[n - 2])                  chisl=chisl-v[0]-v[1]-v[2]-v[4]-v[5]                           #tempV[i] = k[i] * (6 * k[i] - 4 * k[i + 1] + k[i + 2] + k[i - 2] - 4 * k[i - 1])                  v[0] = d[0] * (6 * d[0] - 4 * d[1] + d[2] + d[n - 1] - 4 * d[n])         v[1] = d[1] * (6 * d[1] - 4 * d[2] + d[3] + d[n] - 4 * d[0])         v[2] = d[2] * (6 * d[2] - 4 * d[3] + d[4] + d[0] - 4 * d[1])   #\u0435\u0433\u043e \u043d\u0435 \u043e\u0431\u044f\u0437\u0430\u0442\u0435\u043b\u044c\u043d\u043e \u0441\u0447\u0438\u0442\u0430\u0442\u044c \u043f\u043e\u0442\u043e\u043c           # v[3] = d[i] * (6 * d[i] - 4 * d[i + 1] + d[i + 2] + d[i - 2] - 4 * d[i - 1])  #\u0435\u0433\u043e \u0443\u0434\u0430\u043b\u044f\u0442\u044c \u043d\u0435 \u043d\u0430\u0434\u043e...\u043d\u043e \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u043d\u0430\u0434\u043e         v[4] = d[n - 1] * (6 * d[n - 1] - 4 * d[n] + d[0] + d[n - 3] - 4 * d[n - 2])         v[5] = d[n] * (6 * d[n] - 4 * d[0] + d[1] + d[n - 2] - 4 * d[n - 1])                  znam=znam-v[0]-v[1]-v[2]-v[4]-v[5]                  d=arr[k:k+deslength]                  v[0] = d[0] * (2 * d[0] - d[2] - d[n-1])         v[1] = d[1] * (2 * d[1] - d[3] - d[n])                  v[3] = d[i] * (2 * d[i] - d[i + 2] - d[i - 2])  #\u0435\u0433\u043e \u0443\u0434\u0430\u043b\u044f\u0442\u044c \u043d\u0435 \u043d\u0430\u0434\u043e...!!! \u043d\u043e \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u043d\u0430\u0434\u043e         v[4] = d[n - 1] * (2 * d[n - 1] - d[0] - d[n - 3])         v[5] = d[n] * (2 * d[n] - d[1] - d[n - 2])                  chisl=chisl+ v[0] + v[1]+v[3]+v[4]+v[5]                           v[0] = d[0] * (6 * d[0] - 4 * d[1] + d[2] + d[n - 1] - 4 * d[n])         v[1] = d[1] * (6 * d[1] - 4 * d[2] + d[3] + d[n] - 4 * d[0])         #v[2] = d[2] * (6 * d[2] - 4 * d[3] + d[4] + d[0] - 4 * d[1])   #\u0435\u0433\u043e \u043d\u0435 \u043e\u0431\u044f\u0437\u0430\u0442\u0435\u043b\u044c\u043d\u043e \u0441\u0447\u0438\u0442\u0430\u0442\u044c \u043f\u043e\u0442\u043e\u043c           v[3] = d[i] * (6 * d[i] - 4 * d[i + 1] + d[i + 2] + d[i - 2] - 4 * d[i - 1])  #\u0435\u0433\u043e \u0443\u0434\u0430\u043b\u044f\u0442\u044c \u043d\u0435 \u043d\u0430\u0434\u043e...\u043d\u043e \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u043d\u0430\u0434\u043e         v[4] = d[n - 1] * (6 * d[n - 1] - 4 * d[n] + d[0] + d[n - 3] - 4 * d[n - 2])         v[5] = d[n] * (6 * d[n] - 4 * d[0] + d[1] + d[n - 2] - 4 * d[n - 1])                  znam=znam + v[0] + v[1]+v[3]+v[4]+v[5]                  vSNR[k]=chisl\/znam                   return vSNR<\/code><\/pre>\n<\/p>\n<\/div>\n<\/div>\n<\/div>\n<p><!----><!----><\/div>\n<p><!----><!----><br \/> \u0441\u0441\u044b\u043b\u043a\u0430 \u043d\u0430 \u043e\u0440\u0438\u0433\u0438\u043d\u0430\u043b \u0441\u0442\u0430\u0442\u044c\u0438 <a href=\"https:\/\/habr.com\/ru\/articles\/806889\/\"> https:\/\/habr.com\/ru\/articles\/806889\/<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<div><!--[--><!--]--><\/div>\n<div id=\"post-content-body\">\n<div>\n<div class=\"article-formatted-body article-formatted-body article-formatted-body_version-2\">\n<div xmlns=\"http:\/\/www.w3.org\/1999\/xhtml\">\n<p>\u0414\u043b\u044f \u0440\u0435\u0448\u0435\u043d\u0438\u044f \u0437\u0430\u0434\u0430\u0447 \u043c\u0430\u0448\u0438\u043d\u043d\u043e\u0433\u043e \u043e\u0431\u0443\u0447\u0435\u043d\u0438\u044f \u043f\u0440\u0435\u0434\u043b\u0430\u0433\u0430\u0435\u0442\u0441\u044f \u043c\u0435\u0442\u0440\u0438\u043a\u0430, \u043e\u0441\u043d\u043e\u0432\u0430\u043d\u043d\u0430\u044f \u043d\u0430 \u0444\u043e\u0440\u043c\u0443\u043b\u0435 \u043e\u0442\u043d\u043e\u0448\u0435\u043d\u0438\u044f \u0441\u0438\u0433\u043d\u0430\u043b\u0430 \u043a \u0448\u0443\u043c\u0443, SNR:<\/p>\n<figure class=\"full-width\"><\/figure>\n<p>\u0414\u043b\u044f \u043f\u0440\u0438\u043b\u043e\u0436\u0435\u043d\u0438\u0439 \u0432 \u0441\u0444\u0435\u0440\u0435 \u0438\u043d\u0444\u043e\u0440\u043c\u0430\u0446\u0438\u043e\u043d\u043d\u044b\u0445 \u0442\u0435\u0445\u043d\u043e\u043b\u043e\u0433\u0438\u0439, \u0438\u0437\u043c\u0435\u0440\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0439 \u0442\u0435\u0445\u043d\u0438\u043a\u0435, \u0433\u0434\u0435 \u0441\u0443\u0449\u0435\u0441\u0442\u0432\u0443\u0435\u0442 \u043f\u0440\u043e\u0431\u043b\u0435\u043c\u0430 \u0432\u044b\u0431\u043e\u0440\u0430 \u043e\u0431\u044a\u0435\u043a\u0442\u0438\u0432\u043d\u043e\u0433\u043e \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u0435\u043d\u043d\u043e\u0433\u043e \u043a\u0440\u0438\u0442\u0435\u0440\u0438\u044f \u043f\u043e\u043c\u0435\u0445\u043e\u0443\u0441\u0442\u043e\u0439\u0447\u0438\u0432\u043e\u0441\u0442\u0438 \u0443\u0441\u0442\u0440\u043e\u0439\u0441\u0442\u0432 \u043f\u0440\u0438\u043c\u0435\u043d\u044f\u0435\u0442\u0441\u044f \u043e\u0442\u043d\u043e\u0448\u0435\u043d\u0438\u0435 \u0441\u0438\u0433\u043d\u0430\u043b\u0430 \u043a \u0448\u0443\u043c\u0443, SNR. <\/p>\n<p>\u041f\u0440\u0435\u0434\u043b\u043e\u0436\u0435\u043d\u043d\u0430\u044f \u0444\u043e\u0440\u043c\u0443\u043b\u0430 SNR \u0438\u043d\u0432\u0430\u0440\u0438\u0430\u043d\u0442\u043d\u0430 \u043e\u0442\u043d\u043e\u0441\u0438\u0442\u0435\u043b\u044c\u043d\u043e \u043b\u044e\u0431\u044b\u0445 \u043b\u0438\u043d\u0435\u0439\u043d\u044b\u0445 \u043f\u0440\u0435\u043e\u0431\u0440\u0430\u0437\u043e\u0432\u0430\u043d\u0438\u0439 \u043d\u0430\u0431\u043e\u0440\u043e\u0432 <strong>a(i)<\/strong> (\u0438\u043d\u0432\u0430\u0440\u0438\u0430\u043d\u0442\u043d\u043e\u0441\u0442\u044c \u043f\u043e \u043c\u0430\u0441\u0448\u0442\u0430\u0431\u0443 \u0438 \u0441\u0434\u0432\u0438\u0433\u0443) \u0438 \u043d\u0435 \u0437\u0430\u0432\u0438\u0441\u0438\u0442 \u043e\u0442 \u0447\u0438\u0441\u043b\u0430 \u0433\u0440\u0430\u043d\u0443\u043b\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u044f <strong>M<\/strong> \u0434\u043b\u044f \u043c\u043d\u043e\u0436\u0435\u0441\u0442\u0432\u0430 \u0442\u043e\u0447\u0435\u043a \u043d\u043e\u0440\u043c\u0430\u043b\u044c\u043d\u043e\u0433\u043e \u0433\u0430\u0443\u0441\u0441\u043e\u0432\u0430 \u0440\u0430\u0441\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u0438\u044f. \u041a\u0432\u0430\u0434\u0440\u0430\u0442\u0438\u0447\u043d\u044b\u0435 \u0444\u043e\u0440\u043c\u044b S \u0438 N \u043f\u043e\u043b\u043e\u0436\u0438\u0442\u0435\u043b\u044c\u043d\u043e \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0435\u043d\u044b \u0434\u043b\u044f \u043f\u0440\u043e\u0438\u0437\u0432\u043e\u043b\u044c\u043d\u044b\u0445 \u043d\u0430\u0431\u043e\u0440\u043e\u0432 <strong>a(i).<\/strong><\/p>\n<p>\u0412\u044b\u0432\u043e\u0434 \u0444\u043e\u0440\u043c\u0443\u043b\u044b \u043e\u0441\u0443\u0449\u0435\u0441\u0442\u0432\u043b\u044f\u0435\u0442\u0441\u044f \u0438\u0437 \u0444\u0440\u0430\u043a\u0442\u0430\u043b\u0430 \u00ab\u043f\u044b\u043b\u044c \u041a\u0430\u043d\u0442\u043e\u0440\u0430\u00bb. \u0418\u0437 \u0432\u044b\u0432\u043e\u0434\u0430 \u0444\u043e\u0440\u043c\u0443\u043b\u044b \u0438 \u043e\u043f\u044b\u0442\u0430 \u043f\u0440\u0438\u043c\u0435\u043d\u0435\u043d\u0438\u044f, \u0443\u043c\u0435\u0441\u0442\u043d\u043e \u043f\u0440\u0438\u043b\u043e\u0436\u0435\u043d\u0438\u0435 \u0432 \u0441\u043b\u043e\u0436\u043d\u043e\u0439 \u0441\u0438\u0441\u0442\u0435\u043c\u0435, \u043a\u043e\u0433\u0434\u0430 \u0443\u043f\u043e\u0440\u044f\u0434\u043e\u0447\u0435\u043d\u043d\u043e\u0441\u0442\u044c \u043f\u0440\u0438\u0447\u0438\u043d\u043d\u043e-\u0441\u043b\u0435\u0434\u0441\u0442\u0432\u0435\u043d\u043d\u044b\u0445 \u0441\u0432\u044f\u0437\u0435\u0439 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u0432\u043d\u0443\u0442\u0440\u0438 \u0441\u043b\u043e\u0436\u043d\u043e\u0439 \u0441\u0438\u0441\u0442\u0435\u043c\u044b \u043d\u0435 \u0437\u043d\u0430\u0447\u0438\u043c\u0430, \u0430 \u0437\u0430\u043c\u0435\u0442\u043d\u044b \u043d\u0435\u043b\u0438\u043d\u0435\u0439\u043d\u044b\u0435 \u0442\u0435\u043d\u0434\u0435\u043d\u0446\u0438\u0438 \u043a \u0441\u0430\u043c\u043e\u043e\u0440\u0433\u0430\u043d\u0438\u0437\u0430\u0446\u0438\u0438 \u0434\u0430\u043d\u043d\u044b\u0445. <\/p>\n<p>\u041e\u0442\u043d\u043e\u0448\u0435\u043d\u0438\u0435 \u0441\u0438\u0433\u043d\u0430\u043b\u0430 \u043a \u0448\u0443\u043c\u0443 \u043e\u0431\u0440\u0430\u0442\u043d\u043e \u043f\u0440\u043e\u043f\u043e\u0440\u0446\u0438\u043e\u043d\u0430\u043b\u044c\u043d\u043e \u0447\u0438\u0441\u043b\u0443 \u0441\u0442\u0435\u043f\u0435\u043d\u0435\u0439 \u0441\u0432\u043e\u0431\u043e\u0434\u044b \u0441\u043b\u043e\u0436\u043d\u043e\u0439 \u0441\u0438\u0441\u0442\u0435\u043c\u044b. \u0425\u0430\u0440\u0430\u043a\u0442\u0435\u0440\u043d\u044b \u0434\u043b\u044f \u0441\u0438\u0441\u0442\u0435\u043c\u044b \u043c\u0430\u043a\u0441\u0438\u043c\u0430\u043b\u044c\u043d\u044b\u0435 \u0438 \u043c\u0438\u043d\u0438\u043c\u0430\u043b\u044c\u043d\u044b\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f SNR, \u043f\u043e\u0441\u043b\u0435\u0434\u043e\u0432\u0430\u0442\u0435\u043b\u044c\u043d\u043e\u0435 \u0438\u0445 \u0447\u0435\u0440\u0435\u0434\u043e\u0432\u0430\u043d\u0438\u0435. <\/p>\n<p>\u0414\u043b\u044f \u043f\u0440\u0438\u043c\u0435\u043d\u0435\u043d\u0438\u044f \u043c\u0435\u0442\u0440\u0438\u043a\u0438 \u0432 \u043c\u0430\u0448\u0438\u043d\u043d\u043e\u043c \u043e\u0431\u0443\u0447\u0435\u043d\u0438\u0438 \u043a\u0440\u0438\u0442\u0438\u0447\u0435\u0441\u043a\u0438 \u0432\u0430\u0436\u043d\u0430 \u043e\u043f\u0442\u0438\u043c\u0438\u0437\u0430\u0446\u0438\u044f \u043f\u043e \u0432\u0440\u0435\u043c\u0435\u043d\u0438 \u0440\u0430\u0441\u0447\u0451\u0442\u0430. \u0418\u0441\u0445\u043e\u0434\u044f \u0438\u0437 \u043e\u043f\u044b\u0442\u0430, \u043f\u0440\u0435\u0434\u043b\u0430\u0433\u0430\u0435\u0442\u0441\u044f \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u0439 \u043f\u0440\u043e\u0433\u0440\u0430\u043c\u043c\u043d\u044b\u0439 \u043f\u0440\u043e\u0434\u0443\u043a\u0442. <\/p>\n<pre><code class=\"python\"># SNR, \u0430\u0432\u0442\u043e\u0440 \u0410\u043b\u0435\u043a\u0441\u0430\u043d\u0434\u0440 \u041c\u0430\u043a\u0430\u0440\u0438\u043a(\u0441) import numpy as np from numba import njit  @njit(fastmath=True) def calcChisl(m):     n=m.size-1     s=0     for i in range (2, n - 1):         s = s + m[i] * (2 * m[i] - m[i + 2] - m[i - 2])         #'s = s + massivOfElements(i] * (massivOfElements2(i] - massivOfElements(i + 2] - massivOfElements(i - 2])              #'\u0434\u043b\u044f \u0437\u0430\u043c\u043a\u043d\u0443\u0442\u043e\u0433\u043e \u043a\u043e\u043d\u0442\u0443\u0440\u0430     s = s + m[0] * (2 * m[0] - m[2] - m[n - 1])     s = s + m[1] * (2 * m[1] - m[3] - m[n])      s = s + m[n - 1] * (2 * m[n - 1] - m[0] - m[n - 3])     s = s + m[n] * (2 * m[n] - m[1] - m[n - 2])     return s      @njit(fastmath=True) def calcZnam(m):     n=m.size-1     s=0     for i in range (2, n - 1):         s = s + m[i] * (6 * m[i] - 4 * m[i + 1] + m[i + 2] + m[i - 2] - 4 * m[i - 1])         # 's = s + massivOfElements(i] * (massivOfElements6(i] - massivOfElements4(i + 1] + massivOfElements(i + 2] + massivOfElements(i - 2] - massivOfElements4(i - 1])              #'\u0434\u043b\u044f \u0437\u0430\u043c\u043a\u043d\u0443\u0442\u043e\u0433\u043e \u043a\u043e\u043d\u0442\u0443\u0440\u0430     s = s + m[0] * (6 * m[0] - 4 * m[1] + m[2] + m[n - 1] - 4 * m[n])     s = s + m[1] * (6 * m[1] - 4 * m[2] + m[3] + m[n] - 4 * m[0])      s = s + m[n - 1] * (6 * m[n - 1] - 4 * m[n] + m[0] + m[n - 3] - 4 * m[n - 2])     s = s + m[n] * (6 * m[n] - 4 * m[0] + m[1] + m[n - 2] - 4 * m[n - 1])      #        '\u0434\u043b\u044f \u043d\u0435\u0437\u0430\u043c\u043a\u043d\u0443\u0442\u043e\u0433\u043e \u043a\u043e\u043d\u0442\u0443\u0440\u0430     #        's = s + m(1] * (3 * m(1] - 4 * m(2] + m(3])     #        's = s + m(2] * (7 * m(2] - 4 * m(1] + -4 * m(3] + m(4])      #        's = s + m(n - 1] * (7 * m(n - 1] - 4 * m(n - 2] - 4 * m(n] + m(n - 3])     #        's = s + m(n] * (3 * m(n] - 4 * m(n - 1] + m(n - 2])     return s  def calcSNR(arr):     #print(arr[0:10])     z=calcZnam(arr)     #if z!=0:     a=calcChisl(arr)\/z     #else:     #    a=0     if a&lt;0:print ('znam=',znam)     if a&lt;0:print ('Chisl=',znam)     return a      def snrVectorSlow(arr,deslength):     vSNR=np.zeros(arr.size-deslength)      for i in range(arr.size-deslength):         b=arr[i:i+deslength]         # vSNR = np.append(vSNR, SNR.calcSNR(b])         vSNR[i]= calcSNR(b)     return vSNR      @njit(fastmath=True) def snrChislVector(m):     chisl=np.zeros(m.size)     n=m.size-1     for i in range (2, n - 1):         chisl[i] = m[i] * (2 * m[i] - m[i + 2] - m[i - 2])         #'s = s + massivOfElements(i] * (massivOfElements2(i] - massivOfElements(i + 2] - massivOfElements(i - 2])              #'\u0434\u043b\u044f \u0437\u0430\u043c\u043a\u043d\u0443\u0442\u043e\u0433\u043e \u043a\u043e\u043d\u0442\u0443\u0440\u0430     chisl[0] = m[0] * (2 * m[0] - m[2] - m[n - 1])     chisl[1] = m[1] * (2 * m[1] - m[3] - m[n])      chisl[n-1] = m[n - 1] * (2 * m[n - 1] - m[0] - m[n - 3])     chisl[n] = m[n] * (2 * m[n] - m[1] - m[n - 2])     return chisl      @njit(fastmath=True) def snrZnamVector(m):          n=m.size-1          k=m          tempV=np.zeros(m.size)     for i in range (2, n - 1):         tempV[i] = k[i] * (6 * k[i] - 4 * k[i + 1] + k[i + 2] + k[i - 2] - 4 * k[i - 1])          tempV[0] = k[0] * (6 * k[0] - 4 * k[1] + k[2] + k[n - 1] - 4 * k[n])     tempV[1] = k[1] * (6 * k[1] - 4 * k[2] + k[3] + k[n] - 4 * k[0])      tempV[n - 1] = k[n - 1] * (6 * k[n - 1] - 4 * k[n] + k[0] + k[n - 3] - 4 * k[n - 2])     tempV[n] = k[n] * (6 * k[n] - 4 * k[0] + k[1] + k[n - 2] - 4 * k[n - 1])          znam=tempV #*9     return znam  @njit(fastmath=True) def snrVector(arr,deslength):     vSNR=np.zeros(arr.size-deslength)     b=arr[0:deslength]     chisl=snrChislVector(b)     znam=snrZnamVector(b)     vSNR[0]=np.sum(chisl)\/np.sum(znam)          n=deslength-1     #d=arr[0:deslength]     for k in range(1,arr.size-deslength):         if (k%1000==0):print(k)         d=arr[k:k+deslength]         #chisl=np.delete(chisl,0)                  chisl=chisl[1:chisl.size]  #\u0442\u0430\u043a \u0431\u044b\u0441\u0442\u0440\u0435\u0435         chisl = np.append(chisl, 0) # \u043f\u0440\u043e\u0441\u0442\u043e \u0432 \u043a\u043e\u043d\u0435\u0446 \u0434\u043e\u0431\u0430\u0432\u043b\u044f\u0435\u043c \u044d\u043b\u0435\u043c\u0435\u043d\u0442                  chisl[1] = d[1] * (2 * d[1] - d[3] - d[n])         i = n - 2         chisl[i] = d[i] * (2 * d[i] - d[i + 2] - d[i - 2])         chisl[n - 1] = d[n - 1] * (2 * d[n - 1] - d[0] - d[n - 3])         chisl[n] = d[n] * (2 * d[n] - d[1] - d[n - 2])                           znam=znam[1:znam.size]  #\u0442\u0430\u043a \u0431\u044b\u0441\u0442\u0440\u0435\u0435         znam = np.append(znam, 0) # \u043f\u0440\u043e\u0441\u0442\u043e \u0432 \u043a\u043e\u043d\u0435\u0446 \u0434\u043e\u0431\u0430\u0432\u043b\u044f\u0435\u043c \u044d\u043b\u0435\u043c\u0435\u043d\u0442                  znam[0] = d[0] * (6 * d[0] - 4 * d[1] + d[2] + d[n - 1] - 4 * d[n])         znam[1] = d[1] * (6 * d[1] - 4 * d[2] + d[3] + d[n] - 4 * d[0])          znam[i] = d[i] * (6 * d[i] - 4 * d[i + 1] + d[i + 2] + d[i - 2] - 4 * d[i - 1])         znam[n - 1] = d[n - 1] * (6 * d[n - 1] - 4 * d[n] + d[0] + d[n - 3] - 4 * d[n - 2])         znam[n] = d[n] * (6 * d[n] - 4 * d[0] + d[1] + d[n - 2] - 4 * d[n - 1])                  vSNR[k]=np.sum(chisl)\/np.sum(znam)                            return vSNR          @njit(fastmath=True) def snrVectorFastest(arr,deslength):     vSNR=np.zeros(arr.size-deslength)     d=arr[0:deslength]     chisl=np.sum(snrChislVector(d))     znam=np.sum(snrZnamVector(d))          kkk=chisl\/znam          vSNR[0]=kkk          n=deslength-1          v=np.zeros(6)          # chisl[i] = m[i] * (2 * m[i] - m[i + 2] - m[i - 2])          for k in range(1,arr.size-deslength):         # if (k%1000==0):print(k)                  v[0] = d[0] * (2 * d[0] - d[2] - d[n-1])         v[1] = d[1] * (2 * d[1] - d[3] - d[n])         v[2] = d[2] * (2 * d[2] - d[4] - d[0])   #\u0435\u0433\u043e \u043d\u0435 \u043e\u0431\u044f\u0437\u0430\u0442\u0435\u043b\u044c\u043d\u043e...c\u0447\u0438\u0442\u0430\u0442\u044c \u043f\u043e\u0442\u043e\u043c                  i = n - 2         # v[3] = d[i] * (2 * d[i] - d[i + 2] - d[i - 2])  #\u0435\u0433\u043e \u0443\u0434\u0430\u043b\u044f\u0442\u044c \u043d\u0435 \u043d\u0430\u0434\u043e...!!! \u043d\u043e \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u043d\u0430\u0434\u043e         v[4] = d[n - 1] * (2 * d[n - 1] - d[0] - d[n - 3])         v[5] = d[n] * (2 * d[n] - d[1] - d[n - 2])                  chisl=chisl-v[0]-v[1]-v[2]-v[4]-v[5]                           #tempV[i] = k[i] * (6 * k[i] - 4 * k[i + 1] + k[i + 2] + k[i - 2] - 4 * k[i - 1])                  v[0] = d[0] * (6 * d[0] - 4 * d[1] + d[2] + d[n - 1] - 4 * d[n])         v[1] = d[1] * (6 * d[1] - 4 * d[2] + d[3] + d[n] - 4 * d[0])         v[2] = d[2] * (6 * d[2] - 4 * d[3] + d[4] + d[0] - 4 * d[1])   #\u0435\u0433\u043e \u043d\u0435 \u043e\u0431\u044f\u0437\u0430\u0442\u0435\u043b\u044c\u043d\u043e \u0441\u0447\u0438\u0442\u0430\u0442\u044c \u043f\u043e\u0442\u043e\u043c           # v[3] = d[i] * (6 * d[i] - 4 * d[i + 1] + d[i + 2] + d[i - 2] - 4 * d[i - 1])  #\u0435\u0433\u043e \u0443\u0434\u0430\u043b\u044f\u0442\u044c \u043d\u0435 \u043d\u0430\u0434\u043e...\u043d\u043e \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u043d\u0430\u0434\u043e         v[4] = d[n - 1] * (6 * d[n - 1] - 4 * d[n] + d[0] + d[n - 3] - 4 * d[n - 2])         v[5] = d[n] * (6 * d[n] - 4 * d[0] + d[1] + d[n - 2] - 4 * d[n - 1])                  znam=znam-v[0]-v[1]-v[2]-v[4]-v[5]                  d=arr[k:k+deslength]                  v[0] = d[0] * (2 * d[0] - d[2] - d[n-1])         v[1] = d[1] * (2 * d[1] - d[3] - d[n])                  v[3] = d[i] * (2 * d[i] - d[i + 2] - d[i - 2])  #\u0435\u0433\u043e \u0443\u0434\u0430\u043b\u044f\u0442\u044c \u043d\u0435 \u043d\u0430\u0434\u043e...!!! \u043d\u043e \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u043d\u0430\u0434\u043e         v[4] = d[n - 1] * (2 * d[n - 1] - d[0] - d[n - 3])         v[5] = d[n] * (2 * d[n] - d[1] - d[n - 2])                  chisl=chisl+ v[0] + v[1]+v[3]+v[4]+v[5]                           v[0] = d[0] * (6 * d[0] - 4 * d[1] + d[2] + d[n - 1] - 4 * d[n])         v[1] = d[1] * (6 * d[1] - 4 * d[2] + d[3] + d[n] - 4 * d[0])         #v[2] = d[2] * (6 * d[2] - 4 * d[3] + d[4] + d[0] - 4 * d[1])   #\u0435\u0433\u043e \u043d\u0435 \u043e\u0431\u044f\u0437\u0430\u0442\u0435\u043b\u044c\u043d\u043e \u0441\u0447\u0438\u0442\u0430\u0442\u044c \u043f\u043e\u0442\u043e\u043c           v[3] = d[i] * (6 * d[i] - 4 * d[i + 1] + d[i + 2] + d[i - 2] - 4 * d[i - 1])  #\u0435\u0433\u043e \u0443\u0434\u0430\u043b\u044f\u0442\u044c \u043d\u0435 \u043d\u0430\u0434\u043e...\u043d\u043e \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u043d\u0430\u0434\u043e         v[4] = d[n - 1] * (6 * d[n - 1] - 4 * d[n] + d[0] + d[n - 3] - 4 * d[n - 2])         v[5] = d[n] * (6 * d[n] - 4 * d[0] + d[1] + d[n - 2] - 4 * d[n - 1])                  znam=znam + v[0] + v[1]+v[3]+v[4]+v[5]                  vSNR[k]=chisl\/znam                   return vSNR<\/code><\/pre>\n<\/p>\n<\/div>\n<\/div>\n<\/div>\n<p><!----><!----><\/div>\n<p><!----><!----><br \/> \u0441\u0441\u044b\u043b\u043a\u0430 \u043d\u0430 \u043e\u0440\u0438\u0433\u0438\u043d\u0430\u043b \u0441\u0442\u0430\u0442\u044c\u0438 <a href=\"https:\/\/habr.com\/ru\/articles\/806889\/\"> https:\/\/habr.com\/ru\/articles\/806889\/<\/a><br \/><\/br><\/br><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[],"tags":[],"class_list":["post-373319","post","type-post","status-publish","format-standard","hentry"],"_links":{"self":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/373319","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=373319"}],"version-history":[{"count":0,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/373319\/revisions"}],"wp:attachment":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=373319"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=373319"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=373319"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}