理论与试验建模第三章.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 2、最小二乘法求解 a,b * function [num,den]=lsqident(u,y,m,n) %u，y输入输出数列，m，n输入输出数列的阶次 M=length(u); % 获得观测序列u,y的长度 Y=y(n+1:M); for i=1: length(Y) A(i,:)=[-y(n+i-1:-1:i),u(m+i:-1:i)]; End % 建立A,Y阵,A*X=Y,X未知 X= A\Y % A 左除 Y num=X(n+1:length(X)) % 返回差分方程系数 bi den=[1 X(1:n)] % 返回差分方程系数 ai 最小二乘法参数辨识的m函数 * 最小二乘法参数辨识仿真实例分析%crasim1.m bi0=[24.1467,-67.7944,63.4768,-19.8209] %设定参数结构 ai0=[1,-3.6193,4.9124,-2.9633,0.6703] %设定参数结构 u=rand(15,1); %产生随机序列输入 y=dlsim(bi0,ai0,u); %计算输出响应 [bi,ai]=lsqident(u,y,3,4) %调最小二乘计算函数 Fa2=normest(bi-bi0) %误差分析，采用2范数估计 Fb2=normest(ai-ai0) * 最小二乘法参数辨识仿真实例分析（加入观测误差）%crasim2.m bi0=[24.1467,-67.7944,63.4768,-19.8209] %设定参数结构 ai0=[1,-3.6193,4.9124,-2.9633,0.6703] %设定参数结构 u=rand(15,1); %产生随机序列输入 y=dlsim(bi0,ai0,u); %计算输出响应 a=rand(1,length(y))*0.0000001-0 % 产生随机误差矩阵 y=y-a; % 制造带偏差的输出响应 [bi,ai]=lsqident(u,y,3,4) %调最小二乘计算函数 Fa2=normest(bi-bi0) %误差分析，采用2范数估计 Fb2=normest(ai-ai0) * 第四节 建模实例分析 一、理论分析建模并构成Simulink模型 二、阶跃响应法试验获取输入输出数据 三、Simulink仿真确定各环节参数 四、分析误差 * 一、理论分析建模并构成Simulink模型 低温热水散热器热工性能实验台概述 * 1.1 热媒系统 1、高位水箱供水温度传递函数 图中 qi—热水系统电加热器的加热 功率，kW； θ0 、θ1—水箱入口、出口水 的温度，℃ * 1、高位水箱供水温度传递函数 根据能量方程得到： Mb—高位水箱中水的质量，kg； c —水的比热，kJ/kg·℃； θ—水箱的温升，℃ GD—水的质量流量(kg/s)； * 对上式进行拉普拉斯变换得 传递函数 * 2、循环管路传递函数 θ