Nol, Pertama, Kedua ... Tahan ke-2

9

Fungsi persegi panjang didefinisikan sebagai:

rect(t)={0if |t|>1212if |t|=121if |t|<12.

Fungsi segitiga didefinisikan sebagai:

tri(t)={1|t|,|t|<10,otherwise
ini merupakan konvolusi dari dua fungsi unit persegi panjang yang identik:
tri(t)=rect(t)rect(t)=rect(τ)rect(tτ) dτ
Nol-order hold dan penggunaan terus Pertama-order fungsi-fungsi ini. Bahkan, ia memiliki:
xZOH(t)=n=x(n)rect(tn) 
untuk Zero-order hold, dan
xFOH(t)=n=x(n)tri(tn) 
untuk Hold orde pertama. Karenatri(t)=rect(t)rect(t) , saya ingin tahu apakah ini hanya kebetulan atau jika, untuk Hold orde kedua, respon impulsnya adalah
tri(t)tri(t)=(rect(t)rect(t))(rect(t)rect(t)).
Apakah benar juga untuk umumk rangka -th terus? Yaitu, letakkan
xKTH(t)=n=x(n)gk(tn) 
managk(tn) adalah respons impuls darik-th order hold, saya ingin tahu apakah respons impulsnya adalah
gk(tn)=(rect(t)rect(t))(rect(t)rect(t)),
k kali.
Menandai
sumber
saya belum melihat referensi untuk th agar ditahan untuk k > 1 . saya akan berharap itu menjadi fungsi rect ( t ) yang dilingkupi dengan dirinya sendiri k - 1 kali. tetapi saya tidak tahu apa definisi itu. kk>1rect(t)k1
robert bristow-johnson
1
@ robertbristow-johnson: Dalam analogi dengan hold-orde-nol (interpolasi polinomial orde nol, yaitu konstanta piecewise), dan hold-orde-pertama (interpolasi polinomial orde-pertama, yaitu linier piecewise), penahanan urutan ke-n adalah interpolasi sedikit demi sedikit oleh polinomial orde ke-n. Disebutkan di sini (hlm. 6).
Matt L.
1
Ini dan apa yang dijelaskan oleh @ robertbristow-johnson dalam jawabannya di bawah ini disebut B-splines.
Olli Niemitalo
Adakah yang bisa menunjukkan dengan matriks gambar dengan faktor 2? Dan, saya tidak begitu jelas tentang faktor di sini.
user30462

Jawaban:

9

Ini bukan kasusnya. Pertama-tama, penahanan urutan kedua akan menggunakan tiga titik sampel untuk menghitung polinomial interpolasi, tetapi respons respons impuls yang Anda sarankan adalah tidak nol dalam interval ukuran 4 (dengan asumsi sampel interval dari T = 1 , seperti yang Anda lakukan pada pertanyaan Anda). Namun, respons impuls yang terkait dengan penahanan orde kedua harus memiliki dukungan panjang 3 .tri(t)tri(t)4T=13

Sekarang Anda bisa menyarankan bahwa orde kedua ditahan bisa memiliki respon impuls yang merupakan konvolusi dari n fungsi persegi panjang. Dalam hal ini Anda akan mendapatkan ukuran dukungan yang benar, tetapi tentu saja itu tidak cukup.nthn

Sebuah orde kedua terus menghitung interpolasi sepotong-bijaksana menggunakan n + 1 titik data berturut-turut. Ini dalam analogi dengan penahanan urutan nol menggunakan titik data tunggal, dan penahanan urutan pertama, yang menggunakan dua titik data. Definisi ini umumnya digunakan dalam literatur (lihat misalnya di sini dan di sini ).nthn+1

Sangat mudah untuk menunjukkan bahwa polinomial orde kedua yang menginterpolasi tiga titik data , y [ 0 ] , dan y [ 1 ] diberikan olehy[1]y[0]y[1]

(1)P(t)=y[1]t(t1)2+y[0](1t2)+y[1]t(t+1)2

Untuk menemukan respon impuls yang mencapai interpolasi yang diberikan oleh , kita harus menyamakan ( 1 ) dengan ekspresi(1)(1)

(2)y[1]h(t+1)+y[0]h(t)+y[1]h(t1)

Jika kita memilih dukungan dari respon impuls sebagai interval [ - 1 , 2 ] , yang setara dengan memilih interval interpolasi [ 0 , 1 ] , menyamakan ( 1 ) dan ( 2 ) menghasilkan impuls berikut respons penahanan orde kedua:h(t)[1,2][0,1](1)(2)

(3)h(t)={12(t+1)(t+2),1<t<01t2,0t112(t1)(t2),1<t<20,otherwise

Respons impuls dari penahanan orde kedua terlihat seperti ini: (3)masukkan deskripsi gambar di sini

Saya menyerahkan kepada Anda untuk menunjukkan bahwa respons impuls ini tidak dapat dihasilkan dengan menggabungkan tiga fungsi persegi panjang satu sama lain.

Matt L.
sumber
Matt, dapatkah Anda memberikan referensi untuk representasi Anda tentang penangguhan pesanan ke-2? Saya 100% yakin plotnya salah.
robert bristow-johnson
saya memperbaiki Persamaan. (1) (dengan asumsi premis sudah benar). saya akan menyerahkannya kepada Anda untuk merefleksikannya menjadi . h(t)
robert bristow-johnson
P(1)=y[1]P(1)=y[1]
(t1)=2
6

nrect(tT/2T)n

x(t)x[n]x(nT)

x(t)=n=x[n] sinc(tnTT)

yang merupakan output dari filter brickwall yang ideal dengan respons frekuensi:

H(f)=rect(fT)={1|f|<12T0|f|>12T

ketika didorong oleh fungsi sampel idealnya

xs(t)=x(t)n=δ(tnTT)=x(t)Tn=δ(tnT)=Tn=x(t)δ(tnT)=Tn=x(nT)δ(tnT)=Tn=x[n]δ(tnT)

xs(t)H(f)x(t)TH(f)1

itu berarti bahwa respons impuls dari brickwall filter yang ideal ini adalah

h(t)=F1{H(f)}=1Tsinc(tT)

x(t)

x(t)=h(t)xs(t)

h(t)

x[n]

xDAC(t)=n=x[n] rect(tnTT2T)

dan dapat dimodelkan sebagai filter dengan respons impuls

hZOH(t)=1Trect(tT2T)

xs(t)

xDAC(t)=hZOH(t)xs(t)

dan respons frekuensi dari filter rekonstruksi tersirat adalah

HZOH(f)=F1{hZOH(t)}=1ej2πfTj2πfT=ejπfTsinc(fT)

catat penundaan setengah sampel yang konstan dalam respons frekuensi ini. dari situlah Zero-order hold berasal.

xs(t)

xDAC(t)xs(t)xDAC(t)

1j2πfTx[n]x[n1]X(z)z1X(z)=X(z)(1z1)

(1z1)(1(ej2πfT)1)=1(ej2πfT)

HFOH(f)=F1{hFOH(t)}=(1ej2πfTj2πfT)2=ej2πfTsinc2(fT)

tanggapan impuls dari ini adalah

hFOH(t)=F{HFOH(f)}=(rect(tT2T))(rect(tT2T))=1Ttri(tTT)

ejπfTsinc(fT)rect(tT2T)

robert bristow-johnson
sumber
Ini akhirnya akan menyatu dengan respon impuls Gaussian, dan saya tidak bisa membuat banyak intuisi tentang hal ini. Saya sangat percaya bahwa penahanan urutan ke-n adalah - dalam analogi lengkap dengan ZOH dan FOH - interpolator polinomial urutan ke-n. Saya membagikan pandangan ini dengan beberapa penulis lain: mis. Yang ini dan yang ini . Saya belum melihat interpretasi Anda tentang pesanan ke-n berlaku di tempat lain.
Matt L.
nn+1n(n1)
2

Pertanyaan lain ditandai sebagai duplikat dari ini. Di sana juga ditanyai apa yang memegang poligon . Itu dan tahan poligon tampaknya sinonim untuk interpolasi linier, di mana "titik-titik terhubung" daripada output yang terlihat seperti gergaji seperti dalam prediksi urutan pertama. Menghubungkan sampel dengan garis memerlukan mengetahui sampel berikutnya terlebih dahulu sehingga garis tersebut dapat diarahkan ke arah yang benar. Dalam konteks sistem kontrol waktu-nyata di mana sampel tidak diketahui sebelumnya, itu berarti bahwa output harus ditunda oleh satu periode pengambilan sampel agar garis terhubung pada sampel.

Polinomial hold (bukan polygonal hold) mencakup penahan orde nol dan penahan orde pertama.

Olli Niemitalo
sumber