iqcoef2imbal

Преобразуйте коэффициент компенсатора в неустойчивость фазы и амплитуду

Синтаксис

[A,P] = iqcoef2imbal(C)

Описание

пример

[A,P] = iqcoef2imbal(C) преобразовывает коэффициент компенсатора C в его эквивалентную амплитуду и неустойчивость фазы.

Примеры

свернуть все

Используйте iqcoef2imbal, чтобы оценить амплитуду и неустойчивость фазы для данного комплексного коэффициента. Коэффициентами является вывод от функции step IQImbalanceCompensator.

Создайте повышенный фильтр передачи косинуса, чтобы сгенерировать 64-QAM сигнал.

M = 64;
txFilt = comm.RaisedCosineTransmitFilter;

Модулируйте и отфильтруйте случайные 64-ary символы.

data = randi([0 M-1],100000,1);
dataMod = qammod(data,M);
txSig = step(txFilt,dataMod);

Задайте неустойчивость фазы и амплитуда.

ampImb = 2; % dB 
phImb = 15; % degrees

Примените заданную неустойчивость I/Q.

gainI = 10.^(0.5*ampImb/20);
gainQ = 10.^(-0.5*ampImb/20);
imbI = real(txSig)*gainI*exp(-0.5i*phImb*pi/180);
imbQ = imag(txSig)*gainQ*exp(1i*(pi/2 + 0.5*phImb*pi/180));
rxSig = imbI + imbQ;

Нормируйте степень полученного сигнала

rxSig = rxSig/std(rxSig);

Удалите неустойчивость I/Q с помощью Системного объекта comm.IQImbalanceCompensator. Установите объект компенсатора, таким образом, что комплексные коэффициенты сделаны доступными как выходной аргумент.

hIQComp = comm.IQImbalanceCompensator('CoefficientOutputPort',true);
[compSig,coef] = step(hIQComp,rxSig);

Оцените неустойчивость от последнего значения коэффициента компенсатора.

[ampImbEst,phImbEst] = iqcoef2imbal(coef(end));

Сравните предполагаемые значения неустойчивости с заданными единицами. Заметьте, что существует хорошее соглашение.

[ampImb phImb; ampImbEst phImbEst]
ans = 2×2

    2.0000   15.0000
    2.0178   14.5740

Входные параметры

свернуть все

Коэффициент раньше компенсировал неустойчивость I/Q, заданную как комплексный вектор.

Пример: 0.4+0.6i

Пример: [0.1+0.2i; 0.3+0.5i]

Типы данных: single | double

Выходные аргументы

свернуть все

Амплитудная неустойчивость в дБ, возвращенном как вектор с действительным знаком с теми же размерностями как C.

Неустойчивость фазы в градусах, возвращенный как вектор с действительным знаком с теми же размерностями как C.

Больше о

свернуть все

Компенсация Неустойчивости I/Q

Функциональный iqcoef2imbal является функцией поддержки для Системы comm.IQImbalanceCompensator object™.

Учитывая фактор масштабирования и вращения, G, коэффициент компенсатора, C, и полученный сигнал, x, компенсированный сигнал, y, имеет форму

y=G[x+Cсоюз(x)].

В матричной форме это может быть переписано как

Y=RX,

где X 2 1 вектор, представляющий неустойчивый сигнал [XI, XQ] и Y 2 1 вектор, представляющий компенсатор вывод [YI, YQ].

Матрица R выражается как

R=[1+Ре{C}\Im{C}\Im{C}1Ре{C}]

Для компенсатора, чтобы отлично удалить неустойчивость I/Q, R = K-1, потому что X=KS, где K является матрицей 2 на 2, значения которой определяются амплитудой и неустойчивостью фазы, и S является идеальным сигналом. Задайте матрицу M с формой

M=[1αα1]

И M и M-1 могут считаться масштабированием и матрицами вращения, которые соответствуют факторному G. Поскольку K = R-1, продуктом M-1 R K M является единичная матрица, где M-1 R представляет компенсатор, вывод и K M представляют неустойчивость I/Q. Коэффициент α выбран таким образом что

KM=L[Igainпотому что(θI)Qgainпотому что(θQ)Igainsin(θI)Qgainsin(θQ)]

где L является константой. От этой формы мы можем получить Igain, Qgain, θI и θQ. Для данной неустойчивости фазы ΦImb, синфазные углы и квадратурные углы могут быть выражены как

θI=(π/2)(ΦImb/180)θQ=π/2+(π/2)(ΦImb/180)

Следовательно, because(θQ) = sin (θI) и sin (θQ) = because(θI) так, чтобы

L[Igainпотому что(θI)Qgainпотому что(θQ)Igainsin(θI)Qgainsin(θQ)]=L[Igainпотому что(θI)Qgainsin(θI)Igainsin(θI)Qgainпотому что(θI)]

Неустойчивость I/Q может быть выражена как

KM=[K11+αK12αK11+K12K21+αK22αK21+K22]=L[Igainпотому что(θI)Qgainsin(θI)Igainsin(θI)Qgainпотому что(θI)]

Поэтому

(K21+αK22)/(K11+αK12)=(αK11+K12)/(αK21+K22)=sin(θI)/потому что(θI)

Уравнение может быть записано как квадратное уравнение, чтобы решить для переменной α, которая является D1α2 + D2α + D3 = 0, где

D1=K11K12+K22K21D2=K122+K212K112K222D3=K11K12K21K22

Когда |C | ≤ 1, квадратное уравнение имеет следующее решение:

α=D2D24D1D32D1

В противном случае, когда |C |> 1, решение имеет следующую форму:

α=D2+D24D1D32D1

Наконец, амплитудная неустойчивость, AImb, и неустойчивость фазы, ΦImb, получена.

K=K[1αα1]AImb=20журнал10(K11/K22)ΦImb=2загар1(K21/K11)(180/π)

Примечание

  • Если C действителен и |C | ≤ 1, неустойчивость фазы 0, и амплитудная неустойчивость 20log10 ((1–C) / (1+C))

  • Если C действителен и |C |> 1, неустойчивость фазы составляет 180 °, и амплитудная неустойчивость 20log10 ((C +1) / (C −1)).

  • Если C является мнимым, AImb = 0.

Расширенные возможности

Генерация кода C/C++
Генерация кода C и C++ с помощью MATLAB® Coder™.

Смотрите также

|

Введенный в R2014b