I don't think there is a magic theorem here which works all the time.
What you want depends a lot on the specific $\phi$, ${\bf q}'$, ${\bf q}_0$ you are working with. What you said sounds a lot like a toy model for the renormalization group, which
is a way to study central limit type theorems. For instance the classical central limit
theorem can be interpreted as a convergence in the basin of attraction of the Gaussian
$\mathcal{N}(0,1)$ law with respect to a map $\phi$ on probability distributions
of centered variables with variance 1 given by
$$
Z\longrightarrow \frac{X+Y}{\sqrt{2}}\ .
$$
Namely given a probability distribution for a random variable $Z$, one
makes two independent copies $X$, $Y$, and looks at the probability
distribution of $\frac{X+Y}{\sqrt{2}}$.

In the absence of a precise description of your setup I can only throw in some
general ideas.
The first thing to do is to identify the fixed points of your map $\phi$.
I assume it is nonlinear, so this is not a trivial question. If you cannot solve it
completely, you need to at least find some easy fixed points. I assume this is
what your ${\bf q}'$ is.
Then you need to analyze the linearization of $\phi$
at these easy fixed points. I am again assuming that you did that and
found that this linearization $A$ has all eigenvalues of modulus less than 1
so it is contractive.
From this it follows easily that if you start near ${\bf q}'$ you will converge
to ${\bf q}'$. If your starting point ${\bf q}_0$ is not in this neighborhood
of ${\bf q}'$, then indeed you might want to run some simulations on the computer
not only for ${\bf q}_0$ but also for many other starting points.
In the optimistic situation where the simulation indicates that all points go to
${\bf q}'$, then there could be a Lyapunov function for your map,
see http://en.wikipedia.org/wiki/Lyapunov_function
That's a function which changes monotonically under $\phi$.
If you have such a function which controls the norm of ${\bf q}$
then it might tell you that after enough iterations you will be in the neighborhood
of ${\bf q}'$ where you have the contraction property.
Of course the difficult thing is to come up with the correct guess for this
Lyapunov function if it exists. This guess depends on the specifics of your example.
For the central limit theorem there is a notion of entropy which does the job.

Some references:

- the book <a href="http://books.google.com/books?id=r5XI8a0lYykC&printsec=frontcover&dq=johnson+central+limit&hl=en&ei=4iPETbufAsK90AGRhOyZDg&sa=X&oi=book_result&ct=result&resnum=1&ved=0CCoQ6AEwAA#v=onepage&q&f=false">"Information theory and the central limit theorem"</a>
 by Oliver Johnson, Imperial College Press, 2004.

- the book <a href="http://books.google.com/books?id=tlWOphOFRgwC&printsec=frontcover&dq=sinai+probability&hl=en&src=bmrr&ei=XynETcCtBefn0QGGwpGiCA&sa=X&oi=book_result&ct=result&resnum=1&ved=0CCoQ6AEwAA#v=onepage&q&f=false">"Theory of probability and random processes"</a>
 by Koralov and Sinai which mentions the renormalization group approach
to the CLT.

- you may also find more ideas in the paper by Takashi, Tetsuya and Watanabe
"Triviality of hierarchical Ising model in four dimensions"
Comm. Math. Phys. 220 (2001), no. 1, 13–40. A preprint version is available at
http://www.ma.utexas.edu/mp_arc-bin/mpa?yn=00-397