MATLAB官方示例：https://ww2.mathworks.cn/help/ident/ug/estimating-transfer-function-models-for-a-heat-exchanger.html

0. 数据准备

扫频数据如下

1. 代码（.m文件）

close all; clear all; clc
%load td_xandy
A=xlsread('E:\209laboratory\7-29\td_xandy7-30.xls',1);
pt = A(:,2);
ct = A(:,3);
data = iddata(pt,ct,1); %输出p 输入c Ts为采样时间
data.InputName  = '\Delta CTemp';
data.InputUnit  = 'C';
data.OutputName = '\Delta PTemp';
data.OutputUnit = 'C';
data.TimeUnit   = 'minutes';
plot(data)

sysTF = tfest(data,3,0,nan)%估计传递函数，已知极、零点个数

其中，iddata函数用于封装输入、输出、采样时间信息。

2. 结果(稍等片刻才能出结果)

3. 如果要求离散传递函数：