matlab大图中画小图的简单方法
利用matlab画中国地图的几种方法
利用matl ab画中国地图的几种方法1、下载 ar c_map.zip,内有/s hape_files/chin a.shp文件,地址如下:http://ven us.un ive.i t/mat did.p hp?ut ente=smagr ini&b ase=e c++ur bana+e+reg ional e+-+a a+2008-2009%2F primo+modu lo%2F labor atori o+mat lab&c md=di r2、代码:f igure;inf o = s hapei nfo('china.shp')S = shap eread('chi na.sh p');mapsh ow(S);文件中可能还有更多地理,统计信息。
只是地图政治不正确,没有宝岛。
南沙在。
如何在MA TLAB中绘制中国地图了呢?下面介绍两种方法。
第一种方法是利用m_map工具箱。
该工具箱需要自己下载安装,下载的地址是htt p://w ww.eo s.ubc.ca/~rich/map.h tml。
m_map是一个功能很强大的工具箱,它提供多种投影方式,而且用法也很简单,容易上手,更重要的是,利用它可以很方便地绘制自己的数据。
有兴趣的读者可以参阅它的使用说明。
就在上述的网址上就有。
但是m_ma p提供的边界数据使用起来很不方便,数据是按照国家单个给出的,如果想画出世界各国的边界,就需要把每个国家的数据都下载下来,很麻烦。
网上有如何利用m_map来绘制行政边界的说明,例如下面的这个地址的作者就提供了一个具体的操作方法:将下载的.shp文件通过ma pinfo转换成.d xf文件,利用已编译好的for tran程序,读取经纬度信息,输出.da t文件。
【matlab】 输出图像大小调整
Matlab图形窗口大小的控制,plot窗口大小,figure大小,axis设置,实用Matlab中保存图像时,图形窗口大小的控制zz首先要了解的是Matlab是面向对象的。
最高等级的对象是screen,它定义了figure可以用的最大szie。
screen下面是figure。
figue就是你画图的时候跳出来的那个新的对话窗口。
如果figure变化,screen是不会跟着变化的。
但screen变化的话,figure就要跟着变化了。
figure下面是axes。
axes是那个窗口里面你要画的东西。
axes的大小和位置取决于figure,如果你放大缩小figure的大小的话,里面的图线也会跟着变化的。
set(gca,'position',[])因此,set (gca,'position',[0.1,0.1,0.9,0.9] );的作用是:设置坐标轴距离画板(图形窗口figure)边距。
[0.1,0.1,0.9,0.9] 分别为axes在figure中的左边界,下边界,宽度,高度,最小为0,最大为1(左边界,下边界为0,上边界,右边界为1)见下面的例子:-----------------------------------------------------------------------------figureset (gca,'position',[0.1,0.1,0.9,0.9] );x=1:0.1:10;y=sin(x);plot(x,y)-----------------------------------------------------------------------------结果见下图:set(gcf,'position',[])一般matlab绘出来图的框架(图形窗口)大都是正方形或者近似正方形的矩形,能不能画一些扁的矩形呢?使用图形的position属性可以做到。
MATLAB程序设计-关于常用绘图函数的运用方法及说明
MATLAB程序设计作业陈杰杰2013090302072014-11-3MATLAB具有强大的图形处理功能。
下面给出了3个m脚本文件,请在MA TLAB环境下运行,观察其输出。
要求根据每个m文件输出的图形(共18个),用中文翻译并解释产生每个输出图形的函数具体是什么?其功能是什么?文件1:clear all %清除工作区间所有的变量clf %清除图形窗口的内容mfilename('fullpath') %返回当前正在运行的函数所在文件的文件名(全部路径)echo on %显示M文件执行的每一条命令subplot(2,3,1) %使(2*3)幅子图中第一个子图成为当前图t = 0:0.1:10; %将以0为起点、以10为终点、以0.1为步长的一维矩阵赋值给t z = impulse(1, [1 1 1], t); %动力系统的脉冲响应数据,以时间t步长0.1为单位stairs(t(1:5:end),z(1:5:end)) %绘制阶梯状图,从第1行开始,间隔5行取1行,到最后1行为止hold on %保持当前坐标轴和图形,并接受下一次绘制plot(t,z,'r') %用红线绘制横轴为t、纵轴为z的二维函数图plot([0 t(end)], [0 0], 'k:') %用黑色虚线绘制函数图像,要求经过原点平行于横轴、取值范围为0到t的最后一个值title('Impulse Response - (STAIRS)') %将此图命名为Impulse Response - (STAIRS)(脉冲响应-(阶梯图))subplot(2,3,2) %使(2*3)幅子图中第二个子图成为当前图theta = 2*pi*(0:74)/75; %将以0为起点、以2*pi*74/75为终点、2*pi/75为步长的一维矩阵赋值给thetax = cos(theta); %计算cos(theta)的值,并赋给xy = sin(theta); %计算sin(theta)的值,并赋给yz = abs(fft(ones(10,1), 75))'; %ones(10,1)生成十行一列的全一矩阵;fft(ones(10,1), 75)进行快速傅里叶变换;z = abs(fft(ones(10,1), 75))'取幅值并转置stem3(x, y, z) %绘制三维杆状图title('Polar FFT - (STEM3)') %将此图命名为Polar FFT - (STEM3)(极坐标下快速傅里叶变换-(三维针状图))subplot(2,3,3) %使(2*3)幅子图中第二个子图成为当前图[X,Y,Z] = peaks(-2:0.25:2); %产生-2为起点、2为终点、0.25为步长的guassian分布矩阵,返回峰函数的三个坐标轴空间上的数值,X表示在x轴,Y表示在y轴,Z表示在z轴,这样每个点就对应一个(X,Y,)[U,V] = gradient(Z, 0.25); %返回二维数值梯度的U、V部分,这里的0.25指定了沿着梯度的方向取点的间隔为0.25contour(X,Y,Z,10); %绘制矩阵Z的等高线,绘制的等高线被限定在由X、Y指定的区域内,等高线条数为10hold on %保持当前坐标轴和图形,并接受下一次绘制quiver(X,Y,U,V); %使用箭头来直观的显示矢量场,该调用格式表示通过在(X, Y)指定的位置绘制小箭头来表示以该点为起点的向量(U,V)title('Surface Gradient - (CONTOUR & QUIVER)') %将此图命名为Surface Gradient - (CONTOUR & QUIVER)(表面梯度-(等高线图和矢量场图))theta = 0:0.1:4*pi; %将以0为起点、以4*pi为终点、0.1为步长的一维矩阵赋值给theta[x,y] = pol2cart(theta(1:5:end), theta(1:5:end)); %把极坐标(theta(1:5:end), theta(1:5:end))转换为对应的二维笛卡尔坐标(x,y),theta(1:5:end)表示从第一行到最后一行,每五行取一行subplot(2,3,4) %使(2*3)幅子图中第四个子图成为当前图polar(theta,theta) %绘制极坐标图像,第一个theta是用弧度制表示的角度,第二个theta是对应的半径axis([-13 13 -12.5 14.5]) %横坐标范围为-13到13,纵坐标范围为-12.5到14.5title('Spiral Plot - (POLAR)') %将此图命名为Spiral Plot - (POLAR)(螺旋图-极坐标图)subplot(2,3,5) %使(2*3)幅子图中第五个子图成为当前图compass(x,y) %绘制罗盘图axis([-13 13 -12.5 14.5]) %横坐标范围为-13到13,纵坐标范围为-12.5到14.5title('Direction Vectors - (COMPASS)') %将此图命名为Direction Vectors - (COMPASS)(方向矢量-罗盘图)subplot(2,3,6) %使(2*3)幅子图中第六个子图成为当前图feather(x(1:19),y(1:19)) %绘制羽状图,其中x、y表示一组向量,x是向量的横坐标(x分量),y是向量的纵坐标(y分量)。
matlab中image量化方法
matlab中image量化方法
在MATLAB中,图像量化是指将图像的像素值映射到一个较小的像素值集合上的过程。
这可以通过不同的方法来实现,以下是一些常见的图像量化方法:
1. 线性量化:
这是最简单的图像量化方法之一,它通过将图像的像素值映射到一个较小范围内的线性间隔来实现。
例如,如果原始像素值范围是0到255,线性量化可以将这个范围映射到0到15,从而将256个灰度级别减少到16个灰度级别。
2. 均匀量化:
均匀量化是一种常见的量化方法,它将原始像素值范围均匀地分割成较小的区间,然后将每个区间映射到一个代表该区间的像素值。
这种方法可以减少图像的细节,但可以有效地减少像素值的数量。
3. 非均匀量化:
与均匀量化相反,非均匀量化将像素值范围划分为不均匀的
区间,以便更好地保留图像的细节。
这种方法可以通过使用更多的
像素值来更好地表示图像的细节,但会增加图像的数据量。
4. 自适应量化:
自适应量化是一种根据图像的内容动态调整量化级别的方法。
它可以根据图像的局部特征来确定每个像素的量化级别,从而更好
地保留图像的细节。
在MATLAB中,可以使用`imquantize`函数来实现图像的量化。
该函数可以根据指定的量化级别将图像的像素值映射到新的像素值
集合上。
另外,也可以使用`histeq`函数来对图像进行直方图均衡化,从而增强图像的对比度并减少量化引入的失真。
总之,MATLAB提供了多种方法来实现图像量化,可以根据具体
的需求和图像特性选择合适的方法来进行图像量化处理。
MATLAB画图之自定义图片大小
MATLAB画图之⾃定义图⽚⼤⼩解决问题:使⽤MATLAB的plot函数画图时弹出图⽚的⼤⼩和位置修改。
MATLAB画图的⼏个概念:screen: 屏幕;figure: 弹出来的对话框;figure变化,screen不会变化;screen变化,figure不会变化;axes: figure 中的图像,axes的⼤⼩取决于figure,figure⼤⼩变化,axes⼤⼩也会变化;set(gcf, ……): 对figure⼤⼩和位置进⾏设置;set(gca, ……): 对axes⼤⼩和位置进⾏设置;⼀、⾃定义figure的⼤⼩和位置程序:clear;clc;close all;t = 0:0.001:10;y1 = sin(t);figure(1);plot(t,y1);程序运⾏结果:默认情况下,figure⼤致在screen的正中间,⼤致为⼀个⽅形(具体参数没有深究),如果我们想要⾃定义figure在screen的位置和⼤⼩,有两种⽅法可以实现。
第⼀种⽅法:在声明figure时定义⼤⼩,将上述程序修改为:clear;clc;close all;t = 0:0.001:10;y1 = sin(t);figure('Units','centimeter','Position',[5 5 7 3.5]);plot(t,y1);其中'Units','centimeter' ⽤来定义单位;'Position',[5 5 7 3.5]⽤来定义位置。
[5 5 7 3.5]表⽰为[x0 y0 width height]X0, y0: 表⽰figure左下⾓在screen中的位置,参考位置是screen左下⾓;Width, height: 表⽰figure的宽和⾼的⼤⼩。
程序运⾏结果:第⼆种⽅法:使⽤set(gcf, ……)命令,将上述程序修改为:clear;clc;close all;t = 0:0.001:10;y1 = sin(t);figure(1);set(gcf,'Units','centimeter','Position',[5 5 7 3.5]);plot(t,y1);单位和位置的定义⽅式同第⼀种⽅法,程序运⾏结果为:⼆、⾃定义figure中图像的⼤⼩和位置使⽤set(gca, ……)指令,该指令设置的是图的坐标线axes的位置及⼤⼩,并不是坐标线标注的位置。
给大家一个非常好用的matlab程序(一个figure中画多幅图,colormap如何设置)
给⼤家⼀个⾮常好⽤的matlab程序(⼀个figure中画多幅图,colormap如何设置)function freezeColors(varargin)% freezeColors Lock colors of plot, enabling multiple colormaps per figure. (v2.3)%% Problem: There is only one colormap per figure. This function provides% an easy solution when plots using different colomaps are desired% in the same figure.%% freezeColors freezes the colors of graphics objects in the current axis so% that subsequent changes to the colormap (or caxis) will not change the% colors of these objects. freezeColors works on any graphics object% with CData in indexed-color mode: surfaces, images, scattergroups,% bargroups, patches, etc. It works by converting CData to true-color rgb% based on the colormap active at the time freezeColors is called.%% The original indexed color data is saved, and can be restored using% unfreezeColors, making the plot once again subject to the colormap and% caxis.%%% Usage:% freezeColors applies to all objects in current axis (gca),% freezeColors(axh) same, but works on axis axh.%% Example:% subplot(2,1,1); imagesc(X); colormap hot; freezeColors% subplot(2,1,2); imagesc(Y); colormap hsv; freezeColors etc...%% Note: colorbars must also be frozen. Due to Matlab 'improvements' this can% no longer be done with freezeColors. Instead, please% use the function CBFREEZE by Carlos Adrian Vargas Aguilera% that can be downloaded from the MATLAB File Exchange% (/matlabcentral/fileexchange/24371)%% h=colorbar; cbfreeze(h), or simply cbfreeze(colorbar)%% For additional examples, see test/test_main.m%% Side effect on render mode: freezeColors does not work with the painters% renderer, because Matlab doesn't support rgb color data in% painters mode. If the current renderer is painters, freezeColors% changes it to zbuffer. This may have unexpected effects on other aspects% of your plots.%% See also unfreezeColors, freezeColors_pub.html, cbfreeze.%%% John Iversen (iversen@) 3/23/05%% Changes:% JRI (iversen@) 4/19/06 Correctly handles scaled integer cdata% JRI 9/1/06 should now handle all objects with cdata: images, surfaces,% scatterplots. (v 2.1)% JRI 11/11/06 Preserves NaN colors. Hidden option (v 2.2, not uploaded)% JRI 3/17/07 Preserve caxis after freezing--maintains colorbar scale (v 2.3)% JRI 4/12/07 Check for painters mode as Matlab doesn't support rgb in it.% JRI 4/9/08 Fix preserving caxis for objects within hggroups (e.g. contourf)% JRI 4/7/10 Change documentation for colorbars% Hidden option for NaN colors:% Missing data are often represented by NaN in the indexed color% data, which renders transparently. This transparency will be preserved% when freezing colors. If instead you wish such gaps to be filled with% a real color, add 'nancolor',[r g b] to the end of the arguments. E.g.% freezeColors('nancolor',[r g b]) or freezeColors(axh,'nancolor',[r g b]),% where [r g b] is a color vector. This works on images & pcolor, but not on% surfaces.% Thanks to Fabiano Busdraghi and Jody Klymak for the suggestions. Bugfixes% attributed in the code.% Free for all uses, but please retain the following:% Original Author:% John Iversen, 2005-10% john_iversen@appdatacode = 'JRI__freezeColorsData';[h, nancolor] = checkArgs(varargin);%gather all children with scaled or indexed CDatacdatah = getCDataHandles(h);%current colormapcmap = colormap;nColors = size(cmap,1);cax = caxis;% convert object color indexes into colormap to true-color data using% current colormapfor hh = cdatah',g = get(hh);%preserve parent axis climparentAx = getParentAxes(hh);originalClim = get(parentAx, 'clim');% Note: Special handling of patches: For some reason, setting% cdata on patches created by bar() yields an error,% so instead we'll set facevertexcdata instead for patches.if ~strcmp(g.Type,'patch'),cdata = g.CData;elsecdata = g.FaceVertexCData;end%get cdata mapping (most objects (except scattergroup) have it)if isfield(g,'CDataMapping'),scalemode = g.CDataMapping;elsescalemode = 'scaled';end%save original indexed data for use with unfreezeColorssiz = size(cdata);setappdata(hh, appdatacode, {cdata scalemode});%convert cdata to indexes into colormapif strcmp(scalemode,'scaled'),%4/19/06 JRI, Accommodate scaled display of integer cdata:% in MATLAB, uint * double = uint, so must coerce cdata to double % Thanks to O Yamashita for pointing this need outidx = ceil( (double(cdata) - cax(1)) / (cax(2)-cax(1)) * nColors);else %direct mappingidx = cdata;/8/09 in case direct data is non-int (e.g. image;freezeColors)% (Floor mimics how matlab converts data into colormap index.)% Thanks to D Armyr for the catchidx = floor(idx);end%clamp to [1, nColors]idx(idx<1) = 1;idx(idx>nColors) = nColors;%handle nans in idxnanmask = isnan(idx);idx(nanmask)=1; %temporarily replace w/ a valid colormap index%make true-color data--using current colormaprealcolor = zeros(siz);for i = 1:3,c = cmap(idx,i);c = reshape(c,siz);c(nanmask) = nancolor(i); %restore Nan (or nancolor if specified)realcolor(:,:,i) = c;end%apply new true-color color data%true-color is not supported in painters renderer, so switch out of that if strcmp(get(gcf,'renderer'), 'painters'),set(gcf,'renderer','zbuffer');end%replace original CData with true-color dataif ~strcmp(g.Type,'patch'),set(hh,'CData',realcolor);elseset(hh,'faceVertexCData',permute(realcolor,[1 3 2]))end%restore clim (so colorbar will show correct limits)if ~isempty(parentAx),set(parentAx,'clim',originalClim)endend %loop on indexed-color objects% ============================================================================ %% Local functions%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% getCDataHandles -- get handles of all descendents with indexed CData %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function hout = getCDataHandles(h)% getCDataHandles Find all objects with indexed CData%recursively descend object tree, finding objects with indexed CData% An exception: don't include children of objects that themselves have CData:% for example, scattergroups are non-standard hggroups, with CData. Changing% such a group's CData automatically changes the CData of its children,% (as well as the children's handles), so there's no need to act on them.error(nargchk(1,1,nargin,'struct'))hout = [];if isempty(h),return;endch = get(h,'children');for hh = ch'g = get(hh);if isfield(g,'CData'), %does object have CData?%is it indexed/scaled?if ~isempty(g.CData) && isnumeric(g.CData) && size(g.CData,3)==1,hout = [hout; hh]; %#ok<AGROW> %yes, add to listendelse %no CData, see if object has any interesting childrenhout = [hout; getCDataHandles(hh)]; %#ok<AGROW>endend%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% getParentAxes -- return handle of axes object to which a given object belongs %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function hAx = getParentAxes(h)% getParentAxes Return enclosing axes of a given object (could be self)error(nargchk(1,1,nargin,'struct'))%object itself may be an axisif strcmp(get(h,'type'),'axes'),hAx = h;returnendparent = get(h,'parent');if (strcmp(get(parent,'type'), 'axes')),hAx = parent;elsehAx = getParentAxes(parent);end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% checkArgs -- Validate input arguments %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function [h, nancolor] = checkArgs(args)% checkArgs Validate input arguments to freezeColorsnargs = length(args);error(nargchk(0,3,nargs,'struct'))%grab handle from first argument if we have an odd number of argumentsif mod(nargs,2),h = args{1};if ~ishandle(h),error('JRI:freezeColors:checkArgs:invalidHandle',...'The first argument must be a valid graphics handle (to an axis)')end% 4/2010 check if object to be frozen is a colorbarif strcmp(get(h,'Tag'),'Colorbar'),if ~exist('cbfreeze.m'),warning('JRI:freezeColors:checkArgs:cannotFreezeColorbar',...['You seem to be attempting to freeze a colorbar. This no longer'...'works. Please read the help for freezeColors for the solution.'])elsecbfreeze(h);returnendendargs{1} = [];nargs = nargs-1;elseh = gca;end%set nancolor if that option was specifiednancolor = [nan nan nan];if nargs == 2,if strcmpi(args{end-1},'nancolor'),nancolor = args{end};if ~all(size(nancolor)==[1 3]),error('JRI:freezeColors:checkArgs:badColorArgument',...'nancolor must be [r g b] vector');endnancolor(nancolor>1) = 1; nancolor(nancolor<0) = 0;elseerror('JRI:freezeColors:checkArgs:unrecognizedOption',...'Unrecognized option (%s). Only ''nancolor'' is valid.',args{end-1})endendfunction CBH = cbfreeze(varargin)�FREEZE Freezes the colormap of a colorbar.%% SYNTAX:% cbfreeze% cbfreeze('off')% cbfreeze(H,...)% CBH = cbfreeze(...);%% INPUT:% H - Handles of colorbars to be freezed, or from figures to search% for them or from peer axes (see COLORBAR).% DEFAULT: gcf (freezes all colorbars from the current figure)% 'off' - Unfreezes the colorbars, other options are:% 'on' Freezes% 'un' same as 'off'% 'del' Deletes the colormap(s).% DEFAULT: 'on' (of course)%% OUTPUT (all optional):% CBH - Color bar handle(s).%% DESCRIPTION:% MATLAB works with a unique COLORMAP by figure which is a big% limitation. Function FREEZECOLORS by John Iversen allows to use % different COLORMAPs in a single figure, but it fails freezing the% COLORBAR. This program handles this problem.%% NOTE:% * Optional inputs use its DEFAULT value when not given or [].% * Optional outputs may or not be called.% * If no colorbar is found, one is created.% * The new frozen colorbar is an axes object and does not behaves% as normally colorbars when resizing the peer axes. Although, some % time the normal behavior is not that good.% * Besides, it does not have the 'Location' property anymore.% * But, it does acts normally: no ZOOM, no PAN, no ROTATE3D and no% mouse selectable.% * No need to say that CAXIS and COLORMAP must be defined before using% this function. Besides, the colorbar location. Anyway, 'off' or% 'del' may help.% * The 'del' functionality may be used whether or not the colorbar(s)% is(are) froozen. The peer axes are resized back. Try:% >> colorbar, cbfreeze del%% EXAMPLE:% surf(peaks(30))% colormap jet% cbfreeze% colormap gray% title('What...?')%% SEE ALSO:% COLORMAP, COLORBAR, CAXIS% and% FREEZECOLORS by John Iversen% at /matlabcentral/fileexchange%%% ---% MFILE: cbfreeze.m% VERSION: 1.1 (Sep 02, 2009) (<a href="matlab:web('/matlabcentral/fileexchange/authors/11258')">download</a>) % MATLAB: 7.7.0.471 (R2008b)% AUTHOR: Carlos Adrian Vargas Aguilera (MEXICO)% CONTACT: nubeobscura@% REVISIONS:% 1.0 Released. (Jun 08, 2009)% 1.1 Fixed BUG with image handle on MATLAB R2009a. Thanks to Sergio% Muniz. (Sep 02, 2009)% DISCLAIMER:% cbfreeze.m is provided "as is" without warranty of any kind, under the% revised BSD license.% Copyright (c) 2009 Carlos Adrian Vargas Aguilera% INPUTS CHECK-IN% -------------------------------------------------------------------------% Parameters:cbappname = 'Frozen'; % Colorbar application data with fields:% 'Location' from colorbar% 'Position' from peer axes befor colorbar% 'pax' handle from peer axes.axappname = 'FrozenColorbar'; % Peer axes application data with frozen% colorbar handle.% Set defaults:S = 'on'; Sopt = {'on','un','off','del'};H = get(0,'CurrentFig');% Check inputs:if nargin==2 && (~isempty(varargin{1}) && all(ishandle(varargin{1})) && ...isempty(varargin{2}))% Check for CallBacks functionalities:% ------------------------------------varargin{1} = double(varargin{1});if strcmp(get(varargin{1},'BeingDelete'),'on')% Working as DeletFcn:if (ishandle(get(varargin{1},'Parent')) && ...~strcmpi(get(get(varargin{1},'Parent'),'BeingDeleted'),'on'))% The handle input is being deleted so do the colorbar:S = 'del';if ~isempty(getappdata(varargin{1},cbappname))% The frozen colorbar is being deleted:H = varargin{1};else% The peer axes is being deleted:H = ancestor(varargin{1},{'figure','uipanel'});endelse% The figure is getting close:returnendelseif (gca==varargin{1} && ...gcf==ancestor(varargin{1},{'figure','uipanel'}))% Working as ButtonDownFcn:cbfreezedata = getappdata(varargin{1},cbappname);if ~isempty(cbfreezedata)if ishandle(cbfreezedata.ax)% Turns the peer axes as current (ignores mouse click-over):set(gcf,'CurrentAxes',cbfreezedata.ax);returnendelse% Clears application data:rmappdata(varargin{1},cbappname)endH = varargin{1};endelse% Checks for normal calling:% --------------------------% Looks for H:if nargin && ~isempty(varargin{1}) && all(ishandle(varargin{1}))H = varargin{1};varargin(1) = [];end% Looks for S:if ~isempty(varargin) && (isempty(varargin{1}) || ischar(varargin{1}))S = varargin{1};endend% Checks S:if isempty(S)S = 'on';endS = lower(S);iS = strmatch(S,Sopt);if isempty(iS)error('CVARGAS:cbfreeze:IncorrectStringOption',...['Unrecognized ''' S ''' argument.' ])elseS = Sopt{iS};end% Looks for CBH:CBH = cbfreeze(H); �H = cbhandle(H);if ~strcmp(S,'del') && isempty(CBH)% Creates a colorbar and peer axes:pax = gca;CBH = colorbar('peer',pax);elsepax = [];end% -------------------------------------------------------------------------% MAIN% -------------------------------------------------------------------------% Note: only CBH and S are necesary, but I use pax to avoid the use of the % "hidden" 'Axes' COLORBAR's property. Why... �% Saves current position:fig = get( 0,'CurrentFigure');cax = get(fig,'CurrentAxes');% Works on every colorbar:for icb = 1:length(CBH)% Colorbar axes handle:h = double(CBH(icb));% This application data:cbfreezedata = getappdata(h,cbappname);% Gets peer axes:if ~isempty(cbfreezedata)pax = cbfreezedata.pax;if ~ishandle(pax) % just in casermappdata(h,cbappname)continueendelseif isempty(pax) % not generatedtrypax = double(get(h,'Axes')); % NEW feature in COLORBARs catchcontinueendend% Choose functionality:switch Scase 'del'% Deletes:if ~isempty(cbfreezedata)% Returns axes to previous size:oldunits = get(pax,'Units');set(pax,'Units','Normalized');set(pax,'Position',cbfreezedata.Position)set(pax,'Units',oldunits)set(pax,'DeleteFcn','')if isappdata(pax,axappname)rmappdata(pax,axappname)endendif strcmp(get(h,'BeingDelete'),'off')delete(h)endcase {'un','off'}% Unfrozes:if ~isempty(cbfreezedata)delete(h);set(pax,'DeleteFcn','')if isappdata(pax,axappname)rmappdata(pax,axappname)endoldunits = get(pax,'Units');set(pax,'Units','Normalized')set(pax,'Position',cbfreezedata.Position)set(pax,'Units',oldunits)CBH(icb) = colorbar(...'peer' ,pax,...'Location',cbfreezedata.Location);endotherwise % 'on'% Freezes:% Gets colorbar axes properties:cb_prop = get(h);% Gets colorbar image handle. Fixed BUG, Sep 2009hi = findobj(h,'Type','image');% Gets image data and transform it in a RGB:CData = get(hi,'CData');if size(CData,3)~=1% It's already frozen:continueend% Gets image tag:Tag = get(hi,'Tag');% Deletes previous colorbar preserving peer axes position: oldunits = get(pax,'Units');set(pax,'Units','Normalized')Position = get(pax,'Position');delete(h)cbfreezedata.Position = get(pax,'Position');set(pax,'Position',Position)set(pax,'Units',oldunits)% Generates new colorbar axes:% NOTE: this is needed because each time COLORMAP or CAXIS is used, % MATLAB generates a new COLORBAR! This eliminates that behaviour % and is the central point on this function.h = axes(...'Parent' ,cb_prop.Parent,...'Units' ,'Normalized',...'Position',cb_prop.Position...);% Save location for future call:cbfreezedata.Location = cb_prop.Location;% Move ticks because IMAGE draws centered pixels:XLim = cb_prop.XLim;YLim = cb_prop.YLim;if isempty(cb_prop.XTick)% Vertical:X = XLim(1) + diff(XLim)/2;Y = YLim + diff(YLim)/(2*length(CData))*[+1 -1];else % isempty(YTick)% Horizontal:Y = YLim(1) + diff(YLim)/2;X = XLim + diff(XLim)/(2*length(CData))*[+1 -1];end% Draws a new RGB image:image(X,Y,ind2rgb(CData,colormap),...'Parent' ,h,...'HitTest' ,'off',...'Interruptible' ,'off',...'SelectionHighlight','off',...'Tag' ,Tag...)% Removes all '...Mode' properties:cb_fields = fieldnames(cb_prop);indmode = strfind(cb_fields,'Mode');for k=1:length(indmode)if ~isempty(indmode{k})cb_prop = rmfield(cb_prop,cb_fields{k});endend% Removes special COLORBARs properties:cb_prop = rmfield(cb_prop,{...'CurrentPoint','TightInset','BeingDeleted','Type',... % read-only'Title','XLabel','YLabel','ZLabel','Parent','Children',... % handles'UIContextMenu','Location',... % colorbars'ButtonDownFcn','DeleteFcn',... % callbacks'CameraPosition','CameraTarget','CameraUpVector','CameraViewAngle',... 'PlotBoxAspectRatio','DataAspectRatio','Position',...'XLim','YLim','ZLim'});% And now, set new axes properties almost equal to the unfrozen% colorbar:set(h,cb_prop)% CallBack features:set(h,...'ActivePositionProperty','position',...'ButtonDownFcn' ,@cbfreeze,... % mhh...'DeleteFcn' ,@cbfreeze) % againset(pax,'DeleteFcn' ,@cbfreeze) % and again!% Do not zoom or pan or rotate:setAllowAxesZoom (zoom ,h,false)setAllowAxesPan (pan ,h,false)setAllowAxesRotate(rotate3d,h,false)% Updates data:CBH(icb) = h;% Saves data for future undo:cbfreezedata.pax = pax;setappdata( h,cbappname,cbfreezedata);setappdata(pax,axappname,h);end % switch functionalityend % MAIN loop% OUTPUTS CHECK-OUT% -------------------------------------------------------------------------% Output?:if ~nargoutclear CBHelseCBH(~ishandle(CBH)) = [];end% Returns current axes:if ishandle(cax)set(fig,'CurrentAxes',cax)end% [EOF] cbfreeze.m⼤家只要把这个.m⽂件matlab⼯作⽂件夹中即可,此函数⽤来在⼀个figure上作出多个⽴体图形。
MATLAB画图——基础篇
MATLAB画图——基础篇MATLAB画图——基础篇在MATLAB使⽤的过程中,学会画图是⼀项必要的技能。
在这⾥,我总结了部分简单的画图函数,同时附上代码(本⽂中的程序为了⽅便给出的数据都很简单,⼤家可以⾃⼰去尝试其他数据)。
这对刚刚开始接触MATLAB的⼩⽩来说,我认为还是很有帮助的。
⽂章⽬录⼀、plot()函数1.⼆维图形(1)绘图选项线型颜⾊标记符号-实线b蓝⾊.点s⽅块:虚线g绿⾊o圆圈d菱形.-点划线r红⾊x叉v朝下三⾓符号-双划线c青⾊+加号^朝上三⾓符号m品红*星号<朝左三⾓符号y黄⾊>朝右三⾓符号p五⾓星k⿊⾊h六⾓星w⽩⾊(2)图形的辅助标注和窗⼝的分割title(图形说明)xlabel(x轴说明)ylabel(y轴说明)text(x,y图形说明)——在x,y轴处添加⽂字说明legend(图例⼀,图例⼆,…)subplot(m,n,p)——将绘图区域分割成m*n个⼦区域,并按照⾏从左⾄ 右,从上⾄下依次编号。
p表⽰第p个绘图⼦区域。
注意:如果是要两个图画到同⼀个坐标⾥⾯,则在两个plot函数之间添加⼀⾏hold on(3)格式plot(x)——缺省⾃变量绘图格式plot(x,y)——基本格式。
以y(x)的函数关系作图。
如果y是n*m的矩 阵,则x为⾃变量,作出m条曲线。
plot(x1,y1,x2,y2,…,xn,yn)——多条曲线绘图格式plot(x1,y1,选项1,x2,y2,选项2,…,xn,yn,选项n)——含选项的绘图格式x1=[1 2 3 4 5 6 7 8 9];x2=[2 4 6 8 10 12 14 16 18];y1=[1 4 9 16 25 36 49 64 81];y2=[18 16 14 12 10 8 6 4 2];subplot(4,1,1);plot(x1);title('例⼀');xlabel('⾃变量');ylabel('因变量');subplot(4,1,2);plot(x1,y1);title('例⼆');xlabel('⾃变量');ylabel('因变量');subplot(4,1,3);plot(x1,y1,x2,y2);title('例三');xlabel('⾃变量');ylabel('因变量');subplot(4,1,4);plot(x1,y1,'m+',x2,y2,'c*');title('例四');xlabel('⾃变量');ylabel('因变量');2.三维图形(1)格式plot3(x1,y1,z1,‘选项⼀’,x2,y2,z1,‘选项⼆’,…)x,y,z是长度相同的向量:⼀条曲线x,y,z是维度相同的矩阵:多条曲线(2)⽹格矩阵⽣成函数:meshgrid[X,Y]=meshgrid(x,y)x,y是给定的向量,X,Y是⽹格划分后得到的⽹格矩阵注意,这个函数⽤来⽣成⽹格矩阵,不是直接⽤来画图的,配合mesh使⽤。
详尽全面的matlab绘图教程
详尽全⾯的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 绘制曲线这是以参数形式给出的曲线⽅程,只要给定参数向量,再分别求出x,y向量即可输出曲线:>> t=-pi:pi/100:pi;>> x=t.*cos(3*t);>> y=t.*sin(t).*sin(t);>> plot(x,y)程序执⾏后,打开⼀个图形窗⼝,在其中绘制出如下曲线以上提到plot函数的⾃变量x,y为长度相同的向量,这是最常见、最基本的⽤法。
如何在Matlab中进行二维和三维绘图
如何在Matlab中进行二维和三维绘图在科学研究和工程领域,数据可视化是一项十分重要的任务,而Matlab作为一种功能强大的数值计算和数据分析软件,自然也提供了丰富的绘图功能。
本文将介绍如何在Matlab中进行二维和三维绘图,并探讨一些常见的绘图技巧和应用。
一、二维绘图Matlab中的二维绘图是最常见和基础的绘图任务之一。
在绘制二维图形时,我们通常会用到plot函数。
这个函数可以接受单个向量作为输入,将这个向量的值作为y轴上的数据点,自动生成与该向量长度相同的x轴坐标。
例如,我们可以用以下代码绘制一个简单的二维折线图:```x = 0:0.1:2*pi;y = sin(x);plot(x, y);```上述代码中,x参量取从0到2π的均匀间隔的值,而y则是根据x计算得到的sin函数值。
plot函数会自动根据输入绘制折线图,并添加相应的轴标签和图例。
在实际应用中,我们经常需要绘制多条曲线在同一个坐标系中进行对比分析。
可以通过在plot函数中传入多个x和y向量实现这一功能。
例如,我们可以通过以下代码绘制一个简单的双曲线图:```x = 0:0.1:2*pi;y1 = sin(x);y2 = cos(x);plot(x, y1, x, y2);```这样,就会在同一个坐标系中同时绘制sin曲线和cos曲线。
除了折线图,Matlab还支持其他常见的二维绘图类型,如散点图、柱状图和面积图等。
这些绘图类型可以通过不同的函数实现,例如scatter、bar和area等。
这里不再一一赘述,读者可以通过Matlab的帮助文档或官方网站了解更多的用法和示例。
二、三维绘图除了二维绘图,Matlab也提供了丰富的三维绘图功能,用于可视化更为复杂的数据和模型。
在绘制三维图形时,我们通常会用到surf函数。
这个函数可以接受两个二维矩阵作为输入,将这两个矩阵的值分别作为x、y轴上的坐标,而将第三个二维矩阵的值作为z轴上的数据点。
Matlab中使用Plot函数动态画图方法总结
Matlab中使用Plot函数动态画图方法总结Matlab除了强大的矩阵运算,仿真分析外,绘图功能也是相当的强大,静态画图没什么问题,由于Matlab本身的多线程编程缺陷,想要动态的画图,并且能够很好的在GUI中得到控制,还不是一件很容易的事情,下面总结几种方法。
一. AXIS 移动坐标系这种方法是最简单的一种方法,适合于数据已经全部生成的场合,先画图,然后移动坐标轴。
实例代码如下:%%%先画好,然后更改坐标系%在命令行中使用Ctrl+C 结束t=0:0.1:100*pi;m=sin(t);plot(t,m);x=-2*pi;axis([x,x+4*pi,-2,2]);grid onwhile 1if x>max(t)break;endx=x+0.1;axis([x,x+4*pi,-2,2]); %移动坐标系pause(0.1);end二. Hold On 模式此种方法比较原始,适合于即时数据,原理是先画上一帧,接着保留原始图像,追加下一幀图像,此种方式比较繁琐,涉及画图细节,并且没有完整并连续的Line对象数据。
例如:%%% Hold On 法% 此种方法只能点,或者分段划线hold offt=0;m=0;t1=[0 0.1]; %要构成序列m1=[sin(t1);cos(t1)];p = plot(t,m,'*',t1,m1(1,:),'-r',t1,m1(2,:),'-b','MarkerSize',5);x=-1.5*pi;grid on;for i=1:100hold ont=0.1*i; %下一个点m=t-floor(t);t1=t1+0.1; %下一段线(组)m1=[sin(t1);cos(t1)];p = plot(t,m,'*',t1,m1(1,:),'-r',t1,m1(2,:),'-b','MarkerSize',5);x=x+0.1;axis([x x+2*pi -1.5 1.5]);pause(0.01);end三. Plot 背景擦除模式这种模式比较适合画动画,效率比较高,刷新闪烁小,适合即时数据,最终的Line结构数据完整。
MATLAB作图(超详细)
2020/5/31
数学建模
3. 对数坐标图
在很多工程问题中,通过对数据进行对数转换可以 更清晰地看出数据的某些特征,在对数坐标系中描绘数 据点的曲线,可以直接地表现对数转换.对数转换有双对 数坐标转换和单轴对数坐标转换两种.用loglog函数 可以实现双对数坐标转换,用semilogx和semilogy 函数可以实现单轴对数坐标转换. loglog(Y) 表示 x、y坐标都是对数坐标系
单击鼠标左键,则在当前图形窗口中,以鼠标点中的点为 中心的图形放大2倍;单击鼠标右键,则缩小2倍.
zoom off 关闭缩放模式
grid on
%标注格栅
MATLAB liti37
例 创建一个简单的半对数坐标图. 解 输入命令:
x=0:.1:10;
semilogy(x,10.^x)
MATLAB liti38
例 绘制y=x3的函数图、对数坐标图、半对数坐标图.
2020/5/31
MATLAB liti22 数学建模
返回
三维图形 1. 空间曲线 2. 空间曲面
semilogx(Y) 表示 x坐标轴是对数坐标系
semilogy(…) 表示y坐标轴是对数坐标系
plotyy 有两个y坐标轴,一个在左边,一个在右边
2020/5/31
数学建模
例 用方形标记创建一个简单的loglog.
解 输入命令:
x=logspace(-1,2);
loglog(x,exp(x),’-s’)
数学建模
返回
2. 定制坐标 Axis([xmin xmax ymin ymax zmin zmax])定制图形坐标
x、y、z的最大、最小值
Axis
将坐标轴返回到自动缺省值
MATLAB 中图形的绘制
plot3(x,y,z) %画出连接点列的蓝色实心线(默认)
plot3(x,y,z, ‘r-’)
%画出连接点列的红色实线
例 画出螺旋线:x=sin(t),y=cos(t),z=t,
t [0,10 ]
上一段曲线。
► t=0:pi/50:10*pi; %生成参数t数组
%保持第二个子窗中绘图
► plot(x,y,'bo',x,z,'k+') %用'o'和'+'标记曲线上分点
► hold off 结果如图
%取消图形保持
二维绘图命令plot的图形效果
例2: 画出 [0,2 ]上正弦、余弦曲线并对线型加粗、 点型加大,重新定置坐标系。
► x=0:0.1*pi:2*pi;
1) 对投影域进行划分:
► x=a:p1:b ► y=c:p2:d
%按步长p1对 [a,b]区间进行划分并生成向量 %按步长p2 对[c,d]区间进行划分并生成向量
2) 按上述划分生成投影域上全部网格节点的坐标矩阵: ► [X,Y]=meshgrid(x,y)
3) 根据函数表达式生成全部网格节点处对应的函数值 (坐标z)矩阵Z:
实现在同一个窗口中同时显示多个图像的命令 subplot。 使用格式为:
subplot(m,n,i) 其含义为 :把图形窗口分割为 m 行 n 列子窗口, 然后选定第 i 个窗口为当前窗口。 subplot 命令不仅用于二维图形,对三维图形一样 适用。其本质是将 figure 窗口分为几个区域,再 在每个区域内分别绘图。
画出顺次连接这些点的曲线。用法: plot(x,y) %画出连接点列的蓝色实心线(默认)
matlab中画图的时各种设置
2)坐标轴控制:MATLAB的缺省方式是在绘图时,将所在的坐标系也画出
来,为隐去坐标系,可用axis off;axis on则显示坐标轴
(缺省值).
3)通常MATLAB的坐标系是长方形,长宽比例大约是4:3,为了得到一个
正方形的坐标系可用:axis square
4)坐标系横纵轴的比例是自动设置的,比例可能不一样,要得到相同比例
'MarkerFaceColor','y',... %设置标记点填充颜色为黄色
'MarkerSize',10) %设置标记点的尺寸为10
绘出图形如下:
00.511.522.533.5
-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0.3
0.4
0.5
34
4.双Y轴绘图:plotyy()函数.
其调用格式为:
plotyy(x1,y1,x2,y2)------绘制由x1,y1和x2,y2确定的两组曲线,其中x1,
y1的坐标轴在图形窗口的左侧,x2,y2的坐标轴在
图形窗口的右侧.
Plotyy(x1,y1,x2,y2, 'function1','function2')------功能同上,function是指那些
5.打开图形窗口的方法有三种:
1)调用绘图函数时自动打开;
2)用File---New---Figure新建;
3)figure命令打开,close命令关闭.
在运行绘图程序前若已打开图形窗口,则绘图函数不再打开,而直接利用已
打开的图形窗口;若运行程序前已存在多个图形窗口,并且没有指定哪个窗
Daubechies小波分析的部分Matlab实现
实验5:Daubechies 小波分析的部分Matlab 实现注:本实验在原本操作过程中在重构部分存在着问题,经思考检验调试之后,排除了隐藏的问题,实验得到成功进行同时达到了预期的结果。
第一部分Daubechies 尺度序列,二进点上的尺度函数图像及小波函数图像第一部分实验的目标有两个:1,给出各阶Daubechies 小波的低通滤波器,即尺度函数两尺度序列。
具体的程序算法由《小波导论》中7.3节中的引理7.16,定理7.17给现,其中定理7.17的证明过程实际上给出了一个显式的计算过程,这是程序实现的重要依据。
2,根据两尺度关系,由第一步骤中得到的两尺度序列逐步计算得到Daubechies 尺度函数及小波函数的二进点值,进而给出近似图像。
第一步骤相关函数程序解释:1、sinj(n),实现将引理7.16中的2(sin)2n ω转化为形如cos k ω,n 即表示(sin)2ω的次数。
2、sinj(n)中涉及函数cosn(n), cosn(n)实现将cos kω转化为cos()kkw ∑,n 表示三角函数次数。
程序编写的思路就是先将2(sin)2n ω转化为cos k ω再转化为cos()kkw ∑,然后根据定理7.17的证明过程计算出两尺度序列,其中用于提取多项式系数的函数利用了三角函数的正交性,使用了积分函数int()。
daucoef(n)给出n 阶daubechies 尺度函数的两尺度序列。
具体的原程序代码见附件,部分程序运行结果如下:各阶(2--12)Daubechies 尺度函数的两尺度序列如下:)第二步骤:给出n 阶daubechies 尺度函数图像,由函数dauinterp(n,max1)实现,即给出n 阶daubechies 尺度函数的max1次迭代的函数值,所得到的点都是函数的二进点。
Dauinterpwa(n,max1)则给出小波函数的二进点上的函数值。
这一部分内容中,比较难以实现的地方在于下标的变换,这常常引起一些数据的混乱或错误,相关的数学推导总是难以避免的,由于过程过于繁复,这里不给出详细的推导过程及解释。
怎样用Matlab绘图
help plothelp axisa1=plot();hlod ona2=plot();legend([a1 a2],'图1 名',‘图2 名')hold offx1=-pi:pi/12:pi;x2=-pi:pi/12:pi;y1=sin(x1);y2=cos(x2);plot(x1,y1,x2,y2);axis([-2*pi 2*pi -2 2]);xlabel('x');ylabel('y');title('sin(x) & cos(x)');MATLAB受到控制界广泛接受的一个重要原因是因为它提供了方便的绘图功能.这里主要介绍2 维图形对象的生成函数及图形控制函数的使用方法,还将简单地介绍一些图形的修饰与标注函数及操作和控制MATLAB 各种图形对象的方法.第一节图形窗口与坐标系一.图形窗口1.MATLAB 在图形窗口中绘制或输出图形,因此图形窗口就像一张绘图纸.2.在MATLAB 下,每一个图形窗口有唯一的一个序号h,称为该图形窗口的句柄.MATLAB 通过管理图形窗口的句柄来管理图形窗口;3.当前窗口句柄可以由MATLAB 函数gcf 获得;4. 在任何时刻, 只有唯一的一个窗口是当前的图形窗口( 活跃窗口);figure(h)----将句柄为h 的窗口设置为当前窗口;5.打开图形窗口的方法有三种:1)调用绘图函数时自动打开;2)用File---New---Figure 新建;3)figure 命令打开,close 命令关闭.在运行绘图程序前若已打开图形窗口,则绘图函数不再打开,而直接利用已打开的图形窗口;若运行程序前已存在多个图形窗口,并且没有指定哪个窗口为当前窗口时,则以最后使用过的窗口为当前窗口输出图形.6.窗口中的图形打印:用图形窗口的File 菜单中的Print 项.7.可以在图形窗口中设置图形对象的参数.具体方法是在图形窗口的Edit 菜单中选择Properties 项,打开图形对象的参数设置窗口,可以设置对象的属性.二.坐标系1.一个图形必须有其定位系统,即坐标系;2.在一个图形窗口中可以有多个坐标系,但只有一个当前的坐标系;3.每个坐标系都有唯一的标识符,即句柄值;4.当前坐标系句柄可以由MATLAB 函数gca 获得;5.使某个句柄标识的坐标系成为当前坐标系,可用如下函数:axes(h) h 为指定坐标系句柄值.6.一些有关坐标轴的函数:1)定义坐标范围:一般MATLAB 自动定义坐标范围,如用户认为设定的不合适,可用:axis([Xmin, Xmax, Ymin, Ymax])重新设定;2)坐标轴控制:MATLAB 的缺省方式是在绘图时,将所在的坐标系也画出来,为隐去坐标系,可用axis off;axis on 则显示坐标轴(缺省值).3)通常MATLAB 的坐标系是长方形,长宽比例大约是4:3,为了得到一个正方形的坐标系可用:axis square4)坐标系横纵轴的比例是自动设置的,比例可能不一样,要得到相同比例的坐标系,可用:axis equal第二节二维图形的绘制一. plot 函数plot 函数是最基本的绘图函数,其基本的调用格式为:1.plot(y)------绘制向量y 对应于其元素序数的二维曲线图,如果y 为复数向量,则绘制虚部对于实部的二维曲线图.例:绘制单矢量曲线图.y=[0 0.6 2.3 5 8.3 11.7 15 17.7 19.4 20];plot(y)由于y 矢量有10 个元素,x 坐标自动定义为[1 2 3 4 5 6 7 8 9 10].2.plot(x,y)------绘制由x,y 所确定的曲线.1)x,y 是两组向量,且它们的长度相等,则plot(x,y)可以直观地绘出以x 为横坐标,y 为纵坐标的图形.如:画正弦曲线:t=0:0.1:2*pi;y=sin(t);plot(t,y)2)当plot(x,y)中,x 是向量,y 是矩阵时,则绘制y 矩阵中各行或列对应于向量x的曲线.如果y 阵中行的长度与x 向量的长度相同,则以y 的行数据作为一组绘图数据;如果y 阵中列的长度与x 向量的长度相同,则以y 的列数据作为一组绘图数据;如果y 阵中行,列均与x 向量的长度相同,则以y 的每列数据作为一组绘图数据.例:下面的程序可同时绘出三条曲线.MATLAB 在绘制多条曲线时,会按照一定的规律自动变化每条曲线的的颜色.x=0:pi/50:2*pi;y(1,:)=sin(x);y(2,:)=0.6*sin(x);y(2,:)=0.3*sin(x);plot(x,y)或者还可以这样用:x=0:pi/50:2*pi;y=[ sin(x); 0.6*sin(x); 0.3*sin(x)];plot(x,y)3) 如果x,y 是同样大小的矩阵,则plot(x,y)绘出y 中各列相应于x 中各列的图形.例: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(2,:));y(3,:)=0.3*sin(x(3,:));plot(x,y)x=x';y=y';figureplot(x,y)在这个例子中,x------3x101,y------3x101,所以第一个plot 按列画出101 条曲线,每条3 个点;而x'------101x3,y'------101x3,所以第二个plot 按列画出3 条曲线,每条101 个点.3.多组变量绘图:plot(x1, y1, 选项1, x2, y2, 选项2, ……)上面的plot 格式中,选项是指为了区分多条画出曲线的颜色,线型及标记点而设定的曲线的属性.MATLAB 在多组变量绘图时,可将曲线以不同的颜色,不同的线型及标记点表示出来.这些选项如下表所示:各种颜色属性选项'r' 红色'm' 粉红'g' 绿色'c' 青色'b' 兰色'w' 白色'y' 黄色'k' 黑色各种线型属性选项'-' 实线'--' 虚线':' 点线'-.' 点划线各种标记点属性选项'.' 用点号绘制各数据点'^' 用上三角绘制各数据点'+' 用'+'号绘制各数据点'v' 用下三角绘制各数据点'*' 用'*'号绘制各数据点'>' 用右三角绘制各数据点' .' 用'.'号绘制各数据点'<' 用左三角绘制各数据点's'或squar 用正方形绘制各数据点'p' 用五角星绘制各数据点'd'或diamond 用菱形绘制各数据点'h' 用六角星绘制各数据点这些选项可以连在一起用,如:'-.g'表示绘制绿色的点划线,'g+'表示用绿色的'+'号绘制曲线.注意:1)表示属性的符号必须放在同一个字符串中;2)可同时指定2~3 个属性;3)与先后顺序无关;4)指定的属性中,同一种属性不能有两个以上.例:t=0:0.1:2*pi;y1=sin(t);y2=cos(t);y3=sin(t).*cos(t);plot(t,y1, '-r',t,y2, ':g',t,y3, '*b')该程序还可以按下面的方式写:t=0:0.1:2*pi;y1=sin(t);y2=cos(t);y3=sin(t).*cos(t);plot(t,y1, '-r')hold onplot(t,y2, ':g')plot(t,y3, '*b')hold off注:在MATLAB 中,如画图前已有打开的图形窗口,则再画图系统将自动擦掉坐标系中已有的图形对象,但设置了hold on 后,可以保持坐标系中已绘出的图形.还可以进一步设置包括线的宽度(LineWidth), 标记点的边缘颜色(MarkerEdgeColor),填充颜色(MarkerFaceColor)及标记点的大小(MarkerSize)等其它绘图属性.例:设置绘图线的线型,颜色,宽度,标记点的颜色及大小.t=0:pi/20:pi;y=sin(4*t).*sin(t)/2;plot(t,y,'-bs','LineWidth',2,... %设置线的宽度为2'MarkerEdgeColor','k',... %设置标记点边缘颜色为黑色'MarkerFaceColor','y',... %设置标记点填充颜色为黄色'MarkerSize',10) %设置标记点的尺寸为104.双Y 轴绘图:plotyy()函数.其调用格式为: plotyy(x1,y1,x2,y2)------绘制由x1,y1 和x2,y2 确定的两组曲线,其中x1,y1 的坐标轴在图形窗口的左侧,x2,y2 的坐标轴在图形窗口的右侧.Plotyy(x1,y1,x2,y2, 'function1','function2')------功能同上,function 是指那些绘图函数如:plot,semilogx,loglog 等.例如:在一个图形窗口中绘制双Y 轴曲线.x=0:0.3:12;y=exp(-0.3*x).*sin(x)+0.5;plotyy(x,y,x,y,'plot','stem')stem:绘制stem 形式的曲线(上端带圈的竖线).绘图结果:两条图线自动用不同的颜色区分,两个坐标的颜色与图线的颜色相对应,左边的Y 轴坐标对应的是plot 形式的曲线,右边的Y 坐标对应的是stem 形式的曲线.二.对数坐标图绘制函数:在对数坐标图的绘制中,有三种绘图函数:semilogx,semilogy 和loglog 函数.1)semilogx( )------绘制以X 轴为对数坐标轴的对数坐标图. 其调用格式为:semilogx(x,y,'属性选项')其中属性选项同plot 函数.该函数只对横坐标进行对数变换,纵坐标仍为线性坐标.2)semilogy( )------绘制以Y 轴为对数坐标轴的对数坐标图. 其调用格式为:semilogy(x,y,'属性选项')该函数只对纵坐标进行对数变换,横坐标仍为线性坐标.3)loglog( )------ 绘制X,Y 轴均为对数坐标轴的图形.其调用格式为:loglog(x,y,'属性选项')该函数分别对横,纵坐标都进行对数变换.例:x=0:0.1:6*pi;y=cos(x/3)+1/9;subplot(221), semilogx(x,y);subplot(222), semilogy(x,y);subplot(223), loglog(x,y);4)MATLAB 还提供了一个实用的函数:logspace( )函数,可按对数等间距地分布来产生一个向量,其调用格式为:x=logspace(x1,x2,n)这里,x1 表示向量的起点;x2 表示向量的终点;n 表示需要产生向量点的个数(一般可以不给出,采用默认值50).在控制系统分析中一般采用这种方法来构成频率向量w.关于它的应用后面还要讲到.三.极坐标图的绘制函数:绘极坐标图可用polar( )函数.其调用格式如下:polar(theta, rho,'属性选项')------theta:角度向量,rho:幅值向量,属性内容与plot 函数基本一致.例如:极坐标模型为:3145/)/)cos((+ =θρ, ],[πθ80∈则绘出极坐标图的程序为:theta=0:0.1:8*pi;p=cos((5*theta)/4)+1/3;polar(theta,p)四.绘制多个子图:subplot( )函数MATLAB 允许在一个图形窗口上绘制多个子图(如对于多变量系统的输出),允许将窗口分成nxm 个部分.分割图形窗口用subplot 函数来实现,其调用格式为:subplot(n,m,k)或subplot(nmk)------n,m 分别表示将窗口分割的行数和列数,k 表示要画图部分的代号,表示第几个图形,nmk 三个数可以连写,中间不用符号分开.例如:将窗口划分成2x2=4 个部分,可以这样写:subplot(2,2,1),plot(……)subplot(2,2,2),……subplot(2,2,3),……subplot(2,2,4),……注:subplot 函数没有画图功能,只是将窗口分割.第三节图形的修饰与标注MATLAB 提供了一些特殊的函数修饰画出的图形,这些函数如下: 1)坐标轴的标题:title 函数其调用格式为:title('字符串')------字符串可以写中文如:title('My own plot')2)坐标轴的说明:xlabel 和ylabel 函数格式:xlabel('字符串')ylabel('字符串')如:xlabel('This is my X axis') ylabel('My Y axis')3)图形说明文字:text 和gtext 函数A.text 函数:按指定位置在坐标系中写出说明文字.格式为:text(x1, y1, '字符串', '选项') x1,y1 为指定点的坐标;'字符串'为要标注的文字;'选项'决定x1,y1 的坐标单位,如没有选项,则x1,y1 的坐标单位和图中一致;如选项为'sc', 则x1,y1 表示规范化窗口的相对坐标,其范围为0到1.如:text(1,2, '正弦曲线')B.gtext 函数:按照鼠标点按位置写出说明文字.格式为:gtext('字符串')当调用这个函数时,在图形窗口中出现一个随鼠标移动的大十字交叉线,移动鼠标将十字线的交叉点移动到适当的位置,点击鼠标左键,gtext 参数中的字符串就标注在该位置上.4)给图形加网格:grid 函数在调用时直接写grid 即可.上面的函数的应用实例:例:在图形中加注坐标轴标识和标题及在图形中的任意位置加入文本.t=0:pi/100:2*pi;y=sin(t);plot(t,y),grid,axis([0 2*pi -1 1])xlabel('0 leq itt rm leq pi','FontSize',16)ylabel('sin(t)','FontSize',20)title('正弦函数图形','FontName','隶书' ,'FontSize',20) text(pi,sin(pi),'leftarrowsin(t)=0','FontSize',16)text(3*pi/4,sin(3*pi/4),'leftarrowsin(t)=0.707','FontSize',16)text(5*pi/4,sin(5*pi/4),' sin(t)=-0.707rightarrow',... 'FontSize',16,'HorizontalAlignment','right')5)在图形中添加图例框:legend 函数其调用格式为:A.legend('字符串1', '字符串2', ……)------以字符串1,字符串2……作为图形标注的图例.B.legend('字符串1', '字符串2', ……, pos)------pos 指定图例框显示的位置.图例框被预定了6 个显示位置:0------取最佳位置;1------右上角(缺省值);2------左上角;3------左下角;4------右下角;-1------图的右侧.例:在图形中添加图例.x=0:pi/10:2*pi;y1=sin(x);y2=0.6*sin(x);y3=0.3*sin(x);plot(x,y1,x,y2,'-o',x,y3,'-*')legend('曲线1','曲线2','曲线3')6)用鼠标点选屏幕上的点:ginput 函数格式为:[x, y, button]=ginput(n)其中:n 为所选择点的个数;x,y 均为向量,x 为所选n 个点的横坐标;y 为所选n个点的纵坐标.button 为n 维向量,是所选n 个点所对应的鼠标键的标号:1------左键;2------中键;3------右键.可用不同的鼠标键来选点,以区别所选的点.此语句可以放在绘图语句之后,它可在绘出的图形上操作,选择你所感兴趣的点,如峰值点,达到稳态值的点等,给出点的坐标,可求出系统的性能指标.第四节MATLAB 下图形对象的修改MATLAB 图形对象是指图形系统中最基本,最底层的单元,这些对象包括:屏幕(Root), 图形窗口(Figures), 坐标轴(Axes), 控件(Uicontrol), 菜单(Uimenu),线(Lines),块(Patches),面(Surface),图像(Images),文本(Text)等等.根据各对象的相互关系,可以构成如下所示的树状层次:RootFiguresAxes Uicontrol Uimenu Uicontextmenu (对象菜单)Images Line Patch Surface Text对各种图形对象进行修改和控制,要使用MATLAB 的图形对象句柄(Handle).在MATLAB 中,每个图形对象创立时,就被赋予了唯一的标识,这个标识就是该对象的句柄.句柄的值可以是一个数,也可以是一个矢量.如每个计算机的根对象只有一个,它的句柄总是0,图形窗口的句柄总是正整数,它标识了图形窗口的序号等.利用句柄可以操纵一个已经存在的图形对象的属性,特别是对指定图形对象句柄的操作不会影响同时存在的其它图形对象,这一点是非常重要的.一.对图形对象的修改可以用下面函数:1)set 函数:用于设置句柄所指的图形对象的属性.Set 函数的格式为:set(句柄, 属性名1, 属性值1, 属性名2, 属性值2, ……) 例:h=plot(x,y)set(h, 'Color', [1,0,0])------将句柄所指曲线的颜色设为红色.2)get 函数:获取指定句柄的图形对象指定属性的当前值.格式为:get(句柄, '属性名')如: get(gca, 'Xcolor')------获得X 轴的当前颜色属性值. 执行后可返回X 轴的当前颜色属性值[0,0,0](黑色).3)如果没有设置句柄,则可以使用下列函数获得:gcf:获得当前图形窗口的句柄;gca:获得当前坐标轴对象的句柄;gco:获得当前对象的句柄.如:A.要对图形窗口的底色进行修改,可用:set(gcf, 'Color', [1,1,1])------将图形窗口底色设为白色B.要把当前X 轴的颜色改为绿色,可用:set(gca, 'Xcolor', [0,1,0])C.还可对坐标轴的显示刻度进行定义:t=-pi:pi/20:pi;y=sin(t);plot(t,y)set(gca,'xtick',[-pi:pi/2:pi],'xticklabel',['-pi','-pi/ 2','0','pi/2','pi'])本例中用'xtick'属性设置x 轴刻度的位置(从-pi~pi,间隔pi/2,共设置5 个点),用'xticklabel'来指定刻度的值,由于通常习惯于用角度度量三角函数,因此重新设置['-pi','-pi/2','0','pi/2','pi']5 个刻度值.二.一些常用的属性如下:1)Box 属性:决定图形坐标轴是否为方框形式,选项为'on'(有方框), 'off'(无方框);2)'ColorOrder'属性:设置多条曲线的颜色顺序,默认值为:[1 1 0;1 0 1;0 1 1;1 0 0;0 1 0;0 0 1]黄色粉色天蓝红色绿色兰色颜色向量还有:[1 1 1]------白色;[0 0 0]------黑色.3)坐标轴方向属性:'Xdir','Ydir','Zdir',其选项为:'normal'------正常'reverse'------反向4)坐标轴颜色和线型属性:'Xcolor','Ycolor','Zcolor'------ 轴颜色, 值为颜色向量如何在画好曲线后再在图上标刻度就是想在一些特定的点边上标上一串30.60.90~7200递增的数据,共有96个点要标!!im = imread(url);imshow(im)然后输入:text(100,100,'\o ','Color','red');matlab,用imread 读入一个图片,我想在图上的一些坐标点上做标记。
matlab中subplot函数的用法
matlab中subplot函数的用法subplot函数是MATLAB中用于绘制多个子图的函数,其主要作用是在一个大图中划分出若干个小的子图,并在每一个子图中进行相应的数据可视化操作。
该函数的基本语法如下:subplot(m,n,p)其中m和n表示子图布局的行数和列数,p表示当前子图所在的位置。
如果总共有m*n个子图,那么第一个子图位于第一行第一列,第二个子图位于第一行第二列,依次类推。
需要注意的是,所有的参数都需要使用逗号分隔,且参数p可以是一个矩阵。
以下将详细介绍subplot函数的用法和一些示例。
1.基本的subplot用法最简单的subplot函数用法是通过subplot(1,1,1)绘制单个子图,其效果与直接使用plot函数是一样的。
例子:x=0:0.1:10;y = sin(x);subplot(1,1,1)plot(x,y)这个例子中,我们在一个1×1的网格中设置一个子图。
然后我们使用plot函数绘制了x和y数据。
2.多个子图的使用a.多个子图梯形显示:为了在一个大图中绘制多个子图,我们需要使用subplot函数多次调用来设置每个子图的位置。
例子:x = -pi:0.1:pi;y1 = sin(x);y2 = cos(x);y3 = tan(x);subplot(2,2,1)plot(x,y1)subplot(2,2,2)plot(x,y2)subplot(2,2,3)plot(x,y3)在这个例子中,我们首先定义了x和y1、y2、y3这三组数据。
然后,我们使用subplot函数在一个2×2的网格中设置了三个子图。
b.多个子图网格显示:我们也可以通过修改subplot函数的参数,来设置子图的网格布局。
例子:t = 0:0.1:4*pi;y1 = sin(t);y2 = cos(t);y3 = tan(t);subplot(3,1,1)plot(t,y1)subplot(3,1,2)plot(t,y2)subplot(3,1,3)plot(t,y3)在这个例子中,我们将参数m设置为3,n设置为1,这样就形成了一个3×1的网格布局。
MATLAB中绘图命令介绍
MATLAB中绘图命令介绍本节将介绍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函数即可:hold on 保持当前图形,以便继续画图到当前坐标窗口hold off 释放当前图形窗口title(’图形名称’)(都放在单引号内)xlabel(’x轴说明’)ylabel(’y轴说明’)text(x,y,’图形说明’)legend(’图例1’,’图例2’,…)plot(x, sin(x), x, cos(x));若要改变颜色,在座标对後面加上相关字串即可:plot(x, sin(x), 'c', x, cos(x), 'g');若要同时改变颜色及图线型态,也是在座标对後面加上相关字串即可:plot(x, sin(x), 'co', x, cos(x), 'g*');小整理:plot绘图函数的叁数字元、颜色元、图线型态,y 黄色 .点k 黑色o 圆w 白色x xb 蓝色++g 绿色* *r 红色- 实线c 亮青色: 点线m锰紫色-. 点虚线-- 虚线plot3 三维曲线作图图形完成后,我们可用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围: axis([0, 6, -1.2, 1.2]);axis函数的功能丰富,其常用的用法有:axis equal :纵横坐标轴采用等长刻度axis square:产生正方形坐标系(默认为矩形)axis auto:使用默认设置axis off:取消坐标轴axis on :显示坐标轴此外,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还有其他各种二维绘图函数,以适合不同的应用,详见下表。
论文中对绘图的某个细节进行局部放大(大图套小图)
论⽂中对绘图的某个细节进⾏局部放⼤(⼤图套⼩图)
最近在准备新的论⽂,在实验结果中想对图⽚的某些细节进⾏放⼤处理,实际上就是⼤图中套⼩图的效果。
在⽹上找了⼀些⽅法(多数推荐使⽤magnify.m),但是都不是⾃⼰想要的结果,于是⾃⼰探索出了⼀种⽅法,供需要者参考:
环境:Matlab R2019b + visio 2016
1.在matlab环境中画图并“复制为向量图”(Matlab R2019b中⾃带),粘贴到visio中。
2.对matlab中绘制出来的图进⾏局部放⼤(matlab⾃带的放⼤功能),然后再将放⼤后的区域也“复制为向量图”,粘贴到visio中。
3.在visio中拖拽两者到合适的位置,并“组合”起来
4.最后保存成pdf,必要的时候使⽤Adobe Acrobat Pro DC对pdf格式的图⽚进⾏裁切。
MATLAB画图之多子图画法(subplot和自己确定大小位置两种方法)
MATLAB画图之多⼦图画法(subplot和⾃⼰确定⼤⼩位置两种⽅法)解决问题:在⼀个图中画多个⼦图,⼜能⾃⼰确定⼦图的⼤⼩和位置。
解决⽅法:有两种解决⽅法可以实现:1. 使⽤subplot命令2. 使⽤axes函数来重新规定⼦图的⼤⼩和位置。
这两种⽅法中,第⼀种⽅法相对简单,但是第⼆种⽅法更加灵活,具体使⽤如下:第⼀种⽅法:使⽤subplot先看⼀下subplot的简单使⽤,程序如下:clear;clc;close all;t = 0:0.001:10;y1 = sin(t);y2 = cos(t);figure (1);subplot(2,2,1)plot(t,y1);subplot(2,2,2)plot(t,y1);subplot(2,2,3)plot(t,y2);subplot(2,2,4)plot(t,y2);程序运⾏结果:使⽤subplot的这种⽅法,如果想要⾃定义⼦图的⼤⼩和位置该怎么设置?程序如下:clear;clc;close all;t = 0:0.001:10;y1 = sin(t);y2 = cos(t);figure(1);subplot('position',[0.2,0.7,0.6,0.2]);plot(t,y1);subplot('position',[0.2,0.2,0.6,0.2]);plot(t,y2);程序运⾏结果:第⼆种⽅法:直接使⽤axes函数程序如下:clear;clc;close all;t = 0:0.001:10;y1 = sin(t);y2 = cos(t);figure (1);axes('position',[0.1 0.6 0.3 0.3]);plot(t,y1);axes('position',[0.6 0.6 0.3 0.3]);plot(t,y1);axes('position',[0.1 0.1 0.3 0.3]);plot(t,y2);axes('position',[0.6 0.1 0.3 0.3]);plot(t,y2);程序运⾏结果:其中,'position',[0.6 0.1 0.3 0.3]的含义可以参考我之前的博客⽂章,有详细说明;如果要⾃定义整个figure的⼤⼩,设置gcf,我之前的博客⽂章也有详细说明。
如何在matlab中绘制图形
4.1 二维图形
二、subplot函数 函数
重新绘制上例4个图形,程序变动后如下:
x=linspace(0,2*pi,60); y=sin(x); z=cos(x); t=sin(x)./(cos(x)+eps); ct=cos(x)./(sin(x)+eps); axis ([0 2*pi -1 1]); H3=figure; 同上 plot(x,t); title('tangent(x)'); axis ([0 2*pi -40 40]); H4=figure; 同上 plot(x,ct); title('cotangent(x)'); axis ([0 2*pi -40 40]);
H1=figure;
%创建新窗口并返回句柄到变量H1
plot(x,y); %绘制图形并设置有关属性 title('sin(x)'); axis ([0 2*pi -1 1]);
H2=figure;
%创建第二个窗口并返回句柄到变量H2
plot(x,z); %绘制图形并设置有关属性 title('cos(x)')
grid on; title 'semilogy-logy=x^3';
4.2特殊坐标图形
(二)单对数坐标
以X轴为对数重新绘制上述曲线,程序为:
semilogx(x,y); 单对数X轴绘图命令
同样,可以以Y轴为对数重新绘制上述曲线,程序为:
semilogy(x,y); 单对数Y轴绘图命令
4.2特殊坐标图形
阅读如下程序:
x=linspace(0,2*pi,60); y=sin(x); z=cos(x); plot(x,y,'b');