求助如题~~matlab新手求大虾们帮助~~

problem：y=f（x，I），其中I服从一个Bernoulli分布，I=1的概率为α。给定一组w，w是pro（I=1）的估计值。通过bayes法则结合密度函数不断update w，以估计I。得到I之后再用MLE估计f（x，I）中的参数。请问这个程序肿么编~~求大虾们的帮助啊！！