-
大小: 2KB文件類型: .m金幣: 1下載: 1 次發(fā)布日期: 2021-06-15
- 語言: Matlab
- 標(biāo)簽:
資源簡介
開環(huán)系統(tǒng)參數(shù)辨識,帶遺忘因子的遞推最小二乘估計法(FFRLS),系統(tǒng)為單入單出的CAR(帶控制量的自回歸模型)模型,三階系統(tǒng)
代碼片段和文件信息
%開環(huán)系統(tǒng)參數(shù)辨識,帶遺忘因子的遞推最小二乘估計法(FFRLS),系統(tǒng)為單入單出的CAR(帶控制量的自回歸模型)模型三階系統(tǒng)
clear?all
clc
a=[1?-1.1?0.6?0.1];b=[1?0.7];d=4;?%實(shí)際模型系數(shù)矩陣與純遲延
L=1000;???????????????????????????%仿真長度
na=length(a)-1;nb=length(b)-1;????%nanb為輸出輸入系數(shù)矩陣AB的階數(shù)
yk=zeros(na1);???????????????????%輸出矩陣初始化
yk_m=zeros(na1);?????????????????%模型輸出
uk=zeros(nb+d1);?????????????????%輸入矩陣初始化
theta_e0=zeros(na+nb+11);????????%theta_e0為估計參數(shù)初值,a1a2....anb0b1...bn共na+nb+1個
phi=zeros(na+nb+11);?????????????%phi為當(dāng)前實(shí)際輸出輸入構(gòu)成的矩陣
P=10^6*eye(na+nb+1);??????????????%修正系數(shù)初值
beta=0.99;????????????????????????%遺忘因子,在0.95到1之間
u=randn(L1);?????????????????????%輸入信號,方差為1的白噪聲序列
omega=sqrt(0.1)*randn(L1);???????%干擾信號,方差為0.1的白噪聲序列
for?i=1:L
????theta(:i)=[a(2:na+1)b]‘;????%系統(tǒng)實(shí)際參數(shù)值
????
?
評論
共有 條評論