传递函数频域辨识及M序列生成指南精编.docxVIP

传递函数频域辨识及M序列生成指南 12自动化 摘要：系统辨识、状态估计和控制理论是现代控制论中相互渗透的三个领域。控制理论的应用离不开系统辨识技术，实际中，许多控制系统的模型在工作中是变化的，为了实现自适应控制，需要系统辨识技术不断更新模型参数。通过学习使用MATLAB软件，初步体验系统辨识方法。 关键字：系统辨识，控制理论，MATLAB。 Abstract: System identification, State estimation and The Principle of Automatic Control are three different disciplines of the modern control theory, which are interpenetrated?with?one?another. In practice, the model of system is changing all the time. To control adaptively, the system model should be update its parameters, by the method of System identification. By learning the using of MATLAB, we are supposed to practice the method of system identification. Key Words: System identification, System identification, MATLAB 目录 一、 引言 2 1.1介绍 2 1.2实验目的 3 二、 实验内容和方法 3 2.1实验内容 3 2.1.1实验一: 3 2.1.2实验二: 3 2.2实验步骤 3 2.2.1实验一： 3 非周期激励信号 5 2.2.2实验二：M序列 6 11 三、 实验分析分析 11 一、 引言 1.1介绍 在自然科学和社会科学的许多领域中，人们越来越重视对系统进行定量的系统分析、系统综合、仿真、控制和预测。将被控对象模型化，是开展工作的前提和基础。 基于传递函数的经典控制理论分析前需要确定模型的各个参数，忽略工作中的变化，适用于精度要求不高的情况下。而实际中许多控制系统内部参数是未知的，而且工作中有可能会变化，这就需要系统辨识来确定模型参数。通常为机理分析法和测试法相结合来辨识，对于系统机理已知的部分采用机理分析法，机理未知的部分采用测试法。 1.2实验目的 本次实验通过使用 MATLAB 仿真系统辨识的过程，采用频域特性拟合的Levy方法，按要求完成传递函数的辨识。描述实验验证的数据准备、基本过程和实验结果。 结合课堂学习的M序列的特性和生成方法，使用matlab语言产生一个M序列并分析器自相关特性。 二、 实验内容和方法 2.1实验内容 2.1.1实验一: 1． 自己设定一个稳定系统，采用周期测试信号，测定系统的频率响应。 2． 对题1中的系统，采用非周期测试信号确定系统的频率响应，并与题1的结果对比。 3． 基于题1或题2 产生的频率响应数据，采用课堂讲授的频域特性拟合方法，辨识传递函数的参数。将辨识结果与Matlab工具库中的等价的功能函数invfreqs产生的结果做对比。 2.1.2实验二: 根据最大长度现行反馈寄存器M序列生成机制，编写M序列生成的生成程序。 1． 自己设定移位寄存器的级数和初值，产生响应的M序列。 2． 绘制题1产生的M序列的自相关函数和功率谱密度图形。 2.2实验步骤 2.2.1实验一： 周期性信号： 周期性信号最典型的就是正弦周期函数，在不同频率下会产生不同的周期，进而有不同的频率响应。将不同频率的正弦波组合起来作为输入信号。 U=Asin(w*t+φ) Y=Asin(w*t+ρ) 所以： Am=A’/A Phi= ρ-φ 设计系统函数： num=[4 3]; den=[4 5 6]; sys=tf(num,den); Transfer function: 4 s + 3 --------------- 4 s^2 + 5 s + 6 根据这个原理，设计离散的正弦周期信号。 1、一个周期内取1000个采样点，采样是个周期。单位采样间隔二0.001*T 2、由于 T=2*pi./w，所以每一次频率改变都黑改变周期，进而有不同的激励和响应图像。 3、在十个周期内，对第五和第六个采样周期进行分析，得到周期内的最大值和最大值的坐标。 4、对比激励和响应的最大值的峰值比和相位差，进而得到不同频率的一系列频谱响应。 Matlab代码如下： w=logspace(-1,2,1000);