Làm cách nào để thêm AWGN vào biểu diễn I và Q của tín hiệu?


16

Tôi có một hệ thống liên lạc không dây mà tôi đang mô phỏng trong Matlab. Tôi đang thực hiện một số hình mờ thông qua điều chỉnh một chút pha của tín hiệu truyền đi. Mô phỏng của tôi lấy các giá trị I (inphase) và Q (bậc hai) ban đầu và thêm vào hình mờ. Sau đó tôi phải mô phỏng tỷ lệ lỗi bit kết quả sau khi được truyền đi. Bây giờ tôi chỉ cần thêm một lượng nhiễu nhiệt khác nhau vào tín hiệu.

Vì tôi có tín hiệu được biểu thị là kênh I và Q của nó, nên dễ dàng nhất để thêm AWGN (nhiễu Gaussian trắng phụ gia) vào I và Q trực tiếp. Một ý nghĩ là thêm nhiễu vào cả hai kênh một cách độc lập, nhưng trực giác của tôi nói với tôi rằng điều này không giống với việc thêm nó vào tín hiệu nói chung.

Vì vậy, làm thế nào tôi có thể thêm tiếng ồn cho nó khi nó ở dạng này?


có lẽ nó sẽ hữu ích hơn nếu bạn có thể cung cấp một số chi tiết về hệ thống liên lạc mà bạn đang mô phỏng.
Rajesh Dachiraju

1
Tôi sẽ giả sử bạn chỉ tạo ra tiếng ồn cho cả I và Q và sau đó thêm chúng. Tôi không thấy lý do tại sao tiếng ồn sẽ tương quan giữa hai.
endolith

1
@endolith, Sự khác biệt tiếng ồn sẽ chỉ xuất hiện trong bộ trộn, bên cạnh đó chúng nên chia sẻ tín hiệu nhiễu của chúng.
Kortuk

Bạn đang nói rằng bạn sẽ thêm nó vào tín hiệu ghép đa phương?
Phonon

@phonon, ý của bạn là gì khi ghép kênh?
Kortuk

Câu trả lời:


14

Có, bạn có thể thêm AWGN phương sai riêng biệt cho mỗi hai thuật ngữ, bởi vì tổng của hai Gaussian cũng là một Gaussian và phương sai của chúng cộng lại . Điều này sẽ có tác dụng tương tự như thêm AWGN có phương sai 2 σ 2 vào tín hiệu ban đầu. Dưới đây là một số giải thích thêm nếu bạn quan tâm.σ22σ2

Một tín hiệu tích có thể được viết bằng các thành phần trong pha và vuông góc của nó như làx(t)=a(t)sin(2πft+φ(t))

x(t)=I(t)sin(2πft)+Q(t)cos(2πft)

trong đó Q ( t ) = a ( t ) sin ( φ ( t ) ) . Nếu bạn muốn thêm AWGN đến tín hiệu ban đầu của bạn là x ( t ) + u ( t ) , nơi u ( t ) ~ N ( μ , σ 2I(t)=a(t)cos(φ(t))Q(t)=a(t)sin(φ(t))x(t)+u(t) , sau đó bạn có thể thêm AWGN cho mỗi điều khoản nhưu(t)N(μ,σ2)

y1(t)=[I(t)sin(2πft)+v(t)]+[Q(t)cos(2πft)+w(t)]

nơi v(t),w(t)N(μ/2,σ2/2)

Cũng lưu ý rằng vì các thuật ngữ cùng pha và cầu phương là phụ gia, AWGN cũng có thể được thêm vào một trong hai thuật ngữ trong biểu diễn của x ( t ) ở trên. Nói cách khác,IQx(t)

y 3 = [ I ( t ) sin ( 2 π f t ) + u ( t ) ] + Q ( t ) cos

y2=I(t)sin(2πft)+[Q(t)cos(2πft)+u(t)]
y3=[I(t)sin(2πft)+u(t)]+Q(t)cos(2πft)

là tương đương thống kê với , mặc dù tôi thích sử dụng y 1 vì tôi không phải theo dõi thành phần nào có tiếng ồn được thêm vào nó.y1y1


Vì tín hiệu có nhiễu, dường như nhiễu sẽ xuất hiện trên cả hai kênh với cường độ ban đầu nhưng bị ảnh hưởng bởi quá trình trộn. Tôi nghĩ rằng quá trình trộn sẽ ảnh hưởng đến nhiễu nhiều hơn sau đó thêm phép trừ bằng cách tách tín hiệu.
Kortuk

1
Tất nhiên, nếu bạn có nhiễu trong tín hiệu để bắt đầu và sau đó tách nó thành các thành phần IQ của nó, thì mỗi tín hiệu sẽ có nhiễu liên quan đến nó. Tuy nhiên, OP đang nói về việc mô phỏng nó trong MATLAB và anh ta có các phần I và Q riêng biệt và muốn biết cách thêm nhiễu vào những phần này để mô phỏng thêm nhiễu vào tín hiệu gốc.
Lorem Ipsum

câu trả lời tốt với nhiều chi tiết, nhưng không trả lời chính xác câu hỏi cơ bản - OP: Bỏ qua trực giác của bạn; thêm WGN trên trục thực với WGN trên trục ảo sẽ dẫn đến WGN phức tạp. Hãy nhớ chia tỷ lệ theo 3dB vì phương sai của tổng là gấp đôi so với các phần (stdv2 = 1.413 stdv1)
Mark Borgerding

2
@Yoda, bạn đã có tất cả dữ liệu, nhưng bạn khiến người đọc đọc qua nhiều phương trình trước khi có câu trả lời. Tôi chỉ đề nghị đặt phần in đậm của bạn trước, sau đó cung cấp các chi tiết hỗ trợ.
Mark Borgerding

1
@yoda, tôi đã mệt khi đọc cái này. Câu trả lời của bạn rất sắc sảo. Cảm ơn bạn đã dành thời gian!
Kortuk

5

Kellenjb đã không trả lời các truy vấn từ Rajesh D và endolith, và không dễ để tìm ra chính xác những gì anh ta cần. Nhưng vì tôi không đồng ý với một số chi tiết của Câu trả lời do yoda và Mohammad đưa ra, tôi sẽ đăng một câu trả lời riêng, trong đó, với lời xin lỗi đến Mark Borgerding, tất cả những thứ hữu ích xuất hiện ở cuối cùng sau tất cả các phương trình nhàm chán.

2BfcB

r(t)=I(t)cos(2πfct)Q(t)sin(2πfct)
I(t)Q(t)B
r(t)=Re{[I(t)+jQ(t)]ej2πfct}
I(t)+jQ(t)

2cos(2πfct+θ)2sin(2πfct+θ)θ=0I(t)Q(t)

r(t)[2cos(2πfct)]=I(t)[2cos2(2πfct)]Q(t)[2sin(2πfct)cos(2πfct)]=I(t)+[I(t)cos(2π(2fc)t)Q(t)sin(2π(2fc)t)]r(t)[2sin(2πfct)]=I(t)[2sin(2πfct)cos(2πfct)]+Q(t)[2sin2(2πfct)]=Q(t)+[I(t)sin(2π(2fc)t)Q(t)cos(2π(2fc)t)]
I(t)Q(t) không bị biến dạng.

Nhiễu băng rộng có ở mặt trước của máy thu và các câu hỏi chính cần được trả lời là những gì xảy ra trong máy thu thực tế và những gì phải được thực hiện để mô phỏng thực tế.

  • x(t)=I(t)+NI(t)y(t)=Q(t)+NQ(t)
    NI(t)NQ(t)
    σ2=N02|H(f)|2df
    H(f)t0NI(t0)NQ(t0)σ2NI(t0)NI(t1) I(t)Q(t)
  • r(t)+ Mfcm
    r[m]=r(m/Mfc)+N[m]=I(m/Mfc)cos(2π(m/M))Q(m/Mfc)sin(2π(m/M))+N[m]
    N[m]Các biến ngẫu nhiên Gaussian trung bình bằng không với phương sai chung có giá trị phụ thuộc vào SNR. Chúng có thể được theo dõi thông qua các bộ trộn và bộ lọc thông thấp trong quá trình mô phỏng chi tiết.
  • Tôi không khuyên bạn nên thêm nhiễu giữa các đầu ra của bộ trộn và các bộ lọc thông thấp. Trong khi có được tiếng ồn giới thiệu ở giai đoạn đó, điều này thường bị choáng ngợp bởi tiếng ồn từ các kết thúc trước đó đang đến qua máy trộn.
  • B1mNI[m]NQ[m]NI[m]NI[m+i] độc lập hoặc không đòi hỏi nhiều suy nghĩ và phân tích, và các chi tiết được biết đến với Kellenjb nhưng không phải với chúng tôi.

Cảm ơn, Dilip. Đẹp chi tiết, thực tế tập trung trả lời.
Jason R

-2

Kellenjb,

Tiếng ồn trong cả I và Q trên thực tế sẽ không phải là gaussian. Trong thực tế, chúng sẽ bắt nguồn từ cùng một vector tiếng ồn ban đầu. Điều này là do chỉ có một vectơ nhiễu bắt đầu với máy thu. Vì vậy, những gì đang xảy ra, là tín hiệu của bạn đi vào máy thu, trong đó AWGN được thêm vào tất nhiên. Tuy nhiên, ngay sau đó, người nhận sẽ chiếu (tín hiệu + nhiễu) lên cơ sở tội lỗi và trên cơ sở cosin, từ đó cung cấp cho bạn các thành phần I và Q của bạn.

Vì vậy, bây giờ nhiễu trong một trong hai nhánh không còn là gaussian nữa, nhưng trên thực tế, là sản phẩm của vectơ cơ bản lần và vectơ nhiễu gốc và sản phẩm của cơ sở cosin nhân với vectơ nhiễu gốc.

Cách tôi khuyên bạn nên mô phỏng điều này, (bạn có đang thực hiện tất cả những điều này trong baseband không?), Chỉ đơn giản là xây dựng cơ sở tội lỗi và cosin, và đơn giản là nhân với (tín hiệu + nhiễu), trong đó 'signal' là tín hiệu ban đầu của bạn Tất nhiên, và sau đó tất nhiên đưa nó xuống baseband sau đó. Trong thực tế một khi bạn lọc để đưa nó xuống dải cơ sở, các vectơ nhiễu của bạn sẽ không phải là màu trắng và không phải là gaussian.

Hi vọng điêu nay co ich! :)

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.