%-----CAR(2)-------- clear; clc; close all; N=4000; %步長 rand('seed',10); %改變種子改變數(shù)據(jù) e=0.6*randn(1,N+1); %白噪聲序列 a=-[-0.5000 0.0400]; %poly([0.1 0.4]) b=[3 1]; rand('seed',4); %改變種子改變數(shù)據(jù) u=2*rand(1,N+1)-1; %u(t)在[-1,1]之間服從均勻分布y(1)=e(1); %初值y(1),y(2),y(3) y(2)=a(1)*y(1)+b(1)*u(1)+e(2);for t=3:N+1y(t)=a(1)*y(t-1)+a(2)*y(t-2)+b(1)*u(t-1)+b(2)*u(t-2)+e(t);%AR(3)的RLS算法實現(xiàn) endn=2; lamda=0.998;%遺忘因子 sita(:,2)=zeros(2*n,1);%初值,均為0 P(:,:,2)=10^5*eye(2*n);%P的初值,選取足夠大的單位陣for t=2:Nfai(:,t+1)=[y(t) y(t-1) u(t) u(t-1)]';%最小二乘格式的fai(t)的表達式%參數(shù)估值sita(:,t+1)=sita(:,t)+P(:,:,t)*fai(:,t+1)/(lamda+fai(:,t+1)'*P(:,:,t)*fai(:,t+1))*...[y(t+1)-fai(:,t+1)'*sita(:,t)];P(:,:,t+1)=1/lamda*(P(:,:,t)-P(:,:,t)*fai(:,t+1)/(lamda+fai(:,t+1)'*P(:,:,t)*fai(:,t+1))*...fai(:,t+1)'*P(:,:,t)); endee(1)=y(1); taoe(1)=ee(1)^2;%y(0)的初值為0 for t=2:N%對噪聲方差的估計ee(t)=y(t)-fai(:,t)'*sita(:,t-1);taoe(t)=taoe(t-1)+1/t*[ee(t)^2-taoe(t-1)]; endt=1:N; subplot(3,2,1); plot(t,sita(1,t),'r'); line([0,N],[a(1),a(1)]) axis([0 N 0 1]); title('a(1)'); subplot(3,2,2); plot(t,sita(2,t),'r'); line([0,N],[a(2),a(2)]) axis([0 N -1 1]); title('a(2)'); subplot(3,2,3); plot(t,sita(3,t),'r'); line([0,N],[b(1),b(1)]) axis([0 N 0 5]); title('b(1)'); subplot(3,2,4); plot(t,sita(4,t),'r'); line([0,N],[b(2),b(2)]) axis([0 N 0 3]); title('a(2)'); subplot(3,2,5); plot(t,taoe,'r'); axis([0 N 0 1]); line([0,N],[0.36,0.36]) title('taoe'); subplot(3,2,6); plot(t,u(t),'r'); title('u(t)');

其中,e代表,為白噪聲序列

lamda代表,為遺忘因子

sita代表,估值參數(shù)之一

fai代表,為估值參數(shù)之一

ee代表

taoe代表

其中一次的測試結(jié)果如下圖所示:

總結(jié)

以上是生活随笔為你收集整理的CAR(2)模型的最小二参数估计的MATALB实现的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。