?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（‘原始數(shù)據(jù)‘）;
?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