نمایش مبدأ برای چگالی طیفی

[[File:Fluorescent lighting spectrum peaks labelled.svg|thumb|شکل بالا چگالی طیفی یک [[لامپ فلورسنت|نور فلورسنت]] به صورت تابعی از طول‌موج نوری است که نشان‌دهنده قله‌هایی در گذارهای اتمی است، که این موضوع توسط فلش‌های عدد دار نشان داده شده‌است.]]

[[File:Voice waveform and spectrum.png|thumb|حالت‌موج صدا روی زمان (سمت چپ) دارای یک طیف توان صدایی گسترده‌است (راست).]]

طیف<ref>{{یادکرد فرهنگستان|مصوب=طیف|بیگانه=spectrum|بیگانه در فارسی=|حوزه=فیزیک|دفتر=نخست|بخش=فارسی|سرواژه=طیف}}</ref> توان<ref>{{یادکرد فرهنگستان|مصوب=توان|بیگانه=power|بیگانه در فارسی=|حوزه=شیمی، فیزیک|دفتر=دوم|بخش=فارسی|سرواژه=توان}}</ref> <math>S_{xx}(f)</math> برای یک [[سری زمانی]] <math>x(t)</math> توصیف کننده توزیع [[توان (فیزیک)|توان]] به مولفه‌های فرکانسی تشکیل دهنده سیگنال است.<ref>{{cite web
 | url = http://user.it.uu.se/~ps/SAS-new.pdf
 | title = Spectral Analysis of Signals
 |author1=P Stoica
 |author-link=Peter Stoica
 |author2=R Moses
  |name-list-style=amp | year = 2005
 }}</ref> بر اساس [[آنالیز فوریه]]، هر سیگنال فیزیکی را می‌توان به تعدادی از فرکانس‌های گسسته، یا طیفی از فرکانس‌ها در یک محدوده پیوسته تجزیه کرد. میانگین آماری یک سیگنال معین یا نوعی از سیگنال (شامل [[نویز (الکترونیک)|نویز]]) وقتیکه از دیدگاه محتوای فرکانس آن تحلیل می‌شوند، [[طیف]] آن سیگنال نامیده می‌شوند.

وقتیکه انرژی سیگنال در حول و حوش یک بازه زمانی محدود متمرکز شود، مخصوصاً اگر کل انرژی آن محدود باشد، می‌توان '''چگالی طیفی انرژی''' آن را محاسبه کرد. چیزی که معمول‌تر است، '''چگالی طیفی توان''' (یا به صورت ساده‌تر '''طیف توان''') است که به سیگنال‌های موجود در ''کل'' زمان، یا به یک دوره زمانی به اندازه کافی بزرگ (مخصوصا در رابطه با مدت زمان اندازه‌گیری) که این آخری هم می‌تواند به صورت مشابه روی یک بازه زمانی نامحدود باشد، اعمال می‌شود. چگالی طیف توان (PSD) آنوقت به توزیع انرژی طیفی که «در هر واحد زمان» یافت می‌شود، ارجاع دارد، زیرا انرژی کلی چنین سیگنالی در تمام زمان معمولاً بینهایت است. [[مجموع‌یابی]] یا انتگرال‌گیری مولفه‌های طیفی منجر به توان کلی (برای یک فرایند فیزیکی) یا واریانس (در یک فرایند آماری) می‌شود، این مشابه آن‌چیزی است که از انتگرال‌گیری <math>x^2(t)</math> روی دامنه زمانی به دست می‌آید، که این موضوع توسط [[:en:Parseval's_theorem|قضیه پارسوال]] اجبار می‌شود.<ref>{{cite web
 | url = http://user.it.uu.se/~ps/SAS-new.pdf
 | title = Spectral Analysis of Signals
 |author1=P Stoica
 |author-link=Peter Stoica
 |author2=R Moses
  |name-list-style=amp | year = 2005
 }}</ref>

طیف یک فرایند فیزیکی <math>x(t)</math> معمولاً شامل اطلاعات اساسی دربارهٔ طبیعت <math>x</math> است. یک مثال آن، [[زیروبمی]] و [[رنگ (صدا)|رنگ]] برای یک وسیله موسیقی است به صورت فوری توسط تحلیل طیفی معین می‌شود. مثال دیگر [[رنگ]] یک منبع نوری است که توسط طیف موج الکترومغناطیسی برای فیلد الکتریکی <math>E(t)</math> موقعی که در فرکانس بسیار بالا نوسان می‌کند، تعیین می‌شود. برای به دست آوردن یک فرم طیفی از سری زمانی، مثل موارد بالا، باید از [[تبدیل فوریه]]، و تعمیم‌های مبتنی بر تحلیل فوریه استفاده کرد. در بسیاری از حالات، دامنه زمانی در عمل اعمال نمی‌گردد، مثل موقعی که از یک [[:en:Dispersive_prism|منشور پراکنده]] استفاده می‌کنیم تا به یک طیف نوری در یک [[طیف‌سنج نوری|طیف‌سنج]] برسیم، یا موقعی که یک صدا از طریق تأثیرگذاری روی گیرنده‌های صوتی در [[گوش درونی]] درک می‌شود، که در این‌جا هر گیرنده به یک فرکانس خاص حساس است.

با این حال این مقاله روی وضعیت‌هایی تمرکز می‌کند که سری زمانی‌شان را می‌دانیم (حداقل در مفهوم آماری) یا سری زمانی به صورت مستقیم اندازه‌گیری شده‌است (مثلا توسط میکروفن و توسط رایانه نمونه‌برداری شده). طیف توانی در [[پردازش سیگنال|پردازش سیگنال آماری]] مهم است، همچنین در مطالعه آماری [[فرایند تصادفی|فرایندهای تصادفی]]، و نیز در خیلی از رشته‌های دیگر در [[فیزیک]] و [[مهندسی]] هم اهمیت دارد. معمولاً این فرایند، تابعی از زمان است، اما می‌توان داده را در زمینه فضایی نیز بررسی کرد که این توسط [[فرکانس فضایی]] تجزیه می‌شود.<ref>{{cite web|url=http://user.it.uu.se/~ps/SAS-new.pdf|title=Spectral Analysis of Signals|author1=P Stoica|author-link=Peter Stoica|author2=R Moses|name-list-style=amp|year=2005}}</ref>

== توضیح ==
هر سیگنالی را که بتوان به صورت متغیری نمایش داد که در زمان تغییر می‌کند، یک طیف فرکانسی متناظر دارد. این موضوع شامل موجودیت‌های آشنایی مثل [[نور|نور قابل مشاهده]] (که به صورت [[رنگ]] درک می‌شود)، نت موسیقی (که به صورت [[زیروبمی]] درک می‌شود)، [[فرکانس رادیویی|رادیو/تلویزیون]] (که توسط فرکانس شان، یا گاهی توسط [[طول‌موج]] تعیین می‌شود) و حتی چرخش منظم زمین می‌شود. موقعی که این سیگنال‌ها به حالت یک طیف فرکانسی دیده شوند، جنبه‌های معینی از سیگنال دریافتی یا فرایندهای زیرین تولیدکننده آن‌ها آشکار می‌شود. در بعضی از حالات طیف فرکانسی دارای یک قله مجزا است که با یک مولفه [[موج سینوسی]] متناظر است. و بعلاوه ممکن است قله‌هایی متناظر با [[هارمونیک]]‌های یک قله اصلی وجود داشته باشد، که نشان‌دهنده یک سیگنال تناوبی است که یک سینوسی ساده نیست. یا یک طیف پیوسته می‌تواند بازه‌های فرکانسی باریک را نشان دهد که متناظر با تشدید (رزونانس) به صورت قوی افزایش می‌یابد، یا بازه‌های فرکانسی شامل توان تقریباً صفر باشند که توسط یک [[فیلتر شکافی]] تولید شده‌اند.

در [[فیزیک]]، سیگنال می‌تواند یک موج باشد، مثل یک [[موج الکترومغناطیسی]]، یک [[موج صوتی]] یا اینکه [[ارتعاش]] یک مکانیزم باشد. چگالی طیف توانی (PSD) برای سیگنال، توصیف‌کننده [[توان (فیزیک)|توان]] موجود در سیگنال به صورت تابعی از فرکانس در واحد فرکانس است. چگالی طیف توانی معمولاً به صورت [[وات]] به ازای [[هرتز]] (W/Hz) بیان می‌شود.<ref>{{cite book|title=VSAT Networks|author=Gérard Maral|publisher=John Wiley and Sons|year=2003|isbn=978-0-470-86684-9|url=https://books.google.com/books?id=CMx5HQ1Mr_UC&q=%22power+spectral+density%22+W/Hz&pg=PR20}}</ref>

موقعی که یک سیگنال مثلا از دیدگاه [[ولتاژ]]، تعریف شده‌است، هیچ توان یکتایی مرتبط با دامنه بیان شده وجود ندارد. در این حالت، «توان» به سادگی از دیدگاه مربع سیگنال حساب می‌شود، همان‌طور که این با توان واقعی تحویل شده توسط آن سیگنال به یک [[امپدانس الکتریکی|امپدانس]] معین متناسب است. از این رو می‌توان از واحدهای V<sup>2</sup> Hz<sup>−1</sup> برای PSD استفاده کنیم و از V<sup>2</sup> s Hz<sup>−1</sup> برای ESD (چگالی طیف انرژی energy spectral density)<ref>{{cite book|title=Fundamentals of Noise and Vibration Analysis for Engineers|author1=Michael Peter Norton|author2=Denis G. Karczub|name-list-style=amp|publisher=[[Cambridge University Press]]|year=2003|isbn=978-0-521-49913-2|url=https://books.google.com/books?id=jDeRCSqtev4C&q=%22power+spectral+density%22+%22energy+spectral+density%22&pg=PA352}}</ref> استفاده کنیم، اگرچه هیچ «توان» یا «انرژی» واقعی تعیین نشده‌است.

بعضی اوقات به یک چگالی طیفی دامنه (ASD) (amplitude spectral density) برخورد می‌کنیم، که جذر دوم PSD است؛ ASD برای سیگنال ولتاژ واحد V Hz<sup>−1/2</sup> دارد.<ref>{{cite web|url=http://www.lumerink.com/courses/ece697/docs/Papers/The%20Fundamentals%20of%20FFT-Based%20Signal%20Analysis%20and%20Measurements.pdf|title=The Fundamentals of FFT-Based Signal Analysis and Measurement|author1=Michael Cerna|author2=Audrey F. Harvey|name-list-style=amp|year=2000}}{{پیوند مرده}}</ref> این موضوع موقعی مفید است که بیشتر شکل طیف ثابت است، زیرا آنوقت تغییرات در ASD با تغییرات در خود مرحله ولتاژ سیگنال متناسب است. اما از نظر ریاضیاتی استفاده از PSD ترجیح دارد، زیرا تنها در آن حالت مساحت زیر منحنی از لحاظ توان واقعی تحت همه فرکانس‌ها یا تحت پهنای‌باند معین معنی‌دار است.

در حالت معمول، واحد PSD همان نسبت واحدهای واریانس در واحد فرکانس است؛ از این رو، مثلا، سری‌های مقادیر جابجایی (در متر) روی زمان (در ثانیه) دارای PSD در واحد m<sup>2</sup>/Hz خواهد بود. برای تحلیل ارتعاش تصادفی، به صورت مکرر از واحدهای ''g''<sup>2</sup> Hz<sup>−1</sup> برای [[شتاب]] PSD استفاده می‌شود. در اینجا ''g'' نشان‌دهنده [[نیروی گرانش]] است.<ref>{{cite book|title=Reliability Engineering|author=Alessandro Birolini|publisher=Springer|year=2007|isbn=978-3-540-49388-4|page=83|url=https://books.google.com/books?id=xPIW3AI9tdAC&q=acceleration-spectral-density+g+hz&pg=PA83}}</ref>

از نظر ریاضیاتی، نیازی به انتساب ابعاد فیزیکی به سیگنال یا متغیر مستقل نیست. و در بحثی که در ادامه می‌آید، معنی ''x(t)'' نامعین می‌ماند، اما فرض می‌شود که متغیر مستقل، زمان است.

== تعریف ==
=== چگالی طیف انرژی ===
چگالی طیف انرژی توصیف کننده روشی است که [[انرژی (پردازش سیگنال)|انرژی]] یک سیگنال یا یک [[سری زمانی]] در طول فرکانس توزیع می‌شود. در اینجا، اصطلاح [[انرژی (پردازش سیگنال)|انرژی]] در مفهوم عمومی‌اش در پردازش سیگنال استفاده شده‌است؛<ref name="oppenheim">{{cite book|last1=Oppenheim|last2=Verghese|title=Signals, Systems, and Inference|pages=32–4}}</ref> یعنی، انرژی <math>E</math> از یک سیگنال <math>x(t)</math> به این صورت است:

:<math> E \triangleq \int_{-\infty}^\infty |x(t)|^2\ dt.</math>

چگالی طیف انرژی برای  موارد گذرا بسیار مناسب است-یعنی، سیگنال‌های شبیه پالس، که در آن موقع انرژی کلی محدودی دارد. چه محدود باشد یا نباشد، [[:en:Parseval's_theorem|قضیه پارسوال]]،<ref name="Stein">{{cite book|last=Stein|first=Jonathan Y.|title=Digital Signal Processing: A Computer Science Perspective|url=https://archive.org/details/digitalsignalpro00jona|page=[https://archive.org/details/digitalsignalpro00jona/page/n132 115]|publisher=Wiley|year=2000}}</ref> (یا قضیه پلانچرل) به ما یک عبارت جایگزین برای انرژی سیگنال می‌دهد:

:<math>\int_{-\infty}^\infty |x(t)|^2\, dt = \int_{-\infty}^\infty |\hat{x}(f)|^2\ df,</math>

که در آن:

:<math>\hat{x}(f) \triangleq\int_{-\infty}^\infty e^{-i 2\pi ft}x(t) \ dt</math>

مقدار [[تبدیل فوریه]] <math>x(t)</math> در [[فرکانس]] <math>f</math> (به [[هرتز]]) است. این قضیه در حالت‌های زمان-گسسته نیز درست است. به این دلیل که [[انتگرال]] سمت راست همان انرژی سیگنال است، عبارت‌زیرانتگرال <math>\left |\hat{x}(f) \right |^2</math> را می‌توان به صورت یک [[تابع چگالی احتمال|تابع چگالی]] تفسیر کرد، که توصیف‌کننده انرژی موجود در سیگنال در فرکانس <math>f</math> است. از این رو، '''چگالی طیفی انرژی''' برای <math>x(t)</math> به این صورت تعریف می‌شود:<ref name="Stein"/>

{{Equation box 1
|indent =
|title=
|equation = {{NumBlk|:|<math> \bar{S}_{xx}(f) \triangleq \left |\hat{x}(f) \right |^2 </math>|{{EquationRef|Eq.1}}}}
|cellpadding= 6
|border
|border colour = #0073CF
|background colour=#F5FFFA}}

تابع <math>\bar{S}_{xx}(f)</math> و [[خودهمبستگی]] برای <math>x(t)</math> یک زوج تبدیل فوریه می‌سازند، این نتیجه‌گیری به نام [[قضیه وینر-خینشین]] شناخته می‌شود ([[:en:Periodogram#Definition|تناوب‌سنج]] را ببیند).

به عنوان یک مثال فیزیکی از اینکه چگونه می‌توان چگالی طیف انرژی یک سیگنال را اندازه‌گرفت، فرض کنید که <math>V(t)</math> نمایش دهنده [[پتانسیل الکتریکی|پتانسیل]] (در [[ولت]]) از یک پالس الکتریکی باشد که در طول یک [[خط انتقال]] انتشار با [[امپدانس الکتریکی|امپدانس]] Z منتشر شده‌است، و فرض کنید که این خط توسط یک مقاومت [[تطبیق امپدانس|متناظر]] خاتمه می‌یابد (از این رو همه انرژی پالس به مقاومت تحویل داده می‌شود، و هیچ انرژی بازگشت نمی‌کند). از طریق [[قانون اهم]]، توان تحویل شده به مقاومت در زمان <math>t</math> برابر <math>V(t)^2/Z</math> است، از این رو انرژی کلی از طریق انتگرال‌گیری <math>V(t)^2/Z</math> نسبت به زمان روی مدت پالس به دست می‌آید. برای یافتن مقدار چگالی طیفی انرژی <math>\bar{S}_{xx}(f)</math> در فرکانس <math>f</math>، می‌توان بین خط انتقال و مقاومت یک [[فیلتر میان‌گذر]] قرار داد، که فقط محدوده باریکی از فرکانس‌ها (به آن بگوییم <math>\Delta f</math>) نزدیک فرکانس مورد علاقه را از خود عبور می‌دهد، و سپس کل انرژی <math>E(f)</math> را که از میان مقاومت از بین می‌رود را باید اندازه‌گیری کنیم. مقدار چگالی طیف انرژی در <math>f</math> آنوقت به مقدار <math>E(f)/\Delta f</math> تخمین زده می‌شود. در این مثال، به این دلیل که توان <math>V(t)^2/Z</math> دارای واحد V<sup>2</sup> Ω<sup>−1</sup> است، انرژی <math>E(f)</math> واحد V<sup>2</sup> s Ω<sup>−1</sup> = J دارد، و از این رو تخمین <math>E(f)/\Delta  f</math> از چگالی طیف انرژی واحد J Hz<sup>−1</sup> را دارد که نیاز داشتیم. در بیشتر حالات، می‌توان گام تقسیم بر <math>Z</math> را فراموش کرد و از این رو چگالی طیف انرژی در عوض واحد V<sup>2</sup> Hz<sup>−1</sup> را خواهد داشت.

این تعریف به روش مستقیم به [[سیگنال گسسته]] با تعداد قابل شمارش نامتناهی مقایر <math>x_n</math> تعمیم می‌یابد، مثل یک سیگنال که در زمان‌های گسسته <math>x_n=x_0 + (n\Delta t)</math> نمونه‌گیری شده‌است:

:<math>\bar{S}_{xx}(f) = \lim_{N\to \infty}(\Delta t)^2 \underbrace{\left|\sum_{n=-N}^N x_n e^{-i 2\pi f n \Delta t}\right|^2}_{|\hat x_d(f)|^2},</math>

در اینجا <math>\hat x_d(f)</math> همان [[تبدیل فوریه زمان-گسسته]] برای <math>x_n.</math> است. بازه نمونه‌گیری <math>\Delta t</math> برای نگهداری واحدهای فیزیکی صحیح، و نیز برای اطمینان از اینکه ما در حد <math>\Delta t\to 0.</math> حالت پیوسته را بازیابی می‌کنیم، نیاز است. اما در علوم ریاضیات این بازه معمولاً به ۱ تنظیم می‌شود، که نتایج را با هزینه تعمیم آن، ساده‌سازی می‌کند. ([[:en:Normalized_frequency_(unit)|فرکانس نرمال‌سازی‌شده]] را ببیند).

=== چگالی طیف توانی ===
تعریف بالا برای چگالی طیف انرژی، برای [[حالت گذرا]] (سیگنال‌های مشابه پالس) مناسب است، که در آن انرژی در حول یک پنجره زمانی متمرکز شده‌است؛ در آنصورت تبدیل‌های فوریه برای سیگنال‌ها معمولاً وجود دارد. برای سیگنال‌های پیوسته روی همه زمان، به صورت جایگزین باید ''چگالی طیفی توان'' (PSD) را تعریف کرد، که این تعریف برای [[فرایند مانا|فرایندهای مانا]] وجود دارد؛ این موضوع توصیف می‌کند که چگونه [[توان (فیزیک)|توان]] یک سیگنال یا سری زمانی روی فرکانس توزیع شده‌است، مثل مثال ساده‌ای که در قبل داشتیم. در اینجا، توان می‌تواند توان فیزیکی واقعی باشد، یا به صورت معمول‌تر، برای سادگی با سیگنال‌های انتزاعی، به سادگی توسط مقدار مجذور سیگنال شناسایی شود. برای مثال، آماردانان [[واریانس]] یک تابع را روی زمان (یا روی متغیر مستقل دیگر) <math>x(t)</math> مطالعه می‌کنند، و به کمک قیاس با سیگنال‌های الکتریکی (از میان دیگر فرایندهای فیزیکی)، معمول است که به آن ''طیف توانی'' بگوییم، حتی وقتیکه هیچ توان فیزیکی درگیر نیست. اگر کسی بخواهد یک منبع [[ولتاژ]] فیزیکی را که از <math>x(t)</math> پیروی می‌کند بسازد، و آن را به پایانه‌های یک [[مقاومت (قطعه الکتریکی)|مقاومت]] ۱ [[اهم|اهمی]] اعمال کند، آنوقت توان آنی هدر رفته در آن مقاومت به اندازه <math>x(t)^2</math> [[وات]] است.

توان میانگین <math>P</math> از یک سیگنال <math>x(t)</math> در تمام زمان از این رو توسط این میانگین زمانی به دست می‌آید، که در آن دوره <math> T</math> در یک زمان اختیاری <math> t=t_{0}</math> مرکزدهی شده‌است:

:<math> P = \lim_{T\to \infty} \frac 1 {T} \int_{t_{0}-T/2}^{t_{0}+T/2} |x(t)|^2\,dt</math>

با این حال، به خاطر تعامل با ریاضیاتی که در ادامه می‌آید، ساده‌تر است که با حدود زمانی موجود در خود سیگنال به جای حدود زمانی موجود در مرزهای انتگرال سروکار داشته باشیم. از این رو، ما یک نمایش جایگزین برای توان میانگین داریم که در آن <math> x_{T}(t)=x(t)w_{T}(t)</math> و <math> w_{T}(t)</math> در یک دوره اختیاری یک و در بقیه محل‌ها صفر اند.

:<math> P = \lim_{T\to \infty} \frac 1 {T} \int_{-\infty}^{\infty} |x_{T}(t)|^2\,dt</math>

به صورت واضح، در حالاتی که عبارت بالا برای P غیرصفر است (حتی وقتیکه T به صورت بدون مرز افزایش می‌یابد) خود انتگرال هم باید بدون مرز افزایش بیابد. این همان دلیلی است که ما نمی‌توانیم از خود چگالی طیف انرژی استفاده کنیم، زیرا انتگرال را در این حالات واگرا ''می‌کند''.

در تحلیل محتوای فرکانس سیگنال <math>x(t)</math>، ممکن است بخواهیم تبدیل فوریه معمول <math>\hat{x}(f)</math> را محاسبه کنیم؛ با این حال، برای بسیاری از سیگنال‌های مورد علاقه تبدیل فوریه به صورت صوری موجود نیست.{{#tag:ref|Some authors (e.g. Risken<ref>{{cite book | title = The Fokker–Planck Equation: Methods of Solution and Applications | edition = 2nd | author = Hannes Risken | publisher = Springer | year = 1996 | isbn = 9783540615309 | page = 30 | url = https://books.google.com/books?id=MG2V9vTgSgEC&pg=PA30 }}</ref>) still use the non-normalized Fourier transform in a formal way to formulate a definition of the power spectral density

:<math> \langle \hat x(\omega) \hat x^\ast(\omega') \rangle = 2\pi f(\omega) \delta(\omega-\omega')</math>,

where <math> \delta(\omega-\omega')</math> is the [[Dirac delta function]]. Such formal statements may sometimes be useful to guide the intuition, but should always be used with utmost care.|group="N"}} صرفنظر از این موضوع، [[:en:Parseval's_theorem|قضیه پارساول]] به ما می‌گوید که ما می‌توانیم توان میانگین را به این صورت بازنویسی کنیم:

:<math> P = \lim_{T\to \infty} \frac 1 {T} \int_{-\infty}^{\infty} |\hat{x}_{T}(f)|^2\,df</math>

آنوقت چگالی طیف توان به صورت ساده به صورت زیرانتگرال بالا تعریف می‌شود.<ref>{{cite book|title=Spikes: Exploring the Neural Code (Computational Neuroscience)|author1=Fred Rieke|author2=William Bialek|author3=David Warland|name-list-style=amp|publisher=[[MIT Press]]|year=1999|isbn=978-0-262-68108-7}}</ref><ref name="millers370">{{cite book|title=Probability and random processes|url=https://archive.org/details/probabilityrando00mill|author1=Scott Millers|author2=Donald Childers|name-list-style=amp|publisher=[[Academic Press]]|year=2012|pages=[https://archive.org/details/probabilityrando00mill/page/n309 370]–5}}</ref>

{{Equation box 1
|indent =
|title=
|equation = {{NumBlk|:|<math> S_{xx}(f) = \lim_{T\to \infty} \frac 1 {T}  |\hat{x}_{T}(f)|^2\,</math>|{{EquationRef|Eq.2}}}}
|cellpadding= 6
|border
|border colour = #0073CF
|background colour=#F5FFFA}}

از اینجا ما می‌توانیم <math> |\hat{x}_{T}(f)|^2</math> را به صورت [[تبدیل فوریه]] از [[هم‌گشت]] زمانی <math> x_{T}^*(-t)</math> و <math> x_{T}(t)</math> ببنیم:

:<math> |\hat{x}_{T}(f)|^2=\mathcal{F}\left \{ x_{T}^*(-t)\mathbf {*}x_{T}(t) \right \}=\int_{-\infty}^\infty \left[\int_{-\infty}^\infty x_{T}^*(t-\tau)x_{T}(t) dt \right]e^{-i 2\pi f\tau} \ d\tau</math>

اکنون، ما همگشت زمانی بالا را بر دوره <math>T</math> تقسیم می‌کنیم، و سپس حد <math>T \rightarrow \infty</math> را می‌گیریم، این تبدیل به تابع [[خودهمبستگی]] برای سیگنال فاقد-پنجره <math> x(t)</math> می‌شود، که توسط <math>R_{xx}(\tau)</math> نشان داده می‌شود، با این شرط که <math>x(t)</math> [[ارگادیسیتی|ارگادیک]] باشد، که این موضوع در بیشتر، اما نه همه، حالات عملی درست است.<ref>The [[Wiener–Khinchin theorem]] makes sense of this formula for any [[wide-sense stationary process]] under weaker hypotheses: <math> R_{xx} </math>  does not need to be absolutely integrable, it only needs to exist. But the integral can no longer be interpreted as usual. The formula also makes sense if interpreted as involving [[Distribution (mathematics)|distributions]] (in the sense of [[Laurent Schwartz]], not in the sense of a statistical [[Cumulative distribution function]]) instead of functions. If <math> R_{xx} </math> is continuous, [[Bochner's theorem]] can be used to prove that its Fourier transform exists as a positive [[Measure (mathematics)|measure]], whose distribution function is F (but not necessarily as a function and not necessarily possessing a probability density).</ref>

:<math> \lim_{T\to \infty} \frac 1 {T} |\hat{x}_{T}(f)|^2 =\int_{-\infty}^\infty \left[\lim_{T\to \infty} \frac 1 {T}\int_{-\infty}^\infty x_{T}^*(t-\tau)x_{T}(t) dt \right]e^{-i 2\pi f\tau} \ d\tau=\int_{-\infty}^\infty R_{xx}(\tau)e^{-i 2\pi f\tau} d\tau </math>

از اینجا، ما می‌بینیم که دوباره با فرض ارگادیک بودن <math>x(t)</math>، «چگالی طیفی توان» به صورت «تبدیل فوریه تابع خودهمبستگی» به دست می‌آید ([[قضیه وینر-خینشین]]).

{{Equation box 1
|indent =
|title=
|equation = {{NumBlk|:|<math>S_{xx}(f)=\int_{-\infty}^\infty R_{xx}(\tau) e^{-i 2 \pi f \tau}\,d \tau=\hat R_{xx}(f)</math>|{{EquationRef|Eq.3}}}}
|cellpadding= 6
|border
|border colour = #0073CF
|background colour=#F5FFFA}}

بیشتر نویسنده‌ها از این معادله برای ''تعریف'' واقعی چگالی طیفی توان استفاده می‌کنند.<ref>{{cite book|title=Echo Signal Processing|author=Dennis Ward Ricker|publisher=Springer|year=2003|isbn=978-1-4020-7395-3|url=https://books.google.com/books?id=NF2Tmty9nugC&q=%22power+spectral+density%22+%22energy+spectral+density%22&pg=PA23}}</ref>

توان سیگنال در یک باند فرکانسی معین <math>[f_1, f_2]</math> که در آن، <math> 0<f_1 < f_2</math> توسط انتگرال‌گیری روی فرکانس محاسبه می‌شود. به این دلیل که <math>S_{xx}(-f) = S_{xx}(f)</math>، یک مقدار برابر از توان را می‌توان به باندهای فرکانسی مثبت و منفی منتسب کرد، که به دلیل فاکتور ۲ در حالتی که در ادامه است محسوب می‌شود (چنین فاکتورهای بدیهی بستگی به رسوم استفاده شده دارند):

:<math> P_\mathsf{bandlimited} = 2 \int_{f_1}^{f_2} S_{xx}(f) \, df</math>

به صورت معمول‌تر، فنون مشابهی را می‌توان برای تخمین چگالی طیفی در زمان-متغیر استفاده کرد. در این حالت، بازه زمانی <math> T</math> محدود است به جای آنکه به سمت بی‌نهایت میل کند. این منجر به پوشش و رزلوشن طیفی کمتری می‌شود زیرا فرکانس‌های کمتر از <math>1/T</math> نمونه برداری نمی‌شوند، و منجر به فرکانس‌هایی می‌شوند که یک ضرب عدد صحیح از <math>1/T</math> مستقل نیستند. فقط به کمک چنین سری زمانی منفردی، طیف توانی تخمین‌زده‌شده بسیار «نویزدار» خواهد بود؛ با این حال، این موضوع اگر بتوانیم مقدار انتظاری را (در معادله بالا) ارزیابی کنیم، این خطاها تسکین می‌یابد، که در این حالت باید از تعداد بزرگی (یا تعداد بی‌نهایتی) از طیف‌های با زمان-اندک متناظر با [[آنسامبل آماری]] که از تحقق‌های <math>x(t)</math> که روی پنجره زمانی معین، ارزیابی شده‌اند، استفاده کرد.

درست مثل چگالی طیفی انرژی، تعریف چگالی طیفی توانی را نیز می‌توان به متغیرهای [[زمان پیوسته و زمان گسسته|زمان گسسته]] <math>x_n</math> تعمیم داد. مثل قبل، ما یک پنجره <math>-N\le n\le N</math> را با سیگنال نمونه‌برداری شده در زمان‌های گسسته <math>x_n=x_0+(n\Delta t)</math> برای دوره اندازه‌گیری کلی <math>T=(2N+1) \Delta t</math> در نظر می‌گیریم.

:<math>S_{xx}(f) = \lim_{N\to \infty}\frac{(\Delta t)^2}{T}\left|\sum_{n=-N}^N x_n e^{-i 2\pi f n \Delta t}\right|^2</math>

توجه کنید که یک تخمین منفرد از PSD از طریق یک تعداد محدود نمونه‌برداری قابل دستیابی است. مثل قبل، PSD واقعی موقعی به دست می‌آید که <math>N</math> (و از این رو <math>T</math>) به سمت بی‌نهایت میل کند، و مقدار انتظاری به صورت صوری اعمال شود. در یک کاربرد جهان واقعی، می‌توان معمولاً PSD های محدود اندازه‌گیری شده را روی چندین سعی میانگین‌گیری می‌کنیم، تا به یک تخمین دقیق تر از PSD نظری از فرایند فیزیکی زیربنای اندازه‌گیری‌های منفرد برسیم. به این PSD محاسبه‌شده گاهی دوره‌سنج (periodogram) گفته می‌شود. این [[:en:Periodogram|دوره‌سنج]] موقعی به PSD درست همگرا می‌شود که تعداد تخمین‌ها و همچنین بازه زمانی میانگین <math>T</math> به سمت بی‌نهایت میل کند (Brown & Hwang).<ref>{{cite book|title=Introduction to Random Signals and Applied Kalman Filtering|author1=Robert Grover Brown|author2=Patrick Y.C. Hwang|name-list-style=amp|publisher=[[John Wiley & Sons]]|year=1997|isbn=978-0-471-12839-7}}</ref>

اگر دو سیگنال هر دو چگالی طیفی توانی را داشته باشند، آنوقت [[:en:Spectral_density#Cross-spectral_density|میان-چگالی طیف]] را به صورت مشابه می‌توان محاسبه کرد؛ همان‌طور که PSD با خودهمبستگی مرتبط است، چگالی طیف-میانی هم با [[میان-همبستگی]] مرتبط است.

==== ویژگی‌های چگالی طیفی توان ====
بعضی از ویژگی‌های PSD شامل:<ref>{{Cite book|publisher=Cambridge University Press|isbn=978-0-521-01230-0|last=Storch|first=H. Von|author2=F. W. Zwiers|title=Statistical analysis in climate research|year=2001}}</ref>
* طیف چگالی همیشه حقیقی و غیرمنفی است، و طیف یک فرایند مقدار حقیقی همچنین یک [[تابع زوج]] از فرکانس است: <math>S_{xx}(-f) = S_{xx}(f)</math>.
* برای یک [[فرایند تصادفی]] پیوسته x(t)، تابع خودهمبستگی R<sub>xx</sub>(t) را می‌توان از طیف توان S<sub>xx</sub>(f) آن با استفاده از [[:en:Fourier_inversion_theorem|وارون تبدیل فوریه]] به دست آورد.
* به کمک [[:en:Parseval's_theorem|قضیه پارساول]]، می‌توان [[واریانس]] (توان میانگین) برای یک فرایند را به کمک انتگرال‌گیری طیف توانی روی همه فرکانس‌ها محاسبه کرد:
:<math>P=\text{Var}(x) = \int_{-\infty}^{\infty}\! S_{xx}(f) \, df</math>
* برای یک فرایند حقیقی x(t) با چگالی طیف توان <math>S_{xx}(f)</math> می‌توان ''طیف انتگرال‌گیری شده'' یا ''توزیع طیف توان'' <math>F(f)</math> را محاسبه کرد، که تعیین‌کننده توان ''باندمحدود'' میانگین موجود در فرکانس‌ها از DC تا f است یعنی به کمک:<ref>An Introduction to the Theory of Random Signals and Noise, Wilbur B. Davenport and Willian L. Root, IEEE Press, New York, 1987, {{isbn|0-87942-235-1}}</ref>
:<math>F(f)=  2 \int _{0}^f S_{xx}(f')\, df'. </math>
:توجه کنید که عبارت قبل برای توان کلی (واریانس سیگنال) حالت خاصی است که در آن <bdi>f→∞</bdi>.

=== چگالی میان-طیف توان ===
اگر دو سیگنال <math>x(t)</math> و <math>y(t)</math> داده شده باشد، که هرکدام دارای چگالی طیف توان <math>S_{xx}(f)</math> و <math>S_{yy}(f)</math> باشند، می‌توان یک ''میان چگالی طبف توان'' (''CPSD'') یا ''میان چگالی طیف'' (''CSD'') تعریف کرد. برای شروع، بیایید توان میانگین برای چنین سیگنال ترکیبی را در نظر بگیریم.

:<math>
\begin{align}
P & = \lim_{T\to \infty} \frac 1 {T} \int_{-\infty}^{\infty} \left[x_{T}(t)+y_{T}(t)\right]^*\left[x_{T}(t)+y_{T}(t)\right]dt
\\
  & = \lim_{T\to \infty} \frac 1 {T} \int_{-\infty}^{\infty}|x_{T}(t)|^2 + x^*_{T}(t)y_{T}(t) + y^*_{T}(t)x_{T}(t)+|y_{T}(t)|^2dt
\\
\end{align}</math>

به کمک نمادگذاری و روش‌های مشابهی با روش‌هایی که برای استخراج چگالی طیف توان استفاده می‌شوند، ما از قضیه پارساول بهره می‌گیریم و این را به دست می‌آوریم:

:<math>S_{xy}(f) = \lim_{T\to\infty} \frac 1 {T} \left[\hat{x}^*_{T}(f)\hat{y}_{T}(f)\right] \ \ \ \ \ \ \ S_{yx}(f) = \lim_{T\to\infty} \frac 1 {T} \left[\hat{y}^*_{T}(f)\hat{x}_{T}(f)\right]</math>

که در آن، دوباره، مشارکت‌های <math> S_{xx}(f)</math> و <math> S_{yy}(f)</math> را از قبل می‌دانستیم. توجه کنید که <math> S^*_{xy}(f)=S_{yx}(f)</math>، از این رو مشارکت کامل توان میانی، به صورت کلی، تشکیل دهنده دوبرابر [[قسمت حقیقی]] از هر '''CPSD''' منفرد است. دوباره مثل قبل، از اینجا ما این ضرب‌ها را به صورت تبدیل فوریه برای همگشت زمانی بازسازی می‌کنیم، که وقتیکه به دوره تقسیم شود، و حد <math> T\to \infty</math>  آن گرفته شود، تبدیل به تبدیل فوریه برای یک تابع میان-همبستگی می‌شود.<ref>{{cite web|author=William D Penny|year=2009|title=Signal Processing Course, chapter 7|url=http://www.fil.ion.ucl.ac.uk/~wpenny/course/course.html}}</ref>

:<math>S_{xy}(f) = \int_{-\infty}^{\infty} \left[\lim_{T\to\infty} \frac 1 {T} \int_{-\infty}^{\infty} x^*_{T}(t-\tau) y_{T}(t) dt \right] e^{-i 2 \pi f \tau} d\tau= \int_{-\infty}^{\infty} R_{xy}(\tau) e^{-i 2 \pi f \tau} d\tau</math>
:<math> S_{yx}(f) = \int_{-\infty}^{\infty} \left[\lim_{T\to\infty} \frac 1 {T} \int_{-\infty}^{\infty} y^*_{T}(t-\tau) x_{T}(t) dt \right] e^{-i 2 \pi f \tau} d\tau= \int_{-\infty}^{\infty} R_{yx}(\tau) e^{-i 2 \pi f \tau} d\tau</math>

که در آن <math>R_{xy}(\tau)</math> برابر [[میان-همبستگی]] برای <math>x(t)</math> با <math>y(t)</math> و <math>R_{yx}(\tau)</math> برابر میان-همبستگی برای <math>y(t)</math> با <math>x(t)</math> است. از این دیدگاه، PSD را می‌توان نوع خاصی از CSD برای <math>x(t) = y(t)</math> دید. برای حالتی که  <math>x(t)</math> و <math>y(t)</math> برابر سیگنال‌های ولتاژ و جریان هستند؛ چگالی طیفی دامنه مرتبط فقط یکی از '''CPSD'''‌ها است که با یک فاکتور دو مقیاس‌دهی شده‌است.

:<math>\operatorname{CPSD}_{\rm Full}=2S_{xy}(f)=2 S_{yx}(f)</math>

برای سیگنال‌های گسسته ''x<sub>n</sub>'' و ''y<sub>n</sub>'' ارتباط بین میان-چگالی طیف و میان-همبستگی به این صورت است:
:<math>S_{xy}(f)=\sum_{n=-\infty}^\infty R_{xy}(\tau_n)e^{-i 2 \pi f \tau_n}\Delta\tau</math>

== تخمین ==
هدف از تخمین چگالی طیف همان [[نظریه تخمین|تخمین]] چگالی طیفی برای یک [[فرایند تصادفی|سیگنال تصادفی]] از ترتیبی از نمونه‌های زمانی است. بسته به آنکه چه چیزی را دربارهٔ سیگنال می‌دانیم، فنون تخمین می‌تواند شامل دیدگاه‌های [[آمار پارامتری|پارامتری]] یا [[آمار ناپارامتری|بدون-پارامتر]] باشد، و همچنین می‌تواند بر اساس تحلیل دامنه-زمان یا دامنه-فرکانس باشد. برای مثال، یک فن پارامتری معمول شامل متناسب‌سازی مشاهدات با یک [[مدل خودهمبسته]] است. یک فن بدون-پارامتر معمول [[:en:Periodogram|دوره‌سنج (periodogram)]] است.

چگالی طیفی معمولاً توسط روش‌های [[تبدیل فوریه]] (مثل [[:en:Welch's_method|روش ولچ]]) تخمین زده می‌شود، اما از دیگر فنون مثل روش [[:en:Maximum_entropy_spectral_estimation|آنتروپی حداکثری]] نیز می‌توان استفاده کرد.

== مفاهیم مرتبط ==
* [[:en:Spectral_centroid|گرانیگاه طیفی]] یک سیگنال برابر نقطه میانی تابع چگالی طیفی آن است، یعنی آن فرکانسی است که توزیع را به دو قسمت مساوی تقسیم می‌کند.
* [[:en:Spectral_edge_frequency|فرکانس حاشیه طیفی]] برای یک سیگنال یک گسترش برای مفهوم قبل به هر نسبت بجای دو قسمت برابر است.
* چگالی طیفی یک تابع از فرکانس، و نه تابعی از زمان است. با این حال، چگالی طیفی از یک پنجره کوچک از یک سیگنال بزرگتر را می‌توان محاسبه کرد، و دربرابر زمان مرتبط با پنجره ترسیم کرد. به این گراف [[طیف‌نگاره]] (spectrogram) گفته می‌شود. این موضوع مبنایی برای تعدادی از فنون تحلیل طیفی مثل [[تبدیل فوریه زمان-کوتاه]] و [[موجک]] است.
* یک «طیف» معمولاً به معنی چگالی طیف توان، به صورت توضیح داده شده در بالا است، که توزیع محتوای سیگنال را روی فرکانس ترسیم می‌کند. این را نباید با [[پاسخ فرکانسی]] یک [[تابع انتقال]] اشتباه کرد که آن هم شامل یک [[فاز (موج)|فاز]] است (یا به صورت معادل، یک بخش حقیقی و موهومی از یک تابع از فرکانس است). برای توابع انتقال (مثل [[نمودار بود]] و [[:en:Chirp#Relation_to_an_impulse_signal|چرپ]]) پاسخ فرکانسی کامل را می‌توان در دو قسمت ترسیم کرد، دامنه دربرابر فرکانس و فاز دربرابر فرکانس- چگالی طیفی فاز، طیف فاز یا فاز طیفی (یا به صورت کمتر معمول، به صورت قسمت‌های حقیقی و موهومی از تابع انتقال). [[پاسخ ضربه]] (در [[حوزه زمان|دامنه زمان]]) <math>h(t)</math> را معمولاً نمی‌توان به صورت یکتا فقط از قسمت چگالی طیفی دامنه و بدون تابع فاز بازیابی کرد. اگرچه این‌ها زوج‌های تبدیل فوریه هستند، تقارنی در آن‌ها وجود ندارند (به آن صورت که برای خودهمبستگی وجود دارد) و این تبدیل فوریه را مجبور نمی‌کند تا مقدار حقیقی داشته باشد. [[:en:Phase_noise|نویز فاز]]، [[:en:Group_delay_and_phase_delay|تاخیر زمانی]] و [[:en:Ultrashort_pulse#Spectral_phase|پالس بسیار کوتاه#فاز طیفی]] را ببینید.

== کاربردها ==
=== مهندسی الکترونیک ===
[[File:Spectrogram-fm-radio.png|thumb|خروجی طیف‌سنج برای یک سیگنال [[پخش رادیویی اف‌ام|رادیوی FM]] که نشان‌دهنده فرکانس (روی محور افقی) و زمان (افزایش یابنده به صورت بالا رونده روی محور عمودی) است.]]
مفهوم و استفاده از طیف توانی یک سیگنال در [[مهندسی برق]] بنیادین و اساسی است، مخصوصاً در [[سامانه ارتباطی|سامانه‌های ارتباطی الکترونیکی]] مثل [[ارتباطات رادیویی]]، [[رادار]]ها، و سامانه‌های مرتبط، بعلاوه فن‌آوری [[سنجش از دور]] پالس منفعل هم طیف توان کاربرد دارد. از ابزار الکترونیکی که [[تحلیل‌گر طیف]] نام دارند برای مشاهده و اندازه‌گیری طیف توانی سیگنال‌ها استفاده می‌شود.

تحلیل‌گر طیف اندازه [[تبدیل فوریه زمان-کوتاه]] (STFT) را برای یک سیگنال ورودی اندازه‌گیری می‌کند. اگر سیگنالی که می‌خواهیم آن را تحلیل کنیم را یک فرایند مانا در نظر بگیریم، آنوقت STFT یک تخمین بخوبی صاف از چگالی طیف توان آن است.

=== کیهان‌شناسی ===
[[نوسانات نخستین]]، یعنی تغییرات چگالی در جهان نخستین، توسط یک «طیف توان» کمی‌سازی می‌شوند، که توان تغییرات را به صورت تابعی از مقیاس فضایی به دست می‌دهد.
== یادداشت ==
{{Reflist|group="N"}}
== پانویس ==
{{پانویس|۲}}

== منابع ==
{{یادکرد-ویکی|پیوند=https://en.wikipedia.org/wiki/Spectral_density|عنوان=Spectral density|زبان=انگلیسی|بازیابی=۲۰ اکتبر ۲۰۲۱}}

{{صوت‌شناسی}}

[[رده:آنالیز فوریه]]
[[رده:امواج]]
[[رده:پراکندگی]]
[[رده:پردازش سیگنال]]
[[رده:تحلیل حوزه فرکانس]]
[[رده:طیف (علوم فیزیکی)]]
[[رده:طیف رادیویی]]
[[رده:طیف‌شناسی]]