线性网络在线性系统辩识中的应用。
1 问题描述
对于有限冲击响应线性系统,有输入信号 x,其周期为 5 秒,每 25 毫秒采
样一次,信号 x 为
time=0:0.025:5;
x=sin(sin(time).*time*10);
假设系统的输出测量值为
t=filter([1 0.5 –1.5],1,x);
系统的输入、输出值如图 14.9 所示.
(a) 信号 x (b)系统输出
图 14.9 系统的输入、输出值
1 网络设计
此例中,应用一个线性网络来进行系统辩识.该网络设计为三个输入和一个输
出,其输入由信号 x 的当前值和前两个时刻的值,可由函数 delaysig 生成
p=delaysig(x,0,2)
然后,可利用函数 solvelin 来设计一个线性网络
[w,b]=solvelin(p,w,b)
2 网络性能检验
网络的权值和阈值确定后,可进一步利用函数 simulin 对网络性能进行检验
并比较网络输出与系统实际输出结果。
以下是相应的 MATLAB 源程序
% Linear system identification.
clc
pause % Strike any key to continue...
clc
time = 0:0.025:5; % from 0 to 6 seconds
X = sin(sin(time).*time*10);
0 1 2 3 4 5
-1
-0.5
0
0.5
1
Time
Input Signal
Input S ignal to the System
0 1 2 3 4 5
-2
-1.5
-1
-0.5
0
0.5
1
1.5
2
Time
Output Signal
Output Signal of the System