Mối liên hệ giữa tín hiệu tương tự với tỷ lệ nhiễu và tỷ lệ tín hiệu với nhiễu trong mặt phẳng IQ trong hệ thống giải điều chế cầu phương là gì?


9

Chúng tôi muốn tính toán mối quan hệ định lượng giữa nhiễu tương tự gần tần số LO và thống kê các điểm được tìm thấy trong mặt phẳng IQ sau khi giải điều chế IQ. Để hiểu hoàn toàn câu hỏi, trước tiên chúng tôi đưa ra một mô tả chi tiết về hệ thống giải điều chế IQ.

Hệ thống giải điều chế IQ

Bộ trộn IQ nhận tín hiệu ở tần số cao và đưa chúng đến tần số thấp hơn để chúng có thể được xử lý dễ dàng hơn. Hình 1 cho thấy một sơ đồ của bộ trộn IQ. Tín hiệu dao động cục bộ (LO)cos(Ωt) được sử dụng để trộn tín hiệu RF xuống tần số thấp hơn.

Hệ thống giải điều chế IQ hoàn chỉnh Hình 1: Chuỗi xử lý tín hiệu hoàn chỉnh. Tín hiệu tần số vi sóng (và nhiễu) đi vào bộ trộn IQ thông qua cổng RF. Tín hiệu này được trộn với bộ tạo dao động cục bộ (LO) để chuyển đổi thành tín hiệu tần số trung gianIQ. Các tín hiệu tần số trung gian sau đó được lọc để loại bỏ thành phần tần số cao còn lại (xem văn bản) và được lấy mẫu kỹ thuật số. Việc phát hiện biên độ và pha của từng thành phần tần số được thực hiện thông qua biến đổi Fourier rời rạc trong logic kỹ thuật số.

Tín hiệu mạch lạc - trường hợp dc

Giả sử trong tín hiệu RF đến là Mcos(Ωt+ϕ). Sau đóIQ tín hiệu sẽ là

I(t)=M2cos(ϕ)+M2cos(2Ωt+ϕ)Q(t)=M2sin(ϕ)M2sin(2Ωt+ϕ).
Chúng tôi chuyển các tín hiệu này qua các bộ lọc thông thấp để loại bỏ 2Ω điều khoản, năng suất
IF(t)=M2cos(ϕ)QF(t)=M2sin(ϕ).
Như chúng ta có thể thấy, các dc IQđiện áp có thể được coi là tọa độ Descartes cho biên độ và pha của tín hiệu gốc. Do đó, bộ trộn đã thực hiện nhiệm vụ của mình là cho phép chúng ta tìm biên độ và pha của tín hiệu tần số cao bằng cách chỉ thực hiện các phép đo tần số thấp.

Tín hiệu mạch lạc - trường hợp ac

Trong thực tế, chúng ta thường không giải điều chế tín hiệu RF thành dc. Cái này có một vài nguyên nhân:

  1. Mật độ phổ nhiễu hầu như luôn luôn giảm mạnh ở tần số thấp.

  2. Nếu chúng ta muốn đo đồng thời biên độ và pha của một số thành phần hình sin ở các tần số khác nhau, chúng ta không thể giải điều chế trực tiếp thành dc trong phần tương tự của hệ thống.

Như một ví dụ liên quan đến # 2, chúng ta có thể có

RF(t)=M1cos([Ω+ω1]t+ϕ1)+M2cos([Ω+ω2]t+ϕ2).
Để tìm biên độ và pha của cả hai thành phần tần số, chúng ta phải sử dụng xử lý tín hiệu phức tạp hơn một chút. CácIFQF dạng sóng trong trường hợp này là
IF(t)=M12cos(ω1t+ϕ1)+M22cos(ω2t+ϕ2)QF(t)=M12sin(ω1t+ϕ1)M22sin(ω2t+ϕ2).
Để tìm cả biên độ và cả hai pha, về cơ bản chúng ta cần thực hiện một phép biến đổi Fourier. Để làm điều này, chúng tôi số hóa các dạng sóng, mang lại
In=M12cos(ω1nδt+ϕ1)+M22cos(ω2nδt+ϕ2)Qn=M12sin(ω1nδt+ϕ1)M22sin(ω2nδt+ϕ2)
Ở đâu δtlà khoảng thời gian lấy mẫu kỹ thuật số. Sau đó, trong logic kỹ thuật số, chúng tôi xây dựng chuỗi phức tạpzn Được định nghĩa bởi znIn+iQn. Đối với các tín hiệu được viết ở trên, đây là
zn=M12exp(i[ω1nδt+ϕ1])+M22exp(i[ω2nδt+ϕ2]).
Nếu chúng ta bây giờ, trong logic kỹ thuật số, hãy tính tổng
Z(ωk)=1Nn=0N1zneiωknδt
chúng tôi phục hồi biên độ và pha cho thành phần ở tần số ωk. Ví dụ, nếu chúng ta tính toánZ(ω1) chúng ta sẽ nhận được (M1/2)exp(iϕ1).

Tiếng ồn

Trong thực tế tín hiệu luôn đi kèm với nhiễu. Ảnh hưởng của tiếng ồn là tạo raZ(ω)một biến ngẫu nhiên thay vì một giá trị xác định. Nói cách khác, cho mỗi ω Z(ω)là ngẫu nhiên và sẽ khác nhau cho mỗi lần thực hiện thí nghiệm.

Chúng ta có thể đoán từ trực giác rằng trong sự hiện diện của tiếng ồn Z(ω) có phân bố đối xứng tròn trong mặt phẳng IQ với giá trị trung bình bằng giá trị xác định (M/2)exp(iϕ). Câu hỏi là chính xác phân phối thống kê củaZ trong sự hiện diện của tiếng ồn?

Câu trả lời:


9

Bởi vì mỗi bước trong chuỗi xử lý là tuyến tính, chúng tôi xem xét một trường hợp chỉ có nhiễu và không có tín hiệu kết hợp. Biểu thị tiếng ồnξ(t). CácIQ tín hiệu là

 I(t)=ξ(t)cos(Ωt)Q(t)=ξ(t)sin(Ωt).
Chúng tôi biểu thị hiệu ứng của bộ lọc dưới dạng tích chập với chức năng đáp ứng thời gian h,
IF(t)=dtξ(t)cos(Ωt)h(tt)
và tương tự cho QF. Lưu ý rằng, vì bộ lọc là nhân quả,h(t)=0 cho t<0. Việc lấy mẫu chỉ đơn giản là chọn giá trị củaIFQF vào những thời điểm {nδt},
In=dtξ(t)cos(Ωt)h(nδtt)
và tương tự cho Qn. Sau khi xây dựng được mô tả ở trên cho phần kỹ thuật số của chuỗi xử lý, chúng tôi có
Z(ω)=1Nn=0N1dtξ(t)eiΩth(nδtt)eiωnδt.
Do đó, vấn đề của chúng tôi là tính toán số liệu thống kê của biểu thức này.

Thay đổi biến nδttt sản xuất

Z(ω)=1Nn=0N1dtξ(nδtt)eiΩ(nδtt)h(t)eiωnδt.
Ở giai đoạn này, chúng ta có thể kiểm tra độ tỉnh táo bằng cách tính giá trị trung bình của Z(ω). Hãy nhớ rằng, đây là một quần thể trung bình. Nói cách khác, chúng tôi đang tính toán giá trị trung bình củaZ(ω)mà chúng ta sẽ tìm thấy bằng cách chuyển đổi nhiều trường hợp nhiễu được giải điều chế thành các điểm IQ và sau đó lấy giá trị trung bình của tất cả các điểm đó. Trong mọi trường hợp, kết quả là
Z(ω)=1Nn=0N1dtξ(nδtt)0eiΩ(nδtt)h(t)eiωnδt=0.
Điều này có ý nghĩa như chúng ta mong đợi tiếng ồn không nên thay đổi giá trị trung bình của điểm IQ được giải điều chế, mà chỉ nên thêm một số tính ngẫu nhiên tập trung vào giá trị xác định.

Tôi không biết làm thế nào để tính toán số liệu thống kê của Z(ω) trực tiếp, vì vậy chúng tôi thực hiện một cách tiếp cận khác bằng cách tính toán thay vì bình phương trung bình của Z(ω). Theo định lý giới hạn trung tâm, phần thực và phần ảo củaZ nên có ít nhất xấp xỉ Guassian phân phối (và, như chúng tôi đã chỉ ra, không tương quan) để tìm mô đun bình phương trung bình của Z thực sự nói với chúng ta tất cả những gì chúng ta cần biết

Chúng tôi tiến hành bằng cách trực tiếp xây dựng |Z(ω)|2 và lấy trung bình thống kê (trung bình thống kê được ký hiệu là ).

|Z(ω)|2=dtdt1N2n,m=0N1eiΩ(tt)h(t)h(t)ξ(nδtt)ξ(mδtt)ei(Ω+ω)(nm)δt.()
Bây giờ chúng ta sử dụng định lý Wiener-Khinchin nói rằng đối với một quá trình ngẫu nhiên đứng yên ξ(t) trung bình thống kê ξ(τ)ξ(0) có liên quan đến mật độ phổ công suất Sξ thông qua phương trình sau:
ξ(τ)ξ(0)=12dω2πSξ(ω)eiωτ.
Sử dụng công thức này cho ξ(nδtt)ξ(mδtt) sản lượng
|Z(ω)|2=12dtdtdω2π1N2n,m=0N1eiΩ(tt)h(t)h(t)Sξ(ω)eiω((nm)δt(tt))ei(Ω+ω)(nm)δt=12dω2π|h(ωΩ)|2Sξ(ω)|1Nn=0N1ei(Ω+ωω)nδt|2=12Ndω2π|h(ωΩ)|2Sξ(ω)1N(sin([Ω+ωω]δtN/2)sin([Ω+ωω]δt/2))2Nth order Fejer kernel=12Ndω2π|h(ωΩ)|2Sξ(ω)FN([Ω+ωω]δt/2)
Ở đâu FNNthđặt hàng nhân Fejer . Thay đổi biếnΩωω chúng tôi nhận được
|Z(ω)|2=12Ndω2π|h(ω)|2Sξ(Ωω)FN([ω+ω]δt/2).
Cho đến nay các kết quả đã được chính xác và kết quả chính xác có thể được tìm thấy bằng cách đánh giá số các tích phân. Bây giờ chúng tôi thực hiện một loạt các giả định tương đối yếu để đi đến một công thức thực tế. Nhân FejerFN(x) có trọng lượng tập trung gần x=0. Do đó, chúng tôi hòa nhập hơnSξ chỉ dành cho tần số gần Ω và vì vậy, trong tích phân này, chúng ta có thể tính gần đúng Sξ như một hằng số S(Ωω)Sξ(Ω), cho
|Z(ω)|2=12NSξ(Ω)dω2π|h(ω)|2FN([ω+ω]δt/2).
Ở đây chúng ta có thể thấy rằng các thống kê nhiễu của điểm IQ được giải điều chế chỉ phụ thuộc vào mật độ phổ RF gần tần số LO. Điều này thật ý nghĩa; bộ trộn IQ được thiết kế để đưa nội dung tín hiệu gần tần số LO và đưa nó xuống mức IF thấp hơn, nơi nó có thể được xử lý. Các bộ lọc khử răng cưa loại bỏ tất cả các thành phần tần số ở quá xa LO.

Null đầu tiên của FN(x) xảy ra tại x=2π/Nvà hầu hết trọng lượng được chứa trong một vài thùy đầu tiên. Các null đầu tiên là do

ωnull2π=ω2π±1Nδt.
Điều này có nghĩa là tích phân trên ω bị chi phối bởi các tần số trong một phạm vi được cho bởi tần số lấy mẫu chia cho N. Trong hầu hết các ứng dụng thực tế, phạm vi này rất nhỏh(ω)gần như không đổi trong phạm vi này. Nếu đó là trường hợp, chúng ta có thể thay thếh(ω) với h(ω) (lưu ý rằng h(ω)=h(ω)) Phát hiện
|Z(ω)|2=12NSξ(Ω)|h(ω)|2dω2πFN([ω+ω]δt/2N)1/δt=Sξ(Ω)2T|h(ω)|2
Ở đâu TNδt là tổng thời gian đo.

Tỷ lệ tín hiệu trên tạp âm

Một cách hợp lý cũng được biết rằng nếu một biến ngẫu nhiên Z có Gaussian và các phần thực và ảo phân phối độc lập, và có mô đun bình phương trung bình R, sau đó các phân phối của phần thực và phần ảo của biến đó có độ lệch chuẩn R/2.[a] Do đó, lấy kết quả của chúng tôi cho |Z(ω)|2, quan sát của chúng tôi rằng các phần thực và tưởng tượng của Z được phân phối theo Gaussian và thực tế là chúng không tương thích,[b] chúng ta biết rằng độ lệch chuẩn của các bản phân phối của phần thực và phần ảo là

σ=Sξ(Ω)|h(ω)|2/4T.
Như đã thảo luận ở đầu, một tín hiệu Mcos([Ω+ω]t+ϕ) trở thành (M/2)eiϕtrong mặt phẳng IQ. Tất nhiên, ở đó chúng tôi đã bỏ qua ảnh hưởng của bộ lọc mà chỉ đơn giản là mở rộng biên độ thành
Z(ω)=M|h(ω)|2eiϕ.
Giả sử, như được minh họa trong Hình 2, chúng tôi đang sử dụng hệ thống giải điều chế IQ để phân biệt giữa hai hoặc nhiều tín hiệu, mỗi tín hiệu có một pha khác nhau nhưng có cùng biên độ M. Do nhiễu, mỗi biên độ / pha có thể dẫn đến một đám mây điểm trong mặt phẳng IQ với khoảng cách xuyên tâmM|h(ω)|/2từ nguồn gốc. Khoảng cách giữa tâm của hai đám mây làg(M/2)|h(ω)| Ở đâu glà một yếu tố hình học phụ thuộc vào các giai đoạn của các đám mây. Nếu góc vòng cung giữa hai đám mây làθ và trung tâm của mỗi đám mây là tương đương từ nguồn gốc sau đó g=2sin(θ/2). Ví dụ: nếu hai giai đoạn là±π/2 sau đó g=2sin(π/2)=2. Về mặt hình học, điều này là do khoảng cách giữa các trung tâm của đám mây lớn hơn gấp đôi so với khoảng cách của một trong hai đám mây từ gốc.

Tỷ lệ tín hiệu / nhiễu (SNR) là

SNRseparation22×(cloud std deviation)2=(gM|h(ω)|/2)22Sξ(Ω)|h(ω)|2/4T=(gM)2T2Sξ(Ω)=g2PTSξ(Ω).
Ở đâu PM2/2là công suất tương tự đến. Lưu ý rằng SNR không phụ thuộc vàoh. Để ghi nhớ kết quả này, lưu ý rằng công suất nhiễu là mật độ phổ nhân với băng thôngB. Đang lấyB=1/T chúng tôi thấy rằng kết quả của chúng tôi chỉ nói rằng SNR trong mặt phẳng IQ chính xác bằng SNR tương tự nhân với yếu tố hình học g2.

nhập mô tả hình ảnh ở đây

Hình 2: Hai đám mây IQ. Sự tách biệt giữa tâm của các đám mây tỷ lệ thuận với cường độ xuyên tâm của chúngM, nhưng được thu nhỏ bởi một yếu tố hình học g. Được chiếu lên đường nối giữa tâm của chúng, mỗi đám mây trở thành phân bố Gaussian với chiều rộngSξ(Ω)|h(ω)|2/4T.

[a]: Tra cứu phân phối vuông chi .

[b]: Chúng ta có thể thấy rằng phần thực và phần ảo của Z trong thực tế không tương quan bằng cách viết tương đương của phương trình () nhưng cho ZZ. Làm điều này, chúng ta sẽ thấy rằng tổng đã biến thành kernel Fejer trong trường hợp cho|Z|2 sẽ về không (ít nhất là xấp xỉ) bởi vì nó sẽ gần như là sự chồng chéo của một sin và cos, là trực giao.

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.