This paper makes further researches on the Gradient Regularization method, and not only establishes a fresh iterative gradient regularization algorithm in a brief way but also gives a practical way to choose regularization parameters. These numerical examples which identify parameters of elliptic operator equations not only show that GR method has a wide scope of application and poses good capability of noise resistance, but also illustrate that GR method poses higher convergence speed and larger convergence scope.