贝叶斯信息准则的matlab 程序

以下是一个计算贝叶斯信息准则（BIC）的 Matlab 程序：

function bic = BIC(data, model, params)
% data: 数据
% model: 模型
% params: 模型参数

% 计算负对数似然函数值
n = length(data);
llh = -n/2*log(2*pi) - n/2*log(model(data, params)) - sum(log(normpdf(data, model(data, params), 1)));

% 计算 BIC
num_params = length(params);
bic = -2*llh + num_params*log(n);
end

其中，`data` 是观测数据，`model` 是一个函数句柄，用于计算数据的概率密度函数值，`params` 是模型的参数。该程序首先计算负对数似然函数值，然后使用 BIC 公式计算 BIC 值。程序中使用了 Matlab 内置的正态分布概率密度函数 `normpdf`。