yincwxa2006的个人博客分享 http://blog.sciencenet.cn/u/yincwxa2006

博文

[转载]【建模与辨识程序】二阶系统一次性完成最小二乘辨识程序

已有 390 次阅读 2019-11-4 14:13 |系统分类:科研笔记|文章来源:转载

二阶系统一次性完成最小二乘辨识程序
%FLch3LSeg1
u=[-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,1]; %
系统辨识的输入信号为一个周期的M序列
z=zeros(1,16); %
定义输出观测值的长度
for k=3:16
z(k)=1.5*z(k-1)-0.7*z(k-2)+u(k-1)+0.5*u(k-2); %
用理想输出值作为观测值
end
subplot(3,1,1) %
画三行一列图形窗口中的第一个图形
stem(u) %
画出输入信号u的经线图形
subplot(3,1,2) %
画三行一列图形窗口中的第二个图形
i=1:1:16; %
横坐标范围是116,步长为1
plot(i,z) %
图形的横坐标是采样时刻i, 纵坐标是输出观测值z, 图形格式为连续曲线
subplot(3,1,3) %
画三行一列图形窗口中的第三个图形
stem(z),grid on%
画出输出观测值z的经线图形,并显示坐标网格
u,z%
显示输入信号和输出观测信号
%L=14%
数据长度
HL=[-z(2) -z(1) u(2) u(1);-z(3) -z(2) u(3) u(2);-z(4) -z(3) u(4) u(3);-z(5) -z(4) u(5) u(4);-z(6) -z(5) u(6) u(5);-z(7) -z(6) u(7) u(6);-z(8) -z(7) u(8) u(7);-z(9) -z(8) u(9) u(8);-z(10) -z(9) u(10) u(9);-z(11) -z(10) u(11) u(10);-z(12) -z(11) u(12) u(11);-z(13) -z(12) u(13) u(12);-z(14) -z(13) u(14) u(13);-z(15) -z(14) u(15) u(14)] %
给样本矩阵HL赋值
ZL=[z(3);z(4);z(5);z(6);z(7);z(8);z(9);z(10);z(11);z(12);z(13);z(14);z(15);z(16)]%
给样本矩阵zL赋值
%calculating parameters%
计算参数
c1=HL'*HL; c2=inv(c1); c3=HL'*ZL; c=c2*c3 %
计算并显示
%DISPLAY PARAMETERS
a1=c(1), a2=c(2), b1=c(3), b2=c(4) %
中分离出并显示a1 a2 b1 b2
%End

 




http://blog.sciencenet.cn/blog-748007-1204778.html

上一篇:[转载]【建模与辨识程序】M序列产生程序
下一篇:[转载]【建模与辨识程序】实际压力系统的最小二乘辨识程序

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2020-4-6 10:31

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部