系统函数为H(s)为1/((s+2)(s-1)+k),（2）使用MATLAB命令分别用两种方式画出k=0,1,2,9/4,3时该系统的零极点分布图，并由图讨论k从0增长时，该系统的稳定性变化情况。 （3）对（2）中的稳定系统，使用MATLAB的freqs函数画出它们的线性坐标下的幅频特性和相频特性图，并画出它们的波特图。给出matlab代码

好的，以下是MATLAB代码：

(2) 画出零极点分布图

k = [0, 1, 2, 9/4, 3];

for i = 1:length(k)
    s = tf('s');
    H = 1/((s+2)*(s-1)+k(i));
    zplane([], pole(H));
    title(['k = ', num2str(k(i))]);
    pause;
end

(3) 画出幅频特性、相频特性和波特图

k = [0, 1, 2, 9/4, 3];

for i = 1:length(k)
    s = tf('s');
    H = 1/((s+2)*(s-1)+k(i));
    
    % 幅频特性
    w = logspace(-2, 2, 1000);
    [mag, ~] = bode(H, w);
    subplot(2, 2, i);
    semilogx(w, 20*log10(mag));
    title(['k = ', num2str(k(i))]);
    xlabel('Frequency (rad/s)');
    ylabel('Magnitude (dB)');
    
    % 相频特性
    [~, phase] = bode(H, w);
    subplot(2, 2, i+5);
    semilogx(w, phase);
    xlabel('Frequency (rad/s)');
    ylabel('Phase (deg)');
    
    % 波特图
    [z, p, ~] = tf2zp(H.num{1}, H.den{1});
    subplot(2, 2, i+10);
    zplane(z, p);
    title(['k = ', num2str(k(i))]);
end

希望对你有帮助！