Field experiments were carried out by quadratic orthogonal rotational combination design with four factors and five levels. Functional models of agronomic characters and biomass of Ophiopogon japonicuswere established forN，P 2O5， K2O，organic fertilizer and Ophiopogon japonicus. The effects ofN，P 2O5，K2O and organic fertilizer application rates onagronomiccharacters，biomass and quality of Ophiopogon japonicuswere studied. Through modelanalysis，N，P 2O5，K2Oand organic fertilizer had significant effects on agronomictraits，biomass and quality of Ophiopogon japonicus. N contributedmost to leafnumber，leaflength，tillernumber，tuber number and water soluble extractcontent；P 2O5 contributed most toleaf width；K 2O contributed most to rootlength，fibrous rootnumber，fresh fibrous rootweight，dry fibrous rootweight，drystem and leafweight，fresh yield and total saponincontent；organic fertilizer contributed most to fresh weight of stem and leaf and dry yield. N-organic fertilizer had significant effect on stem and leaf dryweight；N -organic fertilizer interaction hadsignificant effect on freshyield；K 2O-organic fertilizer interaction had extremely significant effect on fresh yield and dryyield； N-K2O interaction had significant effect on total saponin content of Ophiopogon japonicus；N-K2O interaction and N-organic fertilizer interaction had significant effect on water-soluble extract content. The interaction of P2O5 -organic fertilizer had a very significant effect on the content of water-soluble extracts. Through frequency analysis and crossingoptimization，the optimum fertilization schemes for high yield and quality of Ophiopogon japonicuswere obtained as follows:N 676.75～ 819.16 kg/ hm²，P2O5183.07～ 206.30 kg/hm²，K2O 331.07～ 437.53 kg/hm²or 467.31～ 568.76 kg/hm²，and commercial organicfertilizer 2 191.20～ 2732.40 kg/hm².