一个
!
"
!
!
!
,使得
#
(
!
"
)
" !
,
$
(
!
"
)
"
!
,这里
$
(
!
)是一
个矩阵多项式,
#
(
!
)是一个实多项式,由此(
#
)式的特
殊情况可定义为:
!
$ #
" %
,
%
$ #
" !
。
(
&
)初始化。从
!
!
任选
!
!
,若对所有
"
(
!
#
"
#
%
)
成立
&
(
!
)
"
" &
(
!
)
!
,则算法以
’
(
!
)
(
!
)
" &
(
!
)
!
结束;否则定
义
(
!
" &
(
!
)
!
,
!
#
"
!
!
) !
!
,
* " #
。
(
&&
)迭代。若对所有
"
(
*
#
"
#
%
),成立
&
(
* $ #
)
"
"
(
* $ #
或
&
(
* $ #
)
"
" %
,则该插值问题无解,从而算法中止;
否则,从
!
*
选择
!
*
使得
&
(
* $ #
)
*
"
(
* $ #
且
&
(
* $ #
)
*
"
%
。定
义
(
*
"
(
!
*
$ !
* $ #
)
)
(
&
(
* $ #
)
*
$ (
* $ #
)。若
* " %
,算法中
止。否则,对
" " * ’ #
,
* ’ (
,…,
%
,定义
&
(
*
)
"
"
(
!
"
$
!
* $ #
)
)
(
&
(
* $ #
)
"
$ (
* $ #
),注意
%
是
&
(
* $ #
)
"
的允许取值。
定义
!
* ’ #
"
!
*
) !
*
,
* ) " * ’ #
(迭代)。
(
&&&
)结束。对
* " #
,
(
,…,
% $ #
,构造
’
(
*
)
(
!
)
+
(
*
,
!
-
!
*
(
*
,
#
,
…
,
!
-
!
%
-
#
(
%
。
若对任意
*
(
#
#
*
#
% $ #
),成立
’
(
*
)
(
!
* $ #
)
" !
,则该插
值问题无解,否则由算法得(
(
)式。
!
例
设插值数据和运算过程如下表,由算法得
’
(
!
)
(
!
)
+
# !
[ ]
# !
,
!
! !
-
[ ]
# #
,
!
-
(
#
*
! +
-
#
-
[ ]
(
,
!
-
+
! +
-
#
-
[ ]
(
+
#
#,!
(
-
-,!
,
-,
.
#,!
(
-
-,!
,
-, +!
(
!
-
#
)(
!
-
(
)
-
(!
+
,
#.!
(
-
-*!
,
-, !
(
!
-
#
)(
+
-
!
[ ]
)
不难发现
’
(
!
)
(
!
"
)
" &
(
!
)
(
!
"
),
" " !
,
#
,
(
,
+
。
表
! ! # ( +
&
(
!
)
(
!
)
# !
[ ]
# !
# !
[ ]
# !
# !
[ ]
! #
# #
[ ]
! !
&
(
#
)
(
!
)
%
! !
[ ]
$ # #
+
(
! #
[ ]
$ # !
&
(
(
)
(
!
)
! !
#
*
! +
[ ]
$ # $ (
&
(
+
)
(
!
)
! +
[ ]
$ # $ (
%
参 考 文 献
# /0102 3
,
45678 9 :
,
;5&<5=> 3? 9 17@A12&B7 C7=>08 D01 20<B&EF AE@0E2=15&E78
=5EF7E=&5< &E=71G0<5=&0E G10H<7C2? IJJJ 315E2 9A=0C5= K0E=10<
,
#...
;
,,
:
,-(
—
,L.
( MA K>A5EN&EF? 3>&7<7O=6G7 5E8 P51F15EF7O=6G7 F7E715<&Q78 &EB7127 15=&0E5< &EO
=71G0<5=&0E D01 17@=5EFA<51 @0CG<7R C5=1&@72? P&E751 9<F7H15 5GG<
,
#...
;
(.-
:
*
—
+!
+ MA K>A5EN&EF? /&B51&5=7 3>&7<7O=6G7 C5=1&R B5<A78 15=&0E5< &E=71G0<5E=2? S
K0CGA= 9GG< T5=>
,
#..*
;
U!
:
*#
—
U(
(!!!
年
+
月
#L
日收到 国家自然科学基金(
#.U*#!-,
)资助
V(##? L, !!!L;!!(
构造
"#$"
新息模型的
%&’&()*+,-.&()
算法
邓自立 齐国元
(黑龙江大学自动化系,哈尔滨
#-!!U!
)
%&’&()*+,-.&() "/0,(1.23 4,( 5,6).(-7.160 "#$" 866,’9.1,6 $,:&/)
W7EF X&<&
,
Y& MA06A5E
(
W7G51=C7E= 0D 9A=0C5=&0E
,
:7&<0EFZ&5EF [E&B712&=6
,
:51H&E #-!!U!
)
摘 要
9\T9
新息模型被广泛应用于解决状态和信号估计问题。建议用
M7B712O]0A=712
算法构造线性离散随机系统的
9\T9
新息模型。该算法具有算法简单、快速收敛和保证新息模型的可逆性的优点。给出了用
T93P9/
工具包编的
M7BO
712O]0A=712
算法程序清单。两个仿真例子说明了其有效性。
关键词
9\T9
新息模型
M7B712O]0A=712
算法
T93P9/
软件
";).(97. 3>7 5A=017F1722&B7 C0B&EF 5B715F7
(
9\T9
)
&EE0B5=&0E C087<2 517 ^&87<6 5GG<&78 =0 20<B7 =>7 2=5=7 5E8 2&FE5< 72=&C5=&0E G10HO
<7C2? I= &2 2AFF72=78 =0 @0E2=1A@= =>7 9\T9 &EE0B5=&0E C087<2 0D <&E751 8&2@17= 2=0@>52=&@ 262=7C2 H6 A2&EF =>7 M7B712O ]0A=712 5<F01&=>C?
I= >52 58B5E=5F72 =>5= =>7 5<F01&=>C &2 2&CG<7
,
5E8 >52 D52= @0EB71F7E@7 G10G71=6
,
5E8 =>7 &EB71=&H&<&=6 0D =>7 &EE0B5=&0E C087<2 &2 7E2A178?
3>7 G10F15C 0D =>7 M7B712O]0A=712 5<F01&=>C H6 A2&EF =>7 T93P9/ =00<H0R &2 F&B7E? 3^0 2&CA<5=&0E 7R5CG<72 2>0^ &=2 7DD7@=&B7E722?
<&= >,(:) 9\T9 &EE0B5=&0E C087< M7B712O]0A=712 5<F01&=>C T93P9/ 20D=^517
?
问题阐述
*(*
(!!! _0<? L
,
‘0? L
中国学术期刊文摘(科技快报)