Matlab编程与作图实验报告

实验二 MATLAB 编程与作图

1、设x 为一个长度为n 的数组，编程求下列均值和标准差

11,1n i i x x s n n ===>∑

并以[81,70,65,51,76,66,90,87,77]x =为例计算。 解： function [m,s]=myf_1(x) n=length(x); m=sum(x)/n;

s=sqrt((sum(x.^2)-n*m^2)/(n-1)); >>x=[81,70,65,51,76,66,90,87,61,77]; >> [m,s]=myf_1(x) m = 72.4000 s =

12.1124

2、求满足0

ln(1)100

m

n n =+>∑的最小m 值。

解：>> s=log(1);n=0; >> while s<=100 n=n+1;

s=s+log(1+n); end >> m=n m =

37

3、用循环语句形成Fibonacci 数列

12121,,3,4,k k k F F F F F k --===+=L

，并验证极限

1

lim

k k k F F →∞-=

。 解：>> clear;

>> F(1)=1;F(2)=1;k=2;x=0; >> e=1e-8;a=(1+sqrt(5))/2; >> while abs(x-a)>e

k=k+1;F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1); end >> a,x,k a =

1.6180 x =

1.6180 k =

21

4、分别用for 和while

循环结构编写程序，求出

6

10

1

2i

i

K

=

=

，并考虑一种避免循环语句的

程序设计，比较不同算法的运行时间。

解：>> clear;tic;s=0;

for i=1:1000000

s=s+sqrt(3)/2^i;

end

s,toc

tic;s=0;i=1;

while i<=1000000

s=s+sqrt(3)/2^i;i=i+1;

end

s,toc

tic;s=0;

i=1:1000000;

s=sqrt(3)*sum(1./2.^i);

s,toc

s =

1.7321

Elapsed time is 0.944317 seconds.

s =

1.7321

Elapsed time is 1.326814 seconds.

s =

1.7321

Elapsed time is 0.461718 seconds.

5

⑵用Matlab指令dlmwrite将上述数据输出到一个文本文件中，第一列是时刻，第二列是温度，要求用空格分割数据；

⑶从工具条File选Import data 导入上述数据文件中的数据。

解：>> t=0:24;

c=[15 14 14 14 14 15 16 18 20 22 23 25 28 ...

31 32 31 29 27 25 24 22 20 18 17 16];

plot(t,c)

14 16 18 20 22 24 26 28 30 32

6、作出下列函数图像：

⑴曲线

22

sin(2),22

y x x x x

=---≤≤（要求分别使用plot和fplot完成）；

解：x=-2:0.1:2;y=x.^2.*sin(x.^2-x-2);plot(x,y) y=inline('x^2*sin(x^2-x-2)');fplot(y,[-2 2])

⑵椭圆

22

1 49

x y

+=

；

解：t=linspace(0,2*pi,100);

x=2*cos(t);y=3*sin(t); plot(x,y)

⑶抛物面

22,3,3 z x y x y

=+<<

；

解：x=-3:0.1:3;y=x; [x,y]=meshgrid(x,y); z=x.^2+y.^2;

surf(x,y,z)

⑷曲面

422232226,3,313

z x x y x y x y x y =++---+<-<<；

解：x=-3:0.1:3;y=-3:0.1:13; [x,y]=meshgrid(x,y);

z=x.^4+3*x.^2+y.^2-2*x-2*y-2*x.^2.*y+6; surf(x,y,z)

⑸空间曲线sin ,cos ,cos(2),02x t y t z t t π===<<； 解：t=0:0.01:2*pi;

x=sin(t);y=cos(t);z=cos(2*t); plot3(x,y,z)

⑹半球面2sin cos ,2sin sin ,2cos ,02,02x y z ?θ?θ?θπ?π===≤≤≤≤； 解：theta=linspace(0,2*pi,50);fai=linspace(0,pi/2,20); [theta,fai]=meshgrid(theta,fai); x=2*sin(fai).*cos(theta);

y=2*sin(fai).*sin(theta);z=2*cos(fai); surf(x,y,z)

⑺三条曲线合成图

123sin ,sin sin(10),sin ,0y x y x x y x x π===-<<。

解：

7、作下列x=linspace(0,pi,100);

y1=sin(x);y2=sin(x).*sin(10*x);y3=-sin(x);

plot(x,y1,x,y2,x,y3)分段函数图 1.1 1.1 1.11.1 -1.1x y x x x >??

=≤??-

matlab基础作图实例

实验三 MATLAB 的绘图 一、实验目的：掌握利用MATLAB 画曲线和曲面。 二、实验内容： 1、 在不同图形中绘制下面三个函数t ∈[0，4π]的图象，3个图形分别是 figure(1),figure(2),figure(3)。 ) sin(41.0321t e y t y t y t -== =π 说明:y 1 线型：红色实线，y 2 线型：黑色虚线，y 3: 线型：兰色点线 分别进行坐标标注,分别向图形中添加标题‘函数1’，‘函数2’, ‘函数3’ 解答： 源程序与图像： t=0:0.1:4*pi; y_1=t; y_2=sqrt(t); y_3=4*pi.*exp(-0.1*t).*sin(t); figure(1) plot(t,y_1,'-r'); title('函数1'); xlabel('t');ylabel('y_1'); figure(2) plot(t,y_2,'--k'); title('函数2'); xlabel('t');ylabel('y_2'); figure(3) plot(t,y_3,':b'); title('函数3'); xlabel('t');ylabel('y_3'); 246 8101214 02468 10 12 14 函数1 t y 1

0246 8101214 0.511.522.533.54函数2 t y 2 2 4 6 8 10 12 14 -8-6-4-2024 681012函数3 t y 3 2、 在同一坐标系下绘制下面三个函数在t ∈[0，4π]的图象。 （用2种方法来画图，其中之一使用hold on ） 使用text 在图形适当的位置标注“函数1”“函数2”,“函数3” 使用gtext 重复上面的标注，注意体会gtext 和text 之间的区别 解答： 方法一： 程序与图形： t=0:0.1:4*pi; y_1=t; y_2=sqrt(t); y_3=4*pi.*exp(-0.1*t).*sin(t); figure(1) plot(t,y_1,'-r'); gtext('函数1');

MATLAB画图入门篇--各种基本图形绘制的函数与实例

MATLAB画图入门篇--各种基本图形绘制的函数与实例【来自网络】 一．二维图形(Two dimensional plotting) 1.基本绘图函数(Basic plotting function)：Plot,semilogx,semilogy,loglog,polar,plotyy (1).单矢量绘图(single vector plotting)：plot(y),矢量y的元素与y元素下标之间在线性坐标下的关系曲线。 例1：单矢量绘图 y=[00.62.358.311.71517.719.420];plot(y) 可以在图形中加标注和网格， 例2：给例1的图形加网格和标注。 y=[00.62.358.311.71517.719.420];plot(y) title('简单绘图举例');xlabel('单元下标');ylabel('给定的矢量');grid (2).双矢量绘图(Double vector plotting)：如x和y是同样长度的矢量,plot(x,y)命令将绘制y元素对应于x元素的xy曲线图。 例：双矢量绘图。 x=0:0.05:4*pi;y=sin(x);plot(x,y) (3).对数坐标绘图(ploting in logarithm coordinate)：x轴对数semilogx,y轴对数semilogy,双对数loglog, 例：绘制数组y的线性坐标图和三种对数坐标图。 y=[00.62.358.311.71517.719.420]; subplot(2,2,1);plot(y);subplot(2,2,2);semilogx(y) subplot(2,2,3);semilogy(y);subplot(2,2,4);loglog(y) （4）极坐标绘图(Plotting in polar coordinate)： polar(theta,rho)theta—角度，rho—半径 例：建立简单的极坐标图形。 t=0:.01:2*pi;polar(t,sin(2*t).*cos(2*t)) 2.多重曲线绘图(Multiple curve plotting) （1）一组变量绘图(A group variable plotting) plot(x,y) (a)x为矢量，y为矩阵时plot(x,y)用不同的颜色绘制y矩阵中各行或列对应于x的曲线。 例1： x=0:pi/50:2*pi;y(1,:)=sin(x);y(2,:)=0.6*sin(x);y(3,:)=0.3*sin(x);plot(x,y) (b)x为矩阵，y为矢量时绘图规则与（a）的类似，只是将x中的每一行或列对应于y进行绘图。。 例2： x(1,:)=0:pi/50:2*pi;x(2,:)=pi/4:pi/50:2*pi+pi/4;x(3,:)=pi/2:pi/50:2*pi+pi/2; y=sin(x(1,:));plot(x,y) (c)x和y是同样大小的矩阵时,plot(x,y)绘制y矩阵中各列对应于x各列的图形。 例3： x(:,1)=[0:pi/50:2*pi]';x(:,2)=[pi/4:pi/50:2*pi+pi/4]';x(:,3)=[pi/2:pi/50:2*pi+pi/2]'; y(:,1)=sin(x(:,1));y(:,2)=0.6*sin(x(:,1));y(:,3)=0.3*sin(x(:,1)); plot(x,y) 这里x和y的尺寸都是101×3，所以画出每条都是101点组成的三条曲线。如行列转置后就会画出101条曲线，每条线

做一个实例,搞定MATLAB界面编程

做一个实例，搞定MATLAB界面编程 2009-04-12 16:55 作者：彭军 带插图版本的pdf格式文件，请下载： https://www.wendangku.net/doc/19578959.html,/source/1144420 一个实例搞定MATLAB界面编程 作者：彭军 邮件：pengjun@https://www.wendangku.net/doc/19578959.html, 博客：https://www.wendangku.net/doc/19578959.html,/pengjun 下面请跟我一步一步做一个图像处理的程序，如果您坚持做完这个实例，我想MATLAB界面编程对您而言，就没有什么难度了。当然，我这里说的是，您首先要有一定的MATLAB编程基础。还有，我的MATLAB版本是2008a。在2008a以前的版本中没有工具栏编辑器，如果需要工具栏要手动写程序，这个我就不多讲了。好了，废话少说，跟我来吧！ 1、在MATLAB的命令窗口(Command Window)中运行guide命令，来打开GUIDE 界面，如下： 2、然后，选择空模板(Blang GUI)，点击OK，即可打开GUIDE的设计界面，如下： 3、点击工具栏上的菜单编辑器(Menu Editor)，打开菜单编辑器，如下： 4、在Menu Bar中新建一个菜单项，名字为“文件”，其他设置请看下图： 5、在“文件”菜单下添加菜单项：“打开”，“保存”，“退出”。见下图：如果需要在菜单项“退出”上面添加一个分割线的话，选中“Separator above this item”就行了。 保存我的界面为pjimage.fig. 保存完毕之后，会自动打开pjimage.m文件，而我们所有的程序都是要写在这个M文件里面的。在编程中，我们的每一个鼠标动作都对应一个Callback函数。那么我们的菜单项也是如此的。 在界面上，单击鼠标右键选择“Property Inspector”，即可打开属性窗口。当我们点击不同的控件时，其对应的属性都会在这里显示，我们可以进行修改。最主要的属性莫过于Tag属性和String属性。 设置当前Figure窗口的Tag属性为：figure_pjimage，窗口的标题(Name属性)为：图像处理实例。如下： 然后，点击工具栏的保存按钮。之后，点击工具栏的运行按钮(Run Figure)。注意，工具栏的图标都会有提示的，像运行按钮的提示就是Run Figure. 我们会看到如下的界面： 那说明，我们保存的.fig文件的目录不是当前目录，但是没关系啊，我们只要

Matlab经典案例

1、三维曲线 >> t=0:pi/50:10*pi; >> plot3(sin(2*t),cos(2*t),t) >> axis square >> grid on 2、一窗口多图形 >> t=-2*pi:0.01:2*pi; >> subplot(3,2,1) >> plot(t,sin(t)) >> subplot(3,2,2) >> plot(t,cos(t)) >> subplot(3,2,3) >> plot(t,tan(t)) >> axis([-pi pi -100 100]) >> subplot(3,2,4) >> plot(t,cot(t)) >> axis([-pi pi -100 100]) >> subplot(3,2,5) >> plot(t,atan(t)) >> subplot(3,2,6) >> plot(t,acot(t)) 3、图形样式、标注、题字(也可以利用菜单直接 Insert) >> x=0:pi/20:2*pi; >> plot(x,sin(x),'b-.') >> hold on >> plot(x,cos(x),'r--') >> hold on >> plot(x,sin(x)-1,'g:')

>> hold on >> plot(x,cos(x)-1) >> xlabel('x'); >> xlabel('x轴'); >> ylabel('y轴'); >> title('图形样式、标注等'); >> text(pi,sin(pi),'x=\pi'); >> legend('sin(x)','cos(x)','sin(x)-1','cos(x)-1'); >> [x1,y1]=ginput(1) %利用鼠标定位查找线上某点的值x1 = 2.0893 y1 = -0.5000 >> gtext('x=2.5') %鼠标定位放置所需的值在线上 4、 >> fplot('[sin(x),cos(x),sqrt(x)-1]',[0 2*pi]) M文件：myfun.m 内容如下： function y=myfun(x) y(:,1)=sin(x); y(:,2)=cos(x); y(:,3)=x^(1/2)-1; 再运行：>> fplot('myfun',[0 2*pi]) 同样可以得到右图 5、 >> [x,y]=fplot('sin',[0 2*pi]); >> [x1,y1]=fplot('cos',[0 2*pi]); >> plot(x,y,'-r',x1,y1,'-.k') >> legend('y=sinx','y=cosx') 6、

matlab中plot基础绘图实例

%plot基本绘图 x=0:0.1:2*pi; y=sin(x) plot(x,y) %两个参数都是矩阵 x1=0:0.1:2*pi; x2=-pi:0.1:pi; y1=sin(x) y2=cos(x) plot(x1,y1,x2,y2)%多条曲线绘制在统一坐标轴上 %plot只有一个参数 x=linspace(0,2*pi,200) y=sin(x) plot(y) y2=cos(x) y3=y+i*y2%横坐标实部为正弦，纵坐标虚部为余弦，构成一个圆形 plot(y3) axis equal%将上述图型的横纵坐标调整为相同，使得椭圆变为正圆 %plot含有多个参数 x1=linspace(0,2*pi,200) x2=linspace(0,2*pi,100) x3=linspace(0,2*pi,50) y1=cos(x1) y2=sin(x2) y3=0.01*exp(x3) plot(x1,y1,x2,y2,x3,y3)%当x1，x2，x3不同维数（点数不同）可用高方法绘制 %线性选项 x=0:0.1:2*pi; y=sin(x) plot(x,y,'r')%颜色，g为绿色，y为黄色，k为黑色，默认b为蓝色 plot(x,y,'*')%形状，*为*状，p为五角星,.为小方块 plot(x,y,'--')%--为短线，：为虚线，-.为点虚线 plot(x,y,'*r--')%可以组合使用，只需用一对单引号把要求全部括起来 %标注 x=0:0.1:2*pi; y=sin(x) plot(x,y) xlabel('x')%横坐标轴名称 ylabel('y')%纵坐标轴名称

Matlab画图示例

Matlab 画图示例 1 文档说明 期刊文献中对图形格式、图中曲线线型，文字大小等都会有要求，而matlab 默认输出的图形格式和显示方式可能不符合期刊要求。本文档简单给出了几种调整matlab 输出图形格式的方法及示例图形，仅供参考。 2 示例 2.1 图形格式调整 下面给出设置图片输出格式及调整图形大小、坐标轴及题注字体的示例。如作图y=sin(x),要求图中线型为2，颜色为白底黑色；坐标轴字体为14，标注字体为12，图形宽度为21cm ，以tiff 格式输出，要求分辨率为600dpi ； 102030 40506070 角度 幅值 示例图片 图1 图形格式调整示例 点击edit--copy figure ，将图片粘贴到word,图片输出如图1。 Matlab 程序及说明见附录1 2.2 多曲线显示设置 MATLAB 在多组变量绘图时,可将曲线以不同的颜色,不同的线型、线宽及标记点表示出来以示区别。常用的选项设置如下：

各种颜色属性 'r' 红色 'm' 粉红 'g' 绿色 'c' 青色 'b' 兰色 'w' 白色 'y' 黄色 'k' 黑色 各种线型属性 '-' 实线 '--' 虚线 ':' 点线 '-.' 点划线 各种标记点属性选项 '.' 用点号绘制各数据点 '^' 用上三角绘制各数据点 '+' 用'+'号绘制各数据点 'v' 用下三角绘制各数据点 '*' 用'*'号绘制各数据点 '>' 用右三角绘制各数据点 ' .' 用'.'号绘制各数据点 '<' 用左三角绘制各数据点 's'或squar 用正方形绘制各数据点 'p' 用五角星绘制各数据点 'd'或diamond 用菱 形绘制各数据点 'h' 用六角星绘制各数据点 这些选项可以连在一起用,如:'-.g'表示绘制绿色的点划线,'g+'表示用绿色的 '+'号绘制曲线。 注意: 1)表示属性的符号必须放在同一个字符串中; 2)可同时指定2~3个属性; 3) 与先后顺序无关; 4)指定的属性中,同一种属性不能有两个以上. 下面给出几个示例。 2.2.1 线宽显示示例 作出四条曲线，线宽分别为1，3，5，7 10 20 30 40 50 60 70 角度 幅值 示例图片 图2 线宽显示示例

最全的MATLAB绘图命令

Matlab绘图 强大的绘图功能是Matlab的特点之一，Matlab提供了一系列的绘图函数，用户不需要过多的考虑绘图的细节，只需要给出一些基本参数就能得到所需图形，这类函数称为高层绘图函数。此外，Matlab还提供了直接对图形句柄进行操作的低层绘图操作。这类操作将图形的每个图形元素（如坐标轴、曲线、文字等）看做一个独立的对象，系统给每个对象分配一个句柄，可以通过句柄对该图形元素进行操作，而不影响其他部分。 本章介绍绘制二维和三维图形的高层绘图函数以及其他图形控制函数的使用方法，在此基础上，再介绍可以操作和控制各种图形对象的低层绘图操作。 一．二维绘图 二维图形是将平面坐标上的数据点连接起来的平面图形。可以采用不同的坐标系，如直角坐标、对数坐标、极坐标等。二维图形的绘制是其他绘图操作的基础。 一．绘制二维曲线的基本函数 在Matlab中，最基本而且应用最为广泛的绘图函数为plot，利用它可以在二维平面上绘制出不同的曲线。 1． plot函数的基本用法 plot函数用于绘制二维平面上的线性坐标曲线图，要提供一组x坐标和对应的y 坐标，可以绘制分别以x和y为横、纵坐标的二维曲线。plot函数的应用格式 plot(x,y) 其中x,y为长度相同的向量，存储x坐标和y坐标。 例51 在[0 , 2pi]区间，绘制曲线 程序如下：在命令窗口中输入以下命令 >> x=0:pi/100:2*pi; >> y=2*exp(-0.5*x).*sin(2*pi*x); >> plot(x,y) 程序执行后，打开一个图形窗口，在其中绘制出如下曲线 注意：指数函数和正弦函数之间要用点乘运算，因为二者是向量。 例52 绘制曲线

Matlab图形绘制经典案例

Matlab图形绘制经典案例---受用无穷1、 三维曲线 >> t=0:pi/50:10*pi; >> plot3(sin(2*t),cos(2*t),t) >> axis square >> grid on 2、一窗口多图形 >> t=‐2*pi:0.01:2*pi; >> subplot(3,2,1) >> plot(t,sin(t)) >> subplot(3,2,2) >> plot(t,cos(t)) >> subplot(3,2,3) >> plot(t,tan(t)) >> axis([‐pi pi ‐100 100]) >> subplot(3,2,4) >> plot(t,cot(t)) >> axis([‐pi pi ‐100 100])

>> subplot(3,2,5) >> plot(t,atan(t)) >> subplot(3,2,6) >> plot(t,acot(t)) 3、 图形样式、标注、题字 (也可以利用菜单直接Insert) >> x=0:pi/20:2*pi; >> plot(x,sin(x),'b‐.') >> hold on >> plot(x,cos(x),'r‐‐') >> hold on >> plot(x,sin(x)‐1,'g:') >> hold on >> plot(x,cos(x)‐1) >> xlabel('x'); >> xlabel('x轴'); >> ylabel('y轴'); >> title('图形样式、标注等'); >> text(pi,sin(pi),'x=\pi'); >> legend('sin(x)','cos(x)','sin(x)‐1','cos(x)‐1');

Matlab作图例子

基本xy平面绘图命令 MATLAB不但擅长於矩阵相关的数值运算，也适合用在各种科学目视表示（Scientific visualization）。 本节将介绍MA TLAB基本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 三维图形绘制实例

三维图形 一． 三维曲线 plot3(x1,y1,z1,选项1,x2,y2,z2,选项2,…,xn,yn,zn,选项n) 其中每一组x,y,z 组成一组曲线的坐标参数，选项的定义和plot 函数相同。当x,y,z 是同维向量时，则x,y,z 对应元素构成一条三维曲线。当x,y,z 是同维矩阵时，则以x,y,z 对应列元素绘制三维曲线，曲线条数等于矩阵列数。 Example1.绘制三维曲线。 程序如下： clf, t=0:pi/100:20*pi; x=sin(t); y=cos(t); z=t.*sin(t).*cos(t); %向量的乘除幂运算前面要加点 plot3(x,y,z); title('Line in 3-D Space'); xlabel('X');ylabel('Y');zlabel('Z'); grid on; 所的图形如下： X Line in 3-D Space Y Z 二． 三维曲面 1. 产生三维数据 在MATLAB 中，利用meshgrid 函数产生平面区域内的网格坐标矩阵。

语句执行后，矩阵X 的每一行都是向量x ，行数等于向量y 的元素的个数，矩阵Y 的每一列都是向量y ，列数等于向量x 的元素的个数。 2. 绘制三维曲面的函数 surf 函数和mesh 函数 example2. 绘制三维曲面图z=sin(x+sin(y))-x/10。 程序如下： clf, [x,y]=meshgrid(0:0.25:4*pi); %产生平面坐标区域内的网格坐标矩阵 z=sin(x+sin(y))-x./10; surf(x,y,z); axis([0 4*pi 0 4*pi -2.5 1]); title('surf 函数所产生的曲面'); figure; mesh(x,y,z); axis([0 4*pi 0 4*pi -2.5 1]); title('mesh 函数所产生的曲面'); -2.5 -2-1.5-1-0.500.51surf 函数所产生的曲面