伯德图

 s=tf('s') 

G=26/(s*(s+5)*(s^2/79^2+s*0.4/79+1));%构造开环传函数

sys=feedback(G,1);%单位负反馈下的闭环传函数

bode(sys);%画闭环伯德图

绘制pressure-flow curve

x = -1;

 while x < 1.1

hold on;

p = -1 : 0.001 : 1;

if x < 0;

q = x * sqrt(1 + p);

else

 q = x * sqrt(1 - p);

end

subplot(121)

plot(p, q, 'b-', 'linewidth', 1.2)

grid on

x = x + 0.2;

end;

 xlabel('p_L/p_s');

ylabel('Q_L/Q_0');

gtext('-1.0', 'fontsize', 10);

gtext('-0.8', 'fontsize', 10);

gtext('-0.6', 'fontsize', 10);

gtext('-0.4', 'fontsize', 10);

gtext('-0.2', 'fontsize', 10);

gtext('0.0', 'fontsize', 10);

gtext('0.2', 'fontsize', 10);

gtext('0.4', 'fontsize', 10);

gtext('0.6', 'fontsize', 10);

gtext('0.8', 'fontsize', 10);

gtext('1.0', 'fontsize', 10);

gtext('x_v/x_{vmax}', 'fontsize', 10)