实验二--用matlab绘制一元函数与二元函数的图象
实验2 Matlab绘图操作
实验2 Matlab 绘图操作实验目的:1、 掌握绘制二维图形的常用函数;2、 掌握绘制三维图形的常用函数;3、 掌握绘制图形的辅助操作。
实验内容:1. 设sin .cos x y x x ⎡⎤=+⎢⎥+⎣⎦23051,在x=0~2π区间取101点,绘制函数的曲线。
2. 已知: y x =21,cos()y x =22,y y y =⨯312,完成下列操作:(1) 在同一坐标系下用不同的颜色和线性绘制三条曲线; (2) 以子图形式绘制三条曲线;(3) 分别用条形图、阶梯图、杆图和填充图绘制三条曲线。
3.已知:ln(x x e y x x ⎧+≤⎪⎪=⎨⎪+>⎪⎩20102,在x -≤≤55区间绘制函数曲线。
4. 绘制极坐标曲线sin()a b n ρθ=+,并分析参数a 、b 、n 对曲线形状的影响。
5.在xy 平面内选择区域[][],,-⨯-8888,绘制函数z =6. 用plot 函数绘制下面分段函数的曲线。
,(),,x x f x x x x ⎧++>⎪==⎨⎪+<⎩235000507. 某工厂2005年度各季度产值(单位:万元)分别为:450.6、395.9、410.2、450.9,试绘制柱形图和饼图,并说明图形的实际意义。
8. 在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点。
(1).y x =-205(2)sin()cos ,sin()sin x t tt y t tπ=⎧≤≤⎨=⎩303详细实验内容:1.设sin .cos x y x x ⎡⎤=+⎢⎥+⎣⎦23051,在x=0~2π区间取101点,绘制函数的曲线。
>> x=(0:2*pi/100:2*pi);>> y=(0.5+3*sin(x)/(1+x.^2))*cos(x); >> plot(x,y)2.已知: y x =21,cos()y x =22,y y y =⨯312,完成下列操作: (1)在同一坐标系下用不同的颜色和线性绘制三条曲线; >> x= linspace(0, 2*pi, 101);>> y1=x.*x; >> y2=cos(2x); >> y3=y1.*y2;plot(x,y1,'r:',x,y2,'b',x,y3, 'ko')(2)以子图形式绘制三条曲线;>> subplot(2,2,1),plot(x,y1)subplot(2,2,2),plot(x,y2)subplot(2,2,3),plot(x,y3)(3)分别用条形图、阶梯图、杆图和填充图绘制三条曲线。
MATLAB绘图&simulink
1.2 绘制多根二维曲线 1.plot函数的输入参数是矩阵形式 (1) 当x是向量,y是有一维与x同维的矩阵时,则绘 制出多根不同颜色的曲线。曲线条数等于y矩阵的 另一维数,x被作为这些曲线共同的横坐标。 (2) 当x,y是同维矩阵时,则以x,y对应列元素为横、 纵坐标分别绘制曲线,曲线条数等于矩阵的列数。
绘图
线型和颜色控制符
. 点 o(字母)小圆圈 x x标记 + 加号 实线 * 星号 : 虚线 -. 点划线 -间断线 y m c r g b w k 黄 棕色 青色 红色 绿色 蓝色 白色 黑色
1.4 图形标注与坐标控制 1.图形标注 有关图形标注函数的调用格式为: title(图形名称) xlabel(x轴说明) ylabel(y轴说明) text(x,y,图形说明) legend(图例1,图例2,…)
例 6 在同一坐标内,分别用不同线型和颜色绘制曲 线y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx),标记两 曲线交叉点。 程序如下: x=linspace(0,2*pi,1000); y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=2*exp(-0.5*x).*cos(pi*x); k=find(abs(y1-y2)<1e-2); %查找y1与y2相等点 (近似相等)的下标 x1=x(k); %取y1与y2相等点的x坐标 y3=0.2*exp(-0.5*x1).*cos(4*pi*x1); %求y1与y2值 相等点的y坐标 plot(x,y1,x,y2,'k:',x1,y3,'bp');
MATLAB提供的统计分析绘图函数还有很多,例如, 用来表示各元素占总和的百分比的饼图、复数的 相量图等等。 例 14 绘制图形: (1) 某企业全年各季度的产值(单位:万元)分别为: 2347,1827,2043,3025,试用饼图作统计分析。 (2) 绘制复数的相量图:7+2.9i、2-3i和-1. 6i。
Matlab中函数图形的三种绘制方法
Matlab中函数图形的三种绘制方法及局部和全局解
绘制函数的图形2
x
fπ区间[-1,2]
=x
)
10
sin(+
1 利用plot绘制
x=linspace(-1,2,1000);
y=x.*sin(10*pi*x)+2;
plot(x,y)
/ 函数的显式表达式,先设置自变量向量,然后根据表达式计算出函数向量/
2 利用fplot绘制
f='x.*sin(10*pi*x)+2';或f='x*sin(10*pi*x)+2';
fplot(f,[-1,2],1e-4)
/fplot函数可以自适应地对函数进行采样,能更好地反映函数的变化规律/
3 利用ezplot绘制
f='x*sin(10*pi*x)+2';
ezplot(f,[-1,2])
/隐函数绘图:如果函数用隐函数形式给出,可以利用ezplot函数绘制隐函数图形/
1和2 3的区别是2 3可以直接按照函数的原形直接写出,而1中变量相乘或除时都以点乘和点除的形式写出来的
尝试用fminbnd fminunc fminsearch及遗传算法求解上述函数在区间[-1,2]中的最小值,看看它们四个有什么不同?。
Mathematica绘图部分
数值 x i 带入求 xidxi ,当然出错.
解决的办法是使用函数Evaluate[f],告
知Mathematica首先求出表达式f的值.
例3:
2 可选参数
绘图函数的可选参数很多,一下介绍Plot 的常用可选参数.
可选参数分为两类:第一类参数能改变输 出图形的外观,但不影响图形自身的质量;第 二类参数则影响图形自身的质量.
可选参数的格式为:可选项名->可选项值, 当不使用可选参数时该参数去默认值.
1>第一类可选参数 第一类可选参数有以下几种:
<1> PlotRang 指定绘图的范围.其可选值是:
➢Automatic 由Mathematica自动选取范围 ➢ 切除无穷值点和尖峰 ➢All 画出所有点 ➢{min,max} 给出y<三维为z>轴方向的取 值 ➢ 范围
➢{x{1,x2}{,y1,y2}分}别给出x,y<三维
加z>轴方向的取值范围
例4:
<2> AspectRatio 指定图形的高宽比.
可选值是:默认值为0.618<即黄金分割>, 准确值是1/GoldenRatio,其中
GoldenoRa1ti 5 2
是一个Mathematica常数.如果取Automatic, 则高宽比为1,还可以取任何正数.
➢RGBColor[r,g,b] 其中r,g,b是0到1之间的数, 分别表示红,绿,蓝色的强度,[1,1,1]为白色,
[0,0,0]为黑色,[1;11> PlotLabel 用于在图形上方居中加注释. 其值为:
➢None 没有注释<默认值> ➢"字符串" 将字符串里的内容作为注释 ➢
matlabplot二元函数
一、介绍Matlab是一款用于科学计算和技术计算的软件,拥有强大的绘图功能,可以用于绘制各种图形、曲线和函数图像。
在Matlab中,使用plot函数可以绘制二元函数的图像,通过调整参数和设置属性,可以实现不同风格和效果的图像展示。
本文将详细介绍在Matlab中如何绘制二元函数的图像,包括基本的绘图方法和常用的设置技巧。
二、绘制二元函数图像的基本方法1. 准备数据在使用Matlab绘制二元函数的图像前,首先需要准备数据。
通常可以通过生成x、y坐标的网格点,然后计算每个点对应的函数值,从而得到二元函数在指定区域内的数据集。
2. 使用plot函数一旦准备好了数据集,就可以使用Matlab的plot函数进行绘图。
plot函数的基本语法为:plot(x, y),其中x和y分别代表要绘制的点的横坐标和纵坐标。
通过调用plot函数,可以将计算得到的数据点连接起来,形成二元函数的图像。
3. 添加标签和标题为了让图像更加清晰和直观,通常需要添加x和y轴的标签,以及整个图像的标题。
在Matlab中,可以使用xlabel、ylabel和title函数来分别添加x轴、y轴和标题标签。
4. 设置图像属性通过设置图像的属性,可以调整图像的风格和效果。
常用的属性包括线型、线宽、颜色和标记符号等。
在Matlab中,可以使用参数-字符串对的形式来设置图像的属性,例如'LineStyle'、'LineWidth'、'Color'和'Marker'等。
三、绘制常见二元函数的图像在Matlab中,可以绘制各种类型的二元函数的图像,包括线性函数、二次函数、三角函数、指数函数、对数函数等。
下面将分别介绍如何绘制这些常见二元函数的图像。
1. 绘制线性函数图像线性函数的一般形式为y=ax+b,其中a和b分别为常数。
在Matlab 中,可以通过设置a和b的值,然后使用plot函数绘制线性函数的图像。
用MATLAB绘制一元函数和二元函数的图象
《MATLAB语言》课程论文用MATLA B绘制一元函数和二元函数的图象姓名:**学号: ********** 5专业:通信工程班级: 2010级通信1班指导老师:***学院:物理电气信息学院完成日期:2011.12.20用MATL A B 绘制一元函数和二元函数的图像(马军 1201024524 5 2010级通信工程1班)【摘要】大学物理力学中涉及许多复杂的数值计算问题,例如非线性问题,对其手工求解较为复杂,而MATL AB 语言正是处理非线性问题的很好工具,既能进行数值求解,又能绘制有关曲线,非常方便实用。
另外,利用其可减少工作量,节约时间,加深理解,同样可以培养应用能力。
【关键词】一元函数 二元函数 MATLA B 图像的绘制一、问题的提出MATLA B 语言是当今国际上科学界(尤其是自动控制领域)最具影响力、也是最有活力的软件。
它提供了强大的科学运算、灵活的程序设计流程、高质量的图形可视化与界面设计、便捷的与其他程序和语言接口的功能。
MATLA B 语言在各国高校与研究单位起着重大的作用.它是一种集数值计算、符号运算、可视化建模、仿真和图形处理等多种功能…二、实验内容1.平面曲线的表示形式对于平面曲线,常见的有三种表示形式,即以直角坐标方程],[),(b a x x f y ∈=,以参数方程],[),(),(b a t t y y t x x ∈==,和以极坐标],[),(b a r r ∈=ϕϕ表示等三种形式。
2.曲线绘图的M ATLA B 命令MATLA B 中主要用p lot,fplot 二种命令绘制不同的曲线。
plot(x,y) 作出以数据(x (i),y(i))为节点的折线图,其中x,y 为同维数的向量。
plot(x1,y1,x2,y2,…) 作出多组数据折线图 fplot (‘fun’,[a,b]) 作出函数f u n 在区间[a,b]上的函数图。
matlab画二元函数(PDF精品)
图4-4 图4-5 数学实验四 用Matlab 软件作二元函数的图象一、空间曲线调用格式:plot3(x1,y1,z1,x2,y2,z2,…,xn,yn,zn)其中xi 、yi 和zi 为分别是长度相同的向量,用于存储xi 坐标、yi 坐标和zi 坐标数据,i=1,2,…,n .例1 在同一坐标系中分别作出三维曲线 ==x x z x x y cos sin 和+=−=2442x z x y 在]6,0[π∈x 上的图象.解 打开M文件编辑窗口,在其中输入下面命令集:x=0:pi/50:6*pi;y1=x.*sin(x);z1=x.*cos(x);y2=2*x-1;z2=4*x-8;plot3(x,y1,z1,x,y2,z2);xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')legend('圆锥螺线','空间直线')取名为exa1保存,再在命令窗口中输入命令exa1,程序运行结果如图4-4.说明:若要作单根曲线,和前面plot 函数一样,只是多了一个立坐标.二、空间曲面1.三维网格图调用格式一:[x,y]=meshgrid(x1:h1:x2,y1:h2:y2)z=f(x,y)plot3(x,y,z)其中x 、y 和z 为长度相同的向量,分别用于存储x 坐标、y 坐标和z 坐标数据. 例2 作出空间平面324=+−z y x 的图象.解 打开M文件编辑窗口,在其中输入下面命令集:[x,y]=meshgrid(-3:0.1:3,-2:0.1:3); %产生一个x,y 平面上51×61的网格z=3-4*x+2*y;plot3(x,y ,z);xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')title('4x-2y+z=3')取名为exa2保存,再在命令窗口中输入命令exa2,程序运行结果如图4-5.例3 作出空间曲面2222sin y x y x z ++=在区域图4-7图4-8}8||,8|||),{(≤≤y x y x 上的图象.解 打开M文件编辑窗口,在其中输入下面命令集:[x,y]=meshgrid(-8:0.1:8);r=sqrt(x.^2+y.^2); %产生对应的x,y 坐标用点幂z=sin(r)./r; %产生对应于x,y 的z 坐标用点除plot3(x,y,z);xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')取名为exa3保存,再在命令窗口中输入命令exa3,程序运行结果如图4-6.调用格式二:[x,y]=meshgrid(x1:h1:x2,y1:h2:y2)z=f(x,y)mesh(x,y,z)例4 在x 、y 平面内选取一个区域,作出空间曲面22y xe z −−=的图象. 解 打开M文件编辑窗口,在其中输入下面命令集:[x,y]=meshgrid(-4:0.1:3,-3:.1:2); %产生一个x,y 平面上51×71的网格z=exp(-x.^2-y.^2);mesh(x,y,z);xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')取名为exa2保存,再在命令窗口中输入命令exa2,程序运行结果如图4-7.2.三维曲面图调用格式:[x,y]=meshgrid(x1:h1:x2,y1:h2:y2)z=f(x,y)surf(x,y,z) 说明:surf 函数可借助shading 函数作平滑和插值处理,而shading 函数有三个参数,分别为flat (作平滑处理)、interp (去掉连接线条,在各片之间使用颜色插值)和faceted (默认值,对前面两种参数之一的作用进行还原).例5 用子图分别作出马鞍面22y x z −=在区域}2||,1|||),{(≤≤y x y x 图和经插值处理的图.解 打开M文件编辑窗口,在其中输入下面命令集:[x,y]=meshgrid(-1:0.1:1,-2:.1:2);z=x.^2-y.^2;subplot(1,2,1),surf(x,y,z);xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')title('z=x^2-y^2');subplot(1,2,2),surf(x,y,z);图4-9 图4-10xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')shading interp %插值处理title('z=x^2-y^2的插值处理图')取名为exa5保存,再在命令窗口中输入命令exa5,程序运行结果如图4-8.3.柱面图调用格式一:[x,y,z]=cylinder(R,N)mesh(x,y,z)其中,R 是一个向量,其坐标分量依次对应柱面各横截面的半径,N是多边形的边数.默认值是R=[1,1],N=20.例6 作出正六棱柱在R=[2,2]的图象.解 打开M文件编辑窗口,在其中输入下面命令集:[x,y,z]=cylinder([2,2],6);mesh(x,y,z);xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')hidden off %显示隐含线,若不想显示隐含线,则取消此命令或设置为hidden on取名为exa6保存,再在命令窗口中输入命令exa6,程序运行结果如图4-9. 调用格式二:[x,y,z]=cylinder(R,N)surf(x,y,z)其中,R 、N的意义和上面一样.例7 作出正十棱台在R=[1,2]的图象.解 打开M文件编辑窗口,在其中输入下面命令集:[x,y,z]=cylinder([1,2],10);surf(x,y,z);xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')取名为exa7保存,再在命令窗口中输入命令exa7,程序运行结果如图4-10.说明:若要作正棱锥,只须取R的某个分量为0即可.4.球面图调用格式:[x,y,z]=sphere(N)surf(x,y,z) 其中,N (N>2)是球面的边数.默认值是N=20.sphere(N)不返回坐标,直接绘出球面.例8 以子图形式分别作出N=5、20、40和40并作平滑处理的球.解 打开M文件编辑窗口,在其中输入下面命令集:subplot(2,2,1)sphere(5);title('子图1')subplot(2,2,2),sphere;title('子图2')subplot(2,2,3)[x,y,z]=sphere(40);surf(x,y,z);title('子图3')subplot(2,2,4)[x,y,z]=sphere(40);surf(x,y,z);shading flat %对球面作平滑处理title('子图4') 取名为exa8保存,再在命令窗口中输入命令exa8,程序运行结果如图4-11.三、上机实验1.用help 命令查看函数plot3,mesh 和surf 等的用法.2.上机验证上面各例.3.作相关小节练习中空间曲线和曲面的图象.图4-11。
Matlab数学实验2 一元函数极限
x 3 由于左右极限不等,我们就可以指出lim 不存在了。 x 3 | x 3 |
二重极限的命令调用格式: limit(limit(f,x,x0),y,y0)表示
例6.计算 lim ( x y) ln( x 2 y 2 )
( x , y ) (0,0)
( x , y )( x 0, y 0 )
如果我们绘制函数的图像,x 3这一点不连续是很明显的, 为了要显示图像的间断点,我们给Matlab指定了定义域:
也可以编写如下程序代码: x=[-5:0.1:5]; y=(x-3)./abs(x-3); plot(x,y) axis([-5 5 -2 2])
结果如下图所示:
由于x 3是间断点,所以当x 3时要求左右极限, 运行命令如下所示:
命令
limit(f) limit(f,x,a)或 limit(f,a) limit(f,x,inf)
功能
命令
limit(f,x,a, 'right') limit(f,x,a, 'left') limit(f,x, +inf) limit(f,x, -inf)
功能
lim f ( x)
x 0
xa
4、 解:输入命令: limit((cot(x))^(1/log(x)),x,0,'right') 得结果:ans=1/e1.实验目的
掌握函数极限的运算常用命令及常用函数
2.实验要求
熟练掌握有关limit命令的操作及应用
3. 实验内容
MATLAB求解极限的命令
☆ MATLAB使用syms这个函数命令来创建和定 义基本的符号对象.其调用格式为: syms Var1 Var2 … ☆一元函数极限的调用格式如下表所示:
matlab基础课程2(绘图)
四、MATLAB三维曲面绘图
meshgrid——生成网格矩阵 调用格式:
[X,Y]=meshgrid(x,y)--------生成小矩形顶点的坐
标值矩阵
[X,Y]=meshgrid(x) 等价于[X,Y]=meshgrid(x,x)
例:x=[1,2,3];y=[1,2,3,4]; [X,Y]=meshgrid(x,y)
xlabel(‘String’)
在当前图形的x轴旁边加入文字内容
ylabel(‘String’)
zlabel(‘String’)
在当前图形的y轴旁边加入文字内容
在当前图形的z轴旁边加入文字内容
图形标注gtext
gtext(‘String’)
在鼠标指定位置上标注
说明:使用gtext指令后,会在当前图形上出现一个十字 叉,等待用户选定位置进行标注。移动鼠标到所需位置按 下鼠标左键,Matlab就在选定位置标上文字。
马鞍面、平面及交线
练习: 二维曲线4-4,4-5 p42 例4-9 p47 例4-11,4-12
三维曲线:
p57 例4-22,4-23,4-24
三维曲面:
p60 例4-26,4-28,4-29,4-30
作业:P79 3,5,9
x=[1 2 3]; y=[3 3 1 473 581 4 4 4] plot(x,y)
plot3与plot的 用法相同
,
x t sin t 例:绘制三维曲线的图像: y t cos t z t
解:matlab命令为:
(0 t 20 )
t=0:pi/10:20*pi; x=t.*sin(t); y=t.*cos(t) ; z=t; plot3(x,y,z)
[整理]MATLAB一元二次函数的画图.
function varargout = huatu(varargin)% HUATU MATLAB code for huatu.fig% HUATU, by itself, creates a new HUATU or raises the existing% singleton*.%% H = HUATU returns the handle to a new HUATU or the handle to% the existing singleton*.%% HUATU('CALLBACK',hObject,eventData,handles,...) calls the local% function named CALLBACK in HUATU.M with the given input arguments. %% HUATU('Property','Value',...) creates a new HUATU or raises the% existing singleton*. Starting from the left, property value pairsare% applied to the GUI before huatu_OpeningFcn gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to huatu_OpeningFcn via varargin.%% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one % instance to run (singleton)".%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help huatu% Last Modified by GUIDE v2.5 15-Jan-2014 13:13:56% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name', mfilename, ...'gui_Singleton', gui_Singleton, ...'gui_OpeningFcn', @huatu_OpeningFcn, ...'gui_OutputFcn', @huatu_OutputFcn, ...'gui_LayoutFcn', [] , ...'gui_Callback', []);if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});endif nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});elsegui_mainfcn(gui_State, varargin{:});end% End initialization code - DO NOT EDIT% --- Executes just before huatu is made visible.function huatu_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to huatu (see VARARGIN)% Choose default command line output for huatuhandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes huatu wait for user response (see UIRESUME)% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line. function varargout = huatu_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT); % hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structurevarargout{1} = handles.output;% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)a=str2num(get(handles.edit1,'string'));b=str2num(get(handles.edit2,'string'));c=str2num(get(handles.edit3,'string'));x=-10:1:10;hold ony=a*x.^2+b*x+c;plot(x,y);% --- Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)close allfunction edit1_Callback(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit1 as text% str2double(get(hObject,'String')) returns contents of edit1 as a double% --- Executes during object creation, after setting all properties. function edit1_CreateFcn(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit2_Callback(hObject, eventdata, handles)% hObject handle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit2 as text% str2double(get(hObject,'String')) returns contents of edit2 as a double% --- Executes during object creation, after setting all properties.function edit2_CreateFcn(hObject, eventdata, handles)% hObject handle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit3_Callback(hObject, eventdata, handles)% hObject handle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit3 as text% str2double(get(hObject,'String')) returns contents of edit3 as a double% --- Executes during object creation, after setting all properties. function edit3_CreateFcn(hObject, eventdata, handles)% hObject handle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');end% --- Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)% axex2cla(handles.axes2)。
MATLAB绘制函数动态图
创新思维训练指导老师:李安完成日期:2015年7月8日题目:利用MATLAB绘制函数动态图形目的:(1)初步了解MATLAB的界面以及简单应用;(2)学会用MATLAB绘制常见图形;(3)通过制图,进一步加深对函数的理解,观察函数的性质。
内容:1、绘制基本的一元函数动态图例:在一个坐标系中绘制函数y=x^2,y=x,y=cos(x)的图形,并通过movie绘制上述函数动画。
步骤:(1)选定题目:三个小组成员一起选定一个要做的创新思维题目;(2)想法集合:选好题目后,每个成员开动脑筋,说说自己的想法以及思路,然后将所有人的思路汇总在一起;(3)分工合作:三人分别负责一个模块,一人负责程序,一人负责调试,一人负责写作,三人分工明确,彼此独立,却又相互依存。
基本程序:clc;clear;x=-2:0.05:2;y1=x.^2;y2=x;y3=cos(x);figure(1)上述程序比较简单,只能实现三个函数的静态表现,而且图像未经任何修饰。
之后通过小组讨论,上网查阅,将程序进一步的修改和完善,最终达到预期效果。
最终程序如下:最终源码:clc;clear;x=-2:0.05:2;xlabel('x'),ylabel('y')title('y=x^2和y=x和y=cos(x)的图像','Color','b');text(-1.35,1.5,'\rightarrow y1=x^2');text(-1,-1,'\leftarrow y2=x');text(1.8,0,'\leftarrow y3=cos(x)');y1=x.^2;y2=x;y3=cos(x);figure(1)n=length(x);grid onhold onfor i=1:nplot(x(i),y1(i),'gp',x(i),y2(i),'r*',x(i),y3(i),'c+');axis([min(x)-1 max(x)+1 min(y1)-1 max(y1)+1]);axis([min(x)-1 max(x)+1 min(y2)-1 max(y2)+1]);axis([min(x)-1 max(x)+1 min(y3)-1 max(y3)+1]);M(i)=getframe;pause(0.1);endmovie(M);legend('\ity=x^2','\ity=x','\ity=cos(x)')运行程序后的结果:三条函数均能以动态点型绘出,其中函数y=x^2表现为绿色五星型,函数y=x表现为红色星号,函数y=cos(x)表现为青色十字型。
实验二MATLAB的基本绘图方法
实验二MATLAB的基本绘图方法一、实验目的1.二维平面图形的绘制2.三维立体图形的绘制3.隐函数作图二、实验地点:A204三、实验日期:四、实验内容(一)二维平面图形的绘制1、Plot的使用方法介绍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,…元素为纵坐标值绘制多条曲线。
例1:画出一条正弦曲线和一条余弦曲线。
x=0:pi/10:2*pi;y1=sin(x);y2=cos(x);plot(x,y1,x,y2)请尝试给plot()加上标记符参数,如:plot(x,y1,'*',x,y2,'^')。
注:在绘制曲线图形时,常常采用多种颜色或线型来区分不同的数据组,MA TLAB 软件专门提供了这方面的参数选项,我们只要在每个坐标后加上相关字符串,就可实现它们的功能。
具体参见课件。
2、图形修饰MATLAB 软件为用户提供了一些特殊的图形函数,用于修饰已经绘制好的图形。
图形修饰函数表如下:函数含义grid on (/off) 给当前图形标记添加(取消)网络xlable(‘string’) 标记横坐标ylabel(‘string’) 标记纵坐标title(‘string’) 给图形添加标题text(x,y,’string’) 在图形的任意位置增加说明性文本信息gtext(‘string’) 利用鼠标添加说明性文本信息axis([xmin xmax ymin ymax]) 设置坐标轴的最小最大值例2、给例1的图形中加入网络和标记。
第4章 MATLAB绘图
4.2 其他二维图形 4.2.1 二维统计分析图
在MATLAB中,二维统计分析图形很多,常 见的有条形图、阶梯图、杆图和填充图等, 所采用的函数分别是:
bar(x,y,选项) stairs(x,y,选项) stem(x,y,选项) fill(x1,y1,选项1,x2,y2,选项2,…)
例4-13 分别以条形图、阶梯图、杆图和填充图形式 绘制曲线y=2sin(x)。
程序如下:
x=linspace(0,2*pi,1000);
y1=0.2*exp(-0.5*x).*cos(4*pi*x);
y2=2*exp(-0.5*x).*cos(pi*x);
k=find(abs(y1-y2)<1e-2); %查找y1与y2相等点
(近似相等)的下标
x1=x(k);
%取y1与y2相等点的x坐标
例4-1 在0≤x≤2区间内,绘制曲线 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)
例4-2 绘制曲线。 程序如下: t=0:0.1:2*pi; x=t.*sin(3*t); y=t.*sin(t).*sin(t); plot(x,y);
例4-7 在0≤x≤2区间内,绘制曲线y1=2e-0.5x和y2=cos(4πx), 并给图形添加图形标注。
程序如下:
x=0:pi/100:2*pi;
y1=2*exp(-0.5*x);
y2=cos(4*pi*x);
plot(x,y1,x,y2)
title('x from 0 to 2{\pi}');
数据点 黑点 小圆圈 差号 十字标号 星号 小方块
Mathematica教程06函数作图
Dashing[{r1,r2,...}]交替使用数r1,r2,…作为线段和空白 的相对长度画虚线(其中r1,r2,…是远远小于1的数,整 个图形宽度为1).
可以对两条曲线分别指定参数值
Plot[Sin[x^2], {x, 0, 3}, Frame -> True] GridLinese用于加网格线。它的值为:
None不加网格线(默认值)。 Automatic 由Mathematics自动加上网格线。 {{x1,x2,...},{y1,y2,...}} 在横轴上的点x1,x2,...和纵轴上 点y1,y2,...处加上网格线
出多条曲线。
这个函数能添加与Plot一样的可选参数。
ParametricPlot[{Cos[t]^3, Sin[t]^3}, {t, 0, 2Pi}, AspectRatio -> Automatic]
Mathematica没有专门画极坐标图形的内部函数(图形扩展 程序包中有),因此应该先转换成参数式.
等值线图
绘制函数z=f(x,y)的等值线图使用函数: ContourPlot[f,{x,xmin,xmax},{y,ymin,ymax}] 其中f 是二元函
数的表达式. ContourPlot[x^2-y^2,{x,-1,1},{y,-1,1}] 说明:用灰度表示函数值的大小,越亮的地方函数值越大。 这个函数有以下几个可选参数: ①ContourShading用于决定是否使用灰度。
AspectRatio指定图形的高宽比。它的可选值是: 默认值为0.618(即黄金分割),准确值是1/GoldenRatio, 如果取Automatic,则高宽比为1,还可以取任何正数. 不设置这个参数,则圆变成椭圆,设此参数值为 Automatic则可解决问题。
教你用MATLAB快速作一维、二维、三维图
精选可编辑ppt
7
• 加图例
• 给图形加图例命令为legend。该命令把图例放置在图形空 白处,用户还可以通过鼠标移动图例,将其放到希望的位 置。
• 格式:legend('图例说明','图例说明');
• 为正弦、余弦曲线增加图例,其程序为:
• x=0:pi/100:2*pi;
• y1=sin(x);
• grid on 加网格线
• text(x,y,‘string’) adds the string in quotes(引号) to the location specified by the point (x,y).
• \bullet ·
• \pi
π
• \rightarrow 右箭头
• EdgeColor -- Color of the rectangle's edge (none by
例 在[-2,2]范围内绘制函数tanh的图形
解 fplot(‘tanh’,[-2,2])
Matlab liti28
例 x 、 y 的 取 值 范 围 都 在 2 [ - 2 , ] ,
画 函 数 t a n h ( x ) , s i n ( x ) , c o s ( x ) 的 图 形
解 输入命令:
--将多条精线选可画编辑在pp一t 起
3
例 在[0,2*pi]用红线画sin(x),用绿圈画cos(x).
解 x=linspace(0,2*pi,30);
y=sin(x); z=cos(x);
Matlab liti1
plot(x,y,'r',x,z,’g0')
1
0.8
MATLAB一元二次函数的画图
y=a*x A 2+b*x+c 画图 z. I■I.;.■■K 卜: ___________________________do £ 0 5fun cti on varargout = huatu(vararg in)% HUATU MATLAB code for huatu.fig% HUATU, by itself, creates a new HUATU or raises the exist ing% sin glet on*.%% H = HUATU returns the han die to a new HUATU or the han die to% the existi ng sin glet on*.%% HUATU('CALLBACK',hObject,eve ntData,ha ndles,...) calls the local% fun ctio n named CALLBACK in HUATU.M with the give n in put argume nts. %% HUATU(' Property','Value',...) creates a new HUATU or raises the% existi ng sin glet on*. Start ing from the left, property value pairsare% applied to the GUI before huatu_Ope nin gFc n gets called. An% un recog ni zed property n ame or inv alid value makes propertya-1 b 二2 c=1 二元一次函数的绘图applicati on% stop. All in puts are passed to huatu_Ope nin gFc n via vararg in.%% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one % in sta nee to run (sin glet on)".%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the resp onse to help huatu% Last Modified by GUIDE v2.5 15-Ja n-2014 13:13:56% Begin initialization code - DO NOT EDITgui_S in glet on = 1;gui_State = struct( 'gui_Name' , mfilename, ...'gui_S in glet on' , gui_S in glet on, …'gui_OpeningFcn' , @huatu_OpeningFcn,…'gui_OutputFc n' , @huatu_OutputFc n, … 'gui_LayoutFc n',[],...'gui_Callback',[]);if nargin && ischar(varargin{1}) gui_State.gui_Callback = str2fu nc(varargi n{1});endif nargout[varargout{1: nargout}] = gui_mai nfcn( gui_State, varargi n{:});elsegui_mai nfcn( gui_State, varargi n{:});end% End initialization code - DO NOT EDIT% --- Executes just before huatu is made visible.fun cti on huatu_Ope nin gFc n( hObject, eve ntdata, han dles, varargi n)% This fun cti on has no output args, see OutputFc n.% hObject han dle to figure% eventdata reserved - to be defined in a future version of MATLAB% han dies structure with han dies and user data (see GUIDATA)% varargin comma nd line argume nts to huatu (see VARARGIN)% Choose default comma nd line output for huatu han dles.output = hObject;% Update han dles structure guidata(hObject, han dles);% UIWAIT makes huatu wait for user respo nse (see UIRESUME)% uiwait(ha ndles.figure1);% --- Outputs from this fun cti on are returned to the comma nd line, fun cti on varargout = huatu_OutputFc n( hObject, eve ntdata, han dles) % varargout cell array for returning output args (see VARARGOUT);% hObject han dle to figure% eventdata reserved - to be defined in a future version of MATLAB % han dles structure with han dles and user data (see GUIDATA)% Get default comma nd line output from han dles structure varargout{1} = han dles.output;% --- Executes on butt on press in pushbutt on1.fun cti on pushbutt on 1_Callback(hObject, eve ntdata, han dles)% hObject handle to pushbutt on1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % han dles structure with han dles and user data (see GUIDATA) a=str2 nu m(get(ha ndles.edit1,'stri ng' ));b=str2 nu m(get(ha ndles.edit2, 'stri ng' ));c=str2 nu m(get(ha ndles.edit3, 'stri ng' ));x=-10:1:10;hold ony=a*x.A2+b*x+c;plot(x,y);% --- Executes on butt on press in pushbutt on2.fun cti on pushbutt on 2_Callback(hObject, eve ntdata, han dles)% hObject handle to pushbutt on2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% han dles structure with han dles and user data (see GUIDATA) close allfunction edit1_Callback(hObject, eventdata, handles)% hObject han dle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% han dles structure with han dles and user data (see GUIDATA)% Hin ts: get(hObject,'Stri ng') retur ns contents of edit1 as text% str2double(get(hObject,'Stri ng')) returns contents of edit1 as adouble% --- Executes during object creation, after setting all properties.function edit1_CreateFcn(hObject, eventdata, handles)% hObject han dle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% han dles empty - han dles not created un til after all CreateFc nscalled % Hint: edit con trols usually have a white backgro und on Win dows.% See ISPC and COMPUTER.if ispc && isequal(get(hObject, 'BackgroundColor'),get(0, 'defaultUicontrolBackgroundColor' ))set(hObject, 'BackgroundColor' , 'white');end function edit2_Callback(hObject, eventdata, handles)% hObject han dle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% han dles structure with han dles and user data (see GUIDATA)% Hin ts: get(hObject,'Stri ng') retur ns contents of edit2 as text% str2double(get(hObject,'Stri ng')) returns contents of edit2 as adouble% --- Executes on butt on press in pushbutt on3.fun cti on pushbutt on 3_Callback(hObject, eve ntdata, han dles)% --- Executes during object creation, after setting all properties. function edit2_CreateFcn(hObject, eventdata, handles)% hObject han dle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% han dles empty - han dles not created un til after all CreateFc ns % Hint: edit con trols usually have a white backgro und on Win dows. % See ISPC and COMPUTER.if ispc && isequal(get(hObject, 'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor' ))set(hObject, 'BackgroundColor' , 'white');end function edit3_Callback(hObject, eventdata, handles)% hObject han dle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % han dles structure with han dles and user data (see GUIDATA)% Hin ts: get(hObject,'Stri ng') retur ns contents of edit3 as text% str2double(get(hObject,'Stri ng')) returns contents of edit3 as a double % --- Executes during object creation, after setting all properties.function edit3_CreateFcn(hObject, eventdata, handles)% hObject han dle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% han dles empty - han dles not created un til after all CreateFc ns % Hint: edit con trols usually have a white backgro und on Win dows. % See ISPC and COMPUTER.if ispc && isequal(get(hObject, 'BackgroundColor'),get(0, 'defaultUicontrolBackgroundColor' ))set(hObject, 'BackgroundColor' , 'white');endcalled called% hObject handle to pushbutt on3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % han dles structure with han dles and user data (see GUIDATA) % axex2cla(ha ndles.axes2)。
实验二 用matlab绘制一元函数与二元函数的图象
实验二 用matlab 绘制一元函数与二元函数的图象1.平面曲线的表示形式对于平面曲线,常见的有三种表示形式,即以直角坐标方程],[),(b a x x f y ∈=,以参数方程],[),(),(b a t t y y t x x ∈==,和以极坐标],[),(b a r r ∈=ϕϕ表示等三种形式。
2.曲线绘图的MATLAB 命令MATLAB 中主要用plot,fplot 二种命令绘制不同的曲线。
可以用help plot, help fplot 查阅有关这些命令的详细信息例16.2.1 作出函数x y x y cos ,sin ==的图形,并观测它们的周期性。
先作函数x y sin =在]4,4[ππ-上的图形,用MA TLAB 作图的程序代码为:>>x=linspace(-4*pi,4*pi,300); %产生300维向量x >>y=sin(x);>>plot(x,y) %二维图形绘图命令结果如图1.1,上述语句中%后面如“%产生300维向量x ”是说明性语句,无需键入。
图1.1 的图形此图也可用fplot 命令,相应的MATLAB 程序代码为: >>clear; close; %clear 清理内存;close 关闭已有窗口。
>>fplot('sin(x)',[-4*pi,4*pi]) 结果如图1.2.图1.2xy sin=的图形如果在同一坐标系下作出两条曲线xy sin=和xy cos=在]2,2[ππ-上的图形,相应的MA TLAB程序代码为:>>x=-2*pi:2*pi/30:2*pi; %产生向量x>>y1=sin(x); y2=cos(x);>>plot(x,y1,x,y2,’:’)%’:’表示绘出的图形是点线结果如图1.3其中实线是xy sin=的图形,点线是xy cos=的图形。
matlab__函数图形绘制实验报告
实验报告课程名称:数学实验学院名称:数学与统计学院班级:姓名:学号:2012-2013—学年第_________ 学期数学与统计学院制Plot(x,y) grid;图像:(二)参数方程作图3x(t)=2cos t x(t)=2(t—sint)例2 :画出星形线{y(t)=2sin3t及旋轮线{y(t)=2(1-cost)的图形解:输入以下命令:%星形线作图t=li nspace(0,2*pi,5000);x=2*(cos(t)).A3;y=2*(si n(切八3;plot(x,y),grid;结果:%旋轮线作图t=li nspace(0,4*pi,5000); x=2*(t-si n( t));y=2*(1-cos(t)); plot(x,y),axis equal; axis(0,8*pi,0,5);grid;结果:(三)极坐标方程图形例3:画出四叶玫瑰线的图形。
知其极坐标方程:p=acos (2二)解:取a=5做图。
在命令窗口输入下命令theta=li nspace(0,2*pi); r=2*cos(2*theta); polar(theta,r)结果:z=(u.A2-v.A2)./4;surf(x,y,z);bix on;结果:(五)空间曲线在坐标平面上的投影曲面和投影柱面x 二lOcosty"Osint, f R在xOz面上的正投影曲线例5:画出螺旋线{z 二2t的图形。
x=10cost解:化为参数方程{z=2t ,运行下列程序t=lin space(-2*pi,2*pi);x=10*cos(t);z=2*t;h=plot(x,z);grid; xlabel('x');ylabel('z');set(h,'linewidth' ,2);结果:o实验分析:(一)在本次实验中我们初步了解了matlab。
(二)学会了一些简单绘图。
(三)在编制中我们要很明确点乘的重要性”。
Matlab绘制函数图像
Matlab绘制函数图像⼀、绘制图像的常⽤思路在通常情况下,⼿动绘制函数图像的基本思路是确定⾃变量的取值范围、选取合适的⾃变量点、通过函数表达式得出对应的应变量的点、将这些点连接起来,即可得到⼤致图像。
使⽤Matlab绘制函数图像也是使⽤这样的思路,当⾃变量的点取得够多时,连接这些点的图形就越接近函数图像。
⼆、Matlab中绘制图形的常⽤语⾔规范和绘图函数在Matlab中,选取⾃变量的操作变得及其简单,⼤致格式如下:X = 下限:取值间隔:上限。
详细情况可以看案例。
Matlab对于函数的表⽰也⼗分的接近⾃然语⾔,可以参考案例学习。
在绘制图像时,Matlab对于⼆维和三维图像有不同的函数。
绘制⼆维图像⼀般使⽤的是“plot”函数,在调⽤它时,只需要给定对应的⾃变量和应变量即可,例如plot(x,y)就是绘制x为⾃变量y为应变量的函数图像;对于三维图像⼀般使⽤“surf”函数,它的使⽤也和“plot”函数⼤致相同。
当然,这两种绘图函数还有⼀些更加深⼊的⽤法,这⾥不⼀⼀介绍。
三、注意事项由于在Matlab中对于矩阵的乘法有两种⽅法,⼀种是矩阵相乘,满⾜矩阵相乘的规律,⽤“*”符号表⽰;另⼀种是矩阵对应元素的相乘,即需要相乘的矩阵的维度相同,得到的矩阵维度不变,每个位置上的值都是对应元素相乘的结果,⽤“.*”符号表⽰。
同时,在绘制⼆维图像时,plot 函数默认没有使⽤⽹格线标注,所以还需要“grid”函数添加⽹格线条。
在含有除法的函数表达式中,还应该避免除数为零的异常情况,通常解决它的办法是在除数部分表达式后⾯加上⼀个eps(⽆穷⼩)。
四、案例图形绘制⼆维函数图像的绘制:在(0,4)的函数图像:x=0:pi/50:4*pi;y=exp(-t/3).*sin(3*t);plot(x,y);grid图4.1 ⼆维函数图像三维函数图像的绘制:的函数图像,x,y的取值范围都是[-8,8]:x=-8:0.5:8;#确定x的取值范围y=x';X=ones(size(y))*x;#由于是三维图像,所以此时的X不再是⼀维数组,⽽是⼀个⼆维数组,部分结果如下Y=y*ones(size(x));#和X同理R=sqrt(X.^2+Y.^2)+eps;#异常值处理Z=sin(R)./R;surf(X,Y,Z)图4.2 三维函数⾃变量X的表⽰图4.3 三维函数图像END编辑 | 王⽂星责编 | 黄章鱼能⼒越强,责任越⼤。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实验二 用matlab 绘制一元函数与二元函数的图象
1.平面曲线的表示形式
对于平面曲线,常见的有三种表示形式,即以直角坐标方程],[),(b a x x f y ∈=,以参数方程],[),(),(b a t t y y t x x ∈==,和以极坐标],[),(b a r r ∈=ϕϕ表示等三种形式。
2.曲线绘图的MATLAB 命令
MATLAB 中主要用plot,fplot 二种命令绘制不同的曲线。
可以用help plot, help fplot 查阅有关这些命令的详细信息
例16.2.1 作出函数x y x y cos ,sin ==的图形,并观测它们的周期性。
先作函数
x y sin =在]4,4[ππ-上的图形,用MA TLAB 作图的程序代码为:
>>x=linspace(-4*pi,4*pi,300); %产生300维向量x >>y=sin(x);
>>plot(x,y) %二维图形绘图命令
结果如图1.1,上述语句中%后面如“%产生300维向量x ”是说明性语句,无需键入。
图1.1 的图形
此图也可用fplot 命令,相应的MATLAB 程序代码为: >>clear; close; %clear 清理内存;close 关闭已有窗口。
>>fplot('sin(x)',[-4*pi,4*pi]) 结果如图1.2.
图1.2
x
y sin
=的图形
如果在同一坐标系下作出两条曲线
x
y sin
=和x
y cos
=在]
2,
2
[π
π
-上的图形,相应
的MA TLAB程序代码为:
>>x=-2*pi:2*pi/30:2*pi; %产生向量x
>>y1=sin(x); y2=cos(x);
>>plot(x,y1,x,y2,’:’)%’:’表示绘出的图形是点线
结果如图1.3其中实线是
x
y sin
=的图形,点线是x
y cos
=的图形。
图1.3的图形
例16.2.2将例1得到的图形用不同的线型及颜色加以绘制。
>>x=-2*pi:2*pi/30:2*pi; %产生向量x >>y1=sin(x); y2=cos(x);
>>plot(x,y1,x,y2,’gp’) %’ gp’表示绘出的图形是绿色五角星线 结果如图1.4
例16.2.3 作出以极坐标方程]2,0[,1),cos 1(πϕϕ∈=+=a a r 表示的心脏线 相应的MA TLAB 程序代码为:
>>clear; close;
>>t=0:2*pi/30:2*pi; >>r=1+cos(t);
>>x=r.*cos(t); y=r.*sin(t); %极坐标转化为直角坐标 >>plot(x,y)
结果如图1.5
3.曲面绘图的MATLAB 命令
MATLAB中主要用mesh,surf命令绘制二元函数图形。
可以用help mesh, help surf查阅有关这些命令的详细信息
例16.2.4画出函数
2
2y
x
z+
=
的图形,不妨将区域限制在
]3,3
[
]3,3
[
)
,
(-
⨯
-
∈
y
x。
用MA TLAB作图的程序代码为:
>>clear;
>>x=-3:0.1:3; %x的范围为[-3,3]
>>y=-3:0.1:3; %y的范围为[-3,3]
>>[X,Y]=meshgrid(x,y); %将向量x,y指定的区域转化为矩阵X,Y
>>Z=sqrt(X.^2+Y.^2); %产生函数值Z
>>mesh(X,Y,Z)
结果如图5.1。
图5.1是网格线图,如果要画完整的曲面图,只需将上述的MATLAB代码mesh(X,Y,Z)改为surf(X,Y,Z),结果如图5.2
图5.1 锥面
图5.2 锥面
要画等高线,需用contour,contour3命令.其中contour为二维等高线, contour3为三维等高线,如画图5.1的三维等高线, MA TLAB代码为:
>>clear;
>>x=-3:0.1:3;
>>y=-3:0.1:3;
>>[X,Y]=meshgrid(x,y);
>>Z=sqrt(X.^2+Y.^2);
>>contour3(X,Y,Z,10) %画10条等高线
>>xlabel('X-axis'),ylabel('Y-axis'),zlabel('Z-axis') %三个坐标轴的标记
>>title('Contour3 of Surface') %标题
>>grid on %画网格线
结果如图5.3.
图5.3 等高线
如画图5.1的二维等高线, MATLAB代码为:
>>clear; x=-3:0.1:3; y=-3:0.1:3;
>>[X,Y]=meshgrid(x,y); Z=sqrt(X.^2+Y.^2);
>> contour(X,Y,Z,10)
>>xlabel('X-axis'),ylabel('Y-axis')
>>title('Contour of Surface')
>>grid on
结果如图5.4.
图5.4 等高线
如果要画1=z 的等高线,则用命令
>>clear; x=-3:0.1:3; y=-3:0.1:3;
>>[X,Y]=meshgrid(x,y); Z=sqrt(X.^2+Y.^2); >> contour(X,Y,Z,[1 1])
结果如图5.5。
图5.5 等高线
习题16-2
1.画出下列各函数的图形.
(1)立方曲线3x y = (2).立方抛物线3
x y = (3).高斯曲线2
x e y -= (4).
sin x y x = (5).
1
(1)x
y x =+ (6).3sin 2cos y x x =+ 2.画出空间曲线
(1)
2sin,cos,4,
x t y t z t
===z =4t,[],
tππ
∈-
(2)
2
2
2
2
1
sin
10
y
x
y
x
z
+
+
+
=
在
30
,
30<
<
-y
x范围内的图形,并画出相应的等高线。
3. 画出函数
22
z x y
=+的图形,不妨将区域限制在(,)[1,1][1,1]
x y∈-⨯-。