I need to do a simulation to evaluate an integral of a 3 parameter function, we say $f$, which has a very complicated formula. It is asked to use MCMC method to compute it and implement the Metropolis-Hastings algorithm to generate the values distributed as $f$, and it was suggested to use a 3 variated normal as proposal distribution. Reading some examples about it, I've seen that some of then use a normal with fixed parameters $N(\mu, \sigma)$ and some use with a variable mean $N(X, \sigma)$, where $X$ is the last accepted value as distributed according $f$. I have some doubts about both approaches:
1) What is the meaning of choosing the last accepted value as the new mean of our proposal distribution? My intuition says it should guarantee that our values will be closer to values distributed as $f$ and the chances of acceptance would be greater. But doesn't it concentrates too much our sample? It is guaranteed that, if I get more samples, the chain will become stationary?
2) Wouldn't choose fixed parameters (since the $f$ is really difficult to analyize) be really difficult and dependent of the first sample we need to choose to start the algorithm? In this case, what would be the best approach to find which one is better?
Is one of those approaches better than the other or this depends of the case?
I hope my doubts are clear and I would be glad if some literature could be given (I've read some papers about the theme, but more is better!)
Thanks in advance!