一阶液位系统根据系统输入和输出求解传递函数

1、内容简介

略
302-可以交流、咨询、答疑

2、内容说明

略

3、仿真分析

clc
close all
clear
%% 上水槽
data1=[
0	0
7.5	6.636
15.3	11.81
22.6	15.756
30	18.963
37.5	21.5
45	23.207
60	26.021
75	27.689
90	28.506
120	29.78
149.7	29.798
240.7	29.99
300	29.999
];
t1 = data1(:,1);
h1 = data1(:,2);
figure
plot(t1,h1)
xlabel t/s
ylabel h1/cm
title 上水槽

t_1=0:t1(end);
x=1;%单位阶跃信号
y_1=interp1(t1,h1,t_1);
error=0.1;%计算精度
%一阶环节参数确定
[k0_1,t0_1]=ordre1(t_1,x,y_1,error);
sys1=tf(k0_1,[t0_1 1])
y0_1=step(sys1,t_1);
figure
plot(t1,h1,'r',t_1,y0_1,'b')
legend('原始传递函数响应','拟合传递函数响应')
title 上水槽
%% 下水槽
data2=[
0	0
11.35	9.954
22.7	17.82
33.3	23.943
45	28.745
56.3	32.108
67.8	35.021
90	38.912
112.5	41.324
135.3	42.759
180	44.376
225.6	44.697
360	44.885
450	44.998
];
t2 = data2(:,1);
h2 = data2(:,2);
figure
plot(t2,h2)
xlabel t/s
ylabel h1/cm
title 下水槽

t_2=0:t2(end);
x=1;%单位阶跃信号
y_2=interp1(t2,h2,t_2);
error=1;%计算精度
%一阶环节参数确定
[k0_2,t0_2]=ordre1(t_2,x,y_2,error);
sys2=tf(k0_2,[t0_2 1])
y0_2=step(sys2,t_2);
figure
plot(t2,h2,'r',t_2,y0_2,'b')
legend('原始传递函数响应','拟合传递函数响应')
title 下水槽