原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
[系统辨识M序列
M序列的产生
因为Np=15，所以应该用4个移位寄存器。
连接图如下:
仿真图如下：
表如下：
用程序实现如下：
X1=0;X2=0;X3=0;X4=0; %定义移位寄存器Xi初态（1000）， Yi为移位寄存器各级输出
m=32; %置M序列总长度
for i=1:m %1#
Y4=X4; Y3=X3; Y2=X2; Y1=X1;
X4=Y3; X3=Y2; X2=Y1;
X1=xor(Y3,Y4); %异或运算
if Y4==0
U(i)=0.03;
else
U(i)=Y4*(-0.03);
end
end
M=U;
figure(1) %绘M序列图
i1=i
k=1:1:i1;
stairs(M,g-)
xlabel(k)
ylabel(M序列)
title(移位寄存器产生的M序列)
2．从图形得到数据
Y=[2.05 1.80 1.84 1.79 1.08 0.68 0.44 0.80 1.91 2.38 2.47 2.51 3.05 2.69 1.94 1.82 1.82 2.03 2.03 1.03 0.68 0.52 0.86 1.78 2.50 2.50 2.32 3.28 2.82 2.04 2.01 1.67 1.70 1.82 1.04 0.59 0.38 0.81 1.91 2.55 2.28 2.56 3.13 2.70 2.01];
3．相关函数
s=0;c=0; %设置中间变量,初始化
R=zeros(1,15);r=zeros(1,32);%初始化相关数组
for i=1:15 %计算相关函数
for j=16:32
r(j)=U(j-i+1)*Y(j)/30;
s=s+r(j)
end
R(i)=s;
s=0;
end
c=(-1)*R(15); %求补偿量c
for i=1:15 %求脉冲响应g（k）
g(i)=15*(R(i)+c)/(4*0.03^2*16)
end
figure(3)
plot(R)
xlim([1,15])
ylabel(Rym)
title(互相关函数);
4.脉冲响应曲线
hold on
figure(4)
plot(g)
xlim([1,15])
ylabel(g)
title(脉冲响应);
5.阶跃响应曲线
hold on
p=0;
for k=1:15 %用积分法求阶跃响应
for l=1:k
p=g(l)+p
end
v(k)=p;
p=0;
end
figure(5)
plot(v)
title(阶跃响应)
xlim([1,15]);
hold on
6.得出加热炉的数学模型
选择采样间隔T=1，于是得g(t)在t=[1：1：15]这15个时刻对应的值为：
0.2813, 0.3177, 0.6094, 3.896，9.021, 12.76，15.97，18.14，19.82, 21.21，21.79, 21.74
21.86, 22.1, 22.1
用Z变换的方法来求解系统的传递函数，设
(1)
由实验曲线阶跃响应曲线可得其Z变换G（z）（采样间隔为T=1）为
0.2813/z +0.3177 /z^2+0.6094/z^3 (2)
要使（1）式与（2）式相等，则应取 b1=0.3201 a1= -1.2158 a2= -1.2369
则可以得到：
文档评论（0）