过程控制课设

matlab simulink仿真加热炉温度控制系统，采用PID、串级控制两种方案。

 x=[0 1 3 4 5 8 10 15 16.6 20 25 30 40 50 60 70 80 90];
 y=[0 0.46 1.7 3.7 9.0 19.0 26.4 36.0 37.5 33.5 27.2 21.0 10.4 5.1 2.8 1.1 0.5 0.2];
 h=0:1:100;
 yh=interp1（xyh‘cubic‘）;                %采用立方插值
% plot（h*60yh‘r‘）;
%hold on
 %求取单位阶跃响应值ys
 ys=1:1:19;                                %定义一个间隔为1的初始值为1，终值为19的向量
 for i=2:1:19
     ys（i）=yh（（i-1）*5+1）+ys（i-1）;          %注意yc的间隔是1，ys的间隔是5
 end
 %进行差值处理
 xs=[0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90];
 m=0:1:100;
 ym=interp1（xsysm‘cubic‘）;
 plot（m*60ym‘r‘）;
%  xlabel（‘时间/min‘）;
%  ylabel（‘温度/。c‘）;
%  title（‘矩形脉冲响应与阶跃脉冲响应曲线‘）;
%  grid on
 hold on
%归一化处理
%  ym=ym/203.16;
%  figure
%  plot（m*60ym‘r‘）;
%  axis（[0 6000 0 1]）;
 
% xlabel（‘时间/min‘）;
% ylabel（‘温度/。c‘）;
% title（‘归一化后的阶跃脉冲响应曲线‘）;
% grid on
% hold on
%gtext（‘x=16.27y=0.4‘）;
%gtext（‘x=32.57y=0.8‘）;
%验证传函
 num=6.77;
 den=conv（[452.4 1]conv（[452.4 1][452.4 1]））;
 G0=tf（numden）;
 step（G0*30）;
 gtext（‘原始数据‘）;
 gtext（‘求出的对象‘）;
 grid on

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1125  2011-12-03 16:23  shujuchuli.m

     文件      26359  2011-12-03 09:02  yanzhengchuanhan.fig

     文件      53699  2011-12-03 10:45  zhendangqiuxian.fig

     文件      27547  2011-12-03 09:10  Z-N.fig

     文件      91194  2011-12-03 16:35  zong.mdl

     文件      70214  2011-12-03 12:00  chuanjiduibi.fig

     文件      52219  2011-12-03 11:37  chuanjiduibi.mdl

     文件      14647  2011-12-03 11:22  fukongduixiang.fig

     文件      29345  2011-12-03 11:23  fukongduixiang.mdl

     文件     178887  2011-12-03 11:43  genzong.fig

     文件      56567  2011-12-03 11:43  genzong.mdl

     文件      12222  2011-12-03 08:43  guiyihua.fig

     文件      15166  2011-12-03 08:52  shuju.fig

----------- ---------  ---------- -----  ----

               629191                    13