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行。

第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二维图形的绘制(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函数，当输入参数是实矩阵时，则按列绘制每列元素值相对其下标的曲线，曲线条数等于输入参数矩阵的列数。

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

实验四 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，参考图例如下。