Increasingly severe crises, such as climate change, water scarcity and environmental pollution, pose significant challenges to global food security and sustainable agricultural development. For efficient and sustainable tomato cultivation management under resource constraints, quantitatively describing the relationship between yield-quality harvest and water-nitrogen application is practically beneficial. Two successive greenhouse experiments with three irrigation levels (1/3 FI, 2/3 FI, and full irrigation (FI)) and four nitrogen fertilizer treatments (0 FN, 1/3 FN, 2/3 FN, and full nitrogen (FN)) were conducted on tomatoes during the whole phenological stage. The tomato evapotranspiration and nitrogen application amount, yield, comprehensive quality, solid–acid ratio, and lycopene content were measured. Based on crop water production functions, three equation forms of water-nitrogen production functions containing 20 models were established and evaluated to predict tomato harvest parameters. The results show that water increased tomato yield while decreasing fruit quality, and the effect of nitrogen was primarily contrary. Water most significantly impacted tomato formation, and the interaction of water and nitrogen changed among different harvest parameters. Tomato yield and quality formation was more sensitive to water and nitrogen at the flowering and fruit maturation stages. Model Singh-2 outweighed other models for yield estimates, with an R2 of 0.71 and an RMSE of 0.11. Singh-Log, Singh-sigmoid and Rao-Root models were effective models for comprehensive quality, solid–acid ratio, and lycopene content prediction, with an R2 of 0.41, 0.62, and 0.42, and an RMSE of 0.33, 0.50, and 0.16, respectively. Finally, models in the form of f(ETi)·f(N) were ideal for tomato harvest prevision and are recommended for water and nitrogen management in tomato cultivation.