x(t+1)=2x(t)+u(t),y(t)=-2x(t-k),其中k是小于4个采样周期的随机延时，用Matlab绘制该系统谱半径随时延k大小变化情况

可以使用Matlab中的for循环和eig函数来实现该系统谱半径随时延k大小变化情况的绘制。具体的代码如下：

% 系统参数设置
T = 1; % 采样周期
N = 100; % 模拟时长
x = zeros(N, 1); % 系统状态
y = zeros(N, 1); % 系统输出
u = randn(N, 1); % 随机输入

% 绘制系统谱半径随时延k大小变化情况
rho = zeros(1, 4); % 存储系统谱半径
for k = 1:4 % 循环变量k
    % 计算系统状态和输出
    for t = k+1:N
        x(t) = 2*x(t-1) + u(t-1);
        y(t) = -2*x(t-k);
    end
    % 计算系统特征值并存储谱半径
    A = [2 0; -2 0];
    eig_val = eig(A);
    rho(k) = max(abs(eig_val));
end
plot(1:4, rho, 'o-');
xlabel('时延k');
ylabel('系统谱半径');
title('系统谱半径随时延大小变化情况');

运行上述代码，可以得到系统谱半径随时延k大小变化情况的图像。