clear all; close all; clc dt = 0.01; tmax=20; t=0:dt:tmax; s=tf('s'); %设定待辨识传递函数 w=3; %自然频率 f=0.5; %欠阻尼系数 H=w^2/(s^2+2*f*w*s + w^2) %传递函数 %设定输入的阶跃函数,并画出输入与输出函数