/*绘制二元正态分布曲面图*/





%MACRO Graph3DNormal(
                        MuX=1,                                 /*X轴正态分布均数*/
                        VarX=2,                         /*X轴正态分布方差*/
                        MuY=0,                                 /*Y轴正态分布均数*/
                        VarY=4,                         /*Y轴正态分布方差*/
                        Rou=0.75,                         /*二元正态分布相关系数*/
                        XStar=,                                /*X横轴起点，可缺省*/
                        XEnd=,                                 /*X横轴终点，可缺省*/
                        YStar=,                                /*Y横轴起点，可缺省*/
                        YEnd=,                                 /*Y横轴终点，可缺省*/
                        Grid=,                                 /*网格密度，可缺省*/
                        RotateValue=,                 /*图形旋转角度，可缺省*/
                        ZmaxValue=,                        /*Z高度，可缺省*/
                        XTICK=,                                /*X轴分组，可缺省*/
                        YTICK=,                                /*Y轴分组，可缺省*/
                        ZTICK=,                                /*Z轴分组，可缺省*/       
                        DropData=Y);                /*删除绘图数据，Y删除*/       


        OPTIONS NOTES=0 SOURCE=0;
        GOPTIONS RESET=GLOBAL GUNIT=PCT BORDER HTITLE=6 HTEXT=3;

        /*参数判断*/
        %IF (%DATATYP(&MuX.) NE  NUMERIC ) OR
                (%DATATYP(&VarX.) NE  NUMERIC ) OR
                (%DATATYP(&MuY.) NE  NUMERIC ) OR
                (%DATATYP(&VarY.) NE  NUMERIC ) %THEN %DO;

                %PUT %STR(ERROR:参数设置不全！);
                %ABORT;
        %END;



        %IF (%DATATYP(&XStar.) NE  NUMERIC ) %THEN %LET  XStar=%SYSFUNC(CEIL(%SYSEVALF(&MuX.-3*%SYSFUNC(SQRT(&VarX.)))));
        %IF (%DATATYP(&XEnd.) NE  NUMERIC ) %THEN %LET  XEnd=%SYSFUNC(CEIL(%SYSEVALF(&MuX.+3*%SYSFUNC(SQRT(&VarX.)))));
        %IF (%DATATYP(&YStar.) NE  NUMERIC ) %THEN %LET  YStar=%SYSFUNC(CEIL(%SYSEVALF(&MuY.-3*%SYSFUNC(SQRT(&VarY.)))));
        %IF (%DATATYP(&YEnd.) NE  NUMERIC ) %THEN %LET  YEnd=%SYSFUNC(CEIL(%SYSEVALF(&MuY.+3*%SYSFUNC(SQRT(&VarY.)))));
        %IF (%DATATYP(&Grid.) NE  NUMERIC ) %THEN %LET  Grid=0.25;
        %IF (%DATATYP(&RotateValue.) NE  NUMERIC ) %THEN %LET  RotateValue=65;
        %IF (%DATATYP(&ZmaxValue.) NE  NUMERIC ) %THEN %LET  ZmaxValue=0.15;
        %IF (%DATATYP(&XTICK.) NE  NUMERIC ) %THEN %DO;
                %LET  XTICK=11;
        %END;
        %ELSE  %LET  XTICK=%EVAL(&XTICK.+1);

        %IF (%DATATYP(&YTICK.) NE  NUMERIC ) %THEN %DO;
                %LET  YTICK=11;
        %END;
        %ELSE %LET  YTICK=%EVAL(&YTICK.+1);

        %IF (%DATATYP(&ZTICK.) NE  NUMERIC ) %THEN %DO;
                %LET  ZTICK=11;
        %END;
        %ELSE %LET  ZTICK=%EVAL(&ZTICK.+1);

        /*生成数据集*/
        DATA _TEMP_;
                DO X=&XStar. TO &XEnd. BY &Grid.;
                        DO Y=&YStar. TO &YEnd. BY &Grid.;
                                Z=1/( 2*3.14159265359*SQRT(&VarX.*&VarY.*(1-&Rou.**2)) )
                                        * EXP( -1/(2*(1-&Rou.**2))*
                                        (((X-&MuX.)/SQRT(&VarX.))**2+((Y-&MuY.)/SQRT(&VarY.))**2-2*&Rou.*((X-&MuX.)/SQRT(&VarX.))*((Y-&MuY.)/SQRT(&VarY.))    )   ) ;
                                OUTPUT;
                        END;

                END;
        RUN;

        /*绘图3D图*/
        PROC G3D data=_TEMP_;
                PLOT Y*X=Z/ROTATE=&RotateValue. XTICKNUM=&XTICK. YTICKNUM=&YTICK. ZTICKNUM=&ZTICK. ZMAX=&ZmaxValue.;
        RUN;
        QUIT;
       
        /*删除数据*/
        %IF %UPCASE(&DropData.) EQ=Y %THEN %DO;
                PROC DATASETS LIB=WORK NOLIST NOWARN ;DELETE _TEMP_ ;RUN;QUIT;
        %END;
        OPTIONS NOTES=1 SOURCE=1;
%MEND;
%Graph3DNormal(
                        MuX=10,                                 /*X轴正态分布均数*/
                        VarX=3,                         /*X轴正态分布方差*/
                        MuY=8,                                 /*Y轴正态分布均数*/
                        VarY=2,                         /*Y轴正态分布方差*/
                        Rou=0.75,                         /*二元正态分布相关系数*/
                        XStar=,                                /*X横轴起点，可缺省*/
                        XEnd=,                                 /*X横轴终点，可缺省*/
                        YStar=,                                /*Y横轴起点，可缺省*/
                        YEnd=,                                 /*Y横轴终点，可缺省*/
                        Grid=,                                 /*网格密度，可缺省*/
                        RotateValue=,                 /*图形旋转角度，可缺省*/
                        ZmaxValue=,                        /*Z高度，可缺省*/
                        XTICK=,                                /*X轴分组，可缺省*/
                        YTICK=,                                /*Y轴分组，可缺省*/
                        ZTICK=,                                /*Z轴分组，可缺省*/       
                        DropData=Y);                /*删除绘图数据，Y删除*/

这程序太复杂了点吧 我下午找到了个比较简单的，不过非常谢谢你 2# zhangruyang