matlab如何绘制二维正态曲面

傅里叶变换img=imread('RADU}4W~M9]09V7Q)ZQ5%~7.png');%img=double(img);f=fft2(img); %傅里叶变换f=fftshift(f); %使图像对称r=real(f); %图像频域实部i=imag(f); %图像频域虚部margin=log(abs(f)); %图像幅度谱，加log便于显示phase=log(angle(f)*180/pi); %图像相位谱l=log(f);subplot(2,2,1),imshow(img),title('源图像');subplot(2,2,2),imshow(l,[]),title('图像频谱');subplot(2,2,3),imshow(margin,[]),title('图像幅度谱');subplot(2,2,4),imshow(phase,[]),title('图像相位谱');/s/blog_1667198560102wmzu.html傅里叶变换I = imread('RADU}4W~M9]09V7Q)ZQ5%~7.png'); %读入数字图像I = rgb2gray(I);%将图像进行灰度处理J = fft2(I);%将图像实行傅里叶变换figure,imshow(I);%这里能得到频谱图J = fftshift(J);figure,imshow(log(abs(J)),[]); %将频谱平移J(abs(J)<5)=0;%不必要的过滤掉figure,imshow(log(abs(J)+eps),[]);J = ifftshift(J);K = ifft2(J);figure,imshow(K,[0 255]);%傅里叶逆变换自己所写的代码I = imread('RADU}4W~M9]09V7Q)ZQ5%~7.png'); %读入数字图像J = fft2(I); %将图像实行傅里叶变换figure,imshow(I); %这里能得到频谱图J = fftshift(J);figure,imshow(log(abs(J)),[]); %将频谱平移J(abs(J)<5)=0; %不必要的过滤掉figure,imshow(log(abs(J)+eps),[]);J = ifftshift(J);K = ifft2(J);ss=real(ifft2(J));sss=uint8(ss);subplot(1,2,2);imshow(sss)figure,imshow(K,[0 255]); %傅里叶逆变换A=imread('12.png'); %载入图片A=rgb2gray(A)B=fftshift(fft2(A)); % 进行傅立叶变换subplot(231)imshow(A);title('原始图像');subplot(232)imshow(abs(B),[ ]);title('原始频谱图');subplot(233)imshow(log(abs(B)),[ ]);title('取对数后的频谱图');subplot(234)imshow(angle(B),[ ]);title('相位图');subplot(235)imshow(real(B),[ ]);title('实部图');subplot(236)imshow(imag(B),[ ]);title('虚部图');colormap(jet(64)) %给图片上色图像处理程序处理将vtx变成txt然后加载load '12吧.txt';提取出高程行均值①a=load('z.txt');mean(a)②a=load('z.txt');lzz=length(a);ra=sum(abs(a))/lzz标准差>> a=load('z.txt');sd=std(a)断面深度均方根：a=load('z.txt');lzz=length(a);rq=sqrt(sum(a.^2)/lzz)断面平均斜率b=load('z.txt');a=b';dz=abs(a(2:end)-a(1:end-1));d=load('x.txt');c=d';dy=abs(c(2:end)-c(1:end-1));sl ope_x=dz./dy;极差a=load('z.txt');rangA=max(a)-min(a);偏度a=load('z.txt');pian_du=skewness(a); % 偏度：>0 称为右偏态，<0，称为左偏态驼峰度a=load('z.txt');feng_du=kurtosis(a); % 峰度：用作衡量偏离正态分布的尺度之一平均构造深度（MTD）圆形平常铺沙法（代码实现）a=load('z.txt');zmin=min(a);zmax=max(a);lz=length(a);zzeng=a-zmin;v=sum(zzeng);x=load('x.txt'); y=load('y.txt');d=abs((max(x)-min(x))+ (max(y)-min(y)))/2;mtd=(4*v)/((d^2)*pi)[暂行实现]第一步求面积e=load('x.txt');r=e(30366:30836,1);f=r';q=abs(f(2:end)-f(1:end-1));pingjunx=mean(q);w=load('y.tx t');t=w(30366:30836,1);f=t';p=abs(f(2:end)-f(1:end-1));pingjuny=mean(p);mianji=pingjunx*pingju ny;yy=sqrt((pingjunx)^2+(pingjuny)^2)第二部求mtda=load('z.txt');zmin=min(a);zmax=max(a);lz=length(a);Zzeng=zmax-(a-zmin);zzeng=Zzeng-zmax*0. 07;v=sum(zzeng)*0.002025;x=load('x.txt');y=load('y.txt');d=abs((max(x)-min(x))+(max(y)-min(y)))/2;mtd=(4*v)/((d^2)*pi);其中0.00025是第一步骤中的mianji或者a=load('z.txt');zmin=min(a);zmax=max(a);lz=length(a);zzeng=zmax-(a-zmin);v=sum(zzeng)*0.0007 3875;x=load('x.txt');y=load('y.txt');d=abs((max(x)-min(x))+(max(y)-min(y)))/2;mtd=(4*v)/((d^2)*pi);>> v=100-80.5;d=(20+19.8+20)/3;h=40*v/(pi*d^2)第三部v=100-80.5;d=(20+19.8+20)/3;h=40*v/(pi*d^2)平均断面深度(MPD)A=load('xy.txt');x=A(28369:28854,1);y=A(28369:28854,2);plot(x,y);其中x程序表示选取矩阵A中第1列28369:28854行；y=A(28369:28854,2)表示选取A中第2列28369至28854行。

1 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.8 -1 0 1 2 3 4 5 6 7
y3=cos(t);y4=cos(t+0.25);y5=cos(t+0.5); plot(t,y3);hold on; plot(t,y4); plot(t,y5);
1 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.8 -1 0 1 2 3 4 5 6 7

0
figure(1) title('\fontsize{16}y(\omega)=\int^{\infty }_{0}y(t)e^{-j\omegat}dt')
二、绘制曲线的一般步骤
步骤 1 表 4.1 绘制二维、三维图形的一般步骤 内容 曲线数据准备： 对于二维曲线，横坐标和纵坐标数据变量； 对于三维曲面，矩阵参变量和对应的函数值。 指定图形窗口和子图位置： 默认时，打开 Figure No.1 窗口或当前窗口、当前子图； 也可以打开指定的图形窗口和子图。 设置曲线的绘制方式： 线型、色彩、数据点形。 设置坐标轴： 坐标的范围、刻度和坐标分格线 图形注释： 图名、坐标名、图例、文字说明 着色、明暗、灯光、材质处理(仅对三维图形使用) 视点、三度(横、纵、高)比(仅对三维图形使用) 图形的精细修饰(图形句柄操作)： 利用对象属性值设置； 利用图形窗工具条进行设置。
x=peaks;plot(x) x=1:length(peaks);y=peaks;plot(x,y)
10 8 6 4 2 0 -2 -4 -6 -8
0
5
10
15
20
25
30
35
40
45
50
3. 单窗口多曲线分图绘图 subplot(1,3,1); plot(t,y) subplot(1,3,2); plot(t,y3) subplot(1,3,3); plot(t,y2)

ans = 8 9
[i,j,v]=find(A) 返回矩阵A中非零元素所在的行i,
列j,和元素的值v(按所在位置先后 顺序输出)
A=[3 2 0; -5 0 7; 0 0 1]; [i,j,v]=find(A)
i= 1 2 1 2 3 j= 1 1 2 3 3 v = 3 -5 2 7 1
[X,Y]=meshgrid(x,y) 3)根据函数表达式生成全部网格节点出对应的函数值矩阵z： z＝f(X,Y) 4）顺序连接已经产生的空间点（x,y,z)绘制相应曲面： mesh（X,Y,Z) surf（X,Y,Z) shading flat ％去除网格线。
例2-7画出矩形域[-1,1]×[-1,1]旋转抛物面:z=x2+y2. x=linspace(-1,1,100); y=x; [X,Y]=meshgrid(x,y); %生成矩形区[-1,1]×[-1,1]的网格坐标矩阵 Z=X.^2+Y.^2; subplot(1,2,1) mesh(X,Y,Z); subplot(1,2,2) surf(X,Y,Z); shading flat; %对曲面z=x2现方式做保护处理对用户上传分享的文档内容本身不做任何修改或编辑并不能对任何下载内容负责
用matlab绘制二维、三维图形
2.1二维图形的绘制
2.1.1 二维绘图的基本命令 matlab中，最常用的二维绘图命令是plot。
使用该命令，软件将开辟一个图形窗口，并 画出连接坐标面上一系列点的连线。
例2-5 采用不同形式（直角坐标、参数、极坐标），画出 单位圆x2＋y2＝1的图形。
分析：对于直角坐标系方程，y＝ 1 x2，对于参数方 程x＝cost,y=sint,t[0,2 pi] ,利用plot(x,y)命令可以实现。 而在极坐标系中单位圆为r＝1(1+0t),利用polar(t,r)命 令实现。

第2章MA TLAB二维绘图 (2)2.1 二维绘图基本流程 (2)2.2 二维图形的基本绘图命令 (4)2.2.1 高级绘图命令 (4)2.2.2低级绘图命令 (6)2.2 二维图形的修饰 (8)2.2.1 坐标轴的调整 (8)2.2.1.1 调整坐标轴的范围 (8)2.2.1.2 调整坐标轴的状态 (9)2.2.1.3 保存坐标轴的范围 (11)2.2.1.4 保存坐标轴的状态 (11)2.2.2画出或取消网格线 (12)2.2.3设置坐标轴的名称 (12)2.2.4设置图形标题 (13)2.2.5在图形中显示文字 (14)2.2.5.1用坐标轴确定文字位置 (14)2.2.5.2用鼠标确定位置显示文字 (15)2.2.6 图形的标定和颜色条 (16)2.2.7 使用绘图工具栏标注图形 (18)2.3 填充图形的绘制 (19)2.4 多坐标系绘图与图形窗口的分割 (20)2.4.1 图形叠印法 (20)2.4.2 子图的绘制 (21)2.5 特殊坐标图形的绘制 (22)2.5.1 绘制极坐标图形 (23)2.5.2对数/半对数坐标系绘图 (23)2.6 特殊二维图形的绘制 (24)2.4.3 直方图 (24)2.4.1 柱状图和面积图 (26)2.4.2 饼图 (28)2.4.4 离散数据绘图 (28)2.4.5 等高线图 (30)2.4.6 向量图 (31)2.7 函数绘图 (34)2.7.1 fplot函数 (34)2.7.2函数function的定义 (35)2.8 工作空间直接绘图 (36)2.9 手工绘图方式 (38)2.10 小结 (41)第2章MATLAB二维绘图数据可视化是MATLAB一项重要功能，它所提供的丰富绘图功能，使得从繁琐的绘图细节中脱离出来，而能够专心于最关心的本质。

通过数据可视化的方法，工程科研人员可以对自己的样本数据的分布、趋势特性有一个直观的了解。

本章将重点介绍MA TLAB二维图形的绘制方式，并按照完整的步骤来说明一个图形产生的流程，以便将数据以图形形式来识别。

实验四
专业：电子信息工程2班姓名：李书杰学号：3121003210
一、实验目的
1.掌握绘制二维图形及三维图形的方法。

2.掌握图形控制与修饰处理的方法。

3.了解图像处理及动画制作的基本方法。

二、实验内容
1.绘制下列图形曲线。

（1）y=x-x^3/3! (2)x^2+2Y^2=64
解：程序如下
2.设y=1/(1+e^-t),-pi<=t<=pi,在同一个图形窗口中采用子图的形式绘制条形图、阶梯图、杆图和对数坐标等不同图形，并对不同图形加标注说明。

解:程序如下
3.绘制下列极坐标图。

（1）ρ=5cosθ+4 (2)γ=5sin^2φ/cosφ,-π/3<φ<π/3 解：程序如下
思考练习：
2.绘制下列曲线
（1）y=1/2πe^(-x^2/2) (2)x=tsint y=tcost
解：程序如下
（1）
结果如下：
（2）
结果如下：
3.在同一坐标中绘制下列两条曲线并标注两曲线交叉点。

（1）y=2x-0.5
(2)x=sin(3t)cost
Y=sin(3t)sint
解：程序如下
4.分别用plot和fplot函数绘制y=sin(1/x)的曲线，分析两曲线的差别。

解：程序如下
结果如下：
5.绘制下列极坐标图：
(1）p=12/sqrt(θ) (2)γ=3asinφcosφ/(sin^3φ+cos^3φ)解：程序如下
结果如下：。

MATLAB二维曲线绘图 二维曲线绘图 • 一、直角坐标系下的绘图 • 二、隐函数与参数方程的绘图 • 三、极坐标系下的绘图
一、在直角坐标系下的二维绘图
• 格式一：plot(x) 格式一： • 在这种情况下 ， 当 x是实向量时 ， 以该向量元素的 在这种情况下， 是实向量时， 是实向量时 下标为横坐标， 下标为横坐标，元素值为纵坐标画出一条连续曲线 ，这实际上是绘制折线图。 这实际上是绘制折线图。 • 格式二： plot(x,y) 格式二： • 作出以数据(xi,yi)为节点的折线图，其中 和y为维 为节点的折线图， 作出以数据 为节点的折线图 其中x和 为维 数相同的向量，分别用于存储x坐标和 坐标数据。 数相同的向量，分别用于存储 坐标和y坐标数据。 坐标和 坐标数据 思考： 的区别？ 思考：plot(x)，plot(x,y)的区别？ ， 的区别
• 格式三：plot(x1,y1,x2,y2,…,xn,yn) 格式三： • 当输入参数都为向量时，x1和y1，x2和y2，…， 当输入参数都为向量时， 和 ， 和 ， ， xn和yn分别组成一组向量对，每一组向量对的长 分别组成一组向量对， 和 分别组成一组向量对 度可以不同。每一向量对可以绘制出一条曲线， 度可以不同。每一向量对可以绘制出一条曲线， 这样可以在同一坐标内绘制出多条曲线。 这样可以在同一坐标内绘制出多条曲线。 • 例如：1、绘制行向量( 1 8 6 5 4)的图形。 例如： 绘制行向量( 4)的图形 的图形。 • 已知点列(xi,yi)坐标如下 点列(xi,yi)坐标如下： 2、已知点列(xi,yi)坐标如下： x=(0,1,2,3,4,5），y=(3,1,4,2,6,-1),试将该点 ），y=(3,1,4,2,6, x=(0,1,2,3,4,5），y=(3,1,4,2,6,-1),试将该点 列连结成折线。 列连结成折线。

matlab二维图形的绘制(2006-11-20 20:38:35)转载▼分类：matlab基础（电子方向）常用的二维图形命令:plot：绘制二维图形loglog：用全对数坐标绘图semilogx：用半对数坐标（X）绘图semilogy：用半对数坐标（Y）绘图fill：绘制二维多边填充图形polar：绘极坐标图bar：画条形图stem：画离散序列数据图stairs：画阶梯图errorbar：画误差条形图hist：画直方图fplot：画函数图title：为图形加标题xlabel：在X轴下做文本标记ylabel：在Y轴下做文本标记zlabel：在Z轴下做文本标记text：文本注释grid：对二维三维图形加格栅绘制单根二维曲线plot函数，基本调用格式为：plot(x,y)其中x和y为长度相同的向量，分别用于存储x坐标和y坐标数据。

例如：在区间内，绘制曲线y=2e-0.5xcos(4πx)程序如下：x=0:pi/100:2*pi;y=2*exp(-0.5*x).*cos(4*pi*x);plot(x,y)plot函数最简单的调用格式是只包含一个输入参数：plot(x)在这种情况下，当x是实向量时，以该向量元素的下标为横坐标，元素值为纵坐标画出一条连续曲线，这实际上是绘制折线图。

p=[22,60,88,95,56,23,9,10,14,81,56,23];plot(p)绘制多根二维曲线1．plot函数的输入参数是矩阵形式(1) 当x是向量，y是有一维与x同维的矩阵时，则绘制出多根不同颜色的曲线。

曲线条数等于y矩阵的另一维数，x被作为这些曲线共同的横坐标。

(2) 当x,y是同维矩阵时，则以x,y对应列元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。

(3) 对只包含一个输入参数的plot函数，当输入参数是实矩阵时，则按列绘制每列元素值相对其下标的曲线，曲线条数等于输入参数矩阵的列数。

当输入参数是复数矩阵时，则按列分别以元素实部和虚部为横、纵坐标绘制多条曲线。

标记符号选项 选 v ^ < > p （pentagram ） h （hexagram ） 项 标 记 符 号 朝下三角符号 朝上三角符号 朝左三角符号 朝右三角符号 五角星符 六角星符
例 在同一坐标内，分别用不同线型和颜色绘制曲线 y1 = 0.2e−0.5xcos(4x)和y2 = 1.5e−0.5x cos(x)。标记两曲 线交叉点。 x=linspace(0,2*pi,1000); y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=1.5*exp(-0.5*x).*cos(pi*x); k=find(abs(y1-y2)<1e-2); x1=x(k); y3=0.2*exp(-0.5*x1).*cos(4*pi*x1); plot(x,y1,x,y2,'k:',x1,y3,'bp');
MATLAB提供了一些绘图选项，用于确定所绘曲线的线型、 颜色和数据点标记符号。 例如，“b-.”表示蓝色点画线，“y:d”表示黄色虚线并用菱 形符标记数据点。当选项省略时，MATLAB规定，线型一 律用实线，颜色将根据曲线的先后顺序依次采用表3.2给 出的前7种颜色。
表 3.1 线型选项 选 项 : --. 线 型 实线（默认值） 虚线 双画线 点画线
【例 3.10 】表 3.5 所示为某公司 3 类产品各季度的销售额（单位：万元） ，分别按季度绘制簇 状柱形图和堆积条形图。
表 3.5 第 一 季 度 产品 A 产品 B 产品 C 51 67 78 产品全年销售额（单位：万元） 第 二 季 度 82 78 85 第 三 季 度 34 68 65 第 四 季 度 47 90 50
在绘制图形的同时，可以对图形加上一些说明，如图形名 称、坐标轴说明、图形某一部分的含义等，这些操作称为 添加图形标注。有关图形标注函数的调用格式如下： title（图形名称） xlabel（x轴说明） ylabel（y轴说明）

二维数组操作函数
repmat(A,m,n) rot90(A,k)
按指定的“行数、列数”铺放模块 数组，以形成更大的数组 逆时针旋转 k×90 度 组的“行数、列数”
reshape(A,m,n) 在总元素数不便的前提下，改变数 tril(A) triu(A)
提取一个矩阵的下三角部分 提取一个矩阵的上三角部分
π y = sin 2t + 6
plot(x,y)
一维数组, （2）若x 是一维数组 y 是二维数组 ） 的行数相等， 当 x 的长度与 y 的行数相等，则将 x 与 y 中的各列相对应，绘制多条平面曲线； 中的各列相对应，绘制多条平面曲线； 否则， 的列数相等， 否则，若 x 的长度与 y 的列数相等，则将 x 与 y 中的各行相对应，绘制多条平面曲线。 中的各行相对应，绘制多条平面曲线。
7、二维数组的与标量 的运算
• • • • • • • A + c ：A的每个元素加c A - c： A的每个元素减c A.*c：点乘， A的每个元素乘c A./c：右点除， A的每个元素除c A.\c：左点除，c除A的每个元素 A.^c：点幂， A的每个元素做幂运算 c.^A：点幂，c做幂运算
8、函数作用在二维数组上的 运算规则
计算A中每个列向量的元素的和 计算A中每个行向量的元素的和
A中每个列向量的累加和，维数与A相同 A中每个行向量的累加和，维数与A相同
计算A中每个列向量的元素的积 计算A中每个行向量的元素的积
A中每个列向量的累乘积，维数与A相同 A中每个行向量的累乘积，维数与A相同
sort(A)
对A中列向量进行升序排序
由A的各列按自左到右的次序，首尾相接而生成 的“一维长列”数组 引用A中由一维数组s指定的元素。s若是行（或 列），则A(s)就是长度相同的行（或列） 全元素赋值方式。结果：保持A的“行宽、列长” 不变。条件：A、D两个数组的总元素相等，但 “行宽、列长”不一定相同 对A的部分元素重新赋值。结果：保持A的“行 宽、列长”不变。条件：s单下标数组的长度必 须与“一维数组”B的长度相等，但是s、B不一 定同是“行数组”或“列数组”

实验四 MATLAB 二维、三维图形的绘制一 实验目的1 掌握二维、三维图形的绘制；2 掌握特殊二维图形的绘制；3 掌握绘图参数的设置；4 了解并学习简单动画的制作。

二 实验内容1 在0-2π区间上画sin(x)和cos(x)，要求在同一个图像中，其中cos(x)图像用红色小圆圈表示，并在函数图上标注“y=sin(x)”，“y=cos(x)”，坐标轴标签为“x 轴”，“y 轴”，标题为“正弦余弦函数图像”。

2 绘制函数x 2/32+y 2/42=1的边界。

3 绘制三维曲线⎪⎩⎪⎨⎧=+=+=u z y x u u 3)sin 21()cos 21(，]10,0[∈u 。

4 使用极坐标绘制]2,0[,2sin πθθρ∈=。

5 绘制函数]2,2[,)cos()sin(21122121-∈⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛x x x x x x y y 在上的曲线，数据点用菱形表示，再绘制其对应的等高线。

6 在同一坐标内，分别用不同线型和颜色绘制曲线)4c o s (5.012.0x x e y π-=和)cos(5.022x x e y π-=，标记两曲线交叉点，]2,0[π∈x 。

7 在同一张图中用子图的方式分别以条形图、阶梯图、杆图和填充图形式绘制曲线y=2sin(x)。

8 连续函数的可视化：用图形表示调制波形y=sin(t)sin(9t)，变量范围以及步长、曲线表示方法均有自己设定，结果图与下图相似，表达意思相同即可。

9 绘制三维曲线图：x=sin(t), y=cos(t), z=cos(2t)，参考图例如下，学习使用view 和box函数。

10 用曲面画图表示z=x2+y2，参考图例如下。

案例28使用MATLAB 绘制专业图形绘制曲线和曲面图案例主要信息提示 ● 案例内容：使用MATLAB 进行专业作图。

● 关键词：MATLAB ，作图，二维作图，三维作图 ● 建议课时：2课时● 适合专业：理工科各专业●光盘中的数字资源：二维曲线绘图M 文件graph_line.m三维曲面绘图M 文件graph_surface.m一、 实验内容实验内容1.绘制曲线图在1个绘图窗口的4个不同子窗口中分别绘制以下曲线： ①在子窗口1中绘制两条二维曲线，分别为：0.510.2cos(4)x y e x π-=和0.522cos()x y e x π-=并显示网格线。

②在子窗口2中绘制两条二维曲线，分别为：0.512x y e -=和2cos(2)y x π=并添加标题、x/y 轴名称、图形说明和图例。

③在子窗口3中绘制一条极坐标曲线，如下：sin()cos()r t t =④在子窗口4中绘制一条极三维螺旋线。

全部4个子窗口绘制效果如图28-1所示。

图28-1曲线绘制效果图实验内容2.绘制曲面图在1个绘图窗口的4个不同子窗口中分别绘制以下曲面：①在子窗口1和子窗口2中以不同着色方式绘制两个球面。

②在子窗口3和子窗口4中以不同视角绘制两个直径相等的圆管相交的图形。

全部4个子窗口绘制效果如图28-2所示。

图28-1曲面绘制效果图二、预备知识2.1 MATLAB绘图功能简介MA TLAB软件提供了丰富的图形表达功能，包括常用的二维图形和三维图形。

其中，各种二维图形近30种，三维图形20多种。

应用MA TLAB除了能作一般的曲线、散点图、条形图等，还能绘制流线图、三维矢量图等工程实用图形。

下面我们介绍一些基本的二维和三维绘图函数。

2.2 二维图形的绘制二维图形是将平面坐标上的数据点连接起来的平面图形。

可以采用不同的坐标系，如直角坐标、对数坐标、极坐标等。

二维图形的绘制是其他绘图操作的基础。

2.2.1绘制二维曲线的plot函数在MATLAB中，plot函数是最简单、最基本而且应用最为广泛的线性绘图函数，利用它可以在二维平面上生成线段、曲线和参数方程曲线等的函数图形。

[转载]⽤MATLAB画曲⾯图原⽂地址：⽤MATLAB画曲⾯图作者：惠州学院数学建模曲⾯图，除了各线条之间的空档（称作补⽚）⽤颜⾊填充以外，和⽹格图看起来是⼀样的。

这种图⼀般使⽤函数 surf来绘制。

⾃然，函数surf使⽤和函数mesh相同的调⽤语法。

⽐如：[X,Y,Z]=peaks(30);surf(X,Y,Z)grid,xlabel( ‘ x-axis ‘ ),ylabel( ‘ y-axis ‘ ),zlabel( ‘ z-axis ‘ )title( ‘ SURF of PEAKS ‘ )MATLAB 还提供了平滑加颜⾊和插值加颜⾊功能。

这可以通过调⽤函数shading 来实现[X,Y,Z]=peaks(30);surf(X,Y,Z) % same plot as abovegrid,xlabel( ‘ x-axis ‘ ),ylabel( ‘ y-axis ‘ ),zlabel( ‘ z-axis ‘ )title( ‘ SURF of PEAKS ‘ )shading flat如上所⽰平滑加⾊彩的例⼦中，每⼀补⽚仍保存着单⼀的颜⾊，但各块连接处的⿊线已去掉。

shading interp输出见下图如上所⽰内插加⾊彩的例⼦中，同样去掉了线条，但各补⽚以插值加颜⾊，即各补⽚的颜⾊根据赋予顶点的⾊值，对其区间进⾏了插值计算。

很明显，插值⾊彩需要⽐分块和平滑更多的计算量。

在⼀些计算机系统中，插值⾊彩会产⽣⾮常长的打印延时或打印错误。

这问题不在于PostScript ⽂件太⼤，⽽是由于在打印机上产⽣沿图形曲⾯连续变化的阴影所需的巨⼤计算量。

通常对这个问题最简单的解决⽅法是使⽤平滑加⾊彩法来打印。

⾊彩对 surf作图的视觉效果有着巨⼤的影响。

对⽹格图也是如此，尽管由于只有线条有颜⾊，对视觉效果的影响相对要⼩⼀些。

因为曲⾯图不能作成透明，但在⼀些情况下可以很⽅便地移⾛⼀部分表⾯以便看到表⾯以下部分，在 MATLAB 中，这是通过在所期望的洞孔的所在位置，将数据置为特定的NaN 来实现。

实验二matlab图形绘制一、实验目的1、学习MATLAB图形绘制的基本方法；2、熟悉和了解MATLAB图形绘制程序编辑的基本指令；3、熟悉掌握利用MATLAB图形编辑窗口编辑和修改图形界面，并添加图形的各种标注；二、实验原理1．二维数据曲线图（1）绘制单根二维曲线plot(x,y);（2）绘制多根二维曲线plot(x,y) 当x是向量，y是有一维与x同维的矩阵时，则绘制多根不同颜色的曲线。

当x，y是同维矩阵时，则以x，y对应列元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。

（3）含有多个输入参数的plot函数plot(x1,y1,x2,y2,…,xn,yn)（4）具有两个纵坐标标度的图形plotyy(x1,y1,x2,y2)2．图形标注与坐标控制1）title （图形名称）2）xlabel（x轴说明）3）ylabel（y轴说明）4）text（x，y图形说明）5）legend（图例1，图例2，…）6）axis （[xmin xmax ymin ymax zmin zmax]）3．图形窗口的分割subplot（m,n,p）4．三维曲线plot3（x1,y1,z1,选项1，x2,y2，选项2，…,xn,yn,zn,选项n）5．三维曲面mesh(x,y,z,c) 与surf(x,y,z,c)。

一般情况下，x ，y ，z 是维数相同的矩阵。

X ，y 是网格坐标矩阵，z 是网格点上的高度矩阵，c 用于指定在不同高度下的颜色范围。

三、实验内容及步骤1．绘制下列曲线： (1) 21100x y +=x=0:0.02:10; y=100./(1+x.^2); plot(x,y)title('my first plot'); xlabel('x'); ylabel('y'); grid on截图：(2) 2221x e y -=πx=0:0.02:10;y=1./(2*pi).*exp(-(x.^2)./2); plot(x,y)title('my first plot'); xlabel('x'); ylabel('y'); grid on截图：(3) 122=+y x 6t=-4:0.02:4; y=4*sin(t); x=4*cos(t);plot(x,y)title('my first plot'); xlabel('x'); ylabel('y'); grid on截图：(4) ⎩⎨⎧==325t y t x t=0:0.02:10; x=t.^2; y=5.*t.^3; plot(x,y)title('my first plot'); xlabel('x');ylabel('y');grid on截图：2．在一个图形窗口绘制正弦和余弦曲线，要求给图形加标题“正弦和余弦曲线”，X轴Y轴分别标注为“时间t”和“正弦、余弦”，添加图例，在图形的某个位置标注“sin(t)”“cos(t)”，显示网格。

Matlab⼊门教程-⼆维绘图[z] MATLAB不但擅长於矩阵相关的数值运算，也适合⽤在各种科学⽬视表⽰（Scientific visualization）。

本节将介绍MATLAB基本xy平⾯及xyz空间的各项绘图命令，包含⼀维曲线及⼆维曲⾯的绘制、列印及存档。

plot是绘制⼀维曲线的基本函数，但在使⽤此函数之前，我们需先定义曲线上每⼀点的x及y座标。

下例可画出⼀条正弦曲线：close all; x=linspace(0, 2*pi, 100); % 100个点的x座标y=sin(x); % 对应的y座标plot(x,y);====================================================⼩整理：MATLAB基本绘图函数plot: x轴和y轴均为线性刻度（Linear scale）loglog: x轴和y轴均为对数刻度（Logarithmic scale）semilogx: x轴为对数刻度，y轴为线性刻度semilogy: x轴为线性刻度，y轴为对数刻度====================================================若要画出多条曲线，只需将座标对依次放⼊plot函数即可：plot(x, sin(x), x, cos(x));若要改变颜⾊，在座标对后⾯加上相关字串即可：plot(x, sin(x), 'c', x, cos(x), 'g');若要同时改变颜⾊及图线型态（Line style），也是在座标对后⾯加上相关字串即可：plot(x, sin(x), 'co', x, cos(x), 'g*');====================================================⼩整理：plot绘图函数的叁数字元颜⾊字元图线型态y 黄⾊ . 点k ⿊⾊ o 圆w ⽩⾊ x xb 蓝⾊ + +g 绿⾊ * *r 红⾊ - 实线c 亮青⾊ : 点线m 锰紫⾊ -. 点虚线-- 虚线====================================================图形完成后，我们可⽤axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围：axis([0, 6, -1.2, 1.2]);此外，MATLAB也可对图形加上各种注解与处理：xlabel('Input Value'); % x轴注解ylabel('Function Value'); % y轴注解title('Two Trigonometric Functions'); % 图形标题legend('y = sin(x)','y = cos(x)'); % 图形注解grid on; % 显⽰格线我们可⽤subplot来同时画出数个⼩图形於同⼀个视窗之中：subplot(2,2,1); plot(x, sin(x));subplot(2,2,2); plot(x, cos(x));subplot(2,2,3); plot(x, sinh(x));subplot(2,2,4); plot(x, cosh(x));MATLAB还有其他各种⼆维绘图函数，以适合不同的应⽤，详见下表。

MATLAB 的二维绘图基础了解了MATLAB 的矩阵和向量概念与输入方法之后，MATLAB 的二维绘图再简单也不过了。

假设有两个同长度的向量 x 和y, 则用plot(x,y) 就可以自动绘制画出二维图来。

如果打开过图形窗口，则在最近打开的图形窗口上绘制此图，如果未打开窗口，则开一个新的窗口绘图。

〖例〗正弦曲线绘制：>> t=0:.1:2*pi;%生成横坐标向量，使其为0,0.1,0.2,...,6.2y=sin(t); % 计算正弦向量plot(t,y) %绘制图形这样立即可以得出如图所示的二维图[4.1(a)]plot() 函数还可以同时绘制出多条曲线，其调用格式和前面不完全一致，但也好理解。

>> y1=cos(t); plot(t,y,t,y1); %或plot(t,[y; y1]), 即输出为两个行向量组成的矩阵。

图形见 4.1(b)。

plot() 函数最完整的调用格式为：>> plot(x1,y1,选项1, x2,y2, 选项2, x3,y3, 选项3, ...)其中所有的选项如表 4.1 所示。

一些选项可以连用，如'-r' 表示红色实线。

由MATLAB 绘制的二维图形可以由下面的一些命令简单地修饰。

如>> xlabel('字符串') % 给横坐标轴加说明>> ylabel('字符串') % 给纵坐标轴加说明，%并自动旋转90度>> title('字符串') % 给整个图形加图题得出的图形如右图所示。

axis() 函数可以手动地设置x,y 坐标轴范围还可以使用plotyy() 函数绘制具有两个纵坐标刻度的图形。

坐标系的分割在MATLAB 图形绘制中是很有特色的，比较规则的分割方式是用subplot() 函数定义的，其标准调用格式为subplot(n,m,k)其中，n 和m 为将图形窗口分成的行数和列数，而k 为相对的编号。

.实验报告（201 /201 学年第学期）课程名称实验名称二维图形与三维图形的绘制实验时间年月日实验室指导教师学生姓名学号班级专业实验报告三、实验内容及原理（包括硬件原理图、算法、逻辑框图，关键代码等，可续页）（一）二维图形的绘制1、绘制二维曲线的基本函数：○1plot函数plot函数的基本调用格式为：plot(x,y);其中x和y为长度相同的向量，分别用于存储x坐标和y坐标数据。

例，绘制参数方程曲线。

程序如下：含多个输入参数的plot函数调用格式为：p lot(x1,y1,x2,y2,…,xn,yn);含选项的plot函数调用格式为：plot(x1,y1,选项1,x2,y2,选项2,…,xn,yn,选项n);例，用不同线型和颜色在同一坐标内绘制曲线y=2e-0.5x sin(2πx)及其包络线。

程序如下：○2双纵坐标函数plotyyplotyy函数是MATLAB 5.X新增的函数。

它能把函数值具有不同量纲、不同数量级的两个函数绘制在同一坐标中。

调用格式为：plotyy(x1,y1,x2,y2);其中x1-y1对应一条直线，x2-y2对应另一条曲线。

横坐标的标度相同，纵坐标有两个，左纵坐标用于x1-y1数据对，右纵坐标用于x2-y2数据对。

2、绘制二维图形的其他函数在线性直角坐标系中，其他形式的图形有条形图、阶梯图、杆图和填充图等，所采用的函数分别是：bar(x,y,选项);stairs(x,y,选项);stem(x,y,选项);fill(x1,y1,选项1,x2,y2,选项2,…);例，分别以条形图、填充图、阶梯图和杆图形式绘制曲线y=2e-0.5x。

程序如下：（二）三维图形的绘制1、绘制三维曲线的基本函数plot3函数与plot函数用法十分相似，其调用格式为：plot3(x1,y1,z1,选项1,x2,y2,z2,选项2,…,xn,yn,zn,选项n); 例，绘制空间曲线。

程序如下：2、绘制三维曲面的函数surf函数和mesh函数的调用格式为：surf(x,y,z,e);mesh(x,y,z,e);例，绘制两个直径相等的圆管的相交图形。

应用MATLAB绘制二次曲面图1、用surf工mesh函数绘图Surf函数绘制的是三维表面图，mesh函数绘制的是三维网格图，当二次曲面方程是标准方程时，原方程式可化为z=f(x,y),x=f(y,z),y=f(z,x)时，我们就用这两种函数完成绘图。

x2y2z2x2y2z2x2y2z++=1②+-=1③+=在区域例1、绘曲面①49164916494-2≤x≤2,-3≤y≤3,-4≤z≤4上的图像。

4y2x2y2x2y22解：以上三个方程化为：z=±4-、z=±4； -+-1、z=x+949492、用plot3或contour3函数绘图plot3函数绘制的是三维直角坐标曲线图，contour3函数绘制的是三维等高曲线图。

x=-2:0.1:2;y=-3:0.1:3;[x,y]=meshgrid(x,y);z1=4.*sqrt(1-(x.^2)./4-(y.^2)./9);z2=-4.*sqrt(1-(x.^2)./4-(y.^2)./9);subplot(2,3,1);plot3(x,y,z1);hold on;plot3(x,y,z2)grid on3、用ezsurf或ezmesh函数绘图Ezsurf函数和ezmesh函数主要针对参数方程的三维作图函数，它们是专业作图函数，ezsurf绘制三维表面图，ezmesh绘制三维网格图，当二次曲面可化为参数方程时，就可以用这两种函数完成绘图。

⎧x=2cosαcosβ⎪椭球方程的参数方程为：⎨y=3sinαcosβ （0≤α≤2*pi,⎪z=4sinβ⎩⎧x=2tcosα⎪双曲方程的参数方程为：⎨y=3tsinα （0≤α≤2*pi,⎪z=±4t2-1⎩-pipi≤β≤）22t≥1或t≤1）⎧x=2tcosα⎪抛物面方程的参数方程为：⎨y=3tsinα （0≤α≤2*pi,⎪z=4t2⎩syms t1 t2;x=2*cos(t1)*cos(t2);y=3*sin(t1)*cos(t2);z=4*sin(t2);ezmesh(x,y,z,[0,2*pi],[-pi/2,pi/2])所以，把二次曲面的参数方程附在下面： -∞<t<+∞）⎧x=acosϕcosθ⎪球面：x2+y2+z2=a参数方程⎨y=asinϕsinθ⎪z=acosϕ⎩⎧x=asinϕcosθxyz⎪椭球面：2+2+2=1参数方程⎨y=bsinϕsinθ abc⎪z=ccosϕ⎩222 y2z2旋转椭球面：2+2=1或x=0,y=acost,z=bsintab旋转椭球面的参数方程为：0≤t≤2π绕z 轴旋转一周的⎧x=x2(t)+y2(t)coθs⎧x=acostcosθ⎪⎪⎪22θ 即⎨y=acostsinθ ⎨y=x(t)+y(t)sin⎪z=z(t)⎪z=bsint⎪⎩⎩ezmesh('cos(t)*cos(theta)','cos(t)*sin(theta)','2*sin(t)',[0,2*pi],[0,2*pi]) ⎧x=acosθxy⎪椭圆柱面：2+2=1的参数方程：⎨y=bsinθ ab⎪z=t⎩22⎧x=acosθ⎪圆柱面:x2+y2=a2的参数方程⎨y=asinθ⎪z=t⎩⎧x=ucost⎪22旋转抛物面z=x+y的参数方程：⎨y=usint⎪z=u2⎩x2y2椭圆抛物面z=2+2的参数方程：ab⎧x=aucost⎪22⎨y=busint或z=ax+by的参数方程⎪z=u2⎩ucots⎧x=⎪a⎪usint⎪⎨y=⎪2⎪z=u⎪⎩⎧x=achθ⎧x=asecθxy⎪⎪双曲柱面2-2=1的参数方程：⎨y=bshθ或⎨y=btanθab⎪z=t⎪z=u⎩⎩22单叶双曲面：直线x=1,y=t,z=2t绕z轴旋转得到的单叶双曲面的参数方程为：⎧x=x2(t)+y2(t)cosθ⎧x=+t2cosθ⎪⎪⎪⎪222⎨y=x(t)+y(t)sinθ 即⎨y=+tsinθ⎪z=z(t)⎪z=2t⎪⎪⎩⎩λ⎧x=a(+t)⎪222⎪λxy⎪双曲抛物面：z=2-2的参数方程⎨y=b(-t) 2ab⎪⎪z=2λt⎪⎩⎧x=a(u+v)xy⎪或2z=2-2的参数方程⎨y=b(u-v) ab⎪z=2uv⎩22⎧x=v⎪或z=xy的参数方程⎨y=u⎪z=uv⎩⎧x=2pu2⎧x=2u2x⎪⎪2抛物柱面:x=y的参数方程⎨y=2pu或z=的参数方程⎨y=v 2⎪z=2u2⎪z=v⎩⎩⎧x=v⎪2 或z=y的参数方程⎨y=u⎪z=u2⎩⎧x=ucost⎪圆锥面：x2+y2=z2的参数方程⎨y=usint ⎪z=u⎩⎧x=aucosθxyz⎪椭圆锥面2+2=2的参数方程⎨y=businθ abc⎪z=cu⎩222θ+ϕ)-cosθ]}⎧x=a{cosθ+t[cos(xyz⎪单叶双曲面2+2-2=1的参数方程⎨y=b{sinθ+t[sin(θ+ϕ)-sinθ]} abc⎪z=c(2t-1)⎩222⎧x=a+u2cosθ⎧x=asecucosv⎧x=acoshucosv⎪⎪⎪⎪2 或⎨y=b+usinθ或⎨y=bsecusinv 或⎨y=bcoshusinv ⎪z=ctanu⎪z=csinhu⎪z=cu⎩⎩⎪⎩⎧x=au2-1cosθ⎧x=atanucosv⎪x2y2z2⎪⎪2双叶双曲面2+2-2=-1的参数方程⎨y=bu-1sinθ或⎨y=btanusinv abc⎪z=csecu⎪z=cu⎩⎪⎩⎧x=asinhucosv⎪或⎨y=bsinhusinv⎪z=ccoshu⎩ezmesh('sqrt(u^2-1)*cos(theta)','sqrt(u^2-1)*sin(theta)','u')⎧x=(R+rcost)sinθ⎪圆环面的参数方程：⎨y=(R+rcost)cosθ⎪z=rsint⎩v⎧x=(2+ucos)cosv⎪2⎪v⎪莫比乌斯带:⎨y=(2+ucos)sinv 2⎪⎪z=usinv⎪2⎩⎧x=ucosv⎪螺旋面的参数方程:⎨y=usinv⎪z=v⎩。