Abstract: This paper describes the coronary arteries of human being macroscopically and locally, and proposes a method named parameter co-restriction of spherical coordinate to model the heart, which is used to control the spatial configuration of coronary arteries. Then combining the optimum principle of vessel system, a rule-base of vessel bifurcations is built. Finally, the coronary arteries can be simulated by choosing one of the bifurcate modes according to the given probability in each iterated course, which is induced from the idea of iterated function system in fractal theory. At the same time, some random variables are also introduced to make the simulating result more vivid, and the experiment indicates the validity of this method.