统计分析软件课程设计---用SAS绘制立体图形及等高线.doc

燕山大学 统计分析软件课程设计 用SAS绘制立体图形及等高线 学 院 理学院 年级专业 04级统计学2 班 姓 名 周春辉 学 号 040108020050 指导教师 张晓冉 教师职称 讲师 完成日期 2007-6-20 燕山大学课程设计评语表 指导教师评语： 成绩： 摘 要 通过对方程的分析，应用g3d 过程绘制曲面图形, 用gcontour过程绘制曲面的等高线图。本文用到的三个方程为 二元正态分布方程 “山”方程 双曲抛物面方程。 关键词：g3d过程 gcontour过程 曲面图形 等高线图 二元正态分布 双曲抛物面 Abstract By the analysis to the equation, draw a curved surface with g3d process,and draw contour map by gcontour process. Three equation the main body of this book used are binary normal distribution equation mountain equation and hyperbolic paraboloid equation. Keywords:　g3dprocess gcontour process curved surface contour map Binary normal distribution hyperbolic paraboloid equation. 目 录 摘 要 I Abstract I 第1节 问题背景 2 1.1问题背景 2 第2节 二维正态分布的图形 2 2.1二维正态分布的方程 2 2.2程序 2 2.3图形输出 3 第3节 山形图案 4 3.1山形图案方程及程序 4 3.2输出图形及解释 5 第4节 双曲抛物面—马鞍面 6 4.1马鞍面方程及绘制程序 6 4.2输出图形 7 第5节 结论 8 参考文献 8 第1节 问题背景 1.1问题背景 通过统计分析软件课程的学习，本人感觉到SAS功能强大，是统计工作者的好帮手。C语言 ，VB， C++， matlab,SAS都是编程语言，就图形功能来说是各有特色。图形的特点是直观，一目了然。我们在解析几何的课程中学习过许多三维图形，如何用SAS软件绘制呢？本文的目的是绘制出给定方程的立体直观图，并给出其等高线，等高线又称等值线，具体说来就是用一族平行平面来截割曲面，进而研究其截口是如何变化的，从这一截曲线的变化情况，我们就能想象出方程所表示的曲面的具体形状，这是一个认识空间图形的重要方法。 对于一个二元函数,我们有了等间隔取值，等间隔取值时的值，就可以调用g3d过程绘制曲面图形,用gcontour过程绘制曲面的等高线图。 第2节 二维正态分布的图形 2.1二维正态分布的方程 例如我们想绘制一个二维正态分布曲面的图形，假设服从联合正态分布，其均值都是，方差都是1，相关系数r=0.5 这时可以得到 的联合密度函数的公式为： 2.2 程序 data normal; /*建立数据集*/ r=0.5; /*相关系数设为5*/ det=1*(1-r*r); do x=-3 to 3 by 0.3; do y=-3 to 3 by 0.3; z=1/(2*3.1415926*det)*exp(-0.5/det*(x*x+y*y-2*r*x*y)); output; end; end; keep x y z; run; proc g3d data=normal; /*调用g3d过程绘制曲面图*/ plot x*y=z; /*是z=f(x,y)的意思*/ run; proc gcontour data=normal; /*绘制等高线图*/ plot x*y=z ; run; 2.3图形输出 第3节 山形图案 3.1山形图案方程及程序 在数学中有一个著名的方程，叫做山方程， 因其形状酷似山形故而得名。应用与前面类似的方法亦可绘制出来，下面是程序。 data hill; do x=-3 to 3 by 0.3; do y=-3 to 3 by 0.3; z=x*exp(-x*x-y*y); /*计算z的值*/ output; end; end; keep x y z; run; proc g3d data=hill; plot x*y=z;