Модель фильтра Калмана 2 порядка - История изменений

Ippolitov в 13:05, 22 июля 2016

2016-07-22T13:05:43Z

2013-04-25T11:08:42Z

2013-04-25T11:07:57Z

2012-05-04T08:00:13Z

‎С моделью фазы

2012-05-04T07:59:54Z

‎С моделью фазы

2012-05-04T07:42:15Z

‎С моделью фазы

2012-01-20T13:02:58Z

‎С моделью фазы

2012-01-20T12:53:53Z

‎С моделью фазы

2011-11-07T06:43:14Z

‎С моделью фазы

2011-10-12T12:34:49Z

‎С моделью фазы

@@ Строка 61: / Строка 61: @@
 for k = 1:K
      Ud = f(Xextr, Xist);  % Дискриминатор
-     Sd = f(A_IQ);      % Критизна дискриминационной характеристики
+     Sd = f(A_IQ);      % Крутизна дискриминационной характеристики
      Xest = Xextr + Ko*Ud/Sd; % Вектор оценок на k-й интервал
      Xextr = F*Xest;         % Экстраполяция на интервал k+1

@@ Строка 14: / Строка 14: @@
 |repository    =
 |category1     = Статистическая радиотехника
-|category2     = Фазовые измерения
+|category2     = Оценивание частоты
-|category3     = Переходные процессы
+|category3     = Оценивание задержки
 |category4     =
 |category5     =
 |category6     =
 |category7     =
@@ Строка 171: / Строка 171: @@
 == См. также ==

← Предыдущая		Версия 11:07, 25 апреля 2013
Строка 171:		Строка 171:

	== См. также ==		== См. также ==
		+
		+
		+	[[Категория:Оценивание частоты]]
		+	[[Категория:Оценивание задержки]]

@@ Строка 147: / Строка 147: @@
          %     Ud = f(I(k), Q(k), I(k-1), Q(k-1), ...);      % Дискриминатор
          Ud = (I(k)*Qold - Q(k)*Iold);
-         %     Sd = f(A_IQ);             % Критизна дискриминационной характеристики
+         %     Sd = f(A_IQ);             % Крутизна дискриминационной характеристики
          Sd = Tc*(A_IQ(k)*Tf/Tc)^2 * 1.3;
          Xest = Xextr + Ko*Ud/Sd;                 % Вектор оценок на очередной интервал фильтра

@@ Строка 148: / Строка 148: @@
          Ud = (I(k)*Qold - Q(k)*Iold);
          %     Sd = f(A_IQ);             % Критизна дискриминационной характеристики
-         Sd = Tf*(A_IQ(k)*Tf/Tc)^2;
+         Sd = Tc*(A_IQ(k)*Tf/Tc)^2 * 1.3;
          Xest = Xextr + Ko*Ud/Sd;                 % Вектор оценок на очередной интервал фильтра
          Xextr = Ff*Xest;                         % Экстраполяция на следующий интервал

@@ Строка 122: / Строка 122: @@
 nQ = stdn_IQ.*randn(1,K); % // Q-comp noise
-w = 0; Isum = 0; Qsum = 0;
+w = 0; Isum = 0; Qsum = 0; Iold = 1; Qold = 0;
 for k = 1:K
@@ Строка 129: / Строка 129: @@
      EpsW(k) = Xist(2) - Xextr(2);
      EpsTau(k) = 0;
-     [A_IQ(k) qcno] = qcno_change(qcno_ist(k), stdn_IQ(k), Tc);
      A_IQ_eff(k) = A_IQ(k)*sinc(EpsW(k)*Tc/2 /pi)*ro(EpsTau(k));
      mI = A_IQ_eff(k) * cos(EpsW(k)*Tc/2 + EpsPhi(k));

@@ Строка 122: / Строка 122: @@
 nQ = stdn_IQ.*randn(1,K); % // Q-comp noise
-w = 0;
+w = 0; Isum = 0; Qsum = 0;
 for k = 1:K
@@ Строка 135: / Строка 135: @@
      I(k) = mI + nI(k);
      Q(k) = mQ + nQ(k);
+     Isum = Isum + I(k);
      Xextr = Fincorr * Xextr; % Набег фазы в корреляторе к концу накопления
@@ Строка 147: / Строка 149: @@
          Xextr = Ff*Xest;                         % Экстраполяция на следующий интервал
          w = 0;
      end
      Xist = Fc*Xist + [0; 0; 1]*nIst(k)*stdIst; % Здесь может быть любая другая модель изменения истинного вектора состояния