The design and minaturization of microwave circuits involves the analysis and optimization of a large varity of structures, which, in general, are based on arbritrary objects embedded in also arbritary media. The most common situation is found when optimizing terminal antennas, whose main constraint is the antenna's volume, which must be fitted into a specified space, which in general is the phone or laptop case. This kind of antennas can also include stratified media, oblique metallizations, dielectric objects, waveguide components, in their definition making their analysis more difficult from the numeric point of view. The complete treatment of this kind of problems has been performed traditionally using Finite Differences or Finite Elements software, while it is more rare to find Method of Moments implementations for the solution of this kind of general problems. The Method of Moments has a certain number of advantages as the inclusion of radiation conditions or stratified media in the Green's function of the problem. Due to this it is possible to mesh only the surfaces or boundary conditions that are not included in the Green's function definition, obtaining, for instance, a more natural radiation boundary condition compared to other fullwave methods. The other big concern is the optimization of this kind of structures. In general the structures we are dealing with in this work can have very irregular shapes, in order to fit within the required volume constraints. The traditional optimizations techniques based on the gradient or conjugated gradient are not efficient when optimizing functions containning multiple minima and a large search space, which is the case in the optimization of terminal antennas. When miniaturizing terminal antennas or just a simple microstrip filter, the number of unknowns can be very large, and the different optimization variables can vary within wide limits in the search space. In this situation the Genetic Algorithms are ideal for seeking the global optimum solution of our problem. This work covers the analysis and optimization of microwave circuits, including antennas and filters in microstrip and waveguide technologies, paying more attention to the implementation of the Method of Moments (MoM) in conjugation with the Mixed Potential Integral Equation (MPIE) and Electric/Magnetic Field Integral Equation (EFIE/MFIE). The work has been split into 3 main goals: The analysis of 3 dimensional structures in stratified media, including arbritary shaped dielectric bodies using the Method of Moments, (MPIE,EFIE,MFIE). The analysis of shielded environments and more precisely the analysis of retangular waveguide cavities with the Method of Moments. The implementation of several optimization techniques, emcompassed in the frame of genetic algorithms, for the design and miniaturization of terminal antennas and microstrip filters. Some relevant examples linked to the efficiency and accuracy of the different methods and techiques explained in this work are included in each chapter, yielding to practical results suitable to be used in real life design problems. The main contributions of this work can be summarized as follows: In Chapter 2 we have developed the interpolation of 3D Green's function to the general solution of dielectric objects embedded in multilayered media. The interpolation method itself is preceeded by a spectral extraction of the quasi-static terms. This technique, very often applied in complex images techniques, has been applied to our specific Green's function problem yielding very good results in terms of accuracy and Green's function computation time. Another contribution, also emcompassed in Chapter 2, is the integration of the static part of the field dyadic G̿EM, performed by splitting the field dyadic into TM and TE components, allowing then an analytical integration of the dyadic terms in stratified media. In the frame of shielded environments in Chapter 3, the most important contribution is the extension of Ewald's acceleration technique to full electric and magnetic problems, allowing the acceleration of the field dyadic G̿EM and thus permit the extension of the MoM to the simulation of arbritary metallic shapes coupled to apertures through Ewald's approach. In Chapter 4 the most important contribution is the implementation of a Bayesian optimiser based on the estimation of probability made by dependencies trees. The dependency tree based method found in [1] and later in [2] has been adapted to the terminal antenna miniaturization yielding to excellent results in optimisation time and quality of the solution provided by the optimiser. Our implementation of the dependency tree method was found to be better than the traditional method used in this kind of problems. Chapter 4 contains also an implementation of the growing cells method presented in [3] for the specific problem of optimising pseudo periodic structures, like microwave filters, yielding to a powerful interpolating method which accelerates the optimisation process of microwaves devices whose fitness response is very time consuming.