对目标函数target_dist = @(x) x(1).*x(2)-1500用M-H生成 MCMC链的代码

### 回答1：
以下是使用Metropolis-Hastings算法生成MCMC链的Matlab代码，目标函数为target_dist = @(x) x(1).*x(2)-1500：

% 目标函数
target_dist = @(x) x(1).*x(2)-1500;

% 初始状态
x0 = [100, 15];

% 提议分布
q = @(x, y) normrnd(x, y);

% MCMC参数
N = 10000; % 迭代次数
burnin = 1000; % 燃烧期
thin = 10; % 抽样间隔

% 存储MCMC链
chain = zeros(N, 2);

% 初始状态为第一个样本
x_current = x0;

% MCMC迭代过程
for i=1:N
    
    % 从提议分布中生成候选样本
    x_proposed = q(x_current, [10, 1]);
    
    % 计算接受率
    alpha = min(1, target_dist(x_proposed)/target_dist(x_current));
    
    % 决定是否接受候选样本
    if rand < alpha
        x_current = x_proposed;
    end
    
    % 存储当前状态
    chain(i, :) = x_current;
    
end

% 剔除燃烧期
chain = chain(burnin+1:end, :);

% 抽样
samples = chain(1:thin:end, :);

% 绘制样本分布
figure;
scatter(samples(:,1), samples(:,2));
xlabel('x_1');
ylabel('x_2');
title('MCMC样本分布');

在上述代码中，我们使用正态分布作为提议分布，根据当前状态生成候选样本。然后计算接受率alpha，决定是否接受候选样本。最后存储当前状态，重复迭代N次，得到MCMC链。我们剔除燃烧期，然后按照抽样间隔thin进行抽样，得到最终的MCMC样本。最后，我们绘制样本在二维平面上的分布。

### 回答2：
M-H算法是一种用于生成马尔科夫链的一种随机采样算法，并且在贝叶斯统计中经常被使用。下面是使用M-H算法生成MCMC链的代码来计算目标函数target_dist = @(x) x(1).*x(2)-1500。

% 设置参数
num_iterations = 10000;  % 迭代次数
proposal_sigma = 0.1;   % 提议分布标准差
initial_x = [100; 15];  % 初始值

% 初始化
x_current = initial_x;
chain = zeros(2, num_iterations);
accept_count = 0;

% 进行迭代
for i = 1:num_iterations
    % 生成提议点
    x_proposed = x_current + proposal_sigma * randn(2, 1);
    
    % 计算接受率
    acceptance_prob = min(1, target_dist(x_proposed) / target_dist(x_current));
    
    % 决定是否接受提议点
    if rand() < acceptance_prob
        x_current = x_proposed;
        accept_count = accept_count + 1;
    end
    
    % 记录链
    chain(:, i) = x_current;
end

% 输出结果
acceptance_rate = accept_count / num_iterations;
disp(['Acceptance rate: ', num2str(acceptance_rate)]);
disp(['Mean: ', num2str(mean(chain, 2)')]);

% 绘制散点图
figure;
scatter(chain(1, :), chain(2, :), 'filled');
xlabel('x1');
ylabel('x2');
title('MCMC Chain');

在这个代码中，我们首先设置了所需的参数，包括迭代次数、提议分布标准差和初始值。然后，我们初始化当前点、链和接受次数。接下来，我们进行了所需的迭代次数，并在每次迭代中生成提议点、计算接受率，并决定是否接受提议点。同时，我们记录了接受的点来构建MCMC链。最后，我们输出了接受率和MCMC链的均值，并绘制了一个散点图来展示MCMC链在参数空间中的演化。

### 回答3：
对目标函数target_dist = @(x) x(1).*x(2)-1500的M-H方法MCMC链的代码如下：

% 定义目标函数
target_dist = @(x) x(1).*x(2)-1500;

% 定义M-H链的参数
burn_in = 1000; % 燃烧期，用于去除初始的不稳定样本
num_samples = 10000; % 采样数量
proposal_std = 0.1; % 提议分布的标准差

% 初始化链的初始样本
current_sample = [0; 0];
samples = zeros(2, num_samples);

% Metropolis-Hastings方法
for i = 1:(burn_in + num_samples)
    % 生成候选样本
    candidate_sample = current_sample + proposal_std * randn(2, 1);
    
    % 计算接受率
    accept_ratio = target_dist(candidate_sample) / target_dist(current_sample);
    
    % 接受或拒绝候选样本
    if rand() < accept_ratio
        current_sample = candidate_sample;
    end
    
    % 燃烧期结束后保存样本
    if i > burn_in
        samples(:, i - burn_in) = current_sample;
    end
end

% 绘制样本分布
scatter(samples(1, :), samples(2, :));
xlabel('x1');
ylabel('x2');
title('M-H MCMC链样本分布');

以上代码通过Metropolis-Hastings方法生成了一个MCMC链，并且保存了在燃烧期之后的样本。最后通过散点图展示了样本分布，其中x轴表示x1，y轴表示x2。