ریاضیات

مسئله خلاقانه: ژیروسکوپ فوتونی برای ناوبری دقیق ماهواره

یک ماهواره کوچک در مدار زمین به یک سامانه ناوبری دقیق نیاز دارد. برای اندازه‌گیری سرعت زاویه‌ای ماهواره حول محور عمود بر صفحه دستگاه، از یک ژیروسکوپ فوتونی استفاده می‌کنیم. این ژیروسکوپ بر پایه اثر ساگناک کار می‌کند؛ یعنی دو پرتو نور در یک مسیر حلقوی، یکی در جهت ساعت‌گرد و دیگری در جهت پادساعت‌گرد حرکت می‌کنند. اگر دستگاه بچرخد، زمان رسیدن دو پرتو برابر نخواهد بود و از اختلاف فاز آن‌ها می‌توان سرعت زاویه‌ای را به دست آورد.

صورت مسئله

فرض کنید یک ژیروسکوپ فوتونی حلقوی داخل ماهواره نصب شده است. مشخصات دستگاه به شکل زیر است:

• مساحت مؤثر حلقه نوری:
  \(A = 0.02 \ \mathrm{m^2}\)
  


• محیط مسیر نوری:
  \(L = 0.50 \ \mathrm{m}\)
  


• طول موج نور لیزر:
  \(\lambda = 1550 \ \mathrm{nm}\)
  


• سرعت نور در خلأ:
  \(c = 3.0 \times 10^8 \ \mathrm{m/s}\)
  


• سرعت زاویه‌ای واقعی ماهواره:
  \(\Omega = 0.01 \ \mathrm{rad/s}\)
  

خواسته‌ها:

1. اختلاف زمان رسیدن دو پرتو را حساب کنید.


2. اختلاف فاز بین دو پرتو را به دست آورید.


3. فرکانس ضربانی خروجی ژیروسکوپ را محاسبه کنید.


4. نشان دهید چگونه از اختلاف فاز می‌توان سرعت زاویه‌ای ماهواره را تعیین کرد.


5. اگر بخواهیم دقت ژیروسکوپ ۱۰ برابر بهتر شود، چه روش‌هایی پیشنهاد می‌شود؟

گام اول: ایده فیزیکی ژیروسکوپ فوتونی

در یک ژیروسکوپ فوتونی، دو پرتو نور از یک منبع لیزری جدا می‌شوند و در یک مسیر بسته حرکت می‌کنند:

• پرتو اول در جهت ساعت‌گرد حرکت می‌کند.


• پرتو دوم در جهت پادساعت‌گرد حرکت می‌کند.

اگر حلقه نوری ساکن باشد، دو پرتو مسیرهای کاملاً یکسانی طی می‌کنند و هم‌زمان به آشکارساز می‌رسند. اما اگر حلقه با سرعت زاویه‌ای
\(\Omega\)
بچرخد، یکی از پرتوها باید کمی مسیر «بیشتر» و دیگری کمی مسیر «کمتر» مؤثر طی کند.

این باعث ایجاد اختلاف زمان می‌شود:

\[
\Delta t = t_+ - t_-
\]

این اختلاف زمان، باعث اختلاف فاز در موج نوری می‌شود و از روی آن می‌توان سرعت زاویه‌ای را محاسبه کرد.

گام دوم: رابطه اختلاف زمان ساگناک

برای یک مسیر بسته با مساحت مؤثر
\(A\)
که با سرعت زاویه‌ای
\(\Omega\)
می‌چرخد، اختلاف زمان بین دو پرتو برابر است با:

\[
\Delta t = \frac{4A\Omega}{c^2}
\]

این رابطه یکی از مهم‌ترین روابط اثر ساگناک است.

تفسیر رابطه

• هرچه مساحت حلقه بزرگ‌تر باشد، اثر ساگناک قوی‌تر است.


• هرچه سرعت زاویه‌ای بیشتر باشد، اختلاف زمان بیشتر می‌شود.


• چون مخرج رابطه شامل
  \(c^2\)
  است، مقدار اختلاف زمان بسیار بسیار کوچک است.
  

گام سوم: محاسبه اختلاف زمان

داده‌ها:

\[
A = 0.02 \ \mathrm{m^2}
\]
\[
\Omega = 0.01 \ \mathrm{rad/s}
\]
\[
c = 3.0 \times 10^8 \ \mathrm{m/s}
\]

رابطه:

\[
\Delta t = \frac{4A\Omega}{c^2}
\]

جای‌گذاری:

\[
\Delta t =
\frac{4(0.02)(0.01)}{(3.0 \times 10^8)^2}
\]

ابتدا صورت کسر را حساب می‌کنیم:

\[
4 \times 0.02 \times 0.01 = 0.0008
\]

و مخرج:

\[
(3.0 \times 10^8)^2 = 9.0 \times 10^{16}
\]

پس:

\[
\Delta t =
\frac{8.0 \times 10^{-4}}{9.0 \times 10^{16}}
\]

\[
\Delta t \approx 8.89 \times 10^{-21} \ \mathrm{s}
\]

بنابراین اختلاف زمان دو پرتو:

\[
\boxed{\Delta t \approx 8.9 \times 10^{-21} \ \mathrm{s}}
\]

این زمان بسیار کوچک است؛ اما چون فرکانس نور بسیار زیاد است، همین اختلاف زمان بسیار ناچیز می‌تواند اختلاف فاز قابل اندازه‌گیری ایجاد کند.

گام چهارم: محاسبه اختلاف فاز

اختلاف فاز موج نوری از رابطه زیر به دست می‌آید:

\[
\Delta \phi = \omega \Delta t
\]

که در آن:

\[
\omega = 2\pi f
\]

و فرکانس نور:

\[
f = \frac{c}{\lambda}
\]

پس:

\[
\omega = \frac{2\pi c}{\lambda}
\]

بنابراین:

\[
\Delta \phi =
\frac{2\pi c}{\lambda} \Delta t
\]

اکنون
\(\Delta t\)
را جای‌گذاری می‌کنیم:

\[
\Delta \phi =
\frac{2\pi c}{\lambda}
\frac{4A\Omega}{c^2}
\]

ساده می‌شود:

\[
\Delta \phi =
\frac{8\pi A\Omega}{\lambda c}
\]

این رابطه مستقیم اختلاف فاز ساگناک است:

\[
\boxed{\Delta \phi =
\frac{8\pi A\Omega}{\lambda c}}
\]

گام پنجم: جای‌گذاری عددی برای اختلاف فاز

طول موج:

\[
\lambda = 1550 \ \mathrm{nm}
\]

تبدیل به متر:

\[
\lambda = 1550 \times 10^{-9} \ \mathrm{m}
\]

\[
\lambda = 1.55 \times 10^{-6} \ \mathrm{m}
\]

اکنون رابطه اختلاف فاز:

\[
\Delta \phi =
\frac{8\pi A\Omega}{\lambda c}
\]

جای‌گذاری:

\[
\Delta \phi =
\frac{8\pi(0.02)(0.01)}
{(1.55 \times 10^{-6})(3.0 \times 10^8)}
\]

صورت:

\[
8\pi(0.02)(0.01)
=
8\pi \times 2.0 \times 10^{-4}
\]

\[
= 1.6\pi \times 10^{-3}
\]

\[
\approx 5.03 \times 10^{-3}
\]

مخرج:

\[
(1.55 \times 10^{-6})(3.0 \times 10^8)
=
4.65 \times 10^2
\]

\[
= 465
\]

پس:

\[
\Delta \phi =
\frac{5.03 \times 10^{-3}}{465}
\]

\[
\Delta \phi \approx 1.08 \times 10^{-5} \ \mathrm{rad}
\]

بنابراین:

\[
\boxed{\Delta \phi \approx 1.1 \times 10^{-5} \ \mathrm{rad}}
\]

این یعنی چرخش ماهواره با سرعت زاویه‌ای فقط
\(0.01 \ \mathrm{rad/s}\)
باعث اختلاف فازی در حد
\(10^{-5}\)
رادیان می‌شود.

گام ششم: محاسبه فرکانس ضربانی خروجی

در ژیروسکوپ‌های حلقه‌ای، اگر دو موج مخالف روی هم بیفتند، اختلاف فاز متغیر باعث ایجاد ضربان نوری می‌شود. برای ژیروسکوپ حلقه‌ای، فرکانس ضربانی ساگناک معمولاً به شکل زیر نوشته می‌شود:

\[
\Delta f = \frac{4A\Omega}{\lambda L}
\]

که در آن:

• 
  \(A\)
  مساحت حلقه است.
  


• 
  \(L\)
  محیط مسیر نوری است.
  


• 
  \(\Omega\)
  سرعت زاویه‌ای است.
  


• 
  \(\lambda\)
  طول موج نور است.
  

جای‌گذاری:

\[
\Delta f =
\frac{4(0.02)(0.01)}
{(1.55 \times 10^{-6})(0.50)}
\]

صورت:

\[
4(0.02)(0.01)=8.0 \times 10^{-4}
\]

مخرج:

\[
(1.55 \times 10^{-6})(0.50)
=
7.75 \times 10^{-7}
\]

پس:

\[
\Delta f =
\frac{8.0 \times 10^{-4}}
{7.75 \times 10^{-7}}
\]

\[
\Delta f \approx 1032 \ \mathrm{Hz}
\]

بنابراین فرکانس ضربانی خروجی:

\[
\boxed{\Delta f \approx 1.03 \times 10^3 \ \mathrm{Hz}}
\]

یعنی آشکارساز نوری، ضربانی در حدود
\(1.03 \ \mathrm{kHz}\)
ثبت می‌کند.

گام هفتم: استخراج سرعت زاویه‌ای از اختلاف فاز

رابطه اختلاف فاز:

\[
\Delta \phi =
\frac{8\pi A\Omega}{\lambda c}
\]

می‌خواهیم
\(\Omega\)
را بر حسب
\(\Delta \phi\)
به دست آوریم. دو طرف را در
\(\lambda c\)
ضرب می‌کنیم:

\[
\Delta \phi \lambda c = 8\pi A\Omega
\]

حال دو طرف را بر
\(8\pi A\)
تقسیم می‌کنیم:

\[
\Omega =
\frac{\Delta \phi \lambda c}{8\pi A}
\]

پس فرمول ناوبری ژیروسکوپ فوتونی:

\[
\boxed{
\Omega =
\frac{\Delta \phi \lambda c}{8\pi A}
}
\]

این رابطه می‌گوید اگر ماهواره اختلاف فاز را اندازه بگیرد، می‌تواند سرعت زاویه‌ای خودش را محاسبه کند.

گام هشتم: بررسی عددی معکوس

برای اطمینان، مقدار به دست آمده برای اختلاف فاز را در رابطه معکوس قرار می‌دهیم:

\[
\Omega =
\frac{\Delta \phi \lambda c}{8\pi A}
\]

با:

\[
\Delta \phi = 1.08 \times 10^{-5}
\]
\[
\lambda = 1.55 \times 10^{-6} \ \mathrm{m}
\]
\[
c = 3.0 \times 10^8 \ \mathrm{m/s}
\]
\[
A = 0.02 \ \mathrm{m^2}
\]

جای‌گذاری:

\[
\Omega =
\frac{(1.08 \times 10^{-5})(1.55 \times 10^{-6})(3.0 \times 10^8)}
{8\pi(0.02)}
\]

حاصل تقریباً:

\[
\Omega \approx 0.01 \ \mathrm{rad/s}
\]

پس رابطه درست عمل می‌کند:

\[
\boxed{\Omega \approx 0.01 \ \mathrm{rad/s}}
\]

گام نهم: تفسیر فیزیکی برای ماهواره

اگر ماهواره حول محور عمود بر صفحه ژیروسکوپ بچرخد، ژیروسکوپ فوتونی این چرخش را از روی اختلاف فاز نور تشخیص می‌دهد. در واقع نور مانند یک «خط‌کش زمانی بسیار دقیق» عمل می‌کند.

چون فرکانس نور بسیار زیاد است، حتی اختلاف زمان‌هایی در حد:

\[
10^{-21} \ \mathrm{s}
\]

می‌توانند اثر قابل اندازه‌گیری در فاز ایجاد کنند.

در سامانه ناوبری ماهواره، این سرعت زاویه‌ای وارد معادلات تعیین وضعیت می‌شود. اگر بردار سرعت زاویه‌ای ماهواره برابر باشد با:

\[
\vec{\Omega} =
\Omega_x \hat{i} +
\Omega_y \hat{j} +
\Omega_z \hat{k}
\]

برای اندازه‌گیری کامل چرخش سه‌بعدی، معمولاً سه ژیروسکوپ فوتونی در سه راستای عمود بر هم نصب می‌شوند.

بنابراین:

• ژیروسکوپ اول،
  \(\Omega_x\)
  را اندازه می‌گیرد.
  


• ژیروسکوپ دوم،
  \(\Omega_y\)
  را اندازه می‌گیرد.
  


• ژیروسکوپ سوم،
  \(\Omega_z\)
  را اندازه می‌گیرد.
  

گام دهم: افزایش دقت ژیروسکوپ

از رابطه اختلاف فاز:

\[
\Delta \phi =
\frac{8\pi A\Omega}{\lambda c}
\]

می‌بینیم که حساسیت ژیروسکوپ به چرخش به عوامل زیر وابسته است:

\[
\Delta \phi \propto \frac{A}{\lambda}
\]

یعنی برای یک سرعت زاویه‌ای مشخص، اگر اختلاف فاز بزرگ‌تر شود، اندازه‌گیری آسان‌تر و دقیق‌تر خواهد بود.

روش اول: افزایش مساحت حلقه

اگر
\(A\)
را ۱۰ برابر کنیم:

\[
A' = 10A
\]

آنگاه:

\[
\Delta \phi' = 10 \Delta \phi
\]

پس دقت اندازه‌گیری چرخش بهتر می‌شود.

روش دوم: استفاده از فیبر نوری چنددور

در ژیروسکوپ فیبر نوری، نور می‌تواند به جای یک بار، چندین بار دور حلقه بگردد. اگر تعداد دورهای فیبر
\(N\)
باشد، مساحت مؤثر تقریباً برابر می‌شود با:

\[
A_{\mathrm{eff}} = N A
\]

در این صورت اختلاف فاز:

\[
\Delta \phi =
\frac{8\pi N A \Omega}{\lambda c}
\]

اگر
\(N = 1000\)
باشد، حساسیت ژیروسکوپ تقریباً هزار برابر افزایش می‌یابد.

روش سوم: کاهش طول موج

چون:

\[
\Delta \phi \propto \frac{1}{\lambda}
\]

اگر از طول موج کوتاه‌تر استفاده شود، اختلاف فاز بیشتر می‌شود. مثلاً اگر طول موج نصف شود:

\[
\lambda' = \frac{\lambda}{2}
\]

آنگاه:

\[
\Delta \phi' = 2\Delta \phi
\]

روش چهارم: استفاده از تشدیدگر نوری

اگر مسیر نوری به صورت یک تشدیدگر با ضریب کیفیت بالا طراحی شود، نور مدت بیشتری در حلقه باقی می‌ماند و اثر چرخش تقویت می‌شود.

در چنین سامانه‌ای، نور بارها در حلقه رفت‌وبرگشت می‌کند و اثر ساگناک چندین بار جمع می‌شود.

روش پنجم: کاهش نویز آشکارسازی

حتی اگر اثر فیزیکی قوی باشد، نویز می‌تواند اندازه‌گیری را خراب کند. منابع مهم نویز عبارت‌اند از:

• نویز شات فوتونی


• نویز حرارتی آشکارساز


• نوسان توان لیزر


• تغییرات دمایی فیبر یا حلقه


• لرزش مکانیکی ماهواره

برای کاهش نویز می‌توان:

• از لیزر پایدارتر استفاده کرد.


• دما را با کنترل‌کننده حرارتی ثابت نگه داشت.


• حلقه نوری را روی پایه مکانیکی پایدار نصب کرد.


• از آشکارسازی تفاضلی استفاده کرد.


• از فیلتر کالمن برای ترکیب داده‌های ژیروسکوپ و حسگر ستاره‌ای استفاده کرد.

گام یازدهم: نسخه چنددوره‌ای برای افزایش حساسیت

اگر به جای یک حلقه ساده، از فیبر نوری با
\(N\)
دور استفاده کنیم، رابطه اختلاف فاز به صورت زیر درمی‌آید:

\[
\Delta \phi_N =
\frac{8\pi N A\Omega}{\lambda c}
\]

فرض کنید:

\[
N = 1000
\]

در حالت تک‌دور داشتیم:

\[
\Delta \phi \approx 1.08 \times 10^{-5} \ \mathrm{rad}
\]

پس در حالت هزار دور:

\[
\Delta \phi_N =
1000 \times 1.08 \times 10^{-5}
\]

\[
\Delta \phi_N =
1.08 \times 10^{-2} \ \mathrm{rad}
\]

یعنی:

\[
\boxed{
\Delta \phi_N \approx 0.0108 \ \mathrm{rad}
}
\]

این مقدار نسبت به حالت تک‌دور بسیار راحت‌تر قابل اندازه‌گیری است.

گام دوازدهم: تعیین حداقل سرعت زاویه‌ای قابل تشخیص

فرض کنید سامانه آشکارساز بتواند حداقل اختلاف فاز زیر را تشخیص دهد:

\[
\Delta \phi_{\min} = 10^{-6} \ \mathrm{rad}
\]

از رابطه:

\[
\Omega =
\frac{\Delta \phi \lambda c}{8\pi A}
\]

برای حالت تک‌دور:

\[
\Omega_{\min} =
\frac{(10^{-6})(1.55 \times 10^{-6})(3.0 \times 10^8)}
{8\pi(0.02)}
\]

صورت:

\[
(10^{-6})(1.55 \times 10^{-6})(3.0 \times 10^8)
=
4.65 \times 10^{-4}
\]

مخرج:

\[
8\pi(0.02) \approx 0.503
\]

پس:

\[
\Omega_{\min}
\approx
\frac{4.65 \times 10^{-4}}{0.503}
\]

\[
\Omega_{\min}
\approx
9.25 \times 10^{-4} \ \mathrm{rad/s}
\]

بنابراین در حالت تک‌دور:

\[
\boxed{
\Omega_{\min}
\approx
9.3 \times 10^{-4} \ \mathrm{rad/s}
}
\]

اگر از فیبر هزار دور استفاده کنیم، مساحت مؤثر هزار برابر می‌شود و حداقل سرعت زاویه‌ای قابل تشخیص هزار برابر کوچک‌تر خواهد شد:

\[
\Omega_{\min,N}
=
\frac{\Omega_{\min}}{N}
\]

برای
\(N=1000\):

\[
\Omega_{\min,N}
\approx
9.3 \times 10^{-7} \ \mathrm{rad/s}
\]

پس:

\[
\boxed{
\Omega_{\min,N}
\approx
9.3 \times 10^{-7} \ \mathrm{rad/s}
}
\]

این مقدار برای ناوبری دقیق ماهواره بسیار بهتر است.

گام سیزدهم: الگوریتم ساده ناوبری با ژیروسکوپ فوتونی

برای استفاده از ژیروسکوپ در ماهواره، می‌توان مراحل زیر را اجرا کرد:

1. لیزر پایدار، دو پرتو هم‌فاز تولید می‌کند.


2. دو پرتو در دو جهت مخالف وارد حلقه نوری می‌شوند.


3. پس از طی مسیر، دو پرتو روی آشکارساز ترکیب می‌شوند.


4. آشکارساز شدت تداخلی را اندازه می‌گیرد.


5. از شدت تداخلی، اختلاف فاز
   \(\Delta \phi\)
   استخراج می‌شود.
   


6. با رابطه زیر سرعت زاویه‌ای محاسبه می‌شود:
   

   
   \[
   \Omega =
   \frac{\Delta \phi \lambda c}{8\pi A_{\mathrm{eff}}}
   \]
   

   
   


7. با انتگرال‌گیری زمانی از سرعت زاویه‌ای، زاویه چرخش ماهواره به دست می‌آید:
   

   
   \[
   \theta(t) =
   \theta(0) + \int_0^t \Omega(t') \ dt'
   \]
   

   
   

اگر داده‌ها گسسته باشند، یعنی ژیروسکوپ در بازه‌های زمانی
\(\Delta t_s\)
نمونه‌برداری کند، رابطه عددی به شکل زیر است:

\[
\theta_{k+1}
=
\theta_k + \Omega_k \Delta t_s
\]

این رابطه، پایه ساده تعیین وضعیت زاویه‌ای ماهواره است.

گام چهاردهم: شدت تداخلی خروجی

وقتی دو پرتو با اختلاف فاز
\(\Delta \phi\)
روی آشکارساز ترکیب می‌شوند، شدت خروجی به صورت زیر است:

\[
I =
I_0 \left[1 + \cos(\Delta \phi)\right]
\]

اگر اختلاف فاز کوچک باشد، می‌توان تقریب زد:

\[
\cos(\Delta \phi)
\approx
1 - \frac{\Delta \phi^2}{2}
\]

پس:

\[
I
\approx
I_0
\left[
1 + 1 - \frac{\Delta \phi^2}{2}
\right]
\]

\[
I
\approx
2I_0 - \frac{I_0}{2}\Delta \phi^2
\]

اما مشکل اینجاست که برای فازهای خیلی کوچک، تغییر شدت مرتبه دوم است و حساسیت کم می‌شود. به همین دلیل در عمل از روش‌هایی مانند:

• مدولاسیون فاز


• آشکارسازی قفل‌شونده


• آشکارسازی تفاضلی


• بایاس فازی نزدیک
  \(\pi/2\)
  

استفاده می‌شود تا پاسخ سامانه نسبت به اختلاف فاز تقریباً خطی شود.

گام پانزدهم: نتیجه نهایی مسئله

برای ژیروسکوپ فوتونی داده‌شده:

جمع‌بندی مفهومی

ژیروسکوپ فوتونی با استفاده از اثر ساگناک می‌تواند چرخش ماهواره را بدون نیاز به قطعات مکانیکی اندازه‌گیری کند. اختلاف زمان بین دو پرتو بسیار کوچک است، اما چون نور فرکانس بسیار بالایی دارد، این اختلاف زمان به اختلاف فاز قابل آشکارسازی تبدیل می‌شود.

رابطه کلیدی مسئله برای اختلاف زمان:

\[
\boxed{
\Delta t =
\frac{4A\Omega}{c^2}
}
\]

رابطه کلیدی برای اختلاف فاز:

\[
\boxed{
\Delta \phi =
\frac{8\pi A\Omega}{\lambda c}
}
\]

و رابطه استخراج سرعت زاویه‌ای:

\[
\boxed{
\Omega =
\frac{\Delta \phi \lambda c}{8\pi A}
}
\]

برای افزایش دقت، بهترین راه‌ها عبارت‌اند از:

• افزایش مساحت مؤثر حلقه


• استفاده از فیبر نوری چنددور


• کاهش طول موج نور


• استفاده از تشدیدگر نوری با کیفیت بالا


• کاهش نویز آشکارسازی و نوسانات حرارتی

در نتیجه، یک ژیروسکوپ فوتونی می‌تواند به عنوان قلب سامانه تعیین وضعیت ماهواره عمل کند و چرخش‌های بسیار کوچک را با دقت بالا تشخیص دهد.