Mengapa sistem linear menunjukkan kesetiaan sinusoidal?


9

Saya mencari bukti kesetiaan sinusoidal. Dalam DSP kami banyak belajar tentang sistem linear. Sistem linear bersifat homogen dan aditif. Satu lagi kondisi yang memuaskan adalah bahwa jika sinyal adalah gelombang sinus atau cos maka output hanya mengubah fase atau amplitudo. Mengapa? Mengapa tidak bisa output menjadi output yang sama sekali berbeda ketika gelombang sinus diberikan sebagai input?


1
Selamat datang di DSP. Pertanyaan bagus!
Phonon

5
Pemahaman Anda tidak lengkap. Sistem linier (yang berarti homogen dan aditif) tidak harus memiliki sifat bahwa input sinusoid menghasilkan sinusoid dengan frekuensi yang sama tetapi amplitudo dan fasa yang mungkin berbeda. Anda perlu menerapkan batasan lebih lanjut bahwa sistem ini juga invarian waktu. Misalnya, jika input menghasilkan output , sistemnya homogen dan aditif, dan karenanya linier, tetapi tidak memenuhi SISO (sinusoid in-sinusoid out) ) Properti. x ( t ) cos ( 2 π 10 9 t )x(t)x(t)cos(2π109t)
Dilip Sarwate

Dilip (atau seseorang) harus menjawab: "Tidak." Hanya sistem linier invarian waktu yang melakukannya.
hotpaw2

2
Sama seperti catatan, cara lain untuk mengungkapkan pertanyaan ini adalah "Mengapa fungsi eigen eksponensial dari sistem linear time-invariant?"
Jason R

Jawaban:


8

Pelengkap yang agak visual untuk jawaban lainnya

Anda berbicara tentang sistem yang linear dan tidak berubah-ubah waktu.

Fungsi eksponensial memiliki satu properti khusus (dan dapat benar-benar didefinisikan olehnya): melakukan terjemahan waktu menghasilkan fungsi yang sama dikalikan dengan konstanta. Begitu

ett0=et0et

Grafik Mathematica

Eksponensial merah juga bisa menjadi yang biru dibagi dengan e atau bergerak 1 detik ke kanan

Secara umum, ini juga berlaku untuk eksponensial yang kompleks

Dapatkah Anda membayangkan plot harmonik yang kompleks seperti x(t)=ej2πt ? Jika demikian, Anda akan melihatnya seperti pegas: berputar di sepanjang bidang kompleks seiring berjalannya waktu.

Grafik Mathematica

Memutar pegas itu (mengalikan dengan angka kompleks dalam lingkaran unit) sama dengan menerjemahkannya. Anda mungkin pernah mengalami efek visual ini dalam hidup Anda

masukkan deskripsi gambar di sini

Ini adalah prinsip dari sekrup standar juga.

Asumsikan kita memasukkan ini dalam sistem invarian waktu linier. Anda mendapatkan output Sekarang, masukkan versi yang diputar dari pegas ini. Karena linearitas, output harus y diputar dengan jumlah yang sama. Tetapi karena rotasi setara dengan terjemahan waktu, dan sistem ini invarian waktu, output juga harus y -diterjemahkan oleh waktu dengan jumlah yang sama. Jadi, y harus memenuhi properti yang sama dengan input: memutarnya harus sama dengan terjemahan waktu tertentu. Ini hanya terjadi ketika output adalah kelipatan pegas aslinya.yyyy

Berapa banyak terjemahan? Yah, itu berbanding lurus dengan rotasi seperti yang akan terjadi dengan pegas. Semakin erat loop pegas (semakin cepat berputar), semakin sedikit waktu yang diterjemahkan untuk rotasi tertentu. Semakin erat loop sekrup, semakin banyak putaran yang harus Anda lakukan agar pas sepenuhnya. Dan, ketika setengah dari putaran dilakukan, sekrup akan menjadi setengah jalan ... Output harus memenuhi hubungan yang sama sehingga, pegas keluaran berputar pada frekuensi yang sama dengan input.y

Akhirnya, pengingat

cos(t)=ejt+ejt2

dosa(t)=ejt-e-jt2j

Jadi, hal yang terjadi dengan eksponensial sebenarnya tidak perlu terjadi dengan cosinus dan sinus dalam kasus yang paling umum. Tetapi jika sistemnya juga nyata, itu cerita yang berbeda ...

Secara umum, dengan alasan yang sama, setiap eksponensial merupakan "fungsi eigen" (output sebanding dengan input) dari sistem invarian waktu linier. Itu sebabnya untuk sistem ini transformasi Z dan transformasi Laplace sangat berguna


Bagaimana / Dari mana Anda mendapatkan animasi itu?
Spacey

@Mohammad mengambilnya dari halaman wikipedia pada sekrup Archimedes
Rojo

Di mana Anda mendapatkan plot pembuka botol itu? :) math.stackexchange.com/q/144268/2206
endolith

@endolith oh saya baru saja melakukannya di Mathematica. Anda lebih baik;)
Rojo

4

Pertimbangkan sistem dengan input dan output y ( t ) . Meminjam notasi dari jawaban Lars1, kami menyatakan hubungan ini x ( t ) y ( t ) . Sistem ini dikatakan sebagai sistem linear time-invariant (LTI) jika memenuhi properti berikut:x(t)y(t)x(t)y(t)

H. Jika , maka α x ( t ) α y ( t ) .x(t)y(t)αx(t)αy(t)

A. Jika dan x 2 ( t ) y 2 ( t ) , maka x 1 ( t ) + x 2 ( t ) y 1 ( t ) + y 2 ( t ) .x1(t)y1(t)x2(t)y2(t)x1(t)+x2(t)y1(t)+y2(t).

T. Jika , maka x ( t - τ ) y ( t - τ ) untuk bilangan real τ .x(t)y(t)x(t-τ)y(t-τ)τ

Properti H dan A sama dengan Properti L

L. Jika dan x 2 ( t ) y 2 ( t ) , maka α x 1 ( t ) + β x 2 ( t ) α y 1 ( t ) + β y 2 ( t ) .x1(t)y1(t)x2(t)y2(t)αx1(t)+βx2(t)αy1(t)+βy2(t)


Input periodik ke sistem invarian-waktu menghasilkan output periodik
Misalkan adalah sinyal periodik dengan periode T , yaitu, x ( t - n T ) = x ( t ) untuk semua bilangan bulat n . Kemudian, dari Properti T , maka segera bahwa y ( t ) juga merupakan sinyal periodik dengan periode T . Dengan demikian, kita dapat mengekspresikan y ( t ) sebagai seri Fourier:x(t)Tx(tnT)=x(t)ny(t)Ty(t)

manaω=2π/Tadalah frekuensi dasar.

y(t)=a02+n=1ancos(nωt)+bnsin(nωt)
ω=2π/T

Karena dan sin ( ω t ) adalah sinyal periodik, kami memilikinya untuk setiap sistem invarian waktu, baik linier atau tidak, cos ( ω t )cos(ωt)sin(ωt) Bahkan, untuklinearsistem waktu-invariant (LTI),semuayangpn,qn,rn,dansnadalah nolkecuali untukp1,q1,r1,s

cos(ωt)p02+n=1pncos(nωt)+qnsin(nωt)sin(ωt)r02+n=1rncos(nωt)+snsin(nωt).
pn,qn,rn,sn . Untuk mengetahui mengapa demikian, mari kita hitung respons sistem LTI terhadap cos ( ω t - θ ) dengan dua cara berbeda dan bandingkan hasilnya.p1,q1,r1,s1cos(ωtθ)

Karena , kita dapatkan dari Properti L dan persamaan di atas yang cos ( ω t - θ )cos(ωtθ)=cos(θ)cos(ωt)+sin(θ)sin(ωt) Di sisi lain, karenacos(ωt-θ)=cos(ω(t-θ/ω)) hanyalah versi tertundacos(ωt), dari Property

cos(ωtθ)p0cos(θ)+q0sin(θ)2+n=1(pncos(θ)+rnsin(θ))cos(nωt)+n=1(qncos(θ)+snsin(θ))sin(nωt).
cos(ωtθ)=cos(ω(tθ/ω))cos(ωt)T kita mendapatkan Dua seri Fourier ini harus sama tidak peduli berapa pun nilaiθ yangkita pilih. Membandingkan koefisien, kita melihat bahwa p0/2bisa tidak sama(p0cos(θ)+r0cos(θ
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).
θp0/2 untuk semua θ kecuali p 0 = r 0 = 0 . Demikian pula, untuk n > 1 , p n cos ( n θ ) - q n sin ( n θ ) tidak dapat sama dengan p n cos ( θ ) + r n sin ( θ ) dll. Untuk semua θ kecuali p n = q(p0cos(θ)+r0cos(θ))/2θp0=r0=0n>1pncos(nθ)qnsin(nθ)pncos(θ)+rnsin(θ)θ . Namun, untuk n = 1 , p 1 cos ( θ ) - q 1 sin ( θ ) = p 1 cos ( θ ) + r 1 sin ( θ ) menyiratkan bahwa r 1 = - q 1 , dan juga, s 1 = p 1 . Dengan kata lain, untuk sistem LTI, pn=qn=rn=sn=0n=1p1cos(θ)q1sin(θ)=p1cos(θ)+r1sin(θ)r1=q1s1=p1 Sekarang,p1cos(ωt)+q1sin(ωt)=Bcos(ωt-ϕ) di manaB=
cos(ωt)p1cos(ωt)+q1sin(ωt)sin(ωt)q1cos(ωt)+p1sin(ωt).
p1cos(ωt)+q1sin(ωt)=Bcos(ωtϕ) danϕ=arctan(q1/p1). Karena itu, PropertiesTdanHmemberi kita Acos(ωt-θ)ABcos(ωt-ϕ-θ). Setiapsinusoid frekuensiωrad / s dapat dinyatakan sebagaiAcos(ωtB=p12+q12ϕ=arctan(q1/p1)
Acos(ωtθ)ABcos(ωtϕθ).
ω untuk pilihan A dan θ yang tepat , dan itulah hasil yang kita butuhkan.Acos(ωtθ)Aθ

Properti SISO dari sistem linear time-invariant: Jika input ke sistem LTI adalah sinusoid, outputnya adalah sinusoid dengan frekuensi yang sama tetapi amplitudo dan fasa yang mungkin berbeda.

Ini bukan hasil yang diinginkan OP - dia ingin bukti bahwa sistem linier (yang mana Properti H dan A (ekuivalen, Properti L ) berlaku tetapi tidak harus Properti T ) memiliki properti SISO, tetapi sebagai pengembangan menunjukkan di atas, Properti T harus tahan untuk membuktikan hasil yang lebih lemah bahwa input periodik menghasilkan output periodik.


Sebagai komentar terakhir, perhatikan bahwa tidak perlu menggunakan bilangan kompleks atau teorema konvolusi atau Transformasi Fourier atau LaPlace, impuls, fungsi eigen dll untuk membuktikan properti SISO. Ini mengikuti dari Properties L dan * T dan identitas trigonometri

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

Apa yang akan terjadi jika tidak periodik (tidak periodik bisa terjadi untuk frekuensi tidak seimbang)? Perlu T terbatas? Bisakah kita mendapatkan sesuatu dalam hal generalitas dengan mensyaratkan x ( t ) harus dapat diintegrasikan dengan kuadrat dalam interval waktu pengamatan? x(t)Tx(t)
Lars1

@ Lars1 Jika input ke sistem LTI tidak periodik , outputnya juga tidak periodik. Sebagai kasus khusus, jika mana ω 1 / ω 2 tidak rasional (dan inputnya tidak periodik), maka dari Properti L kami memiliki A 1 cos ( ω 1 t ) + \ A 2 cos (x(t)=SEBUAH1cos(ω1t)+\SEBUAH2cos(ω2t)ω1/ω2 yang outputnya juga tidak periodik. Jadi tidak ada masalah.
SEBUAH1cos(ω1t)+\SEBUAH2cos(ω2t)SEBUAH1B1cos(ω1t-ϕ1)+\SEBUAH2B2cos(ω2t-ϕ2)
Dilip Sarwate

@Sarwate: Tidak cukup apa yang saya katakan, maaf. Bertanya-tanya apakah misalnya akan ditangani oleh kasus di atas. Jika kita memerlukan interval waktu pengamatan yang terbatas dengantT=[0; T]sinyal persegi yang dapat diintegrasikan dapat ditulis sebagai deret Fourier dalam interval pengamatan. UntukTterbatas,ini mungkin pendekatan yang paling umum dan derivasi Anda masih berlaku sejauh yang saya bisa lihat. Jelas seri Fourier mendekati pasukan periodisitas luarTtetapi jika kita hanya peduli tentang sinyalt\ padatini tidak terlalu penting. x(t)=cos(πt)+cos(2t)tT=[0;T]TTt\dit
Lars1

@ Lars1 Saya tidak setuju dengan komentar Anda bahwa periodisitas yang dipaksakan di luar tidak masalah. Jika input x ( t ) menghasilkan output y ( t ) dalam sistem LTI, maka menerapkan properti SISO ke seri Fourier tidak memberikan y ( t ) terbatas pada [ 0 , T ] . Sebaliknya, apa yang diperoleh adalah satu periode dari respon periodik y ( t ) dengan sinyal periodik x ([0,T]x(t)y(t)y(t)[0,T]y^(t) di mana untuk setiap kali instan t , - < t < , x ( t ) = x ( t mod T ) . Dengan kata lain,segmen T-kedua dari x ( t ) diulang secara berkala (dengan periode T ) di sepanjang sumbu waktu. x^(t)t-<t<
x^(t)=x(tmodT).
Tx(t)T
Dilip Sarwate

Misalnya dalam sistem RF nonlinear kita sering memilih jumlah sinusoidal yang tidak seimbang untuk memastikan pemetaan frekuensi yang unik dari input ke output. Ini menghasilkan sinyal non-periodik, dan saya hanya ingin tahu mengapa Anda harus menganggap periodisitas di atas yang bagi saya tampaknya mengecualikan sinyal yang paling relevan secara praktis. Integable persegi dan y ( τ ) dalam interval observasi terbatas dapat ditulis sebagai deret Fourier. Saya tidak (berniat) mengklaim bahwa t didefinisikan pada interval yang sama untuk x dan y BTW dan y bisa menjadi versi offset waktu. Saya akan berhenti di sini untuk menghindari kebingungan lebih lanjut. x(t)y(τ)txyy
Lars1

3

Inilah gagasan buktinya. Mari kita asumsikan kita dapat menggambarkan output suatu sistem dengan konvolusi,

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

kt(t)tkt(t)kt t

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

f(t)f(t)=esayaωt

y(t)=k(t-τ)esayaωτdτ=k(τ)esayaω(t-τ)dτ=esayaωtk(τ)e-sayaωτdτ

tK(ω): =k(τ)e-sayaωτdτ

Jadi, kami telah menemukan itu

y(t)=K(ω)esayaωt

y(t)K(ω)t

s

Sekarang, ambil Transformasi Laplace, untuk mengakhiri (karena Transformasi Laplace membutuhkan konvolusi menjadi multiplikasi),

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

ff(t)=eiωtωF(s)=δw(s)

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

K(ω)y(t)

Kebetulan, saya baru saja memperhatikan Anda dapat menemukan ide yang sama ditulis dalam domain waktu di Wikipedia . Penjelasan tingkat yang lebih tinggi (yang dapat Anda abaikan jika terlalu matematis) adalah bahwa teori sistem linear didefinisikan melalui operasi konvolusi, yang didiagonalisasi oleh transformasi Fourier. Dengan demikian, sistem yang inputnya merupakan vektor eigen dari operator transformasi Fourier hanya akan menghasilkan versi skala inputnya.


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

@DilipSarwate Saya menduga dia menggunakan notasi Transformasi Laplace bukan notasi Fourier.
Jim Clay

@sydeulissie Masalahnya adalah Anda menyatakan bahwa K (w) adalah "hanya bilangan kompleks", tetapi Anda belum mengatakan mengapa itu hanya bilangan kompleks pada setiap frekuensi. Itulah jantung buktinya.
Jim Clay

3
Ini memiliki garis besar yang benar tetapi banyak masalah dalam detailnya. Tidak downvoting, tetapi harus diperbaiki.
Telepon

1

x1(t)y1(t)=G(x1(t))x2(t)y2(t)=G(x1(t))

Sebuahx1(t)+bx2(t)y(t)=G(Sebuahx1(t)+bx2(t))=SebuahG(x1(t))+bG(x2(t))=Sebuahy1(t)+by2(t)

Sebuahb

Dari definisi linearitas dan lebih lanjut yang memerlukan sistem invarian waktu, kita dapat langsung melihat bahwa dua (atau lebih sinyal) tidak dapat mengganggu dan menghasilkan komponen frekuensi baru sambil tetap mematuhi persyaratan linearitas. Prinsip superposisi juga mengikuti langsung dari definisi linearitas.

Juga dari definisi linearitas, konsep konvolusi untuk sistem invarian waktu linier mengikuti. Untuk sistem nonlinier kita misalnya memiliki seri Volterra yang merupakan integral konvolusi multi-dimensi - integral konvolusi 1-dimensi adalah kasus khusus dari seri Volterra. Ini jauh lebih rumit daripada teknik linier sekalipun. Tetapi berdasarkan integral konvolusi untuk sistem linear, derivasi mengikuti yang ditunjukkan oleh @sydeulissie.

G:y(t)=x2(t)x1(t)y1(t)=x12(t)x2(t)y2(t)=x22(t)y(t)

y(t)={Sebuahx1(t)+bx2(t)}2=Sebuah2x12(t)+b2x22(t)+2Sebuahbx1(t)x2(t)

atau:

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

x2x(t)=SEBUAHcos(2πf0t+ϕ0) ke sistem G kami memiliki output:

y(t)=x2(t)=SEBUAH2cos2(2πf0t+ϕ0)=SEBUAH22+SEBUAH22cos(2π2f0t+2ϕ0)

Output di sini berisi komponen DC dan komponen lain di frekuensi 2f0. Fungsi nonlinearx2 sehingga menghasilkan komponen frekuensi baru.

Kesimpulannya dapat diamati bahwa sistem linear dapat menghasilkan komponen frekuensi yang tidak ada dalam input (jika sistem adalah varian waktu). Jika sistem invarian waktu linier, output tidak dapat memasukkan komponen frekuensi yang tidak ada dalam input.

Terima kasih kepada @Sarwate untuk komentar paling relevan.


Kamu benar. Saya lupa menyebutkan bahwa saya merujuk ke sistem invarian waktu. Contoh yang Anda berikan adalah sistem yang bervariasi waktu di mana contoh Anda tidak berlaku. Biasanya sinyal seperti itucos(t)diterapkan pada port eksternal sebagai sinyal di mana linieritas tidak terpenuhi. Saya telah mencatat bagian waktu yang tidak berubah dalam jawaban di atas.
Lars1

@DilipSarwate Jadi, apakah hanya sistem LTI yang memiliki properti itu?
Telepon

Hanya memeriksa beberapa buku untuk berada di sisi yang aman. Sebenarnya sepertinya ada beberapa perbedaan dalam detailnya. Satu definisi dalam buku Yang dan Lee tentang sistem sirkuit dari 2007 mengatakan: "Suatu sistem dikatakan linier jika prinsip superposisi berlaku, yaitu outputnya ke kombinasi linear dari beberapa input acak adalah sama dengan kombinasi linear dari output ke input individu ". Dalam hal itu, contoh Sarwate adalah linier - tetapi tidak invarian waktu. Wasit lain kurang tepat. Terima kasih kepada @Sarwate.
Lars1

1
Komentar yang dirujuk oleh Lars1 dengan kesalahan ketik diperbaiki: Pertimbangkan sistem yang menghasilkan outputx(t)cos(t) dari input x(t). Kemudian,Sebuahx1(t)+bx2(t) menghasilkan output
(Sebuahx1(t)+bx2(t))cos(t)=Sebuahx1(t)cos(t)+bx2(t)cos(t)
sehingga sistemnya linear tetapi tanpa properti yang diklaim.
Dilip Sarwate

@Sarwate Bagaimana sistem yang menghasilkan output x (t) cos (t) waktu bervariasi? Saya seorang pemula di DSP
Hobyist

1

Seperti yang Dilip Sarwate tunjukkan, hanya sistem linear shift-invariant (LSIV) yang memiliki properti SISO (sinusoid in-sinusoid out).

Jawaban singkat untuk pertanyaan Anda adalah eksponensial kompleks eȷωtAre adalah fungsi eigen dari sistem LSIV. Dengan definisi fungsi eigen, jika inputnya adalah fungsi eigen (sinus / cos dapat diwakili oleh eksponensial kompleks sesuai dengan rumus Euler), output hanyalah produk dari input dan nilai eigen yang sesuai, yang bisa berupa bilangan kompleks, dan itu dari mana perubahan fase / amplitudo berasal.

Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.