IMRR（镜像抑制比）计算和仿真 (2)

最新推荐文章于 2024-05-15 23:34:34 发布

笨牛慢耕

最新推荐文章于 2024-05-15 23:34:34 发布

阅读量4.5k

点赞数 6

分类专栏：信号处理与数字通信文章标签： matlab 镜像抑制比 IQ imbalance 射频信号处理

本文链接：https://blog.csdn.net/chenxy_bwave/article/details/120054885

版权

信号处理与数字通信专栏收录该内容

14 篇文章 104 订阅

订阅专栏

本文通过MATLAB代码展示了如何进行镜像抑制比(IMRR)的计算和仿真。作者首先介绍了仿真背景，然后提供了详细的MATLAB代码实现，包括信号生成、IQ不平衡模拟和IMRR的估计、理论计算以及近似计算。仿真结果显示理论计算与估计结果高度一致，证实了代码的有效性。最后，给出了不同增益和相位不平衡下IMRR的二维仿真结果。

摘要由CSDN通过智能技术生成

1. 前言

接着上一篇（IMRR（镜像抑制比）计算和仿真 (1)），我们将做一个简单的matlab仿真，看看仿真结果与上述理论推导结果是否会保持一致。

2. 代码(matlab)

% 改编自：http://www.dsplog.com/2013/01/31/imrr-transmit-iq-gain-phase-imbalance/

clear; close all; clc

debug = 1; % Enable the debug mode to visualize the spectrum, etc.

N    = 256; % Simulation data block length
Fsig = 20;  % Modulation signal frequency, i.e, the frequency of Complex sinusoid
gammadB_array = [-3:.1:3]; % Q/I amplitude difference in unit of dB
phiDeg_array  = [-6:.2:6]; % Q/I phase difference in unit of degree
[~, gammadB_zeroIdx ] = min(abs((gammadB_array-0)));
[~, phiDeg_zeroIdx  ] = min(abs((phiDeg_array-0)));

if(debug == 1)
    gammadB_array = -3;
    phiDeg_array  = -6;
end

for (ii = 1:length(gammadB_array))
    for (jj = 1:length(phiDeg_array))
        gammadB    = gammadB_array(ii); 
        phiDeg     = phiDeg_array(jj);
        gammaLin   = 10^(gammadB/20); % Convert dB to linear domain
        phiRad     = phiDeg*pi/180;   % Convert degree to radian, because sin/cos accept rad as input
        alphaLin   = gammaLin -1 ;
        
        % transmitted signal
        xt        = exp(sqrt(-1)*2*pi*Fsig*[0:N-1]/N); % Complex sinusoid
        % received signal with IQ imbalance, after quadrature demodulation
        xht_re    = real(xt) - gammaLin*sin(phiRad)*imag(xt);
        xht_im    = gammaLin*cos(phiRad/2)*imag(xt);
        xht       = complex(xht_re,xht_im);  % "h" mean hat.
        
        if(debug)
            fvtool(xt); title('Transmitter baseband signal spectrum');
            fvtool(xht);title('Receiverer baseband signal spectrum');
        end
        
        % Take fft() to find the +Fsig(desired signal) and -Fsig(image due to I/Q mismatch) components
        yF        = fft(xht,N);
        y_pfm     = yF(Fsig+1);
        y_nfm     = yF(N-Fsig+1);
        
        est_imrr_lin    = (abs(y_pfm)./abs(y_nfm))^2;
        theory_imrr_lin = (gammaLin^2 + 1  + 2*gammaLin*cos(phiRad))/(gammaLin^2 + 1  - 2*gammaLin*cos(phiRad));
        approx_imrr_lin = 4/(alphaLin^2 + phiRad^2);
        
        est_imrr_dB(ii,jj)    = 10*log10(est_imrr_lin);
        theory_imrr_dB(ii,jj) = 10*log10(theory_imrr_lin);
        approx_imrr_dB(ii,jj) = 10*log10(approx_imrr_lin);
    end
end

figure;
plot(gammadB_array,theory_imrr_dB(:,phiDeg_zeroIdx),'bs-'); hold on
plot(gammadB_array,est_imrr_dB(:,phiDeg_zeroIdx),'md-');
plot(gammadB_array,approx_imrr_dB(:,phiDeg_zeroIdx),'gx-');
xlabel('gain imbalance, dB'); ylabel('image rejection, dB'); grid on;
legend('theory','estimated','approx')
title('Image Rejection Ratio with gain imbalance alone');
axis([-3 3 10 50]);

figure;
plot(phiDeg_array,theory_imrr_dB(gammadB_zeroIdx,:),'bs-'); hold on
plot(phiDeg_array,est_imrr_dB(gammadB_zeroIdx,:),'md-');
plot(phiDeg_array,approx_imrr_dB(gammadB_zeroIdx,:),'gx-');
xlabel('phase imbalance, degree'); ylabel('image rejection, dB'); grid on;
legend('theory','estimated','approx')
title('Image Rejection Ratio with phase imbalance alone');
axis([-6 6 10 50]);

figure;
[x,y] = meshgrid(gammadB_array,phiDeg_array);
mesh(x,y,est_imrr_dB); title('IMRR vs {gammadB, phiDeg}');
axis([-3 3 -6 6 10 60]);