y=randn(500,1);
theta=y*pi;
rose(theta);
title('再极坐标下的直方图');