下载文档原格式
利用Matlab进行系统辨识的基本步骤系统辨识是在工程领域中广泛应用的一项技术，它可以通过实验数据来推断系统的特性和行为。
Matlab作为一种强大的数值计算与数据可视化软件，在系统辨识领域也扮演着重要的角色。
本文将介绍利用Matlab进行系统辨识的基本步骤，以帮助读者更好地了解和应用这一技术。
系统辨识的基本步骤可以分为数据采集、预处理、模型选择、参数估计和模型评估五个阶段。
下面将依次介绍每个阶段的具体步骤。
一、数据采集数据采集是系统辨识中至关重要的一步，它需要准备好合适的实验设备，并进行数据记录。
在Matlab中，可以利用数据采集卡或传感器等设备进行数据采集。
一般来说，需要采集输入信号和输出响应信号。
二、预处理预处理是为了提高辨识结果的准确性而对采集到的数据进行处理。
首先要检查数据是否存在异常值或噪声，如果有则进行滤波处理。
另外，还可以进行数据归一化、降采样和数值平滑等处理，以适应不同的辨识方法。
三、模型选择模型选择是在系统辨识过程中非常关键的一步。
在Matlab中，可以选择多种模型结构进行辨识，包括线性模型、非线性模型、时变模型和非参数模型等。
选择合适的模型结构要根据具体问题的要求和实际情况进行综合考量。
四、参数估计参数估计是确定系统辨识结果的过程，它可以通过最小二乘法或极大似然法等统计方法来估计模型的参数。
在Matlab中，可以使用现成的函数或自定义算法进行参数估计。
参数估计的结果可以用来表示系统的传递函数、状态空间模型或差分方程等。
五、模型评估模型评估是对辨识结果进行验证和优化的过程。
在Matlab中，可以利用辨识误差、残差分析和模型预测误差等指标来评估模型的质量。
如果评估结果不理想，还可以尝试调整模型结构或参数估计方法，进行迭代优化。
在实际应用中，系统辨识常常涉及到复杂的数学理论和算法。
Matlab作为一种集成了各类工具箱和函数的软件，可以大大简化系统辨识的步骤，提高辨识的效率和准确性。
利用Matlab进行系统辨识的技术方法一、引言系统辨识是研究系统动态特性的一个重要方法，它广泛应用于控制系统、信号处理、通信等领域。
利用Matlab进行系统辨识能够实现快速、准确的模型建立和参数估计。
本文将介绍在Matlab环境下常用的系统辨识技术方法及其应用。
二、系统辨识的基本概念系统辨识是通过对系统的输入和输出信号进行观测和分析，以推断系统的结构和参数。
一般来说，系统辨识包括建立数学模型、估计系统参数和进行模型验证三个步骤。
1. 建立数学模型建立数学模型是系统辨识的第一步，它是描述系统行为的数学表达式。
常用的数学模型包括线性模型、非线性模型和时变模型等。
2. 估计系统参数在建立了数学模型之后，需要通过对实验数据的分析，估计出系统的参数。
参数估计可以通过最小二乘法、极大似然估计法等方法实现。
3. 模型验证模型验证是为了确定估计得到的系统模型是否准确。
常用的方法有经验验证、残差分析、模型检验等。
三、常用的系统辨识技术方法1. 线性参数模型线性参数模型是最常用的系统辨识方法之一。
它假设系统具有线性特性，并通过估计线性模型的参数来描述系统。
在Matlab中，可以使用函数"arx"进行线性参数模型的辨识。
2. 神经网络模型神经网络模型是一种非线性模型，它通过人工神经元的连接权值来描述系统行为。
在Matlab中，可以使用"nlarx"函数进行神经网络模型的辨识。
3. 系统辨识工具箱Matlab提供了丰富的系统辨识工具箱，包括System Identification Toolbox和Neural Network Toolbox等。
这些工具箱提供了各种方法和函数，方便用户进行系统辨识分析。
四、利用Matlab进行系统辨识的应用案例1. 系统辨识在控制系统中的应用系统辨识在控制系统中具有广泛的应用，如无人机控制、机器人控制等。
通过对系统进行辨识，可以建立准确的数学模型，并用于控制器设计和系统优化。
利用Matlab进行系统辨识的技术方法在Matlab中进行系统辨识的技术方法主要有参数估计法和非参数估计法两种。
1.参数估计法：参数估计法是通过拟合已知输入和输出数据的数学模型来估计系统的参数。
常用的参数估计方法包括最小二乘法（OLS）、最小二乘法（LSE）、最小二乘法（MLE）和极大似然估计法（MLE）等。
a) 最小二乘法（OLS）：OLS方法通过最小化实际输出与模型预测输出之间的误差平方和来估计系统参数。
在Matlab中，可以使用lsqcurvefit函数来实现最小二乘法的系统辨识。
b) 最小二乘法（LSE）：LSE方法是通过最小化实际输出与模型预测输出之间的误差平方和来估计系统参数。
在Matlab中，可以使用lsqnonlin函数来实现最小二乘法的系统辨识。
c) 最小二乘法（MLE）：MLE方法是通过最大化似然函数来估计系统参数。
在Matlab中，可以使用mle函数来实现最大似然估计法的系统辨识。
2.非参数估计法：非参数估计法不需要事先指定系统的数学模型，而是直接根据输入和输出数据的统计特性进行系统辨识。
常用的非参数估计方法包括频域方法、时域方法和时频域方法等。
a) 频域方法：频域方法通过对输入和输出数据进行频谱分析来估计系统的频率响应。
常用的频域方法包括傅里叶变换、功率谱密度估计和频率响应函数估计等。
在Matlab中，可以使用fft函数和pwelch函数来实现频域方法的系统辨识。
b) 时域方法：时域方法通过对输入和输出数据进行时间序列分析来估计系统的时域特性。
常用的时域方法包括自相关函数估计和互相关函数估计等。
在Matlab中，可以使用xcorr函数来实现时域方法的系统辨识。
c) 时频域方法：时频域方法结合了频域方法和时域方法的优势，可以同时估计系统的频率响应和时域特性。
常用的时频域方法包括短时傅里叶变换和小波变换等。
在Matlab中，可以使用spectrogram函数和cwt函数来实现时频域方法的系统辨识。
Matlab系统辨识尝试之详细过程1前面介绍了Matlab系统辨识工具箱的一些用法，这里拿一个直观的例子来尝试工具箱的具体用法。
比较长，给个简单目录吧：1.辨识的准备2.辨识数据结构的构造3.GUI辨识4.辨识效果5.对固有频率的辨识6.结构化辨识7.灰箱辨识8.加入kalman滤波的灰箱辨识1.辨识的准备在辨识前，首先要根据自己辨识的情况，确定要辨识的状态空间模型的一些特点，如连续还是离散的；有无直通分量（即从输入直通到输出的分量）；输入延迟；初始状态等。
了解了这些情况就可以更快速的配置辨识时的一些设置选项。
2.辨识数据结构的构造使用原始数据构造iddata结构：data=iddata(y,u,Ts);这里以一个弹簧质量系统的仿真为例代码如下，其中用到了函数MDOFSolve，这在之前的博文介绍过（/?p=183），拿来用即可。
如果发现运行有错误，可以将MDOFSolve函数开头的一句omega2=real(eval(omega2));注释掉。
%弹簧质量系统建模clcclearclose allm=200;k=980*1000;c=1.5*1000;m1=1*m;m2=1.5*m;k1=1*k;k2=2*k;k3=k1;%%由振动力学知识求固有频率M=[m10;0m2];K=[k1+k2-k2;-k2k3+k2];[omega,phi,phin]=MDOFSolve(M,K);fprintf('固有频率:%fHz\n',subs(omega/2/pi));%%转化到状态空间innum=2;outnum=2;statenum=4;A=[0100;-(k1+k2)/m10k2/m10;0001;k2/m20-(k3+k2)/m20];B=[00;1/m10;00;01/m2];C=[1000;0010];D=zeros(outnum,innum);K=zeros(statenum,innum);mcon=idss(A,B,C,D,K,'Ts',0);%连续时间模型figureimpulse(mcon)%%信号仿真，构造数据供辨识n=511;%输入信号长度Ts=0.001;t=0:Ts:(n-1)*Ts;u1=idinput(n,'prbs');%输入1为伪随机信号u2=zeros(n,1);%输入2为空u=[u1u2];simdat=iddata([],u,Ts);%形成输入数据对象e=randn(n,2)*1e-7;simopt=simOptions('AddNoise',true,'NoiseData',e);%添加噪声yn=sim(mcon,simdat,simopt);%加噪声仿真y=sim(mcon,simdat);%无噪声仿真figurefor i=1:outnumsubplot(outnum,1,i)plot(t,y.OutputData(:,i))hold onplot(t,yn.OutputData(:,i),'r')axis tighttitle(sprintf('输出%d',i))legend({'无噪声仿真','含噪声仿真'})end%保存输入输出数据，供后续辨识data=iddata(y.OutputData,simdat.InputData,Ts);datan=iddata(yn.OutputData,simdat.InputData,Ts);运行后，变量data中保存了无噪声的系统仿真输入输出数据，datan中为含噪声的仿真数据。
Matlab最⼩⼆乘系统辨识原理还是⽐较简单的，不赘述，程序⾥⾯的注释也写的⽐较清楚了%仿真对象:y(k)-1.5y(k-1)+0.7y(k-2)=v(k)+u(k)-0.8u(k-1)%辨识模型:y(k)+a1 y(k-1)+a2 y(k-2)=v(k)+b1 u(k)+b2 u(k-1)%数据长度取n=20000,加权矩阵为I，v(k)是服从正态分布的⽩噪声N（0,1）,u(k)=sin(k)%待估计参数K=[a1 a2 a3 a4]'；准则函数J(K)=(Yn-HnK)'(Yn-HnK);%将辨识模型写为：y(k)=v(k)+a1 y(k-1)+a2 y(k-2)+b1 u(k)+b2 u(k-1)% =v(k)+KHn%Hn=|y(2) y(1)|% |y(3) y(2)|% |.........|clearclose alldata_length=20002;%% 产⽣⽩噪声和输⼊v=randn(1,data_length);v=v./max(v);u=zeros(1,data_length);for k=1:data_lengthu(k)=sin(k);end%% 获得观测值y=zeros(1,data_length);for k=3:data_lengthy(k)=1.5*y(k-1)-0.7*y(k-2)+v(k)+u(k)-0.8*u(k-1);end%% 构造Hn和Y矩阵Hn=zeros(data_length-2,2);count=1;for k=1:10000Hn(k,2)=y(count);count=count+1;Hn(k,1)=y(count);Hn(k,4)=u(count);Hn(k,3)=u(count+1);end%% 求解参数Y=y(3:data_length)';c1=Hn'*Hn;c2=inv(c1);c3=Hn'*Y;K=c2*c3%% 将辨识得到的参数代⼊，得估计输出y_e=zeros(1,data_length);for k=3:data_lengthy_e(k)=K(1)*y_e(k-1)+K(2)*y_e(k-2)+v(k)+K(3)*u(k)+K(4)*u(k-1);end%% 画出实际输出和辨识输出，进⾏对⽐plot((1:data_length),y');title('实际输出')hold onplot((1:data_length),y_e');title('辨识输出')figuresubplot(2,1,1)plot((1:data_length),y');title('实际输出')subplot(2,1,2)plot((1:data_length),y_e');title('辨识输出')。
使用Matlab技术进行系统辨识的基本方法概述：系统辨识是指通过对已知输入输出数据的分析和处理，推断出系统的动态性质和数学模型的过程。
在科学研究、工程设计和控制应用中，系统辨识扮演着重要的角色。
而Matlab作为一种强大的数值计算和数据分析软件，为系统辨识提供了便利且高效的工具。
本文将介绍使用Matlab进行系统辨识的基本方法，并结合实例进行讲解。
一、数据采集与准备在进行系统辨识之前，首先需要采集到对应的输入输出数据。
一般来说，输入信号是已知的，可以通过外部激励或者系统自身的变动来获取；而输出信号则是根据输入信号通过系统响应得到的。
在采集数据时，需要注意数据的质量和采样频率的选择。
二、数据预处理在进行系统辨识之前，数据通常需要进行一些预处理，以去除噪声、平滑数据和调整时间序列等。
这可以通过Matlab中的数据处理函数和滤波器实现。
例如，可以使用高斯滤波器对数据进行平滑处理，或者使用降噪算法去除不必要的噪声。
三、参数估计参数估计是系统辨识的核心步骤之一，它通过对已知数据进行分析和处理，推断出系统的数学模型和参数。
在Matlab中，有多种方法和工具可供选择，如最小二乘法、最大似然法、系统辨识工具箱等。
这些工具可以根据不同的模型和数据类型灵活选择，并提供相应的算法和函数。
四、模型验证与优化根据估计得到的系统模型和参数，可以使用Matlab进行模型验证和优化。
模型验证是指将估计得到的模型与真实系统进行对比，检验其拟合程度和预测能力。
如果模型的拟合程度较差，则需要对参数进行调整和优化，以提高模型的准确性和稳定性。
五、模型预测与应用在系统辨识完成之后，可以使用得到的模型进行系统预测和应用。
通过对未知输入信号进行预测，可以得到相应的输出响应，进而实现对系统动态性质的分析和控制。
Matlab提供了丰富的预测和应用函数，例如时域模拟、频域分析、控制系统设计等，可以满足不同应用场景的需求。
六、案例分析为了更好地理解和掌握使用Matlab进行系统辨识的基本方法，下面通过一个简单的案例进行分析。
Matlab中的系统辨识和参数标识技巧引言：在工程和科学领域中，系统辨识是一项重要的任务。
通过对系统进行辨识和参数标定，我们可以建立数学模型来描述系统的行为。
在Matlab中，有许多工具和技巧可以用来进行系统辨识和参数标定。
本文将介绍一些常用的工具和技巧，以帮助读者更好地理解和应用这些方法。
一、准备工作在进行系统辨识和参数标定之前，我们首先需要准备一些基本的工作。
首先，我们需要收集系统的输入和输出数据。
输入数据通常是对系统施加的激励信号，例如阶跃信号或随机信号。
输出数据是系统对输入信号的响应。
收集足够量的数据对于准确辨识系统非常重要。
其次，我们需要确保数据的质量。
在收集数据时，我们需要注意采样频率和信噪比。
采样频率要足够高，以捕捉系统的快速动态响应；信噪比要足够高，以避免噪声对数据的影响。
如果数据质量不高，将会使得系统辨识和参数标定的结果不准确。
二、线性系统辨识线性系统辨识是系统辨识领域中的一项基本任务。
在Matlab中，可以使用System Identification Toolbox来进行线性系统辨识。
首先，我们可以使用中自相关函数（cross-correlation）来计算输入和输出数据之间的相关性。
这可以通过Matlab中的xcorr函数实现。
相关性分析可以帮助我们了解系统的输入和输出之间的关系，为后续的系统辨识提供参考。
接下来，我们可以使用频域分析方法来对系统进行辨识。
其中一个常用的方法是基于频率响应函数的辨识方法，例如最小二乘法（Least Squares）等。
这些方法可以通过Matlab中的tfest函数来实现。
另外，我们还可以使用时域辨识方法，例如最小均方误差法（Least Mean Squares）等。
时域辨识方法通常通过求解线性方程组来确定系统的参数。
在Matlab中，我们可以使用lsim和inv函数来实现这些方法。
三、非线性系统辨识与线性系统不同，非线性系统的辨识更加复杂。
MATLAB系统辨识工具箱学习详细教程MATLAB系统辨识工具箱是MATLAB软件中的一个工具箱，用于进行系统辨识和模型建模的分析。
该工具箱提供了多种辨识算法和工具，可以对线性和非线性系统进行辨识，并生成对应的数学模型。
下面将为您详细介绍MATLAB系统辨识工具箱的学习过程。
首先，在使用MATLAB系统辨识工具箱前，需要安装MATLAB软件并具备一定的MATLAB编程基础。
如果您还没有安装MATLAB或者对MATLAB不够熟悉，建议您先进行相关的学习和了解。
1.学习基本概念：在开始学习MATLAB系统辨识工具箱之前，需要了解一些基本概念，例如系统辨识、模型建模、参数估计等。
可以通过阅读相关的系统辨识的教材或者进行在线，对相关概念有一个基本的了解。
2.熟悉MATLAB系统辨识工具箱界面：3.数据导入：在进行系统辨识之前，首先需要准备好系统辨识所需的数据。
数据可以是实验数据或者仿真数据，可以是时域数据或者频域数据。
在系统辨识工具箱界面的“数据导入”区域，可以将数据导入到MATLAB中进行后续的辨识分析。
4.选择模型类型：在进行系统辨识之前，需要选择适合的数学模型类型。
MATLAB系统辨识工具箱提供了多种常见的模型类型，包括ARX模型、ARMAX模型、OE模型、TFE模型等。
选择合适的模型类型对辨识结果的精度和准确性有重要的影响。
5.选择辨识算法：在选择模型类型后，需要选择合适的辨识算法进行参数估计和模型建模。
MATLAB系统辨识工具箱提供了多种常用的辨识算法，例如最小二乘法、极大似然法、递推最小二乘法等。
选择合适的辨识算法也对辨识结果的精度和准确性有重要的影响。
6.进行系统辨识：在选择了合适的模型类型和辨识算法后，可以在系统辨识工具箱界面中点击“辨识”按钮，开始进行系统辨识分析。
系统辨识工具箱会根据所选的模型类型和辨识算法，对输入的数据进行参数估计和模型建模，并生成相应的辨识结果。
7.结果分析和评估：在系统辨识完成后，可以在系统辨识工具箱界面中查看辨识结果和模型质量评估。
如何使用MATLAB进行系统辨识与模型建模引言：近年来，随着科学技术的飞速发展，各行各业都在努力寻求更高效、更智能的解决方案。
系统辨识与模型建模作为一种重要方法和工具，被广泛应用于控制系统、信号处理、机器学习等领域。
在这些领域中，MATLAB作为一款功能强大的数值计算软件，为我们提供了丰富的工具和函数，可用于进行系统辨识与模型建模的分析和实现。
本文将详细介绍如何使用MATLAB进行系统辨识与模型建模，并探讨其在实际应用中的意义和局限性。
一、系统辨识的基本原理1.1 系统辨识的概念及意义系统辨识是指通过对已有数据的分析和处理，建立描述该系统行为的数学模型的过程。
在实际应用中，系统辨识可以帮助我们了解系统的结构和特性，预测系统的行为，并为系统控制、优化提供依据。
1.2 系统辨识的方法系统辨识的方法主要包括参数辨识和结构辨识两种。
参数辨识是指通过拟合已知数据，确定数学模型中的参数值的过程。
常用的参数辨识方法有最小二乘法、极大似然估计法等。
结构辨识是指通过选择适当的模型结构和参数化形式，使用已知数据确定模型结构的过程。
常用的结构辨识方法有ARX模型、ARMA模型等。
二、MATLAB在系统辨识中的应用2.1 数据准备与预处理在进行系统辨识之前，我们首先需要准备好相关的数据。
数据的质量和数量对系统辨识的结果有着重要的影响，因此在数据准备阶段应尽量确保数据的准确性和完整性。
MATLAB提供了丰富的数据处理和分析函数，可用于数据预处理、数据清洗、数据归一化等操作，以提高数据的质量和可用性。
2.2 参数辨识的实现参数辨识是系统辨识的重要步骤之一，其主要目标是通过适当的数学模型拟合已知数据，确定模型中的参数值。
在MATLAB中，我们可以使用curve fitting工具箱中的函数，如fit、cftool等，来进行参数辨识的实现。
同时，MATLAB还提供了最小二乘法等常用的参数辨识算法，方便我们根据实际需求进行选择和应用。
使用MATLAB进行系统辨识与参数估计的基本原理近年来，随着人工智能和机器学习的发展，系统辨识和参数估计变得越来越重要。
在工程和科学领域，系统辨识与参数估计可以帮助我们理解和预测复杂系统的行为，从而为决策和控制提供有力支持。
而MATLAB作为一种强大的科学计算软件，在系统辨识与参数估计方面提供了丰富的工具和功能。
本文将介绍MATLAB 中进行系统辨识与参数估计的基本原理。
一、系统辨识的概念系统辨识是指通过一系列的实验和数据分析，确定出系统的数学模型或特性。
在实际工程和科学问题中，我们经常遇到许多系统，如电子电路、生化反应、飞行控制系统等。
通过系统辨识，我们可以了解系统的行为规律，预测未来状态，从而进行优化和控制。
在MATLAB中，可以使用系统辨识工具箱（System Identification Toolbox）进行系统辨识。
该工具箱提供了一系列的函数和算法，可以帮助我们建立和分析系统模型。
例如，使用arx函数可以基于自回归模型建立离散时间系统的模型，使用tfest函数可以进行连续时间系统的模型辨识。
二、参数估计的基本原理参数估计是系统辨识的一个重要部分，它是指通过已知的输入输出数据，估计系统模型中的参数。
在实际应用中，我们通常只能通过实验数据来获得系统的输入输出信息，而无法直接观测到系统内部的参数。
因此，参数估计成为了一种重要的技术，用于从数据中推断出系统的模型参数。
在MATLAB中，参数估计的基本原理是最小二乘估计。
最小二乘估计是指寻找能够最小化实际输出与模型输出之间的误差平方和的参数值。
在MATLAB中，可以使用lsqcurvefit函数进行最小二乘估计，该函数可以用来拟合非线性模型或者线性模型。
此外，还可以使用最大似然估计（MLE，Maximum Likelihood Estimation）进行参数估计，MATLAB通过提供相应的函数，如mle函数和mlecov 函数，支持最大似然估计的使用。
最近在做一个项目的方案设计，应各位老总的要求，只有系统框图和器件选型可不行，为了凸显方案设计的高大上，必须上理论分析，炫一下“技术富”，至于具体有多大实际指导意义，那就不得而知了！本人也是网上一顿百度，再加几日探索，现在对用matlab 实现系统辨识有了一些初步的浅薄的经验，在此略做一小节。
必须要指出的是，本文研究对象是经典控制论理最简单最常用的线性时不变的siso 系统，而且是2阶的哦，也就是具有如下形式的传递函数：121)(22++=Ts s T s G ξ 本文要做的就是，对于有这样传递函数的一个系统，要辨识得到其中的未知数T ， ξ!!这可是控制系统设计分析的基础哦，没有系统模型，啥理论、算法都是白扯，在实际工程中非常重要哦！ 经过总结研究，在得到系统阶跃响应实验数据之后（当然如果是其他响应，也有办法可以辨识，在此还是只讨论最简单的阶跃响应实验曲线，谁让你我是菜鸟呢），利用matlab 至少可以有两种方法实现实现（目前我只会两种，呵呵）！一、函数法二、GUI 系统辨识工具箱下面分别作详细介绍！一、 函数法看官别着急，先来做一段分析（请看下面两排红*之间部分），这段分析是网上找来的，看看活跃一下脑细胞吧，如果不研读一下，对于下面matlab 程序，恐怕真的就是一头雾水咯！*******************************************************************************G(s)可以分解为：))((1)(212ωω++=s s T s G其中， [][]11112221--=-+=ξξωξξωTT1ω、2ω都是实数且均大于零。
则有：211ωω=T ，21212ωωωωξ+=传递函数进一步化为：))(()(2121ωωωω++=s s s G 因此，辨识传递函数就转化为求解1ω、2ω。
当输入为单位阶跃函数时，对上式进行拉普拉斯反变换，得系统时域下的单位阶跃响应为：tteet y 212111221)(ωωωωωωωω---+--=即 tteet y 21211122)(1ωωωωωωωω-----=-令1ω=2ωk )1(>k，得tk t ek e k k t y 22111)(1ωω-----=-⎥⎦⎤⎢⎣⎡--=---t k t e k e k k 2)1(2111ωω 对上式两边取以e 为底的对数得[]⎥⎦⎤⎢⎣⎡-+--=---t k e k t k k t y 2)1(211ln 1ln )(1ln ωω 当∞→t 时，⎥⎦⎤⎢⎣⎡---t k e k 2)1(11ln ω0→，则上式化简为[]t k k t y 21ln )(1ln ω--=-该式的形式满足直线方程b at t y +=)(*其中，)(*t y =[])(1ln t y -，1ln ,2-=-=k kb a ω)1(>k通过最小二乘算法实现直线的拟合，得到a ，b 的值，即可得到1ω、2ω的值，进而可得系统的传递函数。
***************************************************Matlab程序代码:clcclose allt=[1 3 5 7 9 11 13 15 17 19];y=[0.149086 0.5890067 0.830617 0.933990 0.973980 0.991095 0.995868 0.998680 0.999490 0.999850];y2=log(1-y);plot(t,y2,'*');grid onpm=polyfit(t,y2,1)value=polyval(pm,t);hold onplot(t,value,'r')title('\fontname{黑体}\fontsize{20}y(t)=at+b')w2=-pm(1)w1=w2/(1-exp(-pm(2)))T=1/sqrt(w1*w2)theta=(w1+w2)/(2*sqrt(w1*w2))z=[];p=[-w1 -w2];k=w1*w2;sys=zpk(z,p,k)figure(2)step(sys,[0:0.5:20]);axis([0 20 0 1.2])hold onplot(t,y,'r*')打开matlab，新建一个Function,把上述程序段拷进去，保存，运行~~~~~~~~~，运行结果：系统的传递函数为)4797.0)(126.1(54034.0)(++=S S S G很顺利吧？先高兴一个！问题接着马上就来了，上面这个例子，这个传递函数的极点刚好都是负实数，因此辨识得很顺利，但是如果系统是欠阻尼系统，也就是如果传递函数的根是复数，那么上述函数段，就无能为力咯，会出现说“matlab 无法处理增益为复数情况之类······” 例如)1)(1(2)(j S j S S G -+++=对于这个系统，若果用simulink 做一下阶跃响应，再把实验数据代入上述函数段，那就不行咯！怎么办呢，只能另辟蹊径了！二、（System Identification Tool）系统辨识工具箱早听说matlab博大精深，神通广大了，于是乎我确定肯定有更简单、直观、强大的工具来完成这小儿科把戏。
查资料琢磨之后，我做了个小实验，在simulink里验证了该种方法。
该方法的大原则是：在确定了系统的输入输出数据（两个列向量N×1形式，如果是1×N，会提示出错！）之后，设计好一定的辨识原则（比如说是2阶？3阶？，传递函数是零极点形式，还是带阻尼形式，等等），然后就交给强大的matlab，得到辨识结果。
Step by step，plz！Step1、建立模型获取系统输入输出数据图1图1系统的输入是阶跃信号，用Scope1监视，并输出到workspace （这步不会的自己百度哦），采样周期是0.1s，得到输入变量u（101×1的矩阵）；本人在系统的阶跃响应上叠加了一白噪声，当然也可以不加噪声，加了噪声就是期望更真实的模拟实际情况，白噪声参数设置见图2图2同样在Scope2监视，也将结果输出到workspace，得到响应数据y（同样也是101×1的矩阵）Step 2、进入辨识工具箱&设置辨识规则直接在command window 输入 ident，回车，进入辨识工具箱图3图3点击import下拉菜单，选时域数据time domain data，见图4图4在下图5红色圈区域输入之前得到的系统输入和输出数据，u和y图5在下图6绿色圈内输入数据的一些信息，因为之前模型中，阶跃起点我是放在0s处的，这里也设置0，如果前面模型仿真是1s，这里应该也是1s；采样时间是0.1s，根据实际情况设置统一哦图6设置完之后，点击import此时界面变成图7图7如果在下图8勾选红框这个选项，就会出现我们刚才设定输入输出数据的曲线，如图9所示，其他勾选项是频域的分析和显示，暂不用它。
图8图9看看与我们实际设置的输入输出是否符合，如果符合，那么我们离成功就不远咯，如果发现异常，那再好好检查一遍，直到确保数据导入没有问题！下面两段红色斜杠之间的内容，对于本实验，可以直接跳过，看一下对后续复杂模型的处理有好处哦，也算全面熟悉一下工具。
/////////////////////////////////////////////////////////// 到这接着选preprocess也就是对数据进行预处理了，下拉菜单中有很多种处理方法和手段，有这个心思的人可以挨个试一下功能。
图10图10预处理的对象是working data中的数据，每进行一种预处理在左边就会有新的数据生成，这时只要将新的数据移动到working data 的那个方框，就可以将working data换成你所想处理的数据了，可以这样多次进行处理，得到你最终想用来辨识的数据和用于验证的数据（不需要的数据可以拖到那个trash里面删除，就是回收站了，也可以从回收站中找回的）接下来就是辨识了，首先把辨识用的数据拖到working data那个方框，再把验证的数据拖到validation data那个方框，这个validation data就是最原始数据稍作处理得到的一个更接近理论模型的对象数据，在这实验里，本人用的就是默认数据，也就是不做任何preprocess处理。
///////////////////////////////////////////////////////////点下拉菜单estimate，选你想要的模型，在弹出的对话框中设定参数,图11图11这里边的模型种类比较多，有线性的、非线性的、状态空间的、经典传递函数形式的等等，我们所选的就是图中红色方框process models，单击。
弹出如下界面图12！图12在本实验中，我们做如下设置，见图13红框标出部分。
图13是否有零点、有延迟、有积分？根据勾选，传函自动调整传递函数极点个数及其他然后点击最下方Estimate，就有模型生成了。
图14图14先勾选上图绿色框选项，看到了什么？给出了拟合率，best fits 98.72有木有？？传递函数具体的数值，双击上图红色框，见图15？图15在右边的数据栏中；也就是model views中了，下面有很多可以选择，每选一个就可以生成一幅对应的图，是由用于验证的数据生成的。
把模型拖到to workspace那个方框，再去看workspace，多的那个变量就是你所辨识出的模型了注意：在某个模型或某组数据上点一下，线变细了就不会在图中显示出来了！！！！！！！！最后再一次提出，上述分析并不一定完全正确，可能有些概念并不清晰可靠，当然还有很多功能并没有被发掘，可以确定的是解决这个问题的两个大方向没有问题，若想做到精益求精，还需要再仔细研究斟酌哦！。
