Abstract: "Flora Reipublicae Popularis Sinicae "and "Flora of China" had different opinions on the subspecies treatment of Begonia grandis Dryander,which could not be applied in practice according to their classification criteria,and the issues still remain unsolved in the subspecies classification of this species. In order to understand intraspecific diversity of B.grandis,a numerical classification study was conducted with 33 phenotypic traits from 39 wild B.grandis populations. By quantity characters variation analysis,quality traits genetic diversity analysis,principal component analysis and R-type cluster analysis,the main traits were screened out from 33 traits for Q-type cluster analysis. The phenotypic diversity of B. grandis was rich,the variation coefficient of quantity traits was 0.077-0.441,and the Shannon-Wiener diversity index of quality characters was 0.287-1.483. Principal component analysis showed that 33 traits could be integrated into 6 principal components,and the cumulative contribution rate was 72.28%. R-type cluster analysis showed that the characters were relatively independent except for the correlation between leaf size and male flower size.Q-type cluster analysis classified 39 B.grandis populations into four groups with male flower size,leaf size,plant size,leaf color,which was inconsistent with three subspecies classified by Flora of China. Therefore,it was necessary to combine the molecular evidence to discuss the rational intraspecific classification of B.grandis in the future. It was interesting that Japanese naturalized populations of B.grandis very possibly originated from Tianmu Mountain population in Zhejiang Province.