编写matlab代码实现：当k取某一数值时，Pin在0到20范围内找到Pout=Pin*(1-2*(1+cos(Pin*k)))的最大值与最小值的差值T，假设k的取值范围为0到pi，画出x轴为k，y轴为T的曲线图。

在MATLAB中，你可以使用`fminbnd`函数来找到最大值和`fmaxbnd`函数来找最小值，同时使用循环遍历k的取值，并存储每个k对应的T值。下面是一个简单的示例代码：

% 定义函数
function f = myFunction(k, Pin)
    Pout = Pin * (1 - 2 * (1 + cos(Pin * k)));
    T = abs(max(Pout) - min(Pout));
    f = T;
end

% 定义初始k的范围和步长
kMin = 0;
kMax = pi;
dk = 0.01; % 可以调整k的精度

% 初始化最大值和最小值，以及对应的k值
maxT = Inf;
minT = -Inf;
kMaxVal = [];
kMinVal = [];

for k = kMin:dk:kMax
    % 调用函数并更新结果
    [T, success] = fminbnd(@(p) -myFunction(p, 20), 0, max(20, round(20/(Pin*dk))));
    if success
        if T > maxT
            maxT = T;
            kMaxVal = k;
        end
        
        [T, success] = fmaxbnd(@(p) -myFunction(p, 0), 0, round(20/(Pin*dk)));
        if success
            if T < minT
                minT = T;
                kMinVal = k;
            end
        end
    end
end

% 画出曲线图
figure;
plot(k, [minT, maxT]);
hold on;
grid on;
xlabel('k');
ylabel('T');
title(['最大值与最小值之差随k的变化曲线 (Pin = 20)']);

% 显示最大值和最小值对应的k值
fprintf('Max T occurs at k = %.4f with a value of %.4f\n', kMaxVal, maxT);
fprintf('Min T occurs at k = %.4f with a value of %.4f\n', kMinVal, minT);