帮我写出一个matlab的代码：把一组离散数据设置一个阈值，高于阈值的保留原先的值；低于阈值的都等于阈值，然后找出并标记出每个峰值，找出最大的峰并输出它两边的波谷的序列

matlab代码示例：
data = rand(1, 20); % 生成随机数据
threshold = 0.5; % 设置阈值
for i = 1:length(data)
if data(i) < threshold
data(i) = threshold; % 设置低于阈值的数据都等于阈值
end
end
[pks, locs] = findpeaks(data); % 找出峰值
[max_pk, max_loc] = max(pks); % 找出最大的峰
% 找出最大峰两边的波谷
valey1 = data(1:max_loc-1);
valey2 = data(max_loc+1:end);