matlab计算结果的可视化

第五讲计算结果的可视化

本节介绍MATLAB 的两种基本绘图功能：二维平面图形和三维立体图形。5.1 二维平面图形

5.1.1 基本图形函数

plot是绘制二维图形的最基本函数，它是针对向量或矩阵的列来绘制曲线的。也就是

说，使用plot 函数之前，必须首先定义好曲线上每一点的x 及y 坐标，常用格式为：

（1）plot(x) 当x 为一向量时，以x 元素的值为纵坐标，x 的序号为横坐标值绘制

曲线。当x 为一实矩阵时，则以其序号为横坐标，按列绘制每列元素值相对于其序号的曲

线，

当x 为m× n 矩阵时，就由n 条曲线。

（2）plot(x,y) 以x 元素为横坐标值，y 元素为纵坐标值绘制曲线。

（3）plot(x,y1,x,y2,…) 以公共的x 元素为横坐标值，以y1,y2,… 元素为纵坐标值绘

制多条曲线。

例5.1.1 画出一条正弦曲线和一条余弦曲线。

>> x=0:pi/10:2*pi;

>> y1=sin(x);

>> y2=cos(x);

>> plot(x,y1,x,y2)

图5.1.1 函数plot 绘制的正弦曲线

在绘制曲线图形时，常常采用多种颜色或线型来区分不同的数据组，MATLAB 软件专

门提供了这方面的参数选项（见表5.1.1），我们只要在每个坐标后加上相关字符串，就可

实

现它们的功能。

表5.1.1 绘图参数表

色彩字符颜色线型字符线型格式标记符号数据点形式标记符号数据点形式

y 黄- 实线. 点<小于号

m 紫：点线o 圆s 正方形

c 青-. 点划线x 叉号

d 菱形

r 红- - 虚线+ 加号h 六角星

g 绿* 星号p 五角星

b 蓝v 向下三角形

w 白^ 向上三角形

k 黑>大于号

例如，在上例中输入

>> plot(x,y1,'r+-',x,y2,'k*:')

则得图5.1.2

图5.1.2 使用不同标记的plot 函数绘制的正弦曲线

5.1.2 图形修饰

MATLAB 软件为用户提供了一些特殊的图形函数，用于修饰已经绘制好的图形。表5.1.2 图形修饰函数表

函数含义

grid on (/off) 给当前图形标记添加（取消）网络

xlable(‘string’) 标记横坐标

ylabel(‘string’) 标记纵坐标

ti tle(‘string’) 给图形添加标题

text(x,y,’string’) 在图形的任意位置增加说明性文本信息

gtext(‘string’) 利用鼠标添加说明性文本信息

axis([xminxmaxyminymax]) 设置坐标轴的最小最大值

- 3 -

例5.1.2 给例5.1.1 的图形中加入网络和标记。（见图5.1.3 和5.1.4）

>> x=0:pi/10:2*pi;

>> y1=sin(x);

>> y2=cos(x);

>> plot(x,y1,x,y2)

>> grid on

>>xlabel('independent variable X')

>>ylabel('Dependent Variable Y1 & Y2')

>> title('Sine and Cosine Curve')

>> text(1.5,0.3,'cos(x)')

>>gtext('sin(x)')

>> axis([0 2*pi -0.9 0.9])

图5.1.3 使用了图形修饰的plot 函数绘制的正弦曲线

5.1.3 图形的比较显示

在一般默认的情况下，MATLAB 每次使用plot 函数进行图形绘制，将重新产

生一个图

形窗口。但有时希望后续的图形能够和前面所绘制的图形进行比较。一般来说

有两种方法

：

一是采用hold on(/off)命令，将新产生的图形曲线叠加到已有的图形上；二是

采用

subplot(m,n,k)函数，将图形窗口分隔成n m×个子图，并选择第k 个子图作

为当前图形

，然

后在同一个视图窗口中画出多个小图形。

例5.1.3 在同一窗口中绘制线段。（见图5.1.5）

>> x=0:pi/10:2*pi;

>> y1=sin(x);

>> y2=cos(x);

- 4 -

>> y4=log(x);

>> plot(x,y1,x,y2)

>> hold on

>> plot(x,y3)

>> plot(x,y4)

>> hold off

例5.1.4 在多个窗口中绘制图形。（见图5.1.6）

>> x=0:pi/10:2*pi;

>> y1=sin(x);

>> y2=cos(x);

>> y3=exp(x);

>> y4=log(x);

>> subplot(2,2,1);

>> plot(x,y1);

>> subplot(2,2,2);

>> plot(x,y2);

>> subplot(2,2,3);

>> plot(x,y3);

>> subplot(2,2,4);

>> plot(x,y4);

[说明] （1）子窗口的序号按行由上往下，按列从左向右编号。

（2）如果不用指令clf清除，以后图形将被绘制在子图形窗口中。

图5.1.4 设置坐标轴最大最小值的正弦曲线

- 5 -

图5.1.5 图形的比较显示（曲线叠加方法）

图5.1.6 图形的比较显示（图形窗口分割方法）

5.2 三维立体图形

5.2.1 三维曲线图

与二维图形相对应，MATLAB 提供了plot3 函数，可以在三维空间中绘制三维曲线，

它的格式类似于plot，不过多了z 方向的数据。plot3 的调用格式为：

plot3(x1,y1,z1,x2,y2,z2,...)

其中x1,y1,z1,x2,y2,z2,…等分别为维数相同的向量，分别存储着曲线的三个坐标值，该

函数

的使用方式和plot 类似，也可以采用多种的颜色或线型（见表5.1.1）来区分不同的数据

组，

只需在每组变量后面加上相关字符串即可实现该功能。

例5.2.1 绘制方程x=t

y=sin(t)

在t=[0,2*pi]上

的空间方程。（见图5.2.1）

>>clf

>> x=0:pi/10:2*pi;

>> y1=sin(x);

>> y2=cos(x);

>> plot3(y1,y2,x,'m:p')

>> grid on

>>xlabel('Dependent Variable Y1')

>>ylabel('Dependent Variable Y2')

>>zlabel('Independent Variable X')

>> title('Sine and Cosine Curve')

图5.2.1 函数plot 绘制的三维曲线图

5.2.2 三维曲面图

如果要画一个三维的曲面，可以使用mesh(X,Y,Z)或surf(X,Y,Z)函数来实现。mesh 函数为数据点绘制网格线，图形中的每一个已知点和其附近的点用直线连接。surf

函数和mesh 的用法类似，但它可以画出着色表面图，图形中的每一个已知点与其相邻点以

平面连接。

为方便测试立体绘图，MATLAB 提供了一个peaks 函数，它可以产生一个的高

斯分布矩阵，其生成方程是

N N ×

z=3*(1-x).^2.*exp(-(x.^2)-(y+1).^2)-10*(x/5-x.^3-y.^5).*exp(-x.^2-y.^2)-1/3*e x

p(-(x+1).^2-y.^2)

对应的图形是一个凹凸有致的曲面，包含了三个局部极大点及三个局部极小点。下面使用peaks 函数来比较一下mesh 和surf 的区别。

例5.2.2 分别用mesh 函数和surf 函数绘制高斯矩阵的曲面。

>> z=peaks(40);

>> mesh(z);

>> surf(z);

图5.2.2 mesh 函数绘制的三维曲面图

图5.2.3 surf 函数绘制的着色表面图

在曲面绘图中，另一个常用的函数是meshgrid函数，其一般引用格式是：[X, Y]=meshgrid (x, y)

其中x 和y 是向量，通过meshgrid函数就可将x 和y 指定的区域转换成为

矩阵X 和Y。

这样我们在绘图时就可以先用meshgrid函数产生在x-y 平面上的二维的网格数据，再以一

组z 轴的数据对应到这个二维的网格，即可画出三维的曲面。

例5.2.3 绘制方程

sin((x^2+y^2)^(1/2))

z = ---------------------

(x^2+y^2)^(1/2)

在x∈[-7.5,7.5];y∈[-7.5,7.5] 的图形。

>> x=-7.5:0.5:7.5;y=x;

>> [X,Y]=meshgrid(x,y);

>> R=sqrt(X.^2+Y.^2)+eps;

>> Z=sin(R)./R;

>> surf(X,Y,Z)

>>xlabel('X 轴方向')

>>ylabel('Y 轴方向')

>>zlabel('Z 轴方向')

(见图5.2.4)

_(x^2+y^2)

例5.2.4 绘制由方程形成的立体图。（见图5.2.5）z=xe

>> clear

>> x=-2:0.1:2;y=x;

>> [X,Y]=meshgrid(x,y);

>> Z=X.*exp(-X.^2-Y.^2);

- 9 -

>> surf(X,Y,Z)

图5.2.4

图5.2.5

5.2.3 观察点

MTALAB 允许用户设置观察点，其指令是：view(azimuth,elevation)

其中方位角azimuth 是观察点和坐标原点连线在x-y 平面的投影和y 轴负方向的夹角，仰

角

elevation 是观察点与坐标原点的连线和x-y 平面的夹角。对于这两个角度，三维图形的

默认

值分别是-37.5 和30，二维图形的默认值是0 和90。

例5.2.5 从不同的角度观察高斯矩阵的曲面。

>> z=peaks(40);

>> subplot(2,2,1);

>> mesh(z);

>> subplot(2,2,2);

>> mesh(z);

>> view(-37.5,-30);

>> subplot(2,2,3);

>> mesh(z);

>> view(180,0);

>> subplot(2,2,4);

>> mesh(z);

>> view(0,90);

图5.2.6 对应不同观察点的三维曲面图

5.3 其他图形函数

除了plot 绘图函数以外，在有些场合对绘制的曲线会有一些特殊要求，这就要其他函

数来实现，常用的几种函数如下（见表5.3.1）

表5.3.1 其他图形函数表

函数含义

loglog使用对数坐标系绘图

semilogx横坐标为对数坐标轴，纵坐标为线性坐标轴

semilogy横坐标为线性坐标轴，纵坐标为对数坐标轴

polar 绘制极坐标图

fill 绘制实心图

bar 绘制直方图

pie 绘制饼图

area 绘制面积图

quiver 绘制向量场图

stairs 绘制阶梯图

sterm绘制火柴杆图

例5.3.1

>> x=0:pi/10:2*pi;

>> y1=sin(x);

>> subplot(2,2,1);

>> plot(x,y1);

>> subplot(2,2,2);

>> bar(x,y1);

>> subplot(2,2,3);

>> fill(x,y1,'g');

>> subplot(2,2,4);

>> stairs(x,y1,'k');

图5.3.1 其他图形函数

5.3.1 直方图

函数bar(x)可以绘制直方图，这对统计或者数据采集非常直观实用。它共有四种形式：

bar,bar3,barh 和bar3h，其中bar 和bar3 分别用来绘制二维和三维竖直方图，barh和b

ar3h 分

别用来绘制二维和三维水平直方图，调用格式是：

bar(x,y) 其中x 必须单调递增或递减，y 为n m×矩阵，可视化结果为m 组，每

组n 个垂直柱，也就是把y 的行画在一起，同一列的数据用相同的颜色表示；bar(x,y,width) （或bar(y,width)）指定每个直方条的宽度，如width>1，则直方条会重

叠，默认值为width=0.8；

bar(…,’grouped’) 使同一组直方条紧紧靠在一起；

bar(…,’stack’) 把同一组数据描述在一个直方条上。

例5.3.2

>> y=[5 3 2 9;4 7 2 7;1 5 7 3];

>> subplot(2,2,1),bar(y)

>> x=[5 9 11];

>> subplot(2,2,2),bar3(x,y)

>> subplot(2,2,3),bar(x,y,'grouped')

>> subplot(2,2,4),bar(rand(2,3),.75,'stack')

图5.3.2 直方图

5.3.2 面积图

函数area 用来绘制面积图，面积图在plot 的基础上填充x 轴和曲线之间的面积，该图

用于查看某个数在该列所有数的总和中所占的比例。

例5.3.3

>> x=-3:3;

>> y=[3 2 5;6 1 8;7 4 9;6 3 7;8 2 9;4 2 9;3 1 7];

>> area(x,y)

- 13 -

图5.3.3 面积图

5.3.3 饼图

函数pie 用来绘制饼图，它可以形象地表示出向量中各元素所占比例。其调用格式是：

pie(x) x 中的元素通过x/sum(x)进行归一化，以确定饼图中的份额；

pie(x,explode) 向量explode 和x 元素数相同，用来指出需要分开的饼片，explode 中

不为零的部分会被分开。

图5.3.4 饼图

例5.3.4 设某班的某课程的考试成绩如下：90 分以上有32 人，81 至90 有58 人，71

至80 分有27 人，60 至70 分为21 人，60 分以下有16 人，画出饼图。（见图5.3.4）

>> x=[32 58 27 21 16];

>> explode0=[1 0 0 0 0];

>> subplot(1,2,1)

>> pie(x,explode0)

>> explode1=[0 0 0 0 1];

>> subplot(1,2,2)

>> pie(x,explode1)

5.3.4 不同坐标系中的绘图

Semilogx,semilogy,loglo,polar(theta,rho)的使用方法和plot 完全类似，不同的只是绘

制到

不同的图形坐标上。函数semilogx绘制x 轴为对数标度的图形，在半对数坐标系中绘图；

函数semilogy绘制y 轴为对数标度的图形；函数loglog绘制两个轴都为对数间隔的图形

；

函数polar(theta,rho)绘制极坐标图形，其中theta 为相角，rho 为其对应的半径。

例5.3.5 绘制ρ=acos(3θ),a=2 的图形。（见图5.3.5）

>> theta=-pi:pi/80:pi;

>> polar(theta,2*cos(3*theta))

图5.3.5 极坐标图

5.4 符号表达式绘图

MATLAB 软件提供了将表达式进行图形显示的功能。完成此功能需调用fplot 函数和

ezplot函数。

函数fplot用来绘制数学函数，其调用格式为：fplot(fun,lims)

其中fun 就是所要绘制的函数，可以是定义函数的M 文件名，也可以是以x 为变量的可计算字符串。例如’diric(x,10)’或’[sin(x),cos(x)]’，对于向量x 的每个元素，函数fun(x)必须返回一个行向量。如果fun 返回[f1(x),f2(x),f3(x)]，输入[x1;x2]，就会返回矩阵f1(x1) f2(x1) f3(x1).

f1(x2) f2(x2) f3(x2)

lims=[XMIN XMAX YMIN YMAX]限定了x,y轴上的绘图空间。

例5.4.1

>> subplot(2,2,1),fplot('humps',[0 1])

>> subplot(2,2,2),fplot('abs(exp(-j*x*(0:9))*ones(10,1))',[0 2*pi])

>> subplot(2,2,3),fplot('[tan(x),sin(x),cos(x)]',2*pi*[-1 1 -1 1])

>> subplot(2,2,4),fplot('sin(1./x)',[0.01 0.1],1e-3)

图5.4.1 fplot函数绘制表达式图形

ezplot函数是简捷绘图指令之一，它无需数据准备，直接画出函数图形，基本调用格式

为ezplot(f)

其中f 是字符串或代表数学函数的符号表达式，只有一个符号变量，可以是x，缺省情况下

x 轴的绘图区域为[?2π, 2π]但我们可以用ezplot(f,xmin,xmax)或ezp

lot(f,[xmin,xmax])来指

定x 的范围。

例5.4.2

>> y='x^2';

>> subplot(1,2,1)

>>ezplot(y)

>> subplot(1,2,2)

>> y='sin(x)';

>>ezplot(y,[-pi,pi])

图5.4.2 ezplot函数绘制表达式图形

5.5 plot 函数

MATLAB 对数据是按列存储和计算的，运用plot(x)时，当x 为一个向量时，以其元

素为纵坐标，其序号为横坐标值绘制曲线。当x 为实矩阵时，则以其序号为横坐标，按列

绘

制每列元素相对于序号的曲线，当x 为[m,n]矩阵时，就有n 条曲线。

如果x,y是同维向量，plot(x,y)指令以x 元素为横坐标值，y 元素为纵坐标值绘制曲线。

如x 是向量，y 是有一维与x 元素数量相等的矩阵，则以x 为共同横坐标，按列绘制y 每

列

元素值，曲线数为y 的另一维的元素数。如果x,y是同维矩阵，则以x,y对应列元素为、

纵坐标分别绘制曲线，数目等于矩阵的列数。

例5.5.1

>> x=[3 5 10 8];

>> subplot(2,2,1)

>> plot(x)

>> x=[3 5 10 8;7 2 9 4;2 7 2 7]';

>> subplot(2,2,2)

>> plot(x)

>> x=[3 5 6 8];

>> y=[1 5 10 4];

>> subplot(2,2,3)

>> plot(x,y)

>> x=[1 3 5 7;2 4 6 8]';

>> y=[6 2 5 10;3 5 2 6]';

>> subplot(2,2,4)

>> plot(x,y,'k:*')

图5.5.1

5.6 交互式图形指令

ginput是一个比较特殊的图形指令，用作获取图上数据，例如指令

>>[x,y]=ginput(6) % 从图形上选取6 个点

此时，ginput指令将把当前图形调入前台，同时光标变为十字叉，移动光标，使交叉点落

在目标点上，单击鼠标，即可获得该点数据。

例5.5.2

>>fplot('humps',[0 1])

>>ginput(6)

x =

0.0449

0.1832

0.3007

0.3813

0.6417

0.8952

y =

7.4561

38.1579

96.3450

57.4561

10.9649

21.1988

基于MATLAB平台的可视化图像处理系统设计方法

[收稿日期]20050228 [作者简介]黄书先(1962),女,1983年大学毕业,硕士,副教授,现主要从事石油勘探开发研究与科研管理工作。 基于MATLAB 平台的可视化图像处理系统 设计方法 黄书先 (长江大学科学技术处,湖北荆州434023) [摘要]以MATLAB 为工作语言和开发环境,开发了一个在M AT LAB 平台下的可视化图像处理系统, 可实现包括对测井图像在内的一般图像的精细处理,并能和用户开发的程序接口。为M AT LAB 的再开发 和可视化系统的设计作了有益的探索。 [关键词]MATLAB;图像处理;可视化GUI 界面 [中图分类号]TP 311111 [文献标识码]A [文章编号]16731409(2005)04015803 MAT LAB 的图像处理工具为自然科学各学科领域的学者、研究人员和工程师提供了一个直观的灵活的环境,用以解决复杂的图像处理问题。用MAT LAB 语言开发的图像处理算法可以在所有支持MAT LAB 的平台上共享。也可以将m 语言算法和现存的C 程序集成在一起或者将MATLAB 开发的m 语言算法和GU Is 编译为C/C++代码,供其他程序调用,或者发布为一个独立的应用程序。下面笔者提出一个基于MATLAB 平台的可视化图像处理系统设计方法,可用于包括测井图像资料在内的一般图像的精细处理。 1 系统总体设计 在MAT LAB 中有个重要的图像处理工具包[1],该工具包是由一系列支持图像处理操作的函数组成的,按功能可以分为以下几类:图像显示;图像文件输入与输出;几何操作;像素值统计;图像增强;图像识别;图像滤波;图像变换;邻域和块操作;二值图像操作;颜色映射和颜色空间转换;图像格式转换等。和其他工具包一样,用户还可以根据需要书写自己的函数,以满足特定的需要;也可以将这个工具包和信号处理工具包或小波工具包等其他工具包联合起来使用。 MAT LAB 提供了交互式的GU I 开发环境[2~4],用户只需要设置各个对象相应的属性,系统自动生成与之对应的界面,大大减少了开发的难度。本设计的MAT LAB 图像处理系统由封面、主界面和各个子功能界面组成,其框图如图1所示。程序总流程图如图2所示。 2 封面界面设计 首先用MAT LAB 编辑封面的脚本文件,生成系统的封面界面,然后再进入处理的主界面。 在设计封面时,要做好封面的总体布局,力求完美。要插入背景,首先要找到所用的函数(这在一般的书中很难找到);其次要注意设置axes 的属性。在显示背景图像时,不能用imshow ()函数,如果用此函数,背景图像只能显示在封面的一部分,不能整屏覆盖;而要用imagesc ()函数。MATLAB 提供修改文本中文字的字体,它支持华文中宋、华文彩云、华文仿宋、楷书、黑体等,功能较强。 在设计封面的过程中,最主要的技术是要解决时间的显示,如果只用MATLAB 中的clock ()函数,则显示的是静态时间,而不会显示和电脑同步的动态时间。为了解决这个问题,需做一个循环判断语句:while find (get (0,'c hildren'))==h0。成立的条件是,只要是当前窗口循环,否则停止,这样可以减少CPU 的负担。同时,要实时提取clock ()函数,可用fix (clock),使提取的时间更美观。#158#长江大学学报(自科版)2005年4月第2卷第4期/理工卷第2卷第2期 Jour nal of Yangtze University (Nat Sci Edit)Apr 12005Vol 12No 14/Sci &Eng V,Vol 12No 12

潮流计算(matlab)实例计算

潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新

的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图

MATLAB精通科学计算_偏微分方程求解

一、Maple V 系统 Maple V是由Waterloo大学开发的数学系统软件，它不但具有精确的数值处理功能，而且具有无以伦比的符号计算功能。Maple V的符号计算能力还是MathCAD和MATLAB等软件的符号处理的核心。Maple提供了2000余种数学函数，涉及范围包括：普通数学、高等数学、线性代数、数论、离散数学、图形学。它还提供了一套内置的编程语言，用户可以开发自己的应用程序，而且Maple自身的2000多种函数，基本上是用此语言开发的。 Maple采用字符行输入方式，输入时需要按照规定的格式输入，虽然与一般常见的数学格式不同，但灵活方便，也很容易理解。输出则可以选择字符方式和图形方式，产生的图形结果可以很方便地剪贴到Windows应用程序内。 二、MATLAB 系统 MATLAB原是矩阵实验室（Matrix Laboratory）在70年代用来提供Linpack和Eispac k软件包的接口程序，采用C语言编写。从80年代出现3.0的DOS版本，逐渐成为科技计算、视图交互系统和程序语言。MATLAB可以运行在十几个操作平台上，比较常见的有基于W indows 9X/NT、OS/2、Macintosh、Sun、Unix、Linux等平台的系统。 MATLAB程序主要由主程序和各种工具包组成，其中主程序包含数百个内部核心函数，工具包则包括复杂系统仿真、信号处理工具包、系统识别工具包、优化工具包、神经网络工具包、控制系统工具包、μ分析和综合工具包、样条工具包、符号数学工具包、图像处理工具包、统计工具包等。而且5.x版本还包含一套几十个的PDF文件，从MATLAB的使用入门到其他专题应用均有详细的介绍。 MATLAB是数值计算的先锋，它以矩阵作为基本数据单位，在应用线性代数、数理统计、自动控制、数字信号处理、动态系统仿真方面已经成为首选工具，同时也是科研工作人员和大学生、研究生进行科学研究的得力工具。MATLAB在输入方面也很方便，可以使用内部的E ditor或者其他任何字符处理器，同时它还可以与Word6.0/7.0结合在一起，在Word的页面里直接调用MATLAB的大部分功能，使Word具有特殊的计算能力。 三、MathCAD 系统 MathCAD是美国Mathsoft公司推出的一个交互式的数学系统软件。从早期的DOS下的1. 0和Windows下的4.0版本，到今日的8.0版本，功能也从简单的数值计算，直至引用Map le强大的符号计算能力，使得它发生了一个质的飞跃。 MathCAD是集文本编辑、数学计算、程序编辑和仿真于一体的软件。MathCAD7.0 Profe ssional（专业版）运行在Win9X/NT下，它的主要特点是输入格式与人们习惯的数学书写格式很近似，采用WYSWYG（所见所得）界面，特别适合一般无须进行复杂编程或要求比较特殊的计算。MathCAD 7.0 Professional 还带有一个程序编辑器，对于一般比较短小，或者要求计算速度比较低时，采用它也是可以的。这个程序编辑器的优点是语法特别简单。 MathCAD可以看作是一个功能强大的计算器，没有很复杂的规则；同时它也可以和Wor d、Lotus、WPS2000等字处理软件很好地配合使用，可以把它当作一个出色的全屏幕数学公式编辑器。 四、Mathematica 系统 Mathematica是由美国物理学家Stephen Wolfram领导的Wolfram Research开发的数学系统软件。它拥有强大的数值计算和符号计算能力，在这一方面与Maple类似，但它的符

实验五_MATLAB计算的可视化

实验五 MATLAB 计算的可视化（一） 实验目的 1. 熟练掌握MATLAB 二维曲线的绘制 2.掌握图形的修饰 3.掌握三维图的绘制 4.了解各种特殊图形的绘制 内容与步骤 1.在同一幅图形窗口中分别绘制y1=sin(t)和y2=cos(t)二条函数曲线，t 的取值范围为[0,10]。y1用红色虚线表示，y2用蓝色实线表示，横坐标轴名称为“时间t ”，纵坐标轴名称为“正弦、余弦”，整个图形的标题为“正弦和余弦曲线”。在坐标(1.7*pi ，-0.3)处添加文字“sin(t)”， 在坐标(1.6*pi ，0.8)处添加文字“cos(t)”，并在右上角添加图例，其运行界面图如下图所示。之后并尝试修改坐标轴刻度。 2.用subplot 命令在同一个窗口的不同子窗口绘制曲线y=sin(t)，y1=sin(t+0.25) y2=sin(t+0.5)，其中t=[0 10]。 3.绘制三维曲线：?? ? ??=≤≤==)cos()sin()200() cos()sin(t t t z t t y t x π （注意：用plot3命令） 4.三维网线图：绘制z=sin(y)cos(x) 三维网线图。 5. 三维曲面图 绘制22y x z +=的三维曲面图，x 在[-5,5]范围，y 在[-5,5]范围。将曲面图颜色用shading 命令连续变化，并用颜色标尺显示色图（使用函数colorbar 生成）。生成的图形如下图所示。

6.请绘制一个饼形图，数据如下表所示 7. 用semilogx命令绘制传递函数为1//(s+1)(0.5s+1)的对数幅频特性曲线，横坐标为w，纵坐标为Lw，w的范围为10-2-103，按对数分布。

科学计算与MATLAB语言(第四课)

第四讲绘图功能

作为一个功能强大的工具软件，Matlab 具有很强的图形处理功能，提供了大量的二维、三维图形函数。由于系统采用面向对象的技术和丰富的矩阵运算，所以在图形处理方面即常方便又高效。

4.1 二维图形 一、plot函数 函数格式：plot(x,y)其中x和y为坐标向量函数功能：以向量x、y为轴，绘制曲线。【例1】在区间0≤X≤2 内，绘制正弦曲线Y=SIN（X），其程序为： x=0:pi/100:2*pi; y=sin(x); plot(x,y)

一、plot函数 【例2】同时绘制正、余弦两条曲线Y1=SIN（X）和Y2=COS（X），其程序为： x=0:pi/100:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,x,y2) plot函数还可以为plot(x,y1,x,y2，x,y3，…)形式，其功能是以公共向量x为X轴，分别以y1，y2，y3，…为Y轴，在同一幅图内绘制出多条曲线。

一、plot函数 （一）线型与颜色 格式：plot(x,y1,’cs’,...) 其中c表示颜色，s表示线型。 【例3】用不同线型和颜色重新绘制例4.2图形，其程序为：x=0:pi/100:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,'go',x,y2,'b-.') 其中参数'go'和'b-.'表示图形的颜色和线型。g表示绿色，o表示图形线型为圆圈；b表示蓝色，-.表示图形线型为点划线。

一、plot函数 （二）图形标记 在绘制图形的同时，可以对图形加上一些说明，如图形名称、图形某一部分的含义、坐标说明等，将这些操作称为添加图形标记。 title(‘加图形标题'); xlabel('加X轴标记'); ylabel('加Y轴标记'); text(X,Y,'添加文本');

流体力学简单计算MATLAB程式

用matlab进行编程计算 第一问： z=30;p1=50*9.8*10^4;p2=2*9.8*10^4;jdc=0.00015;gama=9800;d=0.257;L=50000 ;mu=6*10^(-6); hf=z+(p1-p2)/(0.86*gama) xdc=2*jdc/d; beta=4.15;m=1; Q=(hf*d^(5-m)/(beta*mu^m*L))^(1/(2-m)); v=4*Q/(pi*d^2); Re=v*d/mu; Re1=59.7/xdc^(8*xdc/7)； Re2=(665-765*log(xdc))/xdc； i=hf/L; if Re<3000 Q=Q; elseif 3000

《Matlab与科学计算》作业 2010010099

《Matlab与科学计算》作业 第一章MATLAB环境 1、MATLAB通用操作界面窗口包括哪些？命令窗口、历史命令窗口、当前目录窗口、工作空间窗口各有哪些功能？ 答：MATLAB通用操作界面窗口包括：命令窗口、历史命令窗口、当前目录浏览器窗口、工作空间窗口、变量编辑器窗口、M文件编辑/调试器窗口、程序性能剖析窗口、MATLAB帮助。 命令窗口是MATLAB命令操作的最主要窗口，可以把命令窗口当做高级的“草稿纸”。在命令窗口中可以输入各种MATLAB的命令、函数和表达式，并显示除图形外的所有运算结果。 历史命令窗口用来记录并显示已经运行过的命令、函数和表达式，并允许用户对它们进行选择、复制和重运行，用户可以方便地输入和修改命令，选择多行命令以产生M文件。 当前目录窗口用来设置当前目录，可以随时显示当前目录下的M、MKL等文件的信息，扬文件类型、文件名、最后个修改时间和文件的说明信息等，并可以复制、编辑和运行M文件及装载MAT数据文件。 工作空间窗口用来显示所有MATLAB工作空间中的变量名、数据结构、类型、大小和字节数。 2、熟悉课本中表格1.4、1.5、1.6、1.7、1.8的内容。 3、如何生成数据文件？如何把数据文件中的相关内容输入到工作空间中，用实例进行操作。 生成数据文件：

把数据文件中的相关内容输入到工作空间中： 结果： 4、在工作空间中可以通过哪些命令管理变量，写出每种语法的具体操作过程。答：（1）把工作空间中的数据存放到MAT数据文件。 语法：save filename 变量1 变量2 ……参数。 （2）从数据文件中取出变量存放到工作空间。 语法：load filename 变量1 变量2 ……。

【原创】MATLAB实验报告-第二次-用MATLAB实现计算数据可视化-北京交通大学

MATLAB 上机实验报告( 2 ) 实验内容: 一、试用如下几种方法来建立向量，观察结果 ( 1) x=1:5, x=(1:5) ' 实验结果：x=1:5是行向量,x=(1:5)是列向量.且1为初始值，5为终止值，默认的步长为 1. >> x=1:5 1 2 3 4 5 >> x=(1:5)' x = 1 2

3 4 5 ( 2) x=0:pi/4:pi 实验结果：x=0:pi/4:pi 指的是x=(0,0.25*pi,0.50*pi,0.75*pi,pi). 其中pi为圆周率初始值为0,终止值为pi,步长为pi/4. >> x=0:pi/4:pi x = 0 0.7854 1.5708 2.3562 3.1416 (3)x=(0：0.2：3) ', y=e-x)p.(*sin(x) 实验结果：x的初始值为0,终止值为3,步长为0.2.而函数y表示将x向量中的每一个数代入函数y=e%x)*sin(x)得到的函数值组成的向量. >> x=(0:0.2:3)', y=exp(-x).*sin(x)

x = 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 2.0000 2.2000 2.4000 2.6000 2.8000 3.0000

0.1627 0.2610 0.3099 0.3223 0.3096 0.2807 0.2430 0.2018 0.1610 0.1231 0.0896 0.0613 0.0383 0.0204 0.0070 (4) k=linspace(-pi,pi,5), k=logspace(-3,-1,5) 实验结果：k=linspace(-pi,pi,5)产生的是初始值为-pi，终止值为 pi,元素总数为5的行向量，即k的步长为pi/2. k=logspace(-3,-1,5产生的是初始值为10八(-3),终止值为10八(-1),元素总数为5的列向量.

《MATLAB与科学计算》期末论文

盐城师范学院《MATLAB与科学计算》期末论文 2016－2017学年度第一学期 用MATLAB解决解析几何的图形问题 学生姓名吴梦成 学院数学与统计学院 专业信息与计算科学 班级数15(5)信计 学号 15213542

用MATLAB 解决解析几何的图形问题 摘 要 将 MATLAB 的图形和动画功能都用于解析几何教学，可使教学形象生动。以图形问题为例，详细给出了实例的程序编写和动画实现过程 。在解析几何教学中有一定的应用价值。 【关键词】： MATLAB ； 解析几何 ；图形 ； 动 画；编程 1 引 言 在解析几何的教学中，使用传统的教学方法。许多曲线及曲面的形成过程与变换过程只通过传统的教师讲授静态图示就很难形象生动地表示出来 。在解析几何教学中使用MATLAB 软件辅助教学，不仅可以很容易绘制出复杂的立体图形。把曲线、曲面的形成和变化过程准确地模拟出来 ，而且还能够对它们进行翻转 、旋转 ，甚 至还能够轻而易举地实现图形的动画效果!这对提高教学效率和培养学生的空间想象能力可起到事半功倍的效果。下面结合实例从几个方面说明MATLAB 在解析几何画图方面的应用。 2 利用 MATLAB 绘制三维曲线 在空间解析几何中，各种曲线和曲面方程的建立都离不开图形 ，而空间曲线和曲面图形既难画又费时。借助MATLAB 的绘图功能 ，可以快捷 、 准确地绘出图形，使教学变得形象 、生动 。有利于学生观察三维空间图形的形状 ， 掌握图形的性质 。 一 般地 ，MATLAB 可用plot3，ezplot3，comet3等函数来各种三维曲线 。 例如画螺旋曲线的图形，其参数方程设为 ：t at cos x =，t b sin t y -=，ct =z 。使用 plot3语句画螺旋曲线图形的方法如下( 设a =2 ,b=4,c=3)： );*3),sin(*.*4),cos(*.*2(3;*10:50/:0t t t t t plot pi pi t -= MATLAB 用两条简单的语句就可以画出螺旋 曲线(图1)，但上述方法是静态的 ，为了体

实验4、matlab的计算可视化和GUI设计

p345 subplot(2,2,1) t1=0:0.1:2; y1=sin(2*pi*t1); plot(t1,y1); title('y=sin(2\pit)') 练习： subplot(2,2,2) t2=0:0.1:2; y2=[exp(-t2);exp(-2*t2);exp(-3*t2)]; plot(t2,y2) axis([0 2 -0.2 1.2]); title('y=e-t,y=e-2t,y=e-3t') 练习： subplot(2,2,3); t3=[0 1 1 2 2 3 4]; y3=[0 0 2 2 0 0 0]; plot(t3,y3); axis([0 4 -0.5 3]); title('脉冲信号') 练习： subplot(2,2,4); t4=0:0.1:2*pi; plot(sin(t4),cos(t4));

axis([-1.2 1.2 -1.2 1.2]); axis equal; title('圆') 练习： P346 x=0:0.1:20; zeta=0 y1=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y1) zeta=0.3; y2=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); hold on plot(x,y2,'r:') zeta=0.5; y3=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y3,'g*') zeta=0.707; y4=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y4,'m-') title('二阶系统曲线') legend('\zeta=0','\zeta=0.3','\zeta=0.5','\zeta=0. 707') grid on gtext('\zeta=0') gtext('\zeta=0.3') gtext('\zeta=0.5') gtext('\zeta=0.707') ginput(3) zeta = ans = 2.6037 0.9035 13.1106 2.0029 4.2166 1.0380 P347 h_fig=gcf h_axis=gca h_line1=gco h_title=get(gca,'title') h_text2=findobj(h_fig,'string','\zeta=0.3') h_fig = 1 h_axis = 151.0018 h_line1 = 1 h_title = 152.0018 h_text2 = Empty matrix: 0-by-1 set(h_line1,'linewidth',5)

Matlab与科学计算样题(加主观题答案)

Matlab 与科学计算考试样题（客观题） 1 下面的MATLAB 语句中正确的有： a) 2a ＝pi 。 b) record_1=3+4i c) a=2.0, d) c=1+6j 2. 已知水的黏度随温度的变化公式如下，其中a=0.03368，b=0.000221，计算温度t 为20，30，40度时的粘度分别是： 2 1at bt μμ=++0μ为0℃水的黏度，值为31.78510-?；a 、b 为常数，分别为0.03368、0.000221。 3. 请补充语句以画出如图所示的图形： [x,y]=meshgrid(-2:0.1:2, -2:0.1:2)。 Z=x.*exp(-x.^2-y.^2)。 。 a) Plot3(x,y,Z) b) plot3(x,y,Z) c) mesh(x,y,Z) d) plot3(x,y,z) 2 a) 0.4900 1.2501 0.8560 b) 0.8560 1.2501 0.4900 c) -0.6341 3.8189 -3.7749 d) 3.8189 -3.7749 2.8533 解释说明：

>> x=0.5:0.5:3.0。 >> y=[1.75,2.45,3.81,4.80,8.00,8.60]。 >> a=polyfit(x,y,2) a = 0.4900 1.2501 0.8560 >> x1=[0.5:0.25:3.0]。 >> y1=a(1)*x1.^2+a(2)*x1+a(3) >> plot(x,y,'*') >> hold on >> plot(x1,y1,'--r') 5. 求方程在 x=0.5附近的根. 21 x x += a) 0.6180 b) -1.1719e-25 c) －1 d) -1.6180 6. 用Newton-Cotes方法计算如下积分 1 5 x? （a）133.6625 (b)23.8600 (c) 87.9027 (d) -1.6180 7. y=ln(1+x)，求x=1时y" a) -0.25 b) 0.5 c) -0.6137 d) -1.6137 8.某公司用3台轧机来生产规格相同的铝合金薄板。取样测量薄板的 厚度，精确至‰厘M。得结果如下： 轧机1：0.236 0.238 0.248 0.245 0.243 轧机2：0.257 0.253 0.255 0.254 0.261 轧机3：0.258 0.264 0.259 0.267 0.262 计算方差分析结果，并判定各台轧机所生产的薄板的厚度有无显著的差异？ a) p＝1.3431e-005，没有显著差异。

matlab计算结果的可视化

第五讲计算结果的可视化 本节介绍MATLAB 的两种基本绘图功能：二维平面图形和三维立体图形。 5.1 二维平面图形 5.1.1 基本图形函数 plot 是绘制二维图形的最基本函数，它是针对向量或矩阵的列来绘制曲线的。也就是说，使用plot 函数之前，必须首先定义好曲线上每一点的x 及y 坐标，常用格式为：（1）plot(x) 当x 为一向量时，以x 元素的值为纵坐标，x 的序号为横坐标值绘制 曲线。当x 为一实矩阵时，则以其序号为横坐标，按列绘制每列元素值相对于其序号的曲线， 当x 为m× n 矩阵时，就由n 条曲线。 （2）plot(x,y) 以x 元素为横坐标值，y 元素为纵坐标值绘制曲线。 （3）plot(x,y1,x,y2,…) 以公共的x 元素为横坐标值，以y1,y2,… 元素为纵坐标值绘 制多条曲线。 例5.1.1 画出一条正弦曲线和一条余弦曲线。 >> x=0:pi/10:2*pi; >> y1=sin(x); >> y2=cos(x); >> plot(x,y1,x,y2) 图5.1.1 函数plot 绘制的正弦曲线 在绘制曲线图形时，常常采用多种颜色或线型来区分不同的数据组，MATLAB 软件专 门提供了这方面的参数选项（见表5.1.1），我们只要在每个坐标后加上相关字符串，就可实 现它们的功能。 - 2 - 表5.1.1 绘图参数表 色彩字符颜色线型字符线型格式标记符号数据点形式标记符号数据点形式 y 黄- 实线. 点<小于号 m 紫：点线o 圆s 正方形 c 青-. 点划线x 叉号 d 菱形 r 红- - 虚线+ 加号h 六角星 g 绿* 星号p 五角星 b 蓝v 向下三角形 w 白^ 向上三角形 k 黑>大于号 例如，在上例中输入 >> plot(x,y1,'r+-',x,y2,'k*:') 则得图5.1.2 图5.1.2 使用不同标记的plot 函数绘制的正弦曲线 5.1.2 图形修饰 MATLAB 软件为用户提供了一些特殊的图形函数，用于修饰已经绘制好的图形。 表5.1.2 图形修饰函数表

中国矿业大学 实验六 MATLAB数据可视化

实验六MATLAB数据可视化 一、实验目的 掌握MATLAB 二维、三维图形绘制，掌握图形属性的设置和图形修饰；掌握图像文件的读取和显示。 二、实验内容 (1) 二维图形绘制。 (2) 三维曲线和三维曲面绘制。 三、实验步骤 1．二维图形绘制 (1) 二维图形绘制主要使用函数plot。 >> clear all; >> x=linspace(0,2*pi,100); >> y1=sin(x); >> plot(x,y) >> hold on >> y2=cos(x) >> plot(x,y) >> hold off

注：hold on 用于保持图形窗口中原有的图形，hold off解除保持。 (2) 函数plot 的参数也可以是矩阵。 >> close all >> x=linspace(0,2*pi,100); >> y1=sin(x); >> y2=cos(x); >> A=[y1 ; y2]'; >> B=[x ; x]' >> plot(B,A)

(3) 选用绘图线形和颜色。>> close all >> plot(x,y1,'g+',x,y2, 'r:') >> grid on

(4) 添加文字标注。 >> title('正弦曲线和余弦曲线') >> ylabel('幅度') >> xlabel('时间') >> legend('sin(x)', 'cos(x)') >> gtext('\leftarrowsinx')

(5) 修改坐标轴范围。 >> axis equal >> axis normal >> axis([0 pi 0 1.5]) 程序如下： x=linspace(0,2*pi,100); y1=sin(x); y2=cos(x); A=[y1 ; y2]'; B=[x ; x]' plot(B,A) plot(x,y1,'g+',x,y2, 'r:') axis equal axis normal axis([0 pi 0 1.5])

matlab实验 数据可视化方法

实验四数据可视化方法 [实验内容] 一．仿照运行，体会数据可视化方法。 1已知n=0,1,……,12，y=，运行下面程序，体会离散数据可视化方法。 说明： ·plot和stem指令均可以实现离散数据的可视化，但通常plot更常用于连续函数中特殊点的标记；而stem广泛运行与数字信号处理中离散点的图示。·用户在运行上面例程时会发现在命令窗口出现警告：Warning: Divide by zero！即警告程序中出现非零数除以0的指令。MATLAB对于这种情况并不中止程序，只是给该项赋值为inf以做标记。 2．下面时用图形表示连续调制波形y=sin（t）sin（9t），仿照运行，分析表现形式不同的原因。 二．编程实现。 1．用图形表示连续调制波形y=sin(t)sin(9t)，过零点及其包络线，如下图所示。

2. 编写函数[x,n]=stepseq(n0,n1,n2),实现： u（n）=, n为整数 并编写脚本文件实现: x(n)=n·[u(n)-u(n-10)]+10[u(n-10)-u(n-20)], 0≤n≤20要求在脚 本文件中调用stepseq 函数，最后绘出序列x(n)在给定区间的波形图。 3．编写一个函数文件[y,n]=sigadd(x1,n1,x2,n2)，实现两个对应样本之间的相加，其中x1是长度为n1的序列,x2是长度为n2的序列，n1、n2分 别是x1、x2的位置信息(n1、n2均为整数)，如: n1={ -3,-2,-1,0,1,2,3,4}，对应的 x1={ 2, 3, 1,4,1,3,1,2}； n2={-4,-3,-2,-1,0,1,2}，对应的 x2={ 1, 3, 2, 5,1,3,4}。 当调用函数[y,n]=sigadd(x1,n1,x2,n2)时，我们应该得到： n={-4,-3,-2,-1,0,1,2,3,4}，对应的 y={ 1, 5, 5, 6,5,4,7,1,2}。 仔细观察 sigadd 函数的功能。编好函数文件程序后，请在命令窗口调用，验证正确性，记录验证结果。 [实验结果] 一．1.

Matlab与工程计算 第二章 Matlab矩阵及其运算

第2章Matlab矩阵及其运算 2.1 Matlab变量 2.2 Matlab数值矩阵 2.3 运算符 2.4 基本数学函数 2.5 稀疏矩阵 2.6 矩阵分析 2.8 字符串 2.9 结构数据 2.10 细胞矩阵

2.1 Matlab变量 1. 变量命名规则 在MATLAB 6.5中，变量名是以字母开头，后接字母、数字或下划线的字符序列，最多63个字符。在MATLAB中，变量名区分字母的大小写。 2．变量赋值 变量=表达式 3．预定义变量 i，j，pi，eps，realmin，realmax，inf，NaN 预定义变量有特定的含义，在使用时，应尽量避免对这些变量重新赋值。

内存变量的管理 1．指令操作法 who whos clear 2. 现场菜单操作法 3. 内存变量文件(.mat) save [文件名] [变量名表] [-append][-ascii] load [文件名] [变量名表] [-ascii] help save help load

数值数据的输出格式 MATLAB用十进制数表示一个常数，具体可采用日常记数法和科学记数法两种表示方法。 在一般情况下，MATLAB内部每一个数值数据元素都是用双精度数来表示和存储的。 数据输出时用户可以用format命令设置或改变数据输出格式。format 命令的格式为： format格式符 其中格式符决定数据的输出格式 help format

2.2 MATLAB数值矩阵 2.2.1 矩阵的建立 1．直接输入法 最简单的建立矩阵的方法是从键盘直接输入矩阵的元素。具体方法如下：将矩阵的元素用方括号括起来，按矩阵行的顺序输入各元素，同一行的各元素之间用空格或逗号分隔，不同行的元素之间用分号分隔。

科学计算与MATLAB

单元测验已完成成绩：100.0分 1 【单选题】 MATLAB一词来自（）的缩写。 ?A、 Mathematica Laboratory ?B、 Matrix Laboratory ?C、 MathWorks Lab ?D、 Matrices Lab 我的答案：B得分：50.0分 2 【单选题】 下列选项中能反应MATLAB特点的是（）。 ?A、 算法最优 ?B、 不需要写程序 ?C、 程序执行效率高 . .

?D、 编程效率高 我的答案：D得分：50.0分 单元测验已完成成绩：96.4分 1 【单选题】 当在命令行窗口执行命令时，如果不想立即在命令行窗口中输出结果，可以在命令后加上（）。 ?A、 冒号（:） ?B、 逗号（,） ?C、 分号（;） ?D、 百分号（%） 我的答案：C得分：7.1分 2 【单选题】 fix(264/100)+mod(264,10)*10的值是（）。 ?A、 86 . .

?B、 62 ?C、 423 ?D、 42 我的答案：D得分：7.1分 3 【单选题】 在命令行窗口输入下列命令后，x的值是（）。 >> clear >> x=i*j ?A、 不确定 ?B、 -1 ?C、 1 ?D、 i*j 我的答案：B得分：7.1分 . .

4 【单选题】 使用语句x=linspace(0,pi,6)生成的是（）个元素的向量。 ?A、 8 ?B、 7 ?C、 6 ?D、 5 我的答案：C得分：7.1分 5 【单选题】 ceil(-2.1)的结果为（）。 ?A、 -2 ?B、 -3 ?C、 1 . .

?D、 2 我的答案：A得分：7.1分 6 【单选题】 eval('sqrt(4)+2')的值是（）。 ?A、 sqrt(4)+2 ?B、 4 ?C、 2 ?D、 2+2 我的答案：B得分：7.1分 7 【单选题】 已知a为3×5矩阵，则执行完a(:,[2,4])=[]后（）。 ?A、 a变成行向量 ?B、 a变为3行2列 . .

实验二MATLAB计算的可视化

课程实验报告 学年学期2011-2012学年第1学期 课程名称MATLAB与科学计算 实验名称实验二MATLAB计算的可视化实验室测量测绘实验中心计算机室专业年级热动113 学生姓名白治朋 学生学号2011012106 提交时间2013年10月23日 成绩 任课教师许景辉、牛亚斌 水利与建筑工程学院

实验二 MATLAB 计算的可视化 1、目的和要求 （1）熟练掌握MATLAB 二维曲线、三维图形的绘制。 （2）熟练掌握各种特殊图形的绘制。 （3）熟练掌握三维图形绘制命令。 （4）了解GUI 设计的一般过程和方法。 2、内容和步骤 参见教材实验四。 3、实验报告提交要求 （1） x=[1 2 3]，y=[1 2;2 3;5 8]，z=[2 6 9;3 8 8;1 5 7]，绘制plot （x ，y ）、plot （x ，z ），说明其各 自绘制的内容。 （2） 绘制如下图形，建立figure （2），绘图同样曲线，但标题为“你的姓名（黑体，16号字）”， 在x 坐标和y 坐标上分别标识学号和班级名称，并将网格线打开。 数组X 的列个数与矩阵y 的行个数相同， plot （ x ，y ）绘制的是x 为横坐标y 的每列为纵坐标的图像。如图1。 图1 数组X 的列个数与方阵z 的行列个数相同，plot （x ，z ）绘制的是x 为横坐标z 的每列为纵坐标的图像。如图2。 图2

（3）演示P133页，例题4.17 。

（4）完成课本P336图S 4.1实验，并用.m文件显示其程序内容。 （5）完成P302第四章例题4.

（6）通过绘制二阶系统阶跃响应，综合演示图形标识，请注释每条命令的含义。 clf; %清除图形窗口 t=6*pi*(0:100)/100;y=1-exp(-0.3*t).*cos(0.7*t); % 数据准备 tt=t(find(abs(y-1)>0.05)); %找出符合条件(y-1)的绝对值>0.05的对应t，赋值给tt ts=max(tt); %ts为tt中最大值ts=9.6133 plot(t,y,'r-','LineWidth',3) %画曲线t-y,红色实线，线粗3磅 axis([-inf,6*pi,0.6,inf]) %设置坐标轴范围。x轴下限自动产生,上限为6*pi；y轴下限0.6,上限自动产生 set(gca,'Xtick',[2*pi,4*pi,6*pi],'Ytick',[0.95,1,1.05,max(y)]) %二维坐标刻度设置。x轴刻度线取2*pi,4*pi,6*pi,y轴取0.95,1,1.05,max(y) grid on %显示坐标刻度线 title('\it y = 1 - e^{ -\alphat}cos{\omegat}') %用斜体1书写图名 text(13.5,1.2,'\fontsize{12}{\alpha}=0.3') %图形标识，添加文字注释。在x=13.5,y=1.2处,字体大小12磅，标注α=0.3 text(13.5,1.1,'\fontsize{12}{\omega}=0.7') %图形标识，添加文字注释。在x=13.5,y=1.1处,字体大小12磅，标注ω=0.7 hold on; %保持原有图形 plot(ts,0.95,'bo','MarkerSize',10); %在x=ts，y=0.95处画蓝色的空心圆圈,大小为

第6章MATLAB计算结果可视化讲解

第六章MATLAB 计算结果可视化 6.1连续函数和离散函数的可视化 【例6-1】用图形表示离散函数1 ) 6(--=n y 。 n=0:12; %产生一组自变量数据 y=1./abs(n-6); %计算相应点的函数值 plot(n,y,'r*','MarkerSize',20) %用红花标出数据点 grid on %画坐标方格 【例6-2】用图形表示连续调制波形)9sin()sin(t t y =。 t1=(0:11)/11*pi; y1=sin(t1).*sin(9*t1); t2=(0:100)/100*pi; y2=sin(t2).*sin(9*t2); subplot(2,2,1),plot(t1,y1,'r.'),axis([0,pi,-1,1]),title('子图 (1)') subplot(2,2,2),plot(t2,y2,'r.'),axis([0,pi,-1,1]),title('子图 (2)') subplot(2,2,3),plot(t1,y1,t1,y1,'r.') axis([0,pi,-1,1]),title('子图 (3)') subplot(2,2,4),plot(t2,y2)

6.2二维曲线绘图的基本操作 6.2.1 plot 的基本调用格式 【例6-3】用图形表示连续调制波形)9sin()sin(t t y 及其包络线。 t=(0:pi/100:pi)'; %长度为101的时间采样列向量 y1=sin(t)*[1,-1]; %包络线函数值，是（101x2）的矩阵 y2=sin(t).*sin(9*t); %长度为101的调制波列向量 t3=pi*(0:9)/9; y3=sin(t3).*sin(9*t3);plot(t,y1,'r:',t,y2,'b',t3,y3,'bo') 6.2.2泛函绘图指令fplot 【例6-4】fplot 与一般绘图指令的绘图效果比较。 [x,y]=fplot('cos(tan(pi*x))',[-0.4,1.4],0.2e-3);n=length(x); subplot(1,2,1),plot(x,y) title('\fontsize{20}\fontname{隶书}泛函绘图指令效果') t=(-0.4:1.8/n:1.4)'; subplot(1,2,2),plot(t,cos(tan(pi*t))) 6.2.3曲线的色彩、线型和数据点形 【例6-5】用图形演示平面上一个方块四个顶点在仿射投影（Affine Projection ）下的位置、 形状变化。 %平面上的四个点和它们构成的方块