matlab画图说明

函数的可视化与Matlab作图（Matlab）实验4),(~ln2NXY0,00,21)(22)(ln21xxexxpx一、问题研究常见初等函数的图形。设某保险公司的赔款额统计表明，若某笔赔款额为随机变量X(元)，且，称X服从对数正态分布，其密度函数为设均值μ=6.012，方差σ2=1.792。试作出对数正态分布密度图形。二、实验目的学会使用Matlab软件进行函数的图形显示。三、预备知识1.常见初等函数的图形2.本实验所用Matlab命令提示：（1）绘制二维的最常用的命令是plot,基本用法为：plot(x,y)：以x为横坐标、y为纵坐标绘制函数或数据图形plot(y):绘制一个向量的图形，相当于x=[1,2,…,length(y)]时的情形在同一图形窗口中绘制多个曲线：可用plot(x,[y1;y2;…])命令，其中x是横坐标向量，[y1;y2;…]是由几个函数的纵坐标拼成的矩阵，这时Matlab将用不同颜色的曲线表示不同的函数。（2）图形文字标注命令：xlabel(‘x轴名称’);ylabel(‘y轴名称’)；zlabel(‘z轴名称’)title(‘图形标题’);text(‘说明文字’);用鼠标在特定位置输入文字。legend命令：可用于说明不同函数所使用的线型。（3）绘制三维图形的命令：meshgrid命令：用于形成二维网格点。然后计算在这些网格上的函数值，为绘制图形做准备。三维绘图命令：plot3(x,y,z)%线条图mesh(x,y,z)%网格图meshc（x,y,z）%具有基本等高线的网格图suef（x,y,z）%表面图suefc（x,y,z）%具有基本等高线的表面图ezplot3(x,y,z,’animate’):Matlab6.x提供了ezplot3命令，产生三维动画图形。（4）条形图bar(x,y)%竖直条形图barh(x,y)%水平条形图bar3(x,y)%三维竖直条形图barh3(x,y)%三维水平条形图（5）线条和颜色控制符:.点y黄。小圆圈m棕色xx标记c青色+加号r红色-实线g绿色*星号b蓝色：虚线w白色-.点划线h黑色0,00,21)(22)(ln21xxexxpx22yxr四、实验内容与要求1.试作出保险经营及风险管理中赔款额分布模型的图形其中对数正态分布密度函数为（通过100笔赔款样本计算得：均值μ=6.012，方差σ2=1.792。）2、在区间[-π,π]内绘制sin(x)、cos(x)的曲线。3、绘制z=sin(r)/r函数的三维图形。(其中：4、绘制由参数表达式x=cos(t),y=tsin(t),z=sqrt(t)描述的图形)22yxz五、思考与练习1.根据索赔额的频数分布画直方图及三维直方图索赔额的频数分布索赔额（元）索赔频数0—20004882000—40001154000—6000926000—8000548000—100003310000—120001912000—140001514000—160001516000—18000718000—200004总额8422、用不同的图形函数画空间曲面的旋转抛物面六、操作提示1、计算过程（1）建立humps的内置函数：humps(x)=1/sqrt(2*3.14*1.792*x.^2)*exp(-(log(x)-6.012).^2/(2*1.792))运行：x=0.1:0.1:4*3.14;plot(x,humps(x))title(‘对数正态分布密度函数曲线图’)xlabel(‘x=0.1:0.1:4*pi’)(2)x=-3.14:0.2:3.14;y=sin(x);plot(x,y,’+r’)holdony=cos(x);plot(x,y,’-b’)（3）x=-8:0.5:8;y=x’;x=ones(size(y))*x;y=y*ones(size(y))’;R=sqrt(x.^2+y.^2)+eps;z=sin(R)./R;Mesh(z)(其中各语句的意义是：首先建立行向量x，列向量y；然后按向量的长度建立1-矩阵；用向量乘以产生的1-矩阵，生成网格矩阵，它们的值对应与x-y坐标平面；接下来计算各网格点的半径；最后计算函数值矩阵z。用mesh函数即可以得到图形。其中eps表示计算机的最小整数。)（4）ezplot3(‘cos(t)’,’t*sin(t)’,’sqrt(t)’,[0,6*pi])(注：在区间[0，6π]上绘制由参数表达式x=cos(t),y=tsin(t),z=sqrt(t)描述的图形)（5）x=0:2000:18000;y=[48811592543319151574];bar(x,y)x=0:2000:18000;y=[48811592543319151574];bar3(x,y)（6）x=-5:0.5:5;y=x;[x,y]=meshgrid(x,y);z=x.^2+y.^2;subplot(2,2,1)mesh(x,y,z);subplot(2,2,2)meshc(x,y,z);subplot(2,2,3)surf(x,y,z);subplot(2,2,4)surfc(x,y,z);2、计算结果图（2）图（1）图（3）图（4）图（5）直方图：三维直方图：图（6）