用MATLAB写一段双样本t检验代码

以下是用MATLAB实现的双样本t检验代码：

% 双样本t检验（两组样本独立）
% H0：两组样本的均值相等
% H1：两组样本的均值不相等

% 输入数据
x1 = [1, 2, 3, 4, 5];
x2 = [6, 7, 8, 9, 10];

% 计算均值和标准差
mean_x1 = mean(x1);
mean_x2 = mean(x2);
std_x1 = std(x1);
std_x2 = std(x2);

% 计算t检验统计量
t = (mean_x1 - mean_x2) / sqrt(std_x1^2/length(x1) + std_x2^2/length(x2));

% 计算p值
df = length(x1) + length(x2) - 2; % 自由度
p = 2 * (1 - tcdf(abs(t), df)); % 双侧检验

% 输出结果
fprintf('t检验统计量为 %.4f\n', t);
fprintf('p值为 %.4f\n', p);
if p < 0.05
    fprintf('p值小于0.05，拒绝原假设H0\n');
else
    fprintf('p值大于等于0.05，不能拒绝原假设H0\n');
end

其中，`x1`和`x2`是需要进行双样本t检验的两组数据，`mean_x1`和`mean_x2`分别是两组数据的均值，`std_x1`和`std_x2`分别是两组数据的标准差。然后，根据t检验的公式，计算出t检验统计量`t`，再根据自由度计算出p值。最后，根据p值是否小于显著性水平（例如0.05）来判断是否拒绝原假设。