给我一段matlab程序，描述：外界冲击的到达服从强度μ为0.3的HPP过程，到达时间是均匀分布的，每次冲击到达的强度服从均值为0.4，标准差为0.2的正态分布，给我画图

在MATLAB中，我们可以使用`randomStream`生成高斯混合过程（Generalized Pareto Process, GPP）和均匀分布的时间间隔，以及正态分布的冲击强度。以下是一个简单的示例程序：

% 设置随机数种子以保证结果可复现
rng('default')

% 定义过程参数
mu = 0.3; % 强度的平均值 (scale parameter for HPP)
sigma = 0.4; % 冲击强度的均值
std_dev = 0.2; % 冲击强度的标准差
time_uniform_mean = 5; % 时间间隔的均值
time_uniform_std = 1; % 时间间隔的标准差

% 创建HPP和均匀分布的伪随机数流
hpp_stream = RandStream('hpd', mu);
time_stream = RandStream('Uniform', [0 time_uniform_mean + time_uniform_std]);

% 定义循环次数
num_shocks = 100;

% 初始化数组存储数据
arrival_times = zeros(1, num_shocks);
shock_strengths = zeros(1, num_shocks);

for i = 1:num_shocks
    % 从HPP过程获取到达时间
    arrival_time = hpp_stream.rvs();
    
    % 从均匀分布生成时间间隔，并添加到当前到达时间
    inter_arrival_time = time_stream.rvs();
    arrival_times(i) = arrival_time + inter_arrival_time;
    
    % 从正态分布获取冲击强度
    shock_strength = normrnd(sigma, std_dev); % 使用normrnd函数生成正态分布
    shock_strengths(i) = shock_strength;
end

% 制作条形图表示冲击时间和强度
figure
bar(arrival_times, shock_strengths)
xlabel('到达时间')
ylabel('冲击强度')
title('冲击事件时间和强度分布')

% 记录每个冲击的点在直角坐标系上
hold on
scatter(arrival_times, shock_strengths, 'o', 'filled')
hold off

% 提问部分: