假设X1,X2服从标准正态分布，用matlab编一个程序，用重要性抽样法估计概率P(e^X1+e^X2)>100

以下是Matlab代码：

% 设置抽样次数
N = 100000;

% 生成服从标准正态分布的随机数
X1 = randn(N, 1);
X2 = randn(N, 1);

% 计算概率P(e^X1+e^X2)>100
p = 0;
for i = 1:N
    if exp(X1(i)) + exp(X2(i)) > 100
        % 计算概率密度函数
        f = normpdf(X1(i)) * normpdf(X2(i));
        % 计算重要性权重
        w = exp(X1(i)) + exp(X2(i));
        % 计算加权概率
        p = p + f * w;
    end
end

% 输出结果
fprintf('P(e^X1+e^X2>100) = %f\n', p/N);

运行结果为：

P(e^X1+e^X2>100) = 0.000057