Sparked by the Willems' fundamental lemma, a class of data-driven control methods has been developed for LTI system. At the same time, the Koopman operator theory attempts cast a nonliner control problem into a standard linear one albeit infinite dimensional. Motivated by these two ideas, a data-driven control scheme for nonlinear systems is proposed in this work. The proposed scheme is compatible with most differential regressor enabling an offline learning. In particular, the model uncertainty is considered, enabling a novel data-driven simulation framework based on Wasserstein distance. Numerical experiments are performed with Bayesian neural networks to show the effectiveness of both the proposed control and simulation scheme.