A procedure based on MATLAB combined with ANSYS is presented and utilized for the multi-objective aerodynamic and structural optimization of horizontal-axis wind turbine (HAWT) blades. In order to minimize the cost of energy (COE) and improve the overall performance of the blades, materials of carbon fiber reinforced plastic (CFRP) combined with glass fiber reinforced plastic (GFRP) are applied. The maximum annual energy production (AEP), the minimum blade mass and the minimum blade cost are taken as three objectives. Main aerodynamic and structural characteristics of the blades are employed as design variables. Various design requirements including strain, deflection, vibration and buckling limits are taken into account as constraints. To evaluate the aerodynamic performances and the structural behaviors, the blade element momentum (BEM) theory and the finite element method (FEM) are applied in the procedure. Moreover, the non-dominated sorting genetic algorithm (NSGA) II, which constitutes the core of the procedure, is adapted for the multi-objective optimization of the blades. To prove the efficiency and reliability of the procedure, a commercial 1.5 MW HAWT blade is used as a case study, and a set of trade-off solutions is obtained. Compared with the original scheme, the optimization results show great improvements for the overall performance of the blade.