You can write a Python code to distribute the fibers with random positions, aspect ratios and also curvature radius inside the matrix. There are plenty of articles involving the same issue in the literature.
There are different strategies depending on the application, computation cost and the level of accuracy you want.
For example if you want to study the interaction between the fibers and the matrix you should go with a full field modeling approach. Or if you are interested in the overall behavior of the composite structure you can go with a RVE homogenization approach.
You can write a Python code to distribute the fibers with random positions, aspect ratios and also curvature radius inside the matrix. There are plenty of articles involving the same issue in the literature.