matlab如何对复杂的函数进行极大似然估计

输入数据θ是一组弧度，然后生成如图的频率直方图(坐标还没改成弧度，是度数)，再用下图的函数进行拟合，如何用matlab对复杂函数进行极大似然估计算出拟合的参数；其中的一个函数为第一类零阶贝塞尔函数，matlab有自带的代码besselj，根据输入数据θ求里面的12个参数

不知道你要求什么？如果是求这个f里面的参数，而且已经有theta和f(theta)的一组数据，且f(theta)服从高斯分布N(f(theta),sqrt(f(theta))的话，极大似然法的步骤是把这一组概率相乘，然后对总概率求最大值。

注意前面有一些假设，包括对于分布的方差的假设。

一般最大似然法需要用到拉格朗日乘子法，这部分可以查找一下有没有相关的现成函数，不然手动实现也是比较麻烦且低效的