帮我用matlab写一段关于蒙特卡洛模拟的求总成本的程序，总成本为几个成本之和，分别为购置成本198900万元，标准差为4500.2，运营成本为120050万元，标准差为398604，故障成本为30710万元，标准差为391.2，报废成本为-1572.5万元，标准差为121.3.报废成本作为负值加入总体成本，用蒙特卡洛模拟几个成本综合，并输出样本直方图

好的，下面是一个使用蒙特卡洛方法求总成本并输出样本直方图的 Matlab 程序：

% 定义购置成本、运营成本、故障成本、报废成本的均值和标准差
mu1 = 198900; sigma1 = 4500.2;
mu2 = 120050; sigma2 = 398604;
mu3 = 30710; sigma3 = 391.2;
mu4 = -1572.5; sigma4 = 121.3;

% 定义模拟次数
N = 10000;

% 生成购置成本、运营成本、故障成本、报废成本的随机数
cost1 = normrnd(mu1, sigma1, N, 1);
cost2 = normrnd(mu2, sigma2, N, 1);
cost3 = normrnd(mu3, sigma3, N, 1);
cost4 = normrnd(mu4, sigma4, N, 1);

% 计算总成本
total_cost = cost1 + cost2 + cost3 + cost4;

% 输出结果
fprintf('总成本的均值为：%f\n', mean(total_cost));
fprintf('总成本的标准差为：%f\n', std(total_cost));

% 绘制样本直方图
histogram(total_cost, 50);
xlabel('总成本');
ylabel('频数');
title('总成本样本直方图');

在这个程序中，我们添加了一个新的成本项：报废成本，并使用 `histogram()` 函数绘制了总成本的样本直方图。由于报废成本是负值，因此我们直接将其加入总成本中。你可以根据自己的需求调整代码中的参数和模拟次数，以得到更准确的结果。