Doğrusal sistemler neden sinüzoidal doğruluk gösterir?


9

Sinüsoidal sadakat için bir kanıt arıyorum. DSP'de lineer sistemler hakkında çok çalıştık. Doğrusal sistemler homojen ve katkı maddesidir. Satisifiye ettiği bir başka koşul, bir sinyal sinüs veya cos dalgasıysa, çıkışın sadece faz veya genliği değiştirmesidir. Neden? Bir sinüs dalgası giriş olarak verildiğinde çıkış neden tamamen farklı bir çıkış olamaz?


1
DSP'ye hoş geldiniz. Harika bir soru!
Phonon

5
Anlayışınız eksik. Doğrusal (homojen ve katkı maddesi anlamına gelen) bir sistem, bir giriş sinüzoidin aynı frekansta, ancak muhtemelen farklı genlik ve fazda bir sinüzoid üretme özelliğine sahip olması gerekmez. Sistemin aynı zamanda zamanla değişmez olduğu konusunda daha fazla kısıtlama getirmeniz gerekir. Örneğin, girişx(t) çıktı üretir x(t)cos(2π109t), sistem homojen ve katkı maddesidir ve bu nedenle doğrusaldır, ancak SISO (sinüzoid insinusoid dışarı) özelliğini karşılamaz.
Dilip Sarwate

Dilip (veya birileri) cevap vermelidir: "Yapmazlar." Sadece zamanla değişmeyen doğrusal sistemler yapar.
hotpaw2

2
Sadece bir not, bu soru olurdu ifadesine başka bir yolu olarak "Neden üstel olan özfonksiyonlar Doğrusal zamanla değişmeyen sistemlerin?"
Jason R

Yanıtlar:


8

Diğer cevaplara biraz görsel bir tamamlayıcı

Doğrusal ve zamanla değişmeyen sistemler hakkında konuşuyorsunuz.

Üstel fonksiyonlar bir tuhaf özelliğe sahiptir (ve aslında onun tarafından tanımlanabilir): zaman çevirisi yapmak aynı fonksiyonda bir sabit ile çarpılır. Yani

ett0=et0et

Mathematica grafikleri

Kırmızı üstel, mavi olanın e veya 1 saniye sağa taşındı

Genel olarak, bu aynı zamanda karmaşık üsleri de kapsar

Aklınızdaki gibi karmaşık bir harmonik planını hayal edebiliyor musunuz? x(t)=ej2πt? Eğer öyleyse, bir yay gibi göreceksiniz: zaman geçtikçe karmaşık düzlem boyunca döner.

Mathematica grafikleri

Bu yayın döndürülmesi (birim çemberdeki karmaşık bir sayıyla çarpılması) çevrilmesi ile aynıdır. Muhtemelen hayatınızda bir süredir bu görsel etkiye girdiniz

resim açıklamasını buraya girin

Herhangi bir standart vidanın prensibi de budur.

Bunu lineer zamanla değişmeyen bir sisteme girdiğimizi varsayalım. Bir çıktı alırsınızy Şimdi, bu yayın döndürülmüş bir versiyonunu girin. Doğrusallık nedeniyle, çıktıyaynı miktarda döndürülmüş. Ancak bir dönüş zaman çevirisine eşdeğer olduğundan ve sistem zamanla değişmediğinden, çıktı dayaynı miktarda zaman çevrilmiş. Yani,ygirdiyle aynı özelliği sağlamalıdır: döndürme işlemi belirli bir zaman çevirisine eşdeğer olmalıdır. Bu, yalnızca çıktı orijinal yayın katları olduğunda gerçekleşir.

Ne kadar çeviri? Bu, bir yayla olduğu gibi rotasyonla doğru orantılıdır. Yayın döngüleri ne kadar sıkı olursa (o kadar hızlı döner), belirli bir dönüş için daha az zaman çevirir. Bir vidanın halkaları ne kadar sıkı olursa, tamamen oturması için ne kadar fazla tur yapmanız gerekir. Ve, mermilerin yarısı tamamlandığında, vida içeri doğru yarısı olacak ... Çıktı aynı ilişkiyi yerine getirmek zorunda, yani çıktı yayıy girişle aynı frekansta döner.

Sonunda bir hatırlatma

cos(t)=ejt+ejt2

sin(t)=ejtejt2j

Üstel olarak gerçekleşen şeyin aslında en genel durumda kosinüs ve sinüslerle gerçekleşmesi gerekmez. Ancak sistem de gerçekse bu farklı bir hikaye ...

Genel olarak, aynı akıl yürütme ile, herhangi bir üstel doğrusal zamanla değişmeyen sistemlerin bir "özfonksiyonudur" (çıktı girdiyle orantılıdır). Bu yüzden bu sistemler için Z-dönüşümleri ve Laplace dönüşümleri çok faydalıdır


Bu animasyonu nasıl / nereden aldınız?
Spacey

@Mohammad, Archimedes vidasındaki wikipedia sayfasından aldı
Rojo

Tirbuşon planını nereden buldun? :) math.stackexchange.com/q/144268/2206
Endolit

@ endolith oh Mathematica'da yaptım. Sevgiler daha güzel;)
Rojo

4

Girişli bir sistem düşünün x(t) ve çıktı y(t). Lars1'in cevabından notasyon ödünç alarak, bu ilişkiyi belirtiyoruz x(t)y(t). Sistemin, aşağıdaki özellikleri karşılaması halinde doğrusal bir zamanla değişmeyen (LTI) sistem olduğu söylenir :

H. Eğerx(t)y(t), sonra αx(t)αy(t).

A. Eğer x1(t)y1(t) ve x2(t)y2(t), sonra x1(t)+x2(t)y1(t)+y2(t).

T. Eğer x(t)y(t), sonra x(tτ)y(tτ) herhangi bir gerçek sayı için τ.

Özellikler , H ve bir araya İşletme eşdeğer L

L. ise x1(t)y1(t) ve x2(t)y2(t), sonra αx1(t)+βx2(t)αy1(t)+βy2(t).


Bir zaman içinde değişmez sistemine periyodik giriş periyodik çıkışı üretir
varsayalımx(t)periyodlu bir periyodik sinyaldirT, yani, x(tnT)=x(t) tüm tamsayılar için n. Ardından, T Mülkünden , hemeny(t) aynı zamanda periyodlu bir periyodik sinyaldir T. Böylece, y(t) Fourier serisi olarak:

y(t)=a02+n=1ancos(nωt)+bnsin(nωt)
nerede ω=2π/T temel frekanstır.

Dan beri cos(ωt) ve sin(ωt) periyodik sinyallerdir, lineer olsun olmasın, herhangi bir zamanla değişmeyen sistem için,

cos(ωt)p02+n=1pncos(nωt)+qnsin(nωt)sin(ωt)r02+n=1rncos(nωt)+snsin(nωt).
Aslında, için doğrusal zaman değişmez (LTI) sistemleri, herpn,qn,rn, ve snsıfır haricinde içinp1,q1,r1,s1. Bunun neden böyle olduğunu görmek için LTI sisteminin yanıtını hesaplayalımcos(ωtθ) iki farklı şekilde ve sonuçları karşılaştırın.

Dan beri cos(ωtθ)=cos(θ)cos(ωt)+sin(θ)sin(ωt), Mülkiyet L ve yukarıdaki denklemlerden

cos(ωtθ)p0cos(θ)+q0sin(θ)2+n=1(pncos(θ)+rnsin(θ))cos(nωt)+n=1(qncos(θ)+snsin(θ))sin(nωt).
Öte yandan, cos(ωtθ)=cos(ω(tθ/ω)) sadece gecikmeli bir versiyonudur cos(ωt), T mülkünden bunu alırız
cos(ωtθ)p02+n=1pncos(nωtnθ)+qnsin(nωtnθ)=p02+n=1(pncos(nθ)qnsin(nθ))cos(nωt)+n=1(qncos(nθ)+pnsin(nθ))sin(nωt).
Bu iki Fourier serisi, değeri ne olursa olsun aynı olmalıdır. θBiz seciyoruz. Katsayıları karşılaştırdığımızda, p0/2 eşit olamaz (p0cos(θ)+r0cos(θ))/2 hepsi için θ olmadıkça p0=r0=0. Benzer şekilde, herhangi biri içinn>1, pncos(nθ)qnsin(nθ) eşit olamaz pncos(θ)+rnsin(θ) vb. herkes için θ olmadıkça pn=qn=rn=sn=0. Ancak,n=1, p1cos(θ)q1sin(θ)=p1cos(θ)+r1sin(θ) ima ediyor ki r1=q1ve benzer şekilde, s1=p1. Başka bir deyişle, bir LTI sistemi için,
cos(ωt)p1cos(ωt)+q1sin(ωt)sin(ωt)q1cos(ωt)+p1sin(ωt).
Şimdi, p1cos(ωt)+q1sin(ωt)=Bcos(ωtϕ) nerede B=p12+q12 ve ϕ=arctan(q1/p1). Bu nedenle, T ve H Özellikleri bize
Acos(ωtθ)ABcos(ωtϕθ).
Herhangi bir sinüsoid frekansıω rad / s şu şekilde ifade edilebilir: Acos(ωtθ) uygun seçim için A ve θve böylece yukarıdaki sonuç ihtiyacımız olan şey.

Doğrusal zamanla değişmeyen sistemlerin SISO özelliği: Bir LTI sistemine giriş bir sinüzoid ise, çıkış aynı frekansta ancak muhtemelen farklı genlik ve fazda bir sinüzoiddir.

Bu OP'nin istediği sonuç değil - doğrusal bir sistemin (Özellik H ve A'nın (eşdeğer olarak, Mülkiyet L ) sahip olduğu ancak mutlaka Özellik T'nin sahip olduğu bir sistem) SISO özelliğine sahip olduğunu kanıtlamak istedi, ancak gelişme olarak Yukarıda gösterdiği gibi, T periyodu, periyodik girdinin periyodik çıktı ile sonuçlandığı daha zayıf sonucu kanıtlamak için beklemek zorundadır.


Son bir yorum olarak, SISO özelliğini kanıtlamak için karmaşık sayılar veya evrişim teoremleri veya Fourier veya LaPlace dönüşümleri, impulslar, özfonksiyonlar vb. Kullanılmasının gerekli olmadığını unutmayın. Bunu L ve * T Özellikleri ve trigonometrik kimlik takip eder.

cos(αβ)=cos(α)cos(β)+sin(α)sin(β).

Ne olsaydı x(t)periyodik değil (orantısız frekanslar için periyodik olmayabilir)? İhtiyaçTsonlu olmak? Genel olarak bir şey talep ederek kazanabilir miyizx(t)zaman aralığında kare ile bütünleşebilir mi?
Lars1

@ Lars1 Bir LTI sistemine giriş periyodik değilse , çıktı da periyodik değildir. Belirli bir örnek olarak,x(t)=A1cos(ω1t)+\A2cos(ω2t) nerede ω1/ω2mantıksızdır (ve dolayısıyla girdi periyodik değildir), o zaman L Mülkünden
A1cos(ω1t)+\A2cos(ω2t)A1B1cos(ω1tϕ1)+\A2B2cos(ω2tϕ2)
hangi çıktı da periyodik değildir. Yani sorun yok.
Dilip Sarwate

@Sarwate: Söylemek istediğim pek değil, üzgünüm. Örneğin;x(t)=cos(πt)+cos(2t)yukarıdaki dava tarafından ele alınacaktır. Eğer sonlu bir gözlem zaman aralığına ihtiyaç duyarsaktT=[0;T]herhangi bir kare entegre sinyal gözlem aralığında bir Fourier serisi olarak yazılabilir. Sonlu içinTbu muhtemelen en genel yaklaşımdır ve türevleriniz hala görebildiğim kadarıyla. Fourier serisi yaklaşımı, dışarıdaki periyodikliği zorluyorT ama sadece sinyali önemsersek t\ontbu gerçekten önemli değil.
Lars1

@ Lars1 Dışarıdaki zorunlu periyodikliğin yorumuna katılmıyorum [0,T]önemli değil. Giriş varsax(t) çıktı üretir y(t)bir LTI sisteminde, SISO özelliğini Fourier serisine uygulamak ,y(t) sınırlı [0,T]. Bunun yerine, elde edilen şey periyodik yanıtın bir dönemidiry^(t) periyodik sinyale x^(t) her seferinde anında t, <t<,
x^(t)=x(tmodT).
Başka bir deyişle, T- ikinci bölümü x(t) periyodik olarak tekrarlanan T) zaman ekseni boyunca.
Dilip Sarwate

Örneğin, doğrusal olmayan RF sistemlerinde, girişten çıkışa benzersiz bir frekans eşlemesi sağlamak için çoğu zaman orantısız sinüzoidallerin toplamını seçeriz. Bunlar periyodik olmayan bir sinyalle sonuçlanır ve neden üstte pratik olarak en alakalı sinyalleri hariç tuttuğum periyodikliği varsaymanız gerektiğini merak ettim. Kare entegre edilebilirx(t) ve y(τ)sonlu gözlem aralıklarında Fourier serileri olarak yazılabilir. Ben iddia etmedim (niyetinde)t için aynı aralıkta tanımlandı x ve y BTW ve yzaman ofseti versiyonu olabilir. Daha fazla karışıklığı önlemek için burada duracağım.
Lars1

3

İşte ispat fikri. Bir sistemin çıktısını bir evrişim ile tanımlayabileceğimizi varsayalım,

y(t)=kt(tτ)f(τ)dτ

Dikkat edin fonksiyonu (aka "çekirdek") kt(t)Bunu yazdım burada olabilir değiştirmek olaraktdeğişir. Ancak, genelliklekt(t)- zamanla değişmediği. Buna "doğrusal zaman değişmezliği" denir ( Toeplitz matrislerindeki Wikipedia sayfasına da bakın ). Sistemimiz doğrusal zamanla değişmezse,ktherhangi biri için aynıdır tve böylece aboneliği görmezden gelip yazacağız

y(t)=k(tτ)f(τ)dτ

Şimdi diyelim f(t) bir sinüzoid, söyle f(t)=eiωt. Böylece sahibiz

y(t)=k(tτ)eiωτdτ=k(τ)eiω(tτ)dτ=eiωtk(τ)eiωτdτ

Son denklemin aşağıdakilere bağımlı olmadığına dikkat edint! Sonuç olarak, tanımlayalımK(ω):=k(τ)eiωτdτ.

Böylece,

y(t)=K(ω)eiωt

veya başka bir deyişle, y(t)girişle aynı frekansta salınan ancak karmaşık bir sayıya göre ağırlıklandırılmış bir sinüzoitK(ω) ile ilgili olan t (ve böylece çıktının genliğini ve fazını girdiye göre kaydırabilir).

EDIT: Yorumlar bu cevabın oldukça gevşek olduğunu kaydetti. Amacım, Fourier dönüşümünün farklı formları gibi detaylardan kaçınmaktı, ancak Fourier ve Laplace dönüşümlerini körükledim . Daha önce Fourier dönüşümü dediğim şey sadece Fourier dönüşümü stamamen hayali idi. Bu rotayı açıklamanın çok fazla gösterim ekleyeceğine karar verdim, bu yüzden italik yazıyorum.

Şimdi, Laplace dönüşümünü al, son olarak (Laplace dönüşümü evrimi çoğalmaya götürdüğü için),

Y(s)=K(s)F(s)

Şimdi eğer f bir sinüzoid, söyle f(t)=eiωt, Laplace dönüşümü bir delta fonksiyonudur. ω. Yani, F(s)=δw(s). Böylece, çıkışın Laplace dönüşümü de bu frekansta bir delta işlevidir:

Y(s)=K(s)δω(s)=K(ω)δω(s)

Dan beri K(ω) sadece giriş frekansına, çıkışa bağlı karmaşık bir sayıdır y(t) girişle aynı frekansta, ancak potansiyel olarak farklı genlik ve faza sahip bir sinüsoid olacaktır.

Bu arada, Wikipedia'da zaman alanında yazılmış aynı fikri bulabileceğinizi fark ettim . Daha üst düzey bir açıklama (çok fazla matiği yoksa görmezden gelebilirsiniz), lineer sistem teorisinin Fourier dönüşümü ile köşegenleştirilen evrişim işlemiyle tanımlanmasıdır. Böylece, girişi Fourier dönüşüm operatörünün özvektörü olan bir sistem, girişinin yalnızca ölçekli bir versiyonunu çıkarır.


-1 nedir s ve bununla nasıl bir ilişkisi var? ω? Ve ne demek istediğini açıklayabilir misinδω(s)? DenkleminizY(s)=K(s)δωs)saçmalık.
Dilip Sarwate

@DilipSarwate Fourier gösterimi yerine Laplace dönüşüm gösterimi kullandığından şüpheleniyorum.
Jim Clay

@sydeulissie Sorun, K (w) 'nin sadece karmaşık bir sayı olduğunu iddia etmenizdir, ancak neden her frekansta karmaşık bir sayı olduğunu söylememişsinizdir. İspatın kalbi bu.
Jim Clay

3
Bunun doğru bir taslağı var, ancak ayrıntılarda birçok sorun var. Aşağı oy değil, düzeltilmelidir.
Fonon

1

Diyelim ki girdili bir sistemimiz var x1(t) çıktı üreten y1(t)=G(x1(t))ve bir girişle x2(t) çıktıyı alıyoruz y2(t)=G(x1(t)). Aşağıdaki durumlarda sistem doğrusaldır:

ax1(t)+bx2(t)y(t)=G(ax1(t)+bx2(t))=aG(x1(t))+bG(x2(t))=ay1(t)+by2(t)

nerede a ve b(gerçek veya karmaşık) sabitlerdir. Yukarıdaki denklemler yerine getirilmezse, sistem doğrusal değildir. Denklem, zaman ve frekans alanlarında gerçek ve karmaşık sinyaller için kullanılabilir. Bu, süperpozisyon ilkesinin geçerli olmasıyla aynıdır. Sarwate'in bir yorumda açıkladığı gibi, bu sistemin yeni frekanslar üretmesini engellemez. Muhtemelen sadece dolaylı olarak zamanla değişmezlik varsaymaya alışkınız. Nedeni, bir veya daha fazla harici kontrol sinyali uygulayarak zamanla değişen bir sistemi zamanla değişmeyen bir sistemle eşlemek genellikle mümkündür.

Doğrusallık tanımından ve ayrıca zamanla değişmeyen bir sisteme ihtiyaç duyulduğunda, doğrudan doğrusallık gereksinimine uyurken iki (veya daha fazla sinyalin) müdahale edemeyeceğini ve yeni frekans bileşenleri oluşturamadığını doğrudan görebiliriz. Süperpozisyon ilkesi doğrudan doğrusallık tanımından da kaynaklanmaktadır.

Ayrıca doğrusallık tanımından doğrusal zamanla değişmeyen sistemler için evrişim kavramı gelir. Doğrusal olmayan sistemler için, örneğin, çok boyutlu bir evrişim integrali olan Volterra serisine sahibiz - 1 boyutlu evrişim integrali, Volterra serisinin özel bir örneğidir. Bu doğrusal tekniklerden çok daha karmaşıktır. Ancak, doğrusal bir sistem için evrişim integraline dayanarak, türetme @sydeulissie tarafından gösterileni izler.

Yeni frekansların üretildiği doğrusal olmayan bir ilişkinin basit bir karşı örneğini göstermek için kullanabiliriz G:y(t)=x2(t). Önce bunun gerçekten doğrusal olmadığını gösterelim. Girişi uygularsakx1(t) çıktıyı alıyoruz y1(t)=x12(t) ve girdiyi uygularsak x2(t) çıktıyı alıyoruz y2(t)=x22(t). Çıktıy(t) o zaman:

y(t)={ax1(t)+bx2(t)}2=a2x12(t)+b2x22(t)+2abx1(t)x2(t)

veya:

y(t)=a2y1(t)+b2y2(t)±2aby1(t)y2(t)ay1(t)+by2(t)

ve böylece kanıtladık x2doğrusal olmamak (şaşırtıcı olmayabilir). Tek bir sinüzoidal sinyal uygularsakx(t)=Acos(2πf0t+ϕ0) sisteme G çıktı var:

y(t)=x2(t)=A2cos2(2πf0t+ϕ0)=A22+A22cos(2π2f0t+2ϕ0)

Buradaki çıktı, bir DC bileşeni ve frekansta başka bir bileşen içerir 2f0. Doğrusal olmayan fonksiyonx2 böylece yeni frekans bileşenleri üretir.

Sonuç olarak, doğrusal bir sistemin girdide bulunmayan frekans bileşenleri üretebileceği gözlenebilir (sistem zaman varyantıysa). Sistem doğrusal zamanla değişmezse, çıkış girişte bulunmayan frekans bileşenlerini içeremez.

En alakalı yorum için @Sarwate'e teşekkürler.


Haklısın. Zamanla değişmeyen sistemlere atıfta bulunduğumu belirtmeyi unuttum. Sağladığınız örnek, örneğinizin elinde olmadığı zamana göre değişen bir sistemdir. Normalde böyle bir sinyalcos(t)harici bir bağlantı noktasına sinyal olarak uygulanır, bu durumda doğrusallık yerine getirilmez. Yukarıdaki cevapta zamanla değişmeyen kısmı not ettim.
Lars1

@DilipSarwate Yani sadece LTI sistemleri bu özelliğe sahip mi?
Phonon

Sadece güvenli tarafta olmak için birkaç kitap kontrol ettim. Aslında ayrıntılarda bazı farklar var gibi görünüyor. Yang ve Lee'nin 2007'den itibaren devre sistemleri kitabındaki bir tanım şöyle diyor: "Bir sistemin, süperpozisyon prensibi geçerliyse doğrusal olduğu söylenir, yani, birkaç rasgele girişin doğrusal bir kombinasyonuna çıkışı, çıkışların doğrusal kombinasyonu ile aynıdır. bireysel girişler ". Bu açıdan Sarwate'in örneği doğrusaldır - ancak zamanla değişmez. Diğer referanslar daha az kesindir. @Sarwate sayesinde.
Lars1

1
Lars1 tarafından tipografik hatalar düzeltilerek atıfta bulunulan yorum: Çıktı üreten sistemi düşününx(t)cos(t) girişten x(t). Sonra,ax1(t)+bx2(t) çıktı üretir
(ax1(t)+bx2(t))cos(t)=ax1(t)cos(t)+bx2(t)cos(t)
böylece sistem doğrusaldır, ancak talep edilen özellik yoktur.
Dilip Sarwate

@Sarwate Çıkış x (t) cos (t) zamanını üreten sistem nasıl değişir? Ben DSP'nin bir acemi değilim
Hobyist

1

Dilip Sarwate'in işaret ettiği gibi, sadece lineer kayma-değişmez (LSIV) sistemleri SISO (sinüsoid insinusoid dışarı) özelliğine sahiptir.

Sorunuzun kısa cevabı, karmaşık üslerin eȷωtare özfonksiyonlar bir LSIV sisteminin. Özfonksiyon tanımıyla, eğer girdi özfonksiyon ise (sinüs / cos, Euler formülüne göre karmaşık üstel ile temsil edilebilir), çıktı sadece karmaşık bir sayı olabilen girişin ve karşılık gelen özdeğerin ürünüdür ve bu burada faz / genlik değişimleri gelir.

Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.