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) |
|
来自: 行者hydraulics > 《MATLAB》