【正文】
x(1,n1+n2+1:end),x(2,n1+n2+1:end),39。r+39。,t,[pt,ft])。最小平方誤差準(zhǔn)則: 設(shè)其中,目標(biāo): min 目標(biāo)函數(shù)的梯度:令 %MSE b=ones(n1+n2,1)。Y=x1(:,1:n1+n2)。Y=Y39。Yplus=(Y39。*Y)^(1)*Y39。ahat=Yplus*b。mt=ahat(2)*t/ahat(3)ahat(1)/ahat(3)。figure(4)plot(t,[pt。ft。mt],x(1,1:n1),x(2,1:n1),39。o39。,x(1,n1+1:n1+n2),x(2,n1+1:n1+n2),39。*39。,x(1,n1+n2+1:end),x(2,n1+n2+1:end),39。r+39。)。axis([,])。display(39。Using MSE principal is:39。)Apf=find(ahat39。*x1(:,n1+n2+1:end)0)+n1+n2,