jk2_5.m

来自「利用matlab仿真系统的响应」· M 代码 · 共 30 行

clc;
clear;
close all;
num=[330]; den=[72,1];
s1=tf(num,den);
[num1 den1]=pade(8,3);
s2=tf(num1,den1);
Kp=0.02
Ki=0.00033
Kd=0.02
num3=[Kd,Kp,Ki]; den3=[1,0];
s3=tf(num3,den3);
s=s1*s2*s3;
G1=feedback(s,1);
step(G1);
hold on;

Kp=0.015
Ki=0.0002
Kd=0.015
num3=[Kd,Kp,Ki]; den3=[1,0];
s3=tf(num3,den3);
s=s1*s2*s3;
G1=feedback(s,1);
step(G1);
hold on;
gtext('Kp=0.02 Ki=0.00033 Kd=0.02');
gtext('Kp=0.015 Ki=0.0002 Kd=0.015');