تباين ألان

أسهل طريقة لاختبار الساعة هي مقارنتها بساعة مرجعية أكثر دقة . خلال فترة زمنية τ ، كما تقيسها الساعة المرجعية، تتقدم الساعة قيد الاختبار بمقدار τy ، حيث y هو متوسط ​​تردد الساعة (النسبي) خلال تلك الفترة. إذا قمنا بقياس فترتين متتاليتين كما هو موضح، يمكننا الحصول على قيمة ( y - y ) ² ، وتشير القيمة الأصغر إلى ساعة أكثر استقرارًا ودقة. إذا كررنا هذه العملية عدة مرات، فإن متوسط ​​قيمة ( y - y ) ² يساوي ضعف تباين ألان (أو مربع انحراف ألان) لفترة الملاحظة τ .

يُعدّ تباين ألان ( AVAR )، المعروف أيضًا باسم تباين عينتين ، مقياسًا لاستقرار التردد في الساعات والمذبذبات والمضخمات . سُمّي هذا التباين نسبةً إلى ديفيد دبليو ألان ، ويُعبّر عنه رياضيًا على النحو التالي :σy2(τ){\displaystyle \sigma _{y}^{2}(\tau )}الانحراف المعياري لآلان ( ADEV )، المعروف أيضًا باسم سيجما تاو ، هو الجذر التربيعي لتباين آلان.σy(τ){\displaystyle \sigma _{y}(\tau )}.

يُعد تباين العينات المتعددة مقياسًا لاستقرار التردد باستخدام M عينة، والوقت T بين القياسات، ووقت الملاحظة.τ{\displaystyle \tau }يُعبّر عن تباين العينة M على النحو التالي :

σy2(م،تي،τ).{\displaystyle \sigma _{y}^{2}(M,T,\tau ).}

يُستخدم تباين ألان لتقدير الاستقرار الناتج عن عمليات التشويش، وليس عن الأخطاء أو العيوب المنهجية مثل انحراف التردد أو تأثيرات درجة الحرارة. يصف كل من تباين ألان وانحراف ألان استقرار التردد. انظر أيضًا قسم " تفسير القيمة" أدناه.

توجد أيضًا تعديلات أو تغييرات مختلفة على تباين ألان، أبرزها تباين ألان المُعدَّل (MAVAR أو MVAR)، والتباين الكلي [ 1 ] ، وتباين هادامارد [ 2 ] . كما توجد متغيرات تتعلق بالاستقرار الزمني، مثل الانحراف الزمني (TDEV) أو التباين الزمني (TVAR). وقد أثبت تباين ألان ومتغيراته فائدته خارج نطاق قياس الوقت ، وهو عبارة عن مجموعة من الأدوات الإحصائية المُحسَّنة التي تُستخدم عندما لا تكون عمليات الضوضاء مستقرة بشكل مطلق، وبالتالي توجد مشتقة له.

يظل تباين العينة M العام مهمًا، لأنه يسمح بفترة التوقف في القياسات، وتسمح دوال الانحياز بالتحويل إلى قيم تباين ألان. ومع ذلك، بالنسبة لمعظم التطبيقات، فإن الحالة الخاصة لعينتين، أو "تباين ألان" معتي=τ{\displaystyle T=\tau }وهو ذو أهمية بالغة.

مثال على رسم بياني لانحراف ألان لساعة. عند فترة رصد قصيرة جدًا (τ) ، يكون انحراف ألان مرتفعًا بسبب التشويش. ومع زيادة فترة الرصد (τ) ، ينخفض ​​الانحراف لأن التشويش يتلاشى. وعند فترات رصد أطول ( τ) ، يبدأ انحراف ألان بالارتفاع مجددًا، مما يشير إلى أن تردد الساعة ينحرف تدريجيًا نتيجة لتغيرات درجة الحرارة، أو تقادم المكونات، أو عوامل أخرى مماثلة. وتزداد أشرطة الخطأ مع زيادة ( τ) نظرًا لأن الحصول على عدد كبير من نقاط البيانات يستغرق وقتًا طويلًا عند فترات رصد طويلة (τ) .
مخطط انحراف ألان كدالة لزمن المتوسط، يوضح الأنظمة الخمسة النموذجية. [ 3 ] 1. ضوضاء تعديل الطور البيضاء/الوميضية (PM): عند أعلى تردد، تسود ضوضاء الطور. وهذا يتوافق معσ(τ)τ-1{\displaystyle \sigma (\tau )\propto \tau ^{-1}}ومع ذلك، فإن رئيس الوزراء الأبيض لديهS[و]=و3{\displaystyle S[f]=f^{3}}لكن برنامج Flicker PM لديهS[و]=و2{\displaystyle S[f]=f^{2}}لا يميز مخطط تباين ألان بينهما، بل يتطلب مخطط تباين ألان معدلاً للتمييز بينهما. ٢. ضوضاء التضمين الترددي الأبيض (FM): عند الترددات المنخفضة، تسود الضوضاء البيضاء في نطاق الترددات المنخفضة. وهذا يتوافق معσ(τ)τ-1/2،S[و]=و0{\displaystyle \sigma (\tau )\propto \tau ^{-1/2},S[f]=f^{0}} 3. فليكر إف إم:σ(τ)τ0،S[و]و-1{\displaystyle \sigma (\tau )\propto \tau ^{0},S[f]\propto f^{-1}}يُطلق على هذا أيضًا اسم "الضوضاء الوردية". 4. راديو المشي العشوائي FM:σ(τ)τ+1/2،S[و]و-2{\displaystyle \sigma (\tau )\propto \tau ^{+1/2},S[f]\propto f^{-2}}يُطلق على هذا أيضًا اسم "الضوضاء البنية" أو "الضوضاء البراونية". في هذا النظام، يتحرك تردد النظام بشكل عشوائي. بعبارة أخرى،دو/دت{\displaystyle df/dt}يتحول إلى ضوضاء بيضاء. 5. انحراف التردد:σ(τ)τ+1،S[و]و-3{\displaystyle \sigma (\tau )\propto \tau ^{+1},S[f]\propto f^{-3}}في هذا النظام، يسلك تردد النظام مسارًا يشبه حركة الضوضاء الوردية. بعبارة أخرى،دو/دت{\displaystyle df/dt}يتحول إلى ضوضاء وردية.

خلفية

عند دراسة استقرار المذبذبات البلورية والساعات الذرية ، تبيّن أنها لا تحتوي على ضوضاء طورية تتكون من ضوضاء بيضاء فقط ، بل أيضًا من ضوضاء تردد الوميض . تُشكّل هذه الأشكال من الضوضاء تحديًا للأدوات الإحصائية التقليدية، مثل الانحراف المعياري ، حيث لا يتقارب المُقدِّر. ولذلك، يُقال إن الضوضاء متباعدة. وقد شملت الجهود المبكرة في تحليل الاستقرار كلاً من التحليل النظري والقياسات العملية. [ 4 ] [ 5 ]

من أهم النتائج الجانبية لوجود هذا النوع من التشويش، عدم توافق طرق القياس المختلفة، مما حال دون تحقيق عنصر أساسي هو قابلية تكرار القياس. وهذا يحدّ من إمكانية مقارنة المصادر ووضع مواصفات دقيقة يمكن طلبها من الموردين. ونتيجةً لذلك، اقتصرت جميع الاستخدامات العلمية والتجارية تقريبًا على قياسات مخصصة، والتي كان من المأمول أن تلبي احتياجات كل تطبيق.

لمعالجة هذه المشكلات، قدّم ديفيد آلان مفهوم تباين M عينة، وبشكل غير مباشر، مفهوم تباين عينتين. [ 6 ] ورغم أن تباين عينتين لم يسمح بتمييز جميع أنواع الضوضاء بشكل كامل، إلا أنه وفّر وسيلةً لفصل العديد من أشكال الضوضاء بشكلٍ فعّال في السلاسل الزمنية لقياسات الطور أو التردد بين مذبذبين أو أكثر. وقدّم آلان طريقةً للتحويل بين أي تباين M عينة وأي تباين N عينة عبر تباين عينتين الشائع، مما جعل جميع تباينات M عينة قابلةً للمقارنة. كما أثبتت آلية التحويل أن تباين M عينة لا يتقارب عند قيم M الكبيرة ، مما يجعله أقل فائدة. وفي وقت لاحق، حدّد معهد مهندسي الكهرباء والإلكترونيات (IEEE) تباين عينتين باعتباره المقياس المُفضّل. [ 7 ]

كان أحد الشواغل المبكرة يتعلق بأجهزة قياس الزمن والتردد التي تتضمن فترة توقف بين القياسات. لم تُشكل هذه السلسلة من القياسات رصدًا متواصلًا للإشارة، مما أدى إلى إدخال انحياز منهجي في القياس. وقد بُذلت عناية فائقة في تقدير هذه الانحيازات. أدى إدخال عدادات بدون فترة توقف إلى الاستغناء عن هذه الانحيازات، ولكن أدوات تحليل الانحياز أثبتت جدواها.

كان أحد الجوانب المبكرة الأخرى التي أثارت القلق يتعلق بكيفية تأثير عرض نطاق جهاز القياس على القياس، مما استدعى أخذ ذلك في الاعتبار. وقد تبين لاحقًا أنه من خلال تغيير الملاحظة خوارزميًاτ{\displaystyle \tau }، منخفض فقطτ{\displaystyle \tau }ستتأثر القيم، بينما لن تتأثر القيم الأعلى. تغييرτ{\displaystyle \tau }يتم ذلك بجعله مضاعفًا صحيحًان{\displaystyle n}قاعدة زمنية للقياسτ0{\displaystyle \tau _{0}}:

τ=نτ0.{\displaystyle \tau =n\tau _{0}.}

قام دي بي ليسون بتحليل فيزياء المذبذبات البلورية [ 5 ] ، وتُعرف النتيجة الآن باسم معادلة ليسون . يؤدي التغذية الراجعة في المذبذب إلى تحويل الضوضاء البيضاء وضوضاء الوميض لمضخم التغذية الراجعة والبلورة إلى ضوضاء تخضع لقانون القوة.و-2{\displaystyle f^{-2}}ضوضاء التردد الأبيض وو-3{\displaystyle f^{-3}}ضوضاء تردد الوميض على التوالي. تؤدي هذه الأشكال من الضوضاء إلى عدم تقارب مُقدِّر التباين القياسي عند معالجة عينات خطأ الوقت. كانت آلية عمل مذبذبات التغذية الراجعة غير معروفة عند بدء العمل على استقرار المذبذب، ولكن قدمها ليسون في نفس الوقت الذي أتاح فيه ديفيد دبليو آلان مجموعة الأدوات الإحصائية . لمزيد من التفاصيل حول تأثير ليسون ، انظر إلى الأدبيات الحديثة حول ضوضاء الطور. [ 8 ]

تفسير القيمة

يُعرَّف تباين ألان بأنه نصف المتوسط ​​الزمني لمربعات الفروق بين القراءات المتتالية لانحراف التردد المأخوذة خلال فترة أخذ العينات. يعتمد تباين ألان على الفترة الزمنية المستخدمة بين العينات، لذا فهو دالة لفترة أخذ العينات، والتي يُرمز لها عادةً بالرمز τ ، وكذلك التوزيع المقاس، ويُعرض بيانيًا بدلًا من رقم واحد. يُعد انخفاض تباين ألان سمةً مميزةً للساعة ذات الاستقرار الجيد خلال الفترة المقاسة.

يُستخدم الانحراف المعياري لآلان على نطاق واسع في الرسوم البيانية (عادةً بصيغة لوغاريتمية ) وعرض الأرقام. وهو مفضل لأنه يُعطي استقرار السعة النسبية، مما يُسهّل المقارنة مع مصادر الأخطاء الأخرى.

يُفسَّر انحراف ألان البالغ 1.3 × 10⁻⁹ عند زمن الرصد 1 ثانية (أي τ = 1 ثانية) على أنه عدم استقرار في التردد بين رصدين يفصل بينهما ثانية واحدة ، بقيمة جذر متوسط ​​تربيعي نسبي (RMS) تبلغ 1.3 × 10⁻⁹ . بالنسبة لساعة بتردد 10 ميجاهرتز، يُعادل هذا حركة جذر متوسط ​​تربيعي قدرها 13 ملي هرتز. إذا لزم تحديد استقرار طور المذبذب، فينبغي الرجوع إلى متغيرات الانحراف الزمني واستخدامها. 

يمكن تحويل تباين ألان وغيره من التباينات في المجال الزمني إلى مقاييس في المجال الترددي للزمن (الطور) واستقرار التردد. [ 9 ]

التركيبات

تباين العينة M

بالنظر إلى سلسلة زمنيةx(ت){\displaystyle x(t)}، لأي أعداد حقيقية موجبةتي،τ{\displaystyle T,\tau }، تعريف متتالية الأعداد الحقيقيةy¯أنا=x(أناتي+τ)-x(أناتي)τأنا=0،1،2،...{\displaystyle {\bar {y}}_{i}={\frac {x(iT+\tau )-x(iT)}{\tau }}\quad i=0,1,2,...}ثمم{\displaystyle M}يُعرَّف تباين العينة [ 6 ] (هنا بصيغة تدوين حديثة) بأنه التباين المصحح وفقًا لـ Bessel للمتتاليةy¯0،...،y¯م-1{\displaystyle {\bar {y}}_{0},...,{\bar {y}}_{M-1}}:σy2(م،تي،τ)=مم-1(1مأنا=0م-1y¯أنا2-[1مأنا=0م-1y¯أنا]2)،{\displaystyle \sigma _{y}^{2}(M,T,\tau )={\frac {M}{M-1}}\left({\frac {1}{M}}\sum _{i=0}^{M-1}{\bar {y}}_{i}^{2}-\left[{\frac {1}{M}}\sum _{i=0}^{M-1}{\bar {y}}_{i}\right]^{2}\right),}تفسير الرموز هو كما يلي:

  • ت{\displaystyle t}هي القراءة على ساعة مرجعية (بوحدات اعتباطية).
  • x(ت){\displaystyle x(t)}هي قراءة الساعة التي نختبرها (بوحدات اعتباطية)، كدالة لقراءة الساعة المرجعية. ويمكن تفسيرها أيضاً على أنها متوسط ​​سلسلة الترددات الكسرية .
  • y¯ن{\displaystyle {\bar {y}}_{n}}هو متوسط ​​التردد الجزئي رقم n خلال فترة الرصدτ{\displaystyle \tau }.
  • م{\displaystyle M}يمثل عدد فترات قراءة الساعة المستخدمة في حسابم{\displaystyle M}تباين العينة،
  • تي{\displaystyle T}هو الوقت الفاصل بين كل عينة تردد،
  • τ{\displaystyle \tau }يمثل طول الفترة الزمنية لكل تقدير للتردد، أو فترة الملاحظة.

يمكن احتساب وقت التوقف عن طريق ترك الوقتتي{\displaystyle T}يختلف عن ذلك الخاص بـτ{\displaystyle \tau }.

تباين ألان

يُعرَّف تباين ألان على النحو التالي:

σy2(τ)=σy2(2،τ،τ)=12(y¯ن+1-y¯ن)2=12τ2(xن+2-2xن+1+xن)2{\displaystyle \sigma _{y}^{2}(\tau )=\left\langle \sigma _{y}^{2}(2,\tau ,\tau )\right\rangle ={\frac {1}{2}}\left\langle \left({\bar {y}}_{n+1}-{\bar {y}}_{n}\right)^{2}\right\rangle ={\frac {1}{2\tau ^{2}}}\left\langle \left(x_{n+2}-2x_{n+1}+x_{n}\right)^{2}\right\rangle }

أينxن:=x(نτ){\displaystyle x_{n}:=x(n\tau )}و{\displaystyle \langle \dotsm \rangle }يشير إلى عامل التوقع.

الحالةتي=τ{\textstyle T=\tau }وهذا يعني أنه يتم أخذ العينات دون وجود فترة توقف بينها.

انحراف ألان

كما هو الحال مع الانحراف المعياري والتباين ، يُعرَّف انحراف ألان بأنه الجذر التربيعي لتباين ألان:

σy(τ)=σy2(τ).{\displaystyle \sigma _{y}(\tau )={\sqrt {\sigma _{y}^{2}(\tau )}}.}

التعريفات الداعمة

نموذج المذبذب

يُفترض أن المذبذب الذي يتم تحليله يتبع النموذج الأساسي لـ

V(ت)=V0الخطيئة(Φ(ت)).{\displaystyle V(t)=V_{0}\sin(\Phi (t)).}

يُفترض أن يكون للمذبذب تردد اسمي قدرهνن{\displaystyle \nu _{\text{n}}}، مُقاسة بالدورات في الثانية (وحدة النظام الدولي للوحدات: هرتز ). التردد الزاوي الاسميωن{\displaystyle \omega _{\text{n}}}(بالراديان لكل ثانية) يُعطى بواسطة

ωن=2πνن.{\displaystyle \omega _{\text{n}}=2\pi \nu _{\text{n}}.}

يمكن فصل المرحلة الكلية إلى مكون دوري تمامًاωنت{\displaystyle \omega _{\text{n}}t}بالإضافة إلى عنصر متذبذبφ(ت){\displaystyle \varphi (t)}:

Φ(ت)=ωنت+φ(ت)=2πνنت+φ(ت).{\displaystyle \Phi (t)=\omega _{\text{n}}t+\varphi (t)=2\pi \nu _{\text{n}}t+\varphi (t).}

خطأ زمني

دالة الخطأ الزمني x ( t ) هي الفرق بين الوقت الاسمي المتوقع والوقت الطبيعي الفعلي:

x(ت)=φ(ت)2πνن=Φ(ت)2πνن-ت=تي(ت)-ت.{\displaystyle x(t)={\frac {\varphi (t)}{2\pi \nu _{\text{n}}}}={\frac {\Phi (t)}{2\pi \nu _{\text{n}}}}-t=T(t)-t.}

بالنسبة للقيم المقاسة ، تُعرَّف سلسلة خطأ الوقت TE( t ) من دالة الوقت المرجعي T ref ( t ) على النحو التالي:

تيهـ(ت)=تي(ت)-تيمرجع(ت).{\displaystyle TE(t)=T(t)-T_{\text{ref}}(t).}

دالة التردد

دالة الترددν(ت){\displaystyle \nu (t)}هو التردد بمرور الوقت، ويُعرَّف على النحو التالي:

ν(ت)=12πدΦ(ت)دت.{\displaystyle \nu (t)={\frac {1}{2\pi }}{\frac {d\Phi (t)}{dt}}.}

التردد الكسري

التردد الكسري y ( t ) هو الفرق المعياري بين الترددν(ت){\displaystyle \nu (t)}والتردد الاسميνن{\displaystyle \nu _{\text{n}}}:

y(ت)=ν(ت)-νنνن=ν(ت)νن-1.{\displaystyle y(t)={\frac {\nu (t)-\nu _{\text{n}}}{\nu _{\text{n}}}}={\frac {\nu (t)}{\nu _{\text{n}}}}-1.}

متوسط ​​التردد الجزئي

يُعرَّف متوسط ​​التردد الجزئي على النحو التالي:

y¯(ت،τ)=1τ0τy(ت+تv)دتv،{\displaystyle {\bar {y}}(t,\tau )={\frac {1}{\tau }}\int _{0}^{\tau }y(t+t_{v})\,dt_{v},}

حيث يتم حساب المتوسط ​​على مدى وقت الملاحظة τ ، و y ( t ) هو خطأ التردد الجزئي في الوقت t ، و τ هو وقت الملاحظة.

بما أن y ( t ) هي مشتقة x ( t )، فيمكننا دون فقدان للعمومية إعادة كتابتها على النحو التالي:

y¯(ت،τ)=x(ت+τ)-x(ت)τ.{\displaystyle {\bar {y}}(t,\tau )={\frac {x(t+\tau )-x(t)}{\tau }}.}

المقدرون

يستند هذا التعريف إلى القيمة المتوقعة الإحصائية ، بالتكامل على مدى زمني غير محدود. لا يسمح الواقع العملي بمثل هذه السلاسل الزمنية، وفي هذه الحالة، يلزم استخدام مُقدِّر إحصائي بديل. سيتم عرض ومناقشة عدد من المُقدِّرات المختلفة.

الاتفاقيات

  • يُرمز إلى عدد عينات التردد في سلسلة التردد الكسري بالرمز M.
  • يُرمز إلى عدد عينات الخطأ الزمني في سلسلة الخطأ الزمني بالرمز N. وتُحدد العلاقة بين عدد عينات التردد الجزئي وسلسلة الخطأ الزمني في العلاقة التالية:
    شمال=م+1.{\displaystyle N=M+1.}
  • بالنسبة لسلسلة عينات الخطأ الزمني ، يرمز xᵢ إلى العينة رقم i من دالة الزمن المستمر x ( t ) كما هو موضح في
    xأنا=x(أناتي)،{\displaystyle x_{i}=x(iT),}

    حيث T هو الوقت بين القياسات. بالنسبة لتباين ألان، يتم ضبط الوقت المستخدم T على وقت الملاحظة τ .

    لنفترض أن N يمثل عدد العينات ( x0 ... xN - 1 ) في سلسلة عينات الخطأ الزمني . ويستخدم الاصطلاح التقليدي الأرقام من 1 إلى N.
  • بالنسبة لسلاسل العينات ذات التردد الجزئي المتوسط ،y¯أنا{\displaystyle {\bar {y}}_{i}}يرمز إلى العينة رقم i من دالة التردد الكسري المستمر المتوسط ​​y ( t ) كما هو موضح في
    y¯أنا=y¯(تيأنا،τ)،{\displaystyle {\bar {y}}_{i}={\bar {y}}(Ti,\tau ),}

    مما يعطي

    y¯أنا=1τ0τy(أناتي+تv)دتv=x(أناتي+τ)-x(أناتي)τ.{\displaystyle {\bar {y}}_{i}={\frac {1}{\tau }}\int _{0}^{\tau }y(iT+t_{v})\,dt_{v}={\frac {x(iT+\tau )-x(iT)}{\tau }}.}

    بافتراض تباين ألان حيث T هو τ، يصبح

    y¯أنا=xأنا+1-xأناτ.{\displaystyle {\bar {y}}_{i}={\frac {x_{i+1}-x_{i}}{\tau }}.}

    لنفترض أن M يمثل عدد العينات في سلسلة العينات ذات التردد الجزئي المتوسط ​​(y¯0...y¯م-1{\displaystyle {\bar {y}}_{0}\ldots {\bar {y}}_{M-1}}) في السلسلة. يستخدم الاصطلاح التقليدي الفهرس من 1 إلى M.

    كاختصار، يُكتب متوسط ​​التردد الكسري غالبًا بدون علامة المتوسط ​​فوقه. مع ذلك، هذا غير صحيح من الناحية الرسمية، لأن التردد الكسري ومتوسط ​​التردد الكسري دالتان مختلفتان. جهاز القياس القادر على إنتاج تقديرات التردد دون أي تأخير سيُنتج في الواقع سلسلة زمنية لمتوسط ​​التردد، والتي لا تحتاج إلا إلى تحويلها إلى متوسط ​​التردد الكسري ، ويمكن استخدامها مباشرةً بعد ذلك.
  • يُشار إلى الوقت بين القياسات بالرمز T ، وهو مجموع وقت الملاحظة τ ووقت التوقف.

مقدرات τ الثابتة

يتمثل أحد أساليب التقدير البسيطة الأولى في ترجمة التعريف مباشرةً إلى

σy2(τ،م)=أفار(τ،م)=12(م-1)أنا=0م-2(y¯أنا+1-y¯أنا)2،{\displaystyle \sigma _{y}^{2}(\tau ,M)=\operatorname {AVAR} (\tau ,M)={\frac {1}{2(M-1)}}\sum _{i=0}^{M-2}({\bar {y}}_{i+1}-{\bar {y}}_{i})^{2},}

أو بالنسبة للسلاسل الزمنية:

σy2(τ،شمال)=أفار(τ،شمال)=12τ2(شمال-2)أنا=0شمال-3(xأنا+2-2xأنا+1+xأنا)2.{\displaystyle \sigma _{y}^{2}(\tau ,N)=\operatorname {AVAR} (\tau ,N)={\frac {1}{2\tau ^{2}(N-2)}}\sum _{i=0}^{N-3}(x_{i+2}-2x_{i+1}+x_{i})^{2}.}

مع ذلك، لا توفر هذه الصيغ سوى الحساب لحالة τ = τ 0. وللحساب لقيمة مختلفة لـ τ ، يلزم توفير سلسلة زمنية جديدة.

مقدرات متغير τ غير المتداخلة

بأخذ السلسلة الزمنية وتجاوز n  1 عينة، ستظهر سلسلة زمنية جديدة (أقصر) بفاصل زمني τ 0 بين العينات المتجاورة، ويمكن حساب تباين ألان لها باستخدام المقدرات البسيطة. يمكن تعديل هذه المقدرات لإدخال المتغير الجديد n بحيث لا يلزم إنشاء سلسلة زمنية جديدة، بل يمكن إعادة استخدام السلسلة الزمنية الأصلية لقيم مختلفة من n . تصبح المقدرات

σy2(نτ0،م)=أفار(نτ0،م)=12م-1نأنا=0م-1ن-1(y¯نأنا+ن-y¯نأنا)2{\displaystyle \sigma _{y}^{2}(n\tau _{0},M)=\operatorname {AVAR} (n\tau _{0},M)={\frac {1}{2{\frac {M-1}{n}}}}\sum _{i=0}^{{\frac {M-1}{n}}-1}\left({\bar {y}}_{ni+n}-{\bar {y}}_{ni}\right)^{2}}

معنم-12{\displaystyle n\leq {\frac {M-1}{2}}}،

وبالنسبة للسلاسل الزمنية:

σy2(نτ0،شمال)=أفار(نτ0،شمال)=12ن2τ02(شمال-1ن-1)أنا=0شمال-1ن-2(xنأنا+2ن-2xنأنا+ن+xنأنا)2{\displaystyle \sigma _{y}^{2}(n\tau _{0},N)=\operatorname {AVAR} (n\tau _{0},N)={\frac {1}{2n^{2}\tau _{0}^{2}\left({\frac {N-1}{n}}-1\right)}}\sum _{i=0}^{{\frac {N-1}{n}}-2}\left(x_{ni+2n}-2x_{ni+n}+x_{ni}\right)^{2}}

معنشمال-12{\displaystyle n\leq {\frac {N-1}{2}}}.

تعاني هذه التقديرات من عيب كبير يتمثل في أنها ستتجاهل كمية كبيرة من بيانات العينة، حيث يتم استخدام 1/ n فقط من العينات المتاحة.

مقدرات متغير τ المتداخلة

قدمت تقنيةٌ من ابتكار جيه جيه سنايدر [ 10 ] أداةً محسّنة، حيث تمّ دمج القياسات في n سلسلة متداخلة من السلسلة الأصلية. وقد قدّم هاو وآلان وبارنز [ 11 ] مُقدِّر تباين ألان المتداخل. ويمكن إثبات أن هذا يُكافئ حساب متوسط ​​عينات الزمن أو التردد المُعَيَّر في مجموعات من n عينة قبل المعالجة. ويصبح المُتنبئ الناتج

σy2(نτ0،م)=أفار(نτ0،م)=12ن2(م-2ن+1)ج=0م-2ن(أنا=جج+ن-1yأنا+ن-yأنا)2=12(م-2ن+1)ج=0م-2ن(y¯ج+ن-y¯ج)2،{\displaystyle {\begin{aligned}\sigma _{y}^{2}(n\tau _{0},M)&=\operatorname {AVAR} (n\tau _{0},M)={\frac {1}{2n^{2}(M-2n+1)}}\sum _{j=0}^{M-2n}\left(\sum _{i=j}^{j+n-1}y_{i+n}-y_{i}\right)^{2}\\[5pt]&={\frac {1}{2(M-2n+1)}}\sum _{j=0}^{M-2n}\left({\bar {y}}_{j+n}-{\bar {y}}_{j}\right)^{2},\end{aligned}}}

أو بالنسبة للسلاسل الزمنية:

σy2(نτ0،شمال)=أفار(نτ0،شمال)=12ن2τ02(شمال-2ن)أنا=0شمال-2ن-1(xأنا+2ن-2xأنا+ن+xأنا)2.{\displaystyle \sigma _{y}^{2}(n\tau _{0},N)=\operatorname {AVAR} (n\tau _{0},N)={\frac {1}{2n^{2}\tau _{0}^{2}(N-2n)}}\sum _{i=0}^{N-2n-1}(x_{i+2n}-2x_{i+n}+x_{i})^{2}.}

تتفوق المقدرات المتداخلة بشكل ملحوظ على المقدرات غير المتداخلة، خاصةً مع ازدياد قيمة n وطول السلسلة الزمنية المتوسط. وقد اعتُمدت المقدرات المتداخلة كمقدرات تباين ألان المفضلة في معايير IEEE [ 7 ] و ITU-T [ 12 ] وETSI [ 13 ] للقياسات المتشابهة، كما هو مطلوب لتأهيل الاتصالات.

تباين ألان المعدل

للتغلب على صعوبة فصل تعديل الطور الأبيض عن تعديل طور الوميض باستخدام مُقدِّرات تباين ألان التقليدية، تُقلِّل عملية ترشيح خوارزمية عرض النطاق الترددي بمقدار n . يُجري هذا الترشيح تعديلًا على التعريف والمُقدِّرات، ويُعرِّف الآن تباينًا مُستقلًا يُسمى تباين ألان المُعدَّل . يُعدّ مقياس تباين ألان المُعدَّل مقياسًا لاستقرار التردد، تمامًا كما هو الحال مع تباين ألان.

مقدرات استقرار الوقت

يمكن حساب مقياس إحصائي لثبات الزمن (σ x )، والذي يُطلق عليه غالبًا الانحراف الزمني (TDEV)، من انحراف ألان المُعدَّل (MDEV). يعتمد TDEV على MDEV بدلًا من انحراف ألان الأصلي، لأن MDEV قادر على التمييز بين تعديل طور الضوء الأبيض وتعديل طور الوميض (PM). فيما يلي تقدير التباين الزمني بناءً على تباين ألان المُعدَّل:

σx2(τ)=τ23تعديلσy2(τ)،{\displaystyle \sigma _{x}^{2}(\tau )={\frac {\tau ^{2}}{3}}{\bmod {\sigma }}_{y}^{2}(\tau ),}

وبالمثل بالنسبة لانحراف ألان المعدل إلى الانحراف الزمني :

σx(τ)=τ3تعديلσy(τ).{\displaystyle \sigma _{x}(\tau )={\frac {\tau }{\sqrt {3}}}{\bmod {\sigma }}_{y}(\tau ).}

يتم تطبيع TDEV بحيث يساوي الانحراف الكلاسيكي للجسيمات البيضاء الدقيقة عند ثابت زمني τ = τ₀ . لفهم عامل مقياس التطبيع بين المقاييس الإحصائية، إليك القاعدة الإحصائية ذات الصلة: بالنسبة للمتغيرين العشوائيين المستقلين X و Y ، فإن تباين ( σz² ) مجموع أو فرق ( z = x - y ) هو مجموع مربعات تبايناتهما (σz² = σx² + σy² ) . تباين مجموع أو فرق ( y = x²τ - ) عينتين مستقلتين من متغير عشوائي يساوي ضعف تباين المتغير العشوائي (σy² = 2σx² ) . MDEV هو الفرق الثاني بين قياسات الطور المستقلة ( x ) التي لها تباين ( σx² ) . بما أن الحساب هو الفرق المزدوج، والذي يتطلب ثلاث قياسات طور مستقلة ( x 2 τ − 2 x τ + x )، فإن تباين ألان المعدل (MVAR) هو ثلاثة أضعاف تباينات قياسات الطور. 

مقدرات أخرى

أدت التطورات اللاحقة إلى ظهور أساليب تقدير محسّنة لنفس مقياس الاستقرار، وهو تباين/انحراف التردد، ولكنها تُعرف بأسماء منفصلة مثل تباين هادامارد ، وتباين هادامارد المعدّل ، والتباين الكلي ، والتباين الكلي المعدّل ، والتباين النظري . [ 14 ] وتتميز هذه الأساليب باستخدامها الأمثل للإحصاءات لتحسين حدود الثقة أو قدرتها على التعامل مع الانحراف الخطي للتردد.

فترات الثقة ودرجات الحرية المكافئة

تقوم المقدرات الإحصائية بحساب قيمة تقديرية لسلسلة العينات المستخدمة. قد تنحرف هذه التقديرات عن القيمة الحقيقية، ويُشار إلى نطاق القيم الذي يحتوي على القيمة الحقيقية باحتمالية معينة بفترة الثقة . تعتمد فترة الثقة على عدد المشاهدات في سلسلة العينات، ونوع الضوضاء السائد، والمقدر المستخدم. كما يعتمد عرضها على اليقين الإحصائي الذي يجعل قيم فترة الثقة تشكل نطاقًا محدودًا، وبالتالي اليقين الإحصائي بأن القيمة الحقيقية تقع ضمن هذا النطاق. بالنسبة لمقدرات τ المتغيرة، يُعد مضاعف τ₀ ، متغيرًا أيضًا.

فاصل الثقة

يمكن تحديد فترة الثقة باستخدام توزيع مربع كاي مع درجات حرية df باستخدام توزيع تباين العينة : [ 7 ] [ 11 ]

χ2=dfs2σ2،{\displaystyle \chi ^{2}={\frac {{\text{df}}\,s^{2}}{\sigma ^{2}}},}

حيث هو تباين العينة لتقديرنا، و σ² هي قيمة التباين الحقيقي، و df هي درجات حرية المُقدِّر، ويتم حساب χ² بناءً على التوزيع التراكمي العكسي لـ χ² بدرجات حرية df. لاحتمالية 90%، والتي تغطي النطاق من 5% إلى 95% على منحنى الاحتمالية، يمكن إيجاد الحدين الأعلى والأدنى باستخدام المتباينة .

χ2(0.05)dfs2σ2χ2(0.95)،{\displaystyle \chi ^{2}(0.05)\leq {\frac {{\text{df}}\,s^{2}}{\sigma ^{2}}}\leq \chi ^{2}(0.95),}

والتي تصبح بعد إعادة ترتيبها للتباين الحقيقي

dfs2χ2(0.95)σ2dfs2χ2(0.05).{\displaystyle {\frac {{\text{df}}\,s^{2}}{\chi ^{2}(0.95)}}\leq \sigma ^{2}\leq {\frac {{\text{df}}\,s^{2}}{\chi ^{2}(0.05)}}.}

درجات الحرية الفعالة

تمثل درجات الحرية عدد المتغيرات الحرة القادرة على المساهمة في التقدير. وتختلف درجات الحرية الفعالة باختلاف المُقدِّر ونوع الضوضاء. وقد تم التوصل تجريبياً إلى صيغ المُقدِّر التي تعتمد على N (عدد نقاط العينة الكلية) و n (مضاعف صحيح لـ τ₀ ) : [ 11 ]

درجات حرية تباين ألان
نوع الضوضاءدرجات الحرية
تعديل الطور الأبيض (WPM)df(شمال+1)(شمال-2ن)2(شمال-ن){\displaystyle {\text{df}}\cong {\frac {(N+1)(N-2n)}{2(N-n)}}}
تعديل طور الوميض (FPM)dfخبرة[(lnشمال-12نln(2ن+1)(شمال-1)4)-1/2]{\displaystyle {\text{df}}\cong \exp \left[\left(\ln {\frac {N-1}{2n}}\ln {\frac {(2n+1)(N-1)}{4}}\right)^{-1/2}\right]}
تعديل التردد الأبيض (WFM)df[3(شمال-1)2ن-2(شمال-2)شمال]4ن24ن2+5{\displaystyle {\text{df}}\cong \left[{\frac {3(N-1)}{2n}}-{\frac {2(N-2)}{N}}\right]{\frac {4n^{2}}{4n^{2}+5}}}
تعديل تردد الوميض (FFM)df{2(شمال-2)2.3شمال-4.9ن=15شمال24ن(شمال+3ن)ن2{\displaystyle {\text{df}}\cong {\begin{cases}{\frac {2(N-2)}{2.3N-4.9}}&n=1\\{\frac {5N^{2}}{4n(N+3n)}}&n\geq 2\end{cases}}}
تعديل التردد العشوائي (RWFM)dfشمال-2ن(شمال-1)2-3ن(شمال-1)+4ن2(شمال-3)2{\displaystyle {\text{df}}\cong {\frac {N-2}{n}}{\frac {(N-1)^{2}-3n(N-1)+4n^{2}}{(N-3)^{2}}}}

ضوضاء قانون القوة

يُعالج تباين ألان أنواع الضوضاء ذات قانون القوة المختلفة بشكلٍ متباين، مما يُسهّل تحديدها وتقدير قوتها. ويُرمز إلى عرض نظام القياس (تردد القطع العالي) بالرمز f<sub> H</sub> .

استجابة قانون القوة لتباين ألان
نوع الضوضاء ذو ​​قانون القوةميل ضوضاء الطورميل ضوضاء الترددمعامل القدرةضوضاء الطور Sx(و){\displaystyle S_{x}(f)}تباين ألان σy2(τ){\displaystyle \sigma _{y}^{2}(\tau )}انحراف ألان σy(τ){\displaystyle \sigma _{y}(\tau )}
تعديل الطور الأبيض (WPM)و0=1{\displaystyle f^{0}=1}و2{\displaystyle f^{2}}ح2{\displaystyle h_{2}}1(2π)2ح2{\displaystyle {\frac {1}{(2\pi )^{2}}}h_{2}}3وح4π2τ2ح2{\displaystyle {\frac {3f_{H}}{4\pi ^{2}\tau ^{2}}}h_{2}}3وح2πτح2{\displaystyle {\frac {\sqrt {3f_{H}}}{2\pi \tau }}{\sqrt {h_{2}}}}
تعديل طور الوميض (FPM)و-1{\displaystyle f^{-1}}و1=و{\displaystyle f^{1}=f}ح1{\displaystyle h_{1}}1(2π)2وح1{\displaystyle {\frac {1}{(2\pi )^{2}f}}h_{1}}3[γ+ln(2πوحτ)]-ln24π2τ2ح1{\displaystyle {\frac {3[\gamma +\ln(2\pi f_{H}\tau )]-\ln 2}{4\pi ^{2}\tau ^{2}}}h_{1}}3[γ+ln(2πوحτ)]-ln22πτح1{\displaystyle {\frac {\sqrt {3[\gamma +\ln(2\pi f_{H}\tau )]-\ln 2}}{2\pi \tau }}{\sqrt {h_{1}}}}
تعديل التردد الأبيض (WFM)و-2{\displaystyle f^{-2}}و0=1{\displaystyle f^{0}=1}ح0{\displaystyle h_{0}}1(2π)2و2ح0{\displaystyle {\frac {1}{(2\pi )^{2}f^{2}}}h_{0}}12τح0{\displaystyle {\frac {1}{2\tau }}h_{0}}12τح0{\displaystyle {\frac {1}{\sqrt {2\tau }}}{\sqrt {h_{0}}}}
تعديل تردد الوميض (FFM)و-3{\displaystyle f^{-3}}و-1{\displaystyle f^{-1}}ح-1{\displaystyle h_{-1}}1(2π)2و3ح-1{\displaystyle {\frac {1}{(2\pi )^{2}f^{3}}}h_{-1}}2ln(2)ح-1{\displaystyle 2\ln(2)h_{-1}}2ln(2)ح-1{\displaystyle {\sqrt {2\ln(2)}}{\sqrt {h_{-1}}}}
تعديل التردد العشوائي (RWFM)و-4{\displaystyle f^{-4}}و-2{\displaystyle f^{-2}}ح-2{\displaystyle h_{-2}}1(2π)2و4ح-2{\displaystyle {\frac {1}{(2\pi )^{2}f^{4}}}h_{-2}}2π2τ3ح-2{\displaystyle {\frac {2\pi ^{2}\tau }{3}}h_{-2}}π2τ3ح-2{\displaystyle {\frac {\pi {\sqrt {2\tau }}}{\sqrt {3}}}{\sqrt {h_{-2}}}}

كما ورد في [ 15 ] [ 16 ] وفي الصيغ الحديثة. [ 17 ] [ 18 ]

لا يستطيع تباين ألان التمييز بين WPM وFPM، ولكنه قادر على تمييز أنواع الضوضاء الأخرى ذات قانون القوة. وللتمييز بين WPM وFPM، يجب استخدام تباين ألان المُعدَّل .

تفترض الصيغ المذكورة أعلاه أن

τ12πوح،{\displaystyle \tau \gg {\frac {1}{2\pi f_{H}}},}

وبالتالي، فإن نطاق تردد وقت الرصد أقل بكثير من نطاق تردد الأجهزة. وعندما لا يتحقق هذا الشرط، فإن جميع أشكال الضوضاء تعتمد على نطاق تردد الجهاز.

رسم الخرائط αμ

رسم خرائط تفصيلية لتعديل الطور من الشكل

Sx(و)=14π2حαوα-2=14π2حαوβ،{\displaystyle S_{x}(f)={\frac {1}{4\pi ^{2}}}h_{\alpha }f^{\alpha -2}={\frac {1}{4\pi ^{2}}}h_{\alpha }f^{\beta },}

أين

βα-2،{\displaystyle \beta \equiv \alpha -2,}

أو تعديل التردد بالشكل

Sy(و)=حαوα{\displaystyle S_{y}(f)=h_{\alpha }f^{\alpha }}

إلى تباين ألان من الشكل

σy2(τ)=كαحατμ{\displaystyle \sigma _{y}^{2}(\tau )=K_{\alpha }h_{\alpha }\tau ^{\mu }}

يمكن تبسيطها بشكل كبير من خلال توفير علاقة بين α و μ . كما يتم تقديم علاقة بين α و K α للتسهيل: [ 7 ]

ألان التباين αμ رسم الخرائط
αβμK α
-2-412π23{\displaystyle {\frac {2\pi ^{2}}{3}}}
-1-302ln2{\displaystyle 2\ln 2}
0-2-112{\displaystyle {\frac {1}{2}}}
1-1-23[γ+ln(2πوحτ)]-ln24π2{\displaystyle {\frac {3[\gamma +\ln(2\pi f_{H}\tau )]-\ln 2}{4\pi ^{2}}}}
20-23وح4π2{\displaystyle {\frac {3f_{H}}{4\pi ^{2}}}}

التحويل العام من ضوضاء الطور

إشارة ذات ضوضاء طور طيفيةSφ{\displaystyle S_{\varphi }}يمكن تحويل الوحدات rad 2 /Hz إلى تباين ألان بواسطة [ 18 ]

σy2(τ)=2ν020وبSφ(و)الخطيئة4(πτو)(πτ)2دو.{\displaystyle \sigma _{y}^{2}(\tau )={\frac {2}{\nu _{0}^{2}}}\int _{0}^{f_{b}}S_{\varphi }(f){\frac {\sin ^{4}(\pi \tau f)}{(\pi \tau )^{2}}}\,df.}

استجابة خطية

على الرغم من أن تباين ألان مصمم لتمييز أشكال الضوضاء، إلا أنه يعتمد على بعض الاستجابات الخطية للزمن وليس كلها. وهي موضحة في الجدول التالي:

استجابة خطية لتباين ألان
التأثير الخطيزمن الاستجابةاستجابة الترددتباين ألانانحراف ألان
إزاحة الطورx0{\displaystyle x_{0}}0{\displaystyle 0}0{\displaystyle 0}0{\displaystyle 0}
إزاحة الترددy0ت{\displaystyle y_{0}t}y0{\displaystyle y_{0}}0{\displaystyle 0}0{\displaystyle 0}
الانجراف الخطيدت22{\displaystyle {\frac {Dt^{2}}{2}}}دت{\displaystyle Dt}د2τ22{\displaystyle {\frac {D^{2}\tau ^{2}}{2}}}دτ2{\displaystyle {\frac {D\tau }{\sqrt {2}}}}

وبالتالي، سيساهم الانحراف الخطي في نتيجة الإخراج. عند قياس نظام حقيقي، قد يلزم تقدير الانحراف الخطي أو أي آلية انحراف أخرى وإزالتها من السلسلة الزمنية قبل حساب تباين ألان. [ 17 ]

خصائص مرشح الوقت والتردد

عند تحليل خصائص تباين ألان وما شابهه، ثبتت فائدة النظر في خصائص المرشح على التردد المعياري. بدءًا من تعريف تباين ألان لـ

σy2(τ)=12(y¯أنا+1-y¯أنا)2،{\displaystyle \sigma _{y}^{2}(\tau )={\frac {1}{2}}\left\langle \left({\bar {y}}_{i+1}-{\bar {y}}_{i}\right)^{2}\right\rangle ,}

أين

y¯أنا=1τ0τy(أناτ+ت)دت.{\displaystyle {\bar {y}}_{i}={\frac {1}{\tau }}\int _{0}^{\tau }y(i\tau +t)\,dt.}

استبدال السلسلة الزمنية لـyأنا{\displaystyle y_{i}}مع المتغير المحوّل بواسطة فورييهSy(و){\displaystyle S_{y}(f)}يمكن التعبير عن تباين ألان في مجال التردد على النحو التالي:

σy2(τ)=0Sy(و)2الخطيئة4πτو(πτو)2دو.{\displaystyle \sigma _{y}^{2}(\tau )=\int _{0}^{\infty }S_{y}(f){\frac {2\sin ^{4}\pi \tau f}{(\pi \tau f)^{2}}}\,df.}

وبالتالي فإن دالة التحويل لتباين ألان هي

|حأ(و)|2=2الخطيئة4πτو(πτو)2.{\displaystyle \left\vert H_{A}(f)\right\vert ^{2}={\frac {2\sin ^{4}\pi \tau f}{(\pi \tau f)^{2}}}.}

دوال الانحياز

سيتعرض تباين العينة M ، وتباين ألان في الحالة الخاصة المحددة، لانحياز منهجي يعتمد على اختلاف عدد العينات M واختلاف العلاقة بين T و τ . ولمعالجة هذه الانحيازات، تم تعريف دالتي الانحياز B1 و B2 [ 19 ] ، مما يسمح بالتحويل بين قيم M و T المختلفة .

لا تكفي دوال الانحياز هذه لمعالجة الانحياز الناتج عن دمج M عينة في زمن الملاحظة 0 خلال MT مع توزيع زمن التوقف بين كتل القياس M بدلاً من توزيعه في نهاية القياس. وقد أدى ذلك إلى الحاجة إلى الانحياز B 3. [ 20 ]

يتم تقييم دوال الانحياز لقيمة μ محددة، لذا يلزم إجراء عملية ربط α–μ لشكل الضوضاء السائد كما تم تحديده باستخدام تقنية تحديد الضوضاء . بدلاً من ذلك، [ 6 ] [ 19 ] يمكن استنتاج قيمة μ لشكل الضوضاء السائد من القياسات باستخدام دوال الانحياز.

دالة الانحياز B 1

تربط دالة الانحياز B1 تباين العينة M بتباين العينة 2، مع الحفاظ على ثبات الفترة الزمنية بين القياسات T والزمن اللازم لكل قياس τ . وقد عُرّفت [ 19 ] على النحو التالي :

ب1(شمال،ر،μ)=σy2(شمال،تي،τ)σy2(2،تي،τ)،{\displaystyle B_{1}(N,r,\mu )={\frac {\left\langle \sigma _{y}^{2}(N,T,\tau )\right\rangle }{\left\langle \sigma _{y}^{2}(2,T,\tau )\right\rangle }},}

أين

ر=تيτ.{\displaystyle r={\frac {T}{\tau }}.}

تصبح دالة الانحياز بعد التحليل

ب1(شمال،ر،μ)=1+ن=1شمال-1شمال-نشمال(شمال-1)[2(رن)μ+2-(رن+1)μ+2-|رن-1|μ+2]1+12[2رμ+2-(ر+1)μ+2-|ر-1|μ+2].{\displaystyle B_{1}(N,r,\mu )={\frac {1+\sum _{n=1}^{N-1}{\frac {N-n}{N(N-1)}}\left[2(rn)^{\mu +2}-(rn+1)^{\mu +2}-|rn-1|^{\mu +2}\right]}{1+{\frac {1}{2}}\left[2r^{\mu +2}-(r+1)^{\mu +2}-|r-1|^{\mu +2}\right]}}.}

دالة الانحياز B 2

تربط دالة الانحياز B2 تباين العينتين عند زمن أخذ العينات T بتباين العينتين (تباين ألان)، مع الحفاظ على عدد العينات N = 2 وزمن الملاحظة τ ثابتًا. وقد عُرّفت [ 19 ] على النحو التالي :

ب2(ر،μ)=σy2(2،تي،τ)σy2(2،τ،τ)،{\displaystyle B_{2}(r,\mu )={\frac {\left\langle \sigma _{y}^{2}(2,T,\tau )\right\rangle }{\left\langle \sigma _{y}^{2}(2,\tau ,\tau )\right\rangle }},}

أين

ر=تيτ.{\displaystyle r={\frac {T}{\tau }}.}

تصبح دالة الانحياز بعد التحليل

ب2(ر،μ)=1+12[2رμ+2-(ر+1)μ+2-|ر-1|μ+2]2(1-2μ).{\displaystyle B_{2}(r,\mu )={\frac {1+{\frac {1}{2}}\left[2r^{\mu +2}-(r+1)^{\mu +2}-|r-1|^{\mu +2}\right]}{2\left(1-2^{\mu }\right)}}.}

دالة الانحياز B3

تربط دالة الانحياز B3 تباين العينتين لوقت أخذ العينة MT 0 ووقت الملاحظة 0 بتباين العينتين (تباين ألان) ، وتُعرَّف [ 20 ] على النحو التالي:

ب3(شمال،م،ر،μ)=σy2(شمال،م،تي،τ)σy2(شمال،تي،τ)،{\displaystyle B_{3}(N,M,r,\mu )={\frac {\left\langle \sigma _{y}^{2}(N,M,T,\tau )\right\rangle }{\left\langle \sigma _{y}^{2}(N,T,\tau )\right\rangle }},}

أين

تي=متي0،{\displaystyle T=MT_{0},}
τ=مτ0.{\displaystyle \tau =M\tau _{0}.}

تُعد دالة الانحياز B 3 مفيدة لضبط قيم تقدير المتغير τ غير المتداخلة والمتداخلة بناءً على قياسات وقت التوقف لوقت الملاحظة τ 0 والوقت بين الملاحظات T 0 إلى تقديرات وقت التوقف العادية.

تصبح دالة الانحياز بعد التحليل (في حالة N  =  2)

ب3(2،م،ر،μ)=2م+مF(مر)-ن=1م-1(م-ن)[2F(نر)-F((م+ن)ر)+F((م-ن)ر)]مμ+2[F(ر)+2]،{\displaystyle B_{3}(2,M,r,\mu )={\frac {2M+MF(Mr)-\sum _{n=1}^{M-1}(M-n)\left[2F(nr)-F{\big (}(M+n)r{\big )}+F{\big (}(M-n)r{\big )}\right]}{M^{\mu +2}[F(r)+2]}},}

أين

F(أ)=2أμ+2-(أ+1)μ+2-|أ-1|μ+2.{\displaystyle F(A)=2A^{\mu +2}-(A+1)^{\mu +2}-|A-1|^{\mu +2}.}

دالة الانحياز τ

على الرغم من عدم صياغتها رسميًا، فقد تم استنتاجها بشكل غير مباشر كنتيجة لربط αμ . عند مقارنة مقياسين لتباين ألان لقيم τ مختلفة ، بافتراض وجود نفس الضوضاء السائدة في شكل نفس معامل μ، يمكن تعريف الانحياز على النحو التالي:

بτ(τ1،τ2،μ)=σy2(2،τ2،τ2)σy2(2،τ1،τ1).{\displaystyle B_{\tau }(\tau _{1},\tau _{2},\mu )={\frac {\left\langle \sigma _{y}^{2}(2,\tau _{2},\tau _{2})\right\rangle }{\left\langle \sigma _{y}^{2}(2,\tau _{1},\tau _{1})\right\rangle }}.}

تصبح دالة الانحياز بعد التحليل

بτ(τ1،τ2،μ)=(τ2τ1)μ.{\displaystyle B_{\tau }(\tau _{1},\tau _{2},\mu )=\left({\frac {\tau _{2}}{\tau _{1}}}\right)^{\mu }.}

التحويل بين القيم

لتحويل مجموعة من القياسات إلى أخرى، يمكن تجميع دوال الانحياز B1 وB2 وτ. أولًا، تحوّل الدالة B1 القيمة ( N1 , T1 , τ1 ) إلى (2, T1 , τ1 ) ، ومنها تحوّل الدالة B2 القيمة (2, τ1 , τ1 ) ، وبالتالي تباين ألان عند τ1 . يمكن تحويل مقياس تباين ألان باستخدام دالة الانحياز τ من τ1 إلى τ2 ، ومنها يتم تحويل (2, T2 , τ2 ) باستخدام B2 ، وأخيرًا باستخدام B1 إلى تباين ( N2 , T2 , τ2 ) . يصبح التحويل الكامل كما يلي :

σy2(شمال2،تي2،τ2)=(τ2τ1)μ[ب1(شمال2،ر2،μ)ب2(ر2،μ)ب1(شمال1،ر1،μ)ب2(ر1،μ)]σy2(شمال1،تي1،τ1)،{\displaystyle \left\langle \sigma _{y}^{2}(N_{2},T_{2},\tau _{2})\right\rangle =\left({\frac {\tau _{2}}{\tau _{1}}}\right)^{\mu }\left[{\frac {B_{1}(N_{2},r_{2},\mu )B_{2}(r_{2},\mu )}{B_{1}(N_{1},r_{1},\mu )B_{2}(r_{1},\mu )}}\right]\left\langle \sigma _{y}^{2}(N_{1},T_{1},\tau _{1})\right\rangle ,}

أين

ر1=تي1ر1،{\displaystyle r_{1}={\frac {T_{1}}{r_{1}}},}
ر2=تي2ر2.{\displaystyle r_{2}={\frac {T_{2}}{r_{2}}}.}

وبالمثل، بالنسبة للقياسات المتسلسلة باستخدام M مقطعًا، يصبح الامتداد المنطقي

σy2(شمال2،م2،تي2،τ2)=(τ2τ1)μ[ب3(شمال2،م2،ر2،μ)ب1(شمال2،ر2،μ)ب2(ر2،μ)ب3(شمال1،م1،ر1،μ)ب1(شمال1،ر1،μ)ب2(ر1،μ)]σy2(شمال1،م1،تي1،τ1).{\displaystyle \left\langle \sigma _{y}^{2}(N_{2},M_{2},T_{2},\tau _{2})\right\rangle =\left({\frac {\tau _{2}}{\tau _{1}}}\right)^{\mu }\left[{\frac {B_{3}(N_{2},M_{2},r_{2},\mu )B_{1}(N_{2},r_{2},\mu )B_{2}(r_{2},\mu )}{B_{3}(N_{1},M_{1},r_{1},\mu )B_{1}(N_{1},r_{1},\mu )B_{2}(r_{1},\mu )}}\right]\left\langle \sigma _{y}^{2}(N_{1},M_{1},T_{1},\tau _{1})\right\rangle .}

مشاكل القياس

عند إجراء القياسات لحساب تباين ألان أو انحراف ألان، قد تؤدي عدة مشكلات إلى تدهور دقة القياسات. نتناول هنا التأثيرات الخاصة بتباين ألان، حيث تكون النتائج متحيزة.

حدود عرض نطاق القياس

من المتوقع أن يكون لنظام القياس نطاق ترددي يساوي أو يقل عن معدل نايكويست ، كما هو موضح في نظرية شانون-هارتلي . وكما هو واضح في معادلات ضوضاء قانون القوة، فإن تعديلات الضوضاء البيضاء وضوضاء الوميض تعتمد على تردد الزاوية العلوي.وح{\displaystyle f_{H}}(يُفترض أن هذه الأنظمة مُرشَّحة بترددات منخفضة فقط). بالنظر إلى خصائص مُرشِّح التردد، يتضح جليًا أن الضوضاء منخفضة التردد لها تأثير أكبر على النتيجة. بالنسبة لأنواع ضوضاء تعديل الطور المسطحة نسبيًا (مثل WPM وFPM)، يكون للترشيح أهمية، بينما بالنسبة لأنواع الضوضاء ذات الميل الأكبر، يصبح الحد الأعلى للتردد أقل أهمية، بافتراض أن عرض نطاق نظام القياس واسع نسبيًا.τ{\displaystyle \tau }كما ورد في

τ12πوح.{\displaystyle \tau \gg {\frac {1}{2\pi f_{H}}}.}

عندما لا يتحقق هذا الافتراض، فإن عرض النطاق الترددي الفعالوح{\displaystyle f_{H}}يجب تدوين ذلك بجانب القياس. ينبغي على المهتمين الرجوع إلى NBS TN394. [ 15 ]

لكن إذا تم تعديل عرض نطاق المُقدِّر باستخدام مضاعفات صحيحة لوقت أخذ العيناتنτ0{\displaystyle n\tau _{0}}عندئذٍ، يمكن تقليل تأثير عرض نطاق النظام إلى مستويات ضئيلة. ولتلبية احتياجات الاتصالات، كانت هذه الأساليب ضرورية لضمان قابلية مقارنة القياسات وإتاحة بعض الحرية للموردين لتنفيذ تطبيقات مختلفة. توصية الاتحاد الدولي للاتصالات G.813 [ 21 ] لقياس TDEV.

يمكن التوصية بأن يكون الأولτ0{\displaystyle \tau _{0}}يتم تجاهل الإشارات المتعددة، بحيث تكون غالبية الضوضاء المكتشفة ضمن نطاق تمرير عرض نطاق أنظمة القياس.

أُجريت تطويرات إضافية على تباين ألان لتقليل عرض النطاق الترددي للأجهزة باستخدام البرمجيات. وقد سمح هذا التطوير بمعالجة التشويش المتبقي، ويُشار إلى هذه الطريقة الآن باسم تباين ألان المُعدَّل . يجب عدم الخلط بين تقنية تقليل عرض النطاق الترددي هذه والنسخة المُحسَّنة من تباين ألان المُعدَّل ، والتي تُغيِّر أيضًا عرض نطاق مرشح التنعيم.

الوقت الضائع في القياسات

تتضمن العديد من أجهزة قياس الزمن والتردد مراحل التجهيز، ووقت قاعدة الوقت، ووقت المعالجة، وقد تعيد تشغيل التجهيز. يبدأ وقت التجهيز من لحظة تشغيله حتى وقوع حدث البدء على قناة البدء. تضمن قاعدة الوقت مرور أقل وقت ممكن قبل قبول حدث على قناة الإيقاف كحدث إيقاف. يتم تسجيل عدد الأحداث والوقت المنقضي بين حدث البدء وحدث الإيقاف وعرضهما خلال وقت المعالجة. عند حدوث المعالجة (المعروفة أيضًا بوقت التوقف )، لا يستطيع الجهاز عادةً إجراء قياس آخر. بعد اكتمال المعالجة، يقوم الجهاز في الوضع المستمر بتشغيل دائرة التجهيز مرة أخرى. يصبح الوقت بين حدث الإيقاف وحدث البدء التالي وقتًا ميتًا ، لا تتم خلاله مراقبة الإشارة. يُدخل هذا الوقت الميت انحيازات قياس منهجية، والتي يجب تعويضها للحصول على نتائج دقيقة. بالنسبة لأنظمة القياس هذه، يُشير الوقت T إلى الوقت بين أحداث البدء المتجاورة (وبالتالي القياسات)، بينماτ{\displaystyle \tau }يشير إلى طول القاعدة الزمنية، أي الطول الاسمي بين حدثي البداية والنهاية لأي قياس.

تؤثر فترات التوقف على القياسات بشكل كبير على النتائج، ما دفع العديد من الدراسات إلى تحديد خصائصها بدقة. وقد أدى ظهور عدادات زمن التوقف الصفري إلى الاستغناء عن هذا التحليل. يتميز هذا النوع من العدادات بأن حدث التوقف في إحدى القياسات يُستخدم أيضًا كحدث بدء للحدث التالي. تُنشئ هذه العدادات سلسلة من أزواج الأحداث والطوابع الزمنية، زوج لكل قناة، يفصل بينها قاعدة زمنية. وقد أثبتت هذه القياسات فائدتها في تحليل السلاسل الزمنية.

يمكن تصحيح القياسات التي تُجرى مع وجود فترة توقف باستخدام دالة الانحياز B1 و B2 و B3 . وبالتالي، فإن فترة التوقف بحد ذاتها لا تمنع الوصول إلى تباين ألان، ولكنها تجعل الأمر أكثر تعقيدًا. يجب معرفة فترة التوقف، بحيث يمكن تحديد الفترة الزمنية بين العينات T.

طول القياس والاستخدام الفعال للعينات

عند دراسة تأثير طول N لسلسلة العينة وتأثير المعلمة المتغيرة τ n على فترات الثقة ، قد تصبح فترات الثقة كبيرة جدًا نظرًا لأن درجة الحرية الفعالة قد تصبح صغيرة لبعض تركيبات N و n لشكل الضوضاء السائد (لـ τ هذا ).

قد يكون التأثير هو أن القيمة المقدرة قد تكون أصغر بكثير أو أكبر بكثير من القيمة الحقيقية، مما قد يؤدي إلى استنتاجات خاطئة للنتيجة.

يوصى بما يلي:

  • يتم رسم فترة الثقة جنبًا إلى جنب مع البيانات، بحيث يعرف قارئ الرسم البياني عدم اليقين الإحصائي للقيم.
  • يجب الحفاظ على طول تسلسل العينة (أي عدد العينات N ) بأكبر قدر ممكن لضمان أن تكون فترة الثقة صغيرة على مدى نطاق τ ذي الأهمية.
  • يمكن استخدام مُقدِّرات توفر قيمًا أفضل لدرجات الحرية كبديل لمُقدِّرات تباين ألان، أو كمُكمِّلة لها عندما تتفوق عليها. ومن بين هذه المُقدِّرات، ينبغي النظر في مُقدِّرات التباين الكلي ومُقدِّرات تباين ثيو .
  • إن نطاق τ الذي يتم مسحه بواسطة مضاعف τ 0 n محدود في الطرف العلوي بالنسبة إلى N ، بحيث لا يختلط الأمر على قارئ الرسم البياني بسبب قيم المقدر غير المستقرة للغاية.

نوع الضوضاء السائد

يعتمد عدد كبير من ثوابت التحويل وتصحيحات الانحياز وفترات الثقة على نوع الضوضاء السائد. وللحصول على تفسير صحيح، يجب تحديد نوع الضوضاء السائد لقيمة τ المحددة محل الاهتمام من خلال تحليل الضوضاء. سيؤدي عدم تحديد نوع الضوضاء السائد إلى قيم منحازة. قد تصل بعض هذه الانحيازات إلى عدة مراتب ، مما قد يكون ذا أهمية بالغة.

الانحراف الخطي

لا يتم إلغاء التأثيرات المنهجية على الإشارة إلا جزئيًا. يتم إلغاء انحراف الطور والتردد، ولكن لا يتم إلغاء الانحراف الخطي أو الأشكال الأخرى ذات الدرجة العالية لمنحنيات الطور متعددة الحدود، مما يشكل قيدًا على القياس. يمكن استخدام مطابقة المنحنى وإزالة الانحراف المنهجي. غالبًا ما تكون إزالة الانحراف الخطي كافية. يمكن أيضًا استخدام مقدرات الانحراف الخطي مثل تباين هادامارد . يمكن استخدام مقدر قائم على العزوم لإزالة الانحراف الخطي.

تحيز تقدير أداة القياس

كانت الأجهزة التقليدية تقتصر على قياس الأحداث المفردة أو أزواج الأحداث. وقد أتاح إدخال الأداة الإحصائية المحسّنة للقياسات المتداخلة، التي طورها جيه جيه سنايدر [ 10 ] ، دقةً محسّنةً بشكل كبير في قراءات التردد، مما أدى إلى اختلال التوازن التقليدي بين الأرقام وقاعدة الزمن. ورغم فائدة هذه الأساليب للغرض المقصود منها، فإن استخدام هذه القياسات المُنعّمة لحسابات تباين ألان قد يُعطي انطباعًا خاطئًا بدقة عالية [ 22 ] [ 23 ] [ 24 ] ، ولكن مع زيادة قيمة τ، يتلاشى هذا التأثير تدريجيًا، وتظهر قيم متحيزة في منطقة τ المنخفضة من القياس. هذا التحيز يُعطي قيمًا أقل مما ينبغي، لذا فهو تحيز مفرط في التفاؤل (بافتراض أن القيم المنخفضة هي المرغوبة)، مما يُقلل من جدوى القياس بدلًا من تحسينه. ويمكن عادةً تعطيل هذه الخوارزميات الذكية أو تجاوزها باستخدام وضع الطابع الزمني، وهو الوضع المُفضّل إن كان متاحًا.

القياسات العملية

على الرغم من إمكانية ابتكار العديد من الطرق لقياس تباين ألان، إلا أن مثالاً بسيطاً قد يوضح كيفية إجراء القياسات.

قياس

تُعتبر جميع قياسات تباين ألان بمثابة مقارنة بين ساعتين مختلفتين. لنفترض وجود ساعة مرجعية وجهاز قيد الاختبار ، وكلاهما بتردد اسمي مشترك يبلغ 10 ميجاهرتز. يُستخدم عداد فترات زمنية لقياس الفترة الزمنية بين الحافة الصاعدة للساعة المرجعية (القناة A) والحافة الصاعدة للجهاز قيد الاختبار . 

لضمان قياسات متباعدة بانتظام، يُقسّم تردد الساعة المرجعية لتحديد معدل القياس، مما يُفعّل عداد الفترات الزمنية (مدخل ARM). يمكن أن يكون هذا المعدل 1  هرتز (باستخدام خرج 1 نبضة في الثانية من الساعة المرجعية)، ولكن يمكن استخدام معدلات أخرى مثل 10  هرتز و100  هرتز. وتُحدد سرعة إنجاز عداد الفترات الزمنية للقياس، وإخراج النتيجة، والاستعداد للنبضة التالية، تردد التشغيل.

ثم يصبح استخدام الحاسوب مفيداً لتسجيل سلسلة الفروق الزمنية التي يتم رصدها.

المعالجة اللاحقة

تتطلب السلاسل الزمنية المسجلة معالجة لاحقة لفك التفاف الطور، بحيث يتم توفير خطأ طور مستمر. عند الضرورة، يجب أيضًا تصحيح أخطاء التسجيل والقياس. ينبغي إجراء تقدير الانحراف وإزالته، مع ضرورة تحديد آلية الانحراف وفهمها للمصادر. قد تكون قيود الانحراف في القياسات شديدة، لذا من الضروري ترك المذبذبات تستقر لفترة كافية من خلال تشغيلها.

يمكن بعد ذلك حساب تباين ألان باستخدام المُقدِّرات المُعطاة، ولأغراض عملية، يُفضَّل استخدام مُقدِّر التداخل نظرًا لاستخدامه المُتفوق للبيانات مقارنةً بمُقدِّر عدم التداخل. كما يُمكن استخدام مُقدِّرات أخرى، مثل مُقدِّرات التباين الكلي أو مُقدِّرات ثيو، إذا طُبِّقت تصحيحات الانحياز بحيث تُوفِّر نتائج مُتوافقة مع تباين ألان.

لتكوين المخططات الكلاسيكية، يتم رسم انحراف ألان (الجذر التربيعي لتباين ألان) بتنسيق لوغاريتمي مقابل فترة الملاحظة τ . 

المعدات والبرامج

عادةً ما يكون عداد الفترات الزمنية عدادًا جاهزًا متوفرًا تجاريًا. تشمل العوامل المحددة دقة القياس الفردي، وتذبذب الإشارة، وسرعة القياسات، واستقرار ساعة المرجع. يمكن إجراء جمع البيانات ومعالجتها باستخدام برامج تجارية أو برامج متاحة للعموم. توجد حلول متطورة للغاية توفر القياس والحساب في جهاز واحد.

تاريخ البحث

لقد خضع مجال استقرار التردد للدراسة لفترة طويلة. ومع ذلك، خلال ستينيات القرن الماضي، تبيّن وجود نقص في التعريفات المتسقة. وقد أسفرت ندوة مشتركة بين وكالة ناسا ومعهد مهندسي الكهرباء والإلكترونيات (IEEE) حول الاستقرار قصير المدى في نوفمبر 1964 [ 25 ] عن إصدار عدد خاص في فبراير 1966 من مجلة وقائع معهد مهندسي الكهرباء والإلكترونيات (IEEE) حول استقرار التردد.

جمعت ندوة ناسا-IEEE العديد من المجالات والاستخدامات المتعلقة بالاستقرار قصير المدى وطويل المدى، بمشاركة باحثين مختلفين في أوراق بحثية. واتفقت المقالات وحلقات النقاش على وجود ضوضاء التذبذب الترددي، وعلى الرغبة في التوصل إلى تعريف موحد للاستقرار قصير المدى وطويل المدى.

ظهرت أوراق بحثية مهمة، بما في ذلك أوراق ديفيد ألان، [ 6 ] وجيمس أ. بارنز، [ 26 ] وإل إس كاتلر وسي إل سيرل [ 4 ] ودي بي ليسون، [ 5 ] في وقائع IEEE حول استقرار التردد وساعدت في تشكيل المجال.

تحلل مقالة ديفيد آلان تباين التردد الكلاسيكي لـ M عينة، متناولةً مسألة زمن التوقف بين القياسات بالإضافة إلى دالة انحياز أولية. [ 6 ] على الرغم من أن دالة الانحياز الأولية لآلان تفترض عدم وجود زمن توقف، إلا أن صيغه تتضمن حسابات زمن التوقف. تحلل مقالته حالة M عينة تردد (يُشار إليها بـ N في المقالة) ومُقدِّرات التباين. كما تُقدم تعيين α–μ القياسي، مستندةً بوضوح إلى عمل جيمس بارنز [ 26 ] في العدد نفسه.

تُعدّ حالة التباين لعينتين حالةً خاصةً من حالة التباين لـ M عينة، والتي تُنتج متوسط ​​مشتق التردد. يستخدم آلان ضمنيًا حالة التباين لعينتين كحالة أساسية، إذ يُمكن نقل القيم، لأي قيمة مختارة لـ M ، عبر حالة التباين لعينتين إلى حالة التباين لـ M عينة. لم يُذكر تفضيل واضح لحالة التباين لعينتين، حتى مع توفير الأدوات اللازمة. مع ذلك، أرست هذه المقالة الأساس لاستخدام حالة التباين لعينتين كوسيلة لمقارنة حالات التباين الأخرى لـ M عينة.

قام جيمس بارنز بتوسيع نطاق العمل المتعلق بدوال الانحياز بشكل ملحوظ، [ 19 ] حيث قدم دالتي الانحياز الحديثتين B1 و B2 . ويشير إلى تباين M- العينة باسم "تباين ألان"، مع الإشارة إلى مقالة ألان "إحصاءات معايير التردد الذري". [ 6 ] وباستخدام دالتي الانحياز الحديثتين هاتين، أصبح من الممكن إجراء تحويل كامل بين مقاييس تباين M- العينة لقيم M و T و τ المختلفة ، وذلك عن طريق التحويل من خلال تباين العينتين.

قام جيمس بارنز وديفيد آلان بتوسيع دوال الانحياز باستخدام الدالة B3 [ 20 ] لمعالجة انحياز مُقدِّر العينات المتسلسلة. وكان هذا ضروريًا للتعامل مع الاستخدام الجديد لملاحظات العينات المتسلسلة مع وجود فترة توقف بينها .

في عام ١٩٧٠، قدمت اللجنة الفنية للتردد والزمن التابعة لمعهد مهندسي الكهرباء والإلكترونيات (IEEE)، ضمن مجموعة الأجهزة والقياسات التابعة للمعهد، ملخصًا للمجال، نُشر تحت اسم الإشعار الفني رقم ٣٩٤ الصادر عن المكتب الوطني للمعايير (NBS). [ ١٥ ] كانت هذه الورقة البحثية الأولى ضمن سلسلة من الأوراق التعليمية والعملية التي ساعدت المهندسين على فهم هذا المجال. أوصت هذه الورقة باستخدام تباين عينتين مع T = τ ، مشيرةً إليه باسم تباين ألان (بدون علامات الاقتباس). يتيح اختيار هذه المعلمة معالجة جيدة لبعض أشكال الضوضاء والحصول على قياسات قابلة للمقارنة؛ فهو في جوهره القاسم المشترك الأصغر بمساعدة دالتي الانحياز B1 و B2 .

اقترح جيه جيه سنايدر طريقة محسّنة لتقدير التكرار أو التباين، باستخدام إحصاءات العينة لعدادات التكرار. [ 10 ] وللحصول على درجات حرية أكثر فعالية من مجموعة البيانات المتاحة، تكمن الحيلة في استخدام فترات رصد متداخلة. يوفر هذا تحسينًا بمقدار √n ، وقد تم دمجه في مُقدِّر تباين ألان المتداخل . [ 11 ] كما تم دمج معالجة برمجيات متغيرة-τ. [ 11 ] وقد حسّن هذا التطوير مُقدِّرات تباين ألان الكلاسيكية، كما وفّر مصدر إلهام مباشر للعمل على تباين ألان المُعدَّل .

قدم كل من هاو وآلان وبارنز تحليل فترات الثقة ودرجات الحرية والمقدرات المعتمدة. [ 11 ]

موارد تعليمية وعملية

يُعدّ مجال الزمن والتردد، واستخدامه لمصطلحات مثل تباين ألان وانحراف ألان وما شابهها، مجالًا متعدد الجوانب، يتطلب فهمًا دقيقًا للمفاهيم والقياسات العملية والمعالجة اللاحقة. ولذا، يتوفر كمٌّ هائل من المواد التعليمية التي تغطي فترة زمنية تمتد لأربعين عامًا تقريبًا. وبما أن هذه المواد تعكس تطورات البحث العلمي في عصرها، فإنها تركز على تدريس جوانب مختلفة عبر الزمن، وفي هذه الحالة، قد يكون استعراض الموارد المتاحة وسيلة مناسبة للعثور على المورد الأمثل.

يُعدّ الملخص الفني رقم 394 الصادر عن المكتب الوطني للمعايير (NBS) بعنوان "توصيف استقرار التردد" أول ملخص ذي مغزى. [ 15 ] وهو نتاج اللجنة الفنية المعنية بالتردد والزمن التابعة لمجموعة IEEE للأجهزة والقياس. يقدم هذا الملخص نظرة عامة أولية على المجال، موضحًا المشكلات، ومحددًا التعريفات الأساسية الداعمة، ومتناولًا تباين ألان، ودوال الانحياز B1 و B2 ، وتحويل مقاييس المجال الزمني. يُعدّ هذا الملخص مفيدًا، إذ إنه من أوائل المراجع التي تُجدول تباين ألان لأنواع الضوضاء الخمسة الأساسية.

يُعدّ المرجع الكلاسيكي هو دراسة NBS رقم 140 [ 27 ] الصادرة عام 1974، والتي تتضمن في الفصل الثامن "إحصاءات تحليل بيانات الزمن والتردد". [ 28 ] وهي نسخة موسعة من المذكرة الفنية رقم 394 الصادرة عن NBS، وتضيف بشكل أساسي تقنيات القياس والمعالجة العملية للقيم.

ستكون إضافة مهمة هي خصائص مصادر الإشارة وطرق القياس . [ 11 ] فهي تغطي الاستخدام الفعال للبيانات، وفترات الثقة، ودرجة الحرية الفعالة، كما تقدم مقدر تباين ألان المتداخل. يُنصح بشدة بقراءة هذا الكتاب في هذه المواضيع.

إن معيار IEEE 1139 التعريفات القياسية للكميات الفيزيائية للتردد الأساسي وقياس الوقت [ 7 ] يتجاوز كونه مجرد معيار، فهو مرجع شامل ومورد تعليمي.

يُعد كتاب "مزامنة شبكات الاتصالات الرقمية" لستيفانو بريني [ 17 ] مرجعًا حديثًا في مجال الاتصالات. فهو يُلخص ليس فقط المجال، بل أيضًا جزءًا كبيرًا من أبحاثه فيه حتى ذلك الحين. ويهدف إلى تضمين كلٍ من المقاييس التقليدية والمقاييس الخاصة بالاتصالات، مثل MTIE. ويُعتبر هذا الكتاب مرجعًا مفيدًا عند دراسة المقاييس المتعلقة بمعايير الاتصالات.

يُعدّ منشور المعهد الوطني للمعايير والتكنولوجيا الخاص رقم 1065 بعنوان "دليل تحليل استقرار التردد" للمؤلف دبليو جيه رايلي [ 18 ] مرجعًا أساسيًا لكل من يرغب في التخصص في هذا المجال. فهو غني بالمراجع، ويغطي نطاقًا واسعًا من المقاييس والانحيازات والوظائف ذات الصلة التي ينبغي أن تكون متاحة للمحلل الحديث. كما يصف أيضًا المعالجة الشاملة اللازمة لأداة حديثة.

الاستخدامات

يُستخدم تباين ألان كمقياس لاستقرار التردد في مجموعة متنوعة من المذبذبات الدقيقة، مثل المذبذبات البلورية والساعات الذرية والليزر ذي التردد المستقر، على مدى ثانية أو أكثر. ويُعبّر عن الاستقرار قصير المدى (أقل من ثانية) عادةً بضوضاء الطور . كما يُستخدم تباين ألان لتوصيف استقرار الانحياز في الجيروسكوبات ، بما في ذلك جيروسكوبات الألياف الضوئية ، وجيروسكوبات الرنان نصف الكروي، وجيروسكوبات ومقاييس التسارع بتقنية الأنظمة الكهروميكانيكية الدقيقة (MEMS) . [ 29 ] [ 30 ]

الذكرى الخمسين

في عام 2016، ستنشر IEEE-UFFC عددًا خاصًا للاحتفال بالذكرى الخمسين لتباين ألان (1966-2016). [ 31 ] وسيكون محرر هذا العدد زميل ديفيد السابق في المعهد الوطني للمعايير والتكنولوجيا ، يهودا ليفين، وهو أحدث الحاصلين على جائزة رابي الثانية .

انظر أيضاً

مراجع

  1. "التفسير الكلي للتباين*" (ملف PDF) . tf.boulder.nist.gov . مؤرشف من النسخة الأصلية (ملف PDF) بتاريخ 22 ديسمبر 2016. تم الاطلاع عليه بتاريخ 5 يونيو 2026 .
  2. "تباين هادامارد" . www.wriley.com . مؤرشف من الأصل بتاريخ 29 أبريل 2026. تم الاطلاع عليه بتاريخ 5 يونيو 2026 .
  3. منشور خاص رقم 1065 صادر عن المعهد الوطني للمعايير والتكنولوجيا، دليل تحليل استقرار التردد. يوليو 2008
  4. 1 2 كتلر، إل إس؛ سيرل، سي إل (فبراير 1966)، "بعض جوانب نظرية وقياسات تقلبات التردد في معايير التردد" (ملف PDF) ، وقائع معهد مهندسي الكهرباء والإلكترونيات ، 54 (2): 136-154 ، doi : 10.1109/proc.1966.4627 ، مؤرشف (ملف PDF) من الأصل في 9 أكتوبر 2022
  5. 1 2 3 ليسون، د. ب. (فبراير 1966)، "نموذج بسيط لطيف ضوضاء مذبذب التغذية الراجعة" ، وقائع معهد مهندسي الكهرباء والإلكترونيات ، 54 (2): 329-330 ، doi : 10.1109/proc.1966.4682 ، مؤرشف من الأصل في 1 فبراير 2014 ، تم استرجاعه في 20 سبتمبر 2012
  6. 1 2 3 4 5 6 آلان، د. إحصائيات معايير التردد الذري ، الصفحات 221-230. وقائع معهد مهندسي الكهرباء والإلكترونيات، المجلد 54، العدد 2، فبراير 1966.
  7. 1 2 3 4 5 "تعريفات معيار IEEE للكميات الفيزيائية للتردد الأساسي وقياس الزمن - عدم الاستقرار العشوائي". IEEE STD 1139-1999 . 1999. doi : 10.1109/IEEESTD.1999.90575 . ISBN 978-0-7381-1753-9.
  8. روبيولا، إنريكو (2008)، ضوضاء الطور واستقرار التردد في المذبذبات ، مطبعة جامعة كامبريدج، ISBN 978-0-521-88677-2
  9. http://www.allanstime.com/Publications/DWA/Conversion_from_Allan_variance_to_Spectral_Densities.pdf . مؤرشف بتاريخ 6 فبراير 2012 في أرشيف الإنترنت (Wayback Machine) .
  10. 1 2 3 سنايدر، جي جي: مقياس تردد فائق الدقة ، الصفحات 464-469، ندوة التحكم في التردد رقم 35، 1981.
  11. 1 2 3 4 5 6 7 D. A. Howe, DW Allan, JA Barnes: خصائص مصادر الإشارة وطرق القياس ، الصفحات 464-469، ندوة التحكم في التردد رقم 35، 1981.
  12. ITU-T Rec. G.810: التعريفات والمصطلحات الخاصة بالتزامن والشبكات ، ITU-T Rec. G.810 (08/96).
  13. ETSI EN 300 462-1-1: التعريفات والمصطلحات لشبكات التزامن ، ETSI EN 300 462-1-1 V1.1.1 (1998–05).
  14. هاو، د.أ. (1 أغسطس 2006). "ثيوه: إحصائية هجينة عالية الثقة تُحسّن من انحراف ألان" . مجلة القياس . 43 (4): ص322- ص331. doi : 10.1088/0026-1394/43/4/S20 . ISSN 0026-1394 . 
  15. 1 2 3 4 J. A. Barnes, AR Chi, LS Cutler, DJ Healey, DB Leeson, TE McGunigal, JA Mullen, WL Smith, R. Sydnor, RFC Vessot, GMR Winkler: Characterization of Frequency Stability , NBS Technical Note 394, 1970.
  16. JA Barnes, AR Chi, LS Cutler, DJ Healey, DB Leeson, TE McGunigal, JA Mullen, Jr., WL Smith, RL Sydnor, RFC Vessot, GMR Winkler: Characterization of Frequency Stability , IEEE Transactions on Instruments and Measurements 20, pp. 105–120, 1971.
  17. 1 2 3 بريني، ستيفانو: مزامنة شبكات الاتصالات الرقمية ، وايلي 2002، ISBN 0-471-61550-1.
  18. 1 2 3 NIST SP 1065: دليل تحليل استقرار التردد .
  19. 1 2 3 4 5 بارنز، جيه إيه: جداول دوال الانحياز، B 1 و B 2 ، للتباينات القائمة على عينات محدودة من العمليات ذات الكثافات الطيفية لقانون القوة ، مذكرة فنية NBS رقم 375، 1969.
  20. 1 2 3 J. A. Barnes, DW Allan: التباينات القائمة على البيانات مع وقت ميت بين القياسات ، NIST Technical Note 1318، 1990.
  21. ITU-T Rec. G.813: خصائص التوقيت لساعة المعدات التابعة SDH (SEC) ، ITU-T Rec.G.813 (03/2003).
  22. روبيولا، إنريكو (2005). "حول قياس التردد وتباين العينة باستخدام عدادات عالية الدقة" (ملف PDF) . مجلة الأدوات العلمية . 76 (5): 054703–054703–6. arXiv : physics/0411227 . Bibcode : 2005RScI...76e4703R . doi : 10.1063/1.1898203 . S2CID 119062268. مؤرشف من النسخة الأصلية (ملف PDF) بتاريخ 20 يوليو 2011. 
  23. روبيولا، إنريكو: حول قياس التردد وتباين العينة باستخدام عدادات عالية الدقة. مؤرشف في 20 يوليو 2011 في Wayback Machine ، وقائع ندوة IEEE الدولية المشتركة للتحكم في التردد واجتماع أنظمة وتطبيقات الوقت الدقيق والفترات الزمنية، الصفحات 46-49، فانكوفر، كندا، 29-31 أغسطس 2005.
  24. روبيولا، إنريكو: عدادات التردد عالية الدقة (نسخة موسعة، 53 شريحة) مؤرشفة في 20 يوليو 2011 في Wayback Machine ، ندوة ألقيت في معهد FEMTO-ST، في جامعة هنري بوانكاريه، وفي مختبر الدفع النفاث، ناسا-كالتك.
  25. ناسا:استقرار التردد على المدى القصير ، ندوة ناسا-IEEE حول استقرار التردد على المدى القصير، مركز غودارد لرحلات الفضاء، 23-24 نوفمبر 1964، منشور ناسا الخاص 80.
  26. 1 2 بارنز، جيه إيه: ضبط الوقت الذري وإحصائيات مولدات الإشارات الدقيقة ، وقائع IEEE حول استقرار التردد، المجلد 54 رقم 2، الصفحات 207-220، 1966.
  27. بلير، بي إي: الزمن والتردد: النظرية والأساسيات ، دراسة NBS رقم 140، مايو 1974.
  28. ديفيد دبليو. ألان، جون إتش. شوف ودونالد هالفورد: إحصاءات تحليل بيانات الوقت والتردد ، دراسة NBS رقم 140، الصفحات 151-204، 1974.
  29. تحليل تباين ألان لخصائص الخطأ في مقياس تسارع MEMS منخفض التكلفة. MMA8451Q afahc.ro 2014
  30. بوز، س.؛ غوبتا، أ.ك.؛ هاندل، ب. (سبتمبر 2017). "حول أداء الضوضاء والطاقة لنظام تحديد المواقع بالقصور الذاتي متعدد الوحدات المثبتة على الحذاء". المؤتمر الدولي لتحديد المواقع والملاحة الداخلية (IPIN) لعام 2017. الصفحات 1-8 . doi : 10.1109/IPIN.2017.8115944 . ISBN  978-1-5090-6299-7. S2CID 19055090 . 
  31. "IEEE UFFC | المنشورات | معاملات UFFC | مقترح لإصدار خاص من معاملات IEEE حول UFFC" . مؤرشف من الأصل في 3 سبتمبر 2014. تم الاطلاع عليه في 28 أغسطس 2014 .