Basics of MATLAB and Beyond

More documents

Recommendations

Info

$Algoritmos de para calcular la dimensiÃ³n fractal de seÃ±ales 1D La ...$

ut it will be high quality when printed.screen, but are lower quality when printed.Bitmaps do appear on theBitmap Import On PCs choose matlab’s -dmeta option when usingthe print command to produce a file with the extension .emf. Thenfrom CorelDraw select “File”→“Import ... ” to get a file finder dialogbox. Select the .emf file you created with matlab and choose the“Windows Metafile (WMF)” option in the “Files of Type” box. Clickthe “Import” button and CorelDraw will return you to your documentwhere you must drag out a rectangle to define the size of the importedgraphic. Once you drag the rectangle, the matlab graphic will appearin bitmap form.Vector Graphic Import Another option with CorelDraw is to importan HPGL file printed using the -dhpgl option from matlab. This is aformat designed for Hewlett-Packard pen plotters, and has the advantagethat the imported image can be edited from within CorelDraw. Followthe steps as for bitmaps above, but use the “HPGL Plotter File (PLT)”option in the “Files of Type” box. You have the option here of changingthe pen colours used by CorelDraw’s emulation of an HP pen plotter.This format works well for graphics that could easily be plotted on apen plotter. Line drawings work well but images do not. Surface plotsare initially imported with no hidden line removal, but if you change thefill colour to white with CorelDraw, hidden lines will be removed. To dothis fill colour change, ungroup the graphic, select the surface, click theright mouse button, and choose “Properties”, then choose white as thefill colour. If you have problems importing with this format, try openingthe .hgl file in a text editor and deleting the last line.36 Irregular GridsIn this section we discuss how to deal with data that is defined overan irregular grid. In Section 7.6 on page 32 we saw how do handlenonrectangular grids, but these were still regular in the sense that thex and y data grids could be defined using rectangular matrices (thoughthe resulting geometrical domains did not have to be rectangles).Some irregular three-dimensional data is supplied with matlab inthe data file seamount.mat. Load the data and plot the points:>> load seamount>> whosName Size ...caption 1x229 ...c○ 2000 by CRC Press LLC
x 294x1 ...y 294x1 ...z 294x1 ...Grand total is 1111 elements>> plot3(x,y,z,’.’)>> axis tightHere we have 294 measurements of x (latitude), y (longitude), and z(height above sea level, which is negative). They represent measurementsof a mountain under the sea. 17 The data are stored as column vectorsof x, y, and z values. Suppose we want to plot a surface and a contourmap representing this seamount data. If you try to type surf(x,y,z)or contour(x,y,z) with this data, you will not get any meaningful plot.Two ways of generating the desired plots are as follows:1. Interpolate the data over a rectangular grid.2. Use triangular gridding instead of rectangular gridding.Let us look at each of these.36.1 Interpolation over a Rectangular GridWe continue the example above and define vectors of uniformly spacedpoints between the minimum and maximum values of x and y:xiv = linspace(min(x),max(x),50);yiv = linspace(min(y),max(y),50);Each of these vectors has 50 elements. We now use the griddata interpolationfunction to do two things: (1) create matrices of the x and ygrids that correspond to a rectangular grid over the vectors xiv and yiv,and (2) interpolate the data over this new rectangular grid. In the callto griddata that follows, we need to transpose the vector yiv becausegriddata expects it to be a column vector in this case:>> [xi,yi,zi] = griddata(x,y,z,xiv,yiv’);>> whosName Size Bytes Classcaption 1x229 458 char array17 The reference can be found by typing the caption variable:>> captioncaption=Parker, R. L., Shure, L. & Hildebrand, J., "The application of inversetheory to seamount magnetism", Reviews of Geophysics vol 25, pp 17-40,1987. x is latitude (degrees), y is longitude (degrees), z is negativedepth (meters).c○ 2000 by CRC Press LLC
Page 2 and 3:
BASICS OFMATLAB®and Beyondc○ 200
Page 4 and 5:
Library of Congress Cataloging-in-P
Page 6 and 7:
typing pwd (print working directory
Page 8 and 9:
ContentsIBasics of MATLAB1 First St
Page 10 and 11:
20 Time-Frequency Analysis21 Line A
Page 12 and 13:
39Answers to Exercises (Part I)40 A
Page 14 and 15:
Set → Getting Started with MATLAB
Page 16 and 17:
1.5 The Colon OperatorTo generate a
Page 18 and 19:
matlab will recall the last command
Page 20 and 21:
a = zeros(2,3)a =0 0 00 0 0>> b = o
Page 22 and 23:
[a a(a)]ans =1 2 3 1 4 74 5 6 2 5 8
Page 24 and 25:
a + bans =11 3344 6677 993.8 Transp
Page 26 and 27:
4.2 Adding PlotsWhen you issue a pl
Page 28 and 29:
x = [[1 2 3 4]’ [2 3 4 5]’ [3 4
Page 30 and 31:
4.7 AxesSo far we have allowed matl
Page 32 and 33:
v =1122>> v*uans =1 2 0 11 2 0 12 4
Page 34 and 35:
max returns a vector containing the
Page 36 and 37:
s = 100*spiral(3)s =700 800 900600
Page 38 and 39:
We can plot the surface z as a func
Page 40 and 41:
The contour function plots the cont
Page 42 and 43:
m = gray(8);colormap(m)imagesc(1:10
Page 44 and 45:
hold onplot3(x(ind),y(ind),z(ind),
Page 46 and 47:
z = 5*ones(3,3)z =5 5 55 5 55 5 5>>
Page 48 and 49:
contour(x,y,z,30);You may have noti
Page 50 and 51:
variables that are local to themsel
Page 52 and 53:
SwitchThe basic form of a switch st
Page 54 and 55:
i =i =123136Vectorised Codematlab i
Page 56 and 57:
9.1 MATLAB FormatTo save all the va
Page 58 and 59:
11 StartupEach time you start matla
Page 60 and 61:
Where p i is the population for yea
Page 62 and 63:
NaNs by leaving them off the plot.
Page 64 and 65:
clfpolar(t,gdb)In this case you mus
Page 66 and 67:
-0.5000-0.3750-0.2500-0.125000.1250
Page 68 and 69:
the periodograms of the blocks, and
Page 70 and 71:
plot(y(1:1000))axis([0 1000 -5 5])z
Page 72 and 73:
comet(x,y)(You can get a three-dime
Page 74 and 75:
By clicking within the panner box a
Page 76 and 77:
x = linspace(-1,1,8192);Fs = 1000;y
Page 78 and 79:
23.3 Graphical Object Hierarchymatl
Page 80 and 81:
Now we need to get the handles of a
Page 82 and 83:
Beyond the BasicsPreludeThis part o
Page 84 and 85:
clfgplot(A,[x y])axis off(Try zoomi
Page 86 and 87:
str =you’re the one>> str = ’
Page 88 and 89:
strrep(str,’go’,’am’)ans =I
Page 90 and 91:
ans =3.14159e+100Some additional te
Page 92 and 93:
eval(str)v =1 2 3 4 5The eval(str)
Page 94 and 95:
Inline objects, like every other ma
Page 96 and 97:
a(1,1) = {[1 2 3]}a =[1x3double]or>
Page 98 and 99:
t = {’help’ spiral(3) ; eye(2)
Page 101 and 102:
meteo.Temperatureans =24ans =19.000
Page 103 and 104:
3 >> a = [1 2 ;4 5 6;7 8 9]a =1 2 3
Page 105 and 106:
[x,y] = meshgrid(1:5,1:3)x =1 2 3 4
Page 107 and 108:
ans(:,:,2) =999And to sum over “p
Page 109 and 110:
An Application of RGB ImagesTo see
Page 111 and 112:
You can see that over the 8 time st
Page 113 and 114:
name>> staff(2,1,2).name(5:9)ans =B
Page 115 and 116:
q =7 8 96 1 25 4 3>> save saved_dat
Page 117 and 118: fid = fopen(’asc.dat’);The fope
Page 119 and 120: display and the enclosed area is th
Page 121 and 122: MarkerEdgeColor = auto(and so on)In
Page 123 and 124: text(-.7,f(-.7),’f(X)’)But we h
Page 125 and 126: plt(x,f(x/2),’--’)Let us try ca
Page 127 and 128: Example: findobjThe findobj command
Page 129 and 130: 3. Add the property you want to set
Page 131 and 132: is numbered from top to bottom. The
Page 133 and 134: axis ijUsually images like this do
Page 135 and 136: 32.2 Tick Marks and Labelsmatlab’
Page 137 and 138: clfaxes(’pos’,[.2 .1 .7 .4])x =
Page 139 and 140: property to [0 top] in the call to
Page 141 and 142: object’s HorizontalAlignment and
Page 143 and 144: the required symbols. The m-file ti
Page 145 and 146: 1. Use graphical elements that serv
Page 147 and 148: uicontrol(’Callback’,’ezplot(
Page 149 and 150: uicontrol(’Pos’,[110 280 60 19]
Page 151 and 152: ’String’,’0.5’,’Style’,
Page 153 and 154: function exradio(action)if nargin =
Page 155 and 156: fail. One place to put variables th
Page 157 and 158: 34.7 Fast DrawingFor fast drawing o
Page 159 and 160: controlled by guide and which are a
Page 161 and 162: InvertHardcopy: [ {on} | off ]Paper
Page 163 and 164: Another option is to use an image f
Page 165 and 166: yur is the y coordinate of the uppe
Page 167: \end{center}\caption{This is the fi
Page 171 and 172: diagram. Such a triangular grid can
Page 173 and 174: xt = x([1 9 10 1]);yt = y([1 9 10 1
Page 175 and 176: x = [0 0 .5 0.5 .5 1 11 .5 .5 .5];y
Page 177 and 178: fac2 = [1 2 61 6 31 4 61 6 5];The r
Page 179 and 180: t = linspace(0,2*pi,20);x = cos(t);
Page 181 and 182: N = 20;dt = 2*pi/N;t = 0:dt:(N-1)*d
Page 183 and 184: colouring of the sphere itself. To
Page 185 and 186: material dullThe default is materia
Page 187 and 188: time taken by this routine to calcu
Page 189 and 190: first row of a by taking the column
Page 191 and 192: 38.3 Debuggingmatlab has a suite of
Page 193 and 194: 10: y = linspace(ycentre - L,ycentr
Page 195 and 196: We change variables:p ′ = C + Bx.
Page 197 and 198: Exercise 6 (Page 59)You can generat
Page 199 and 200: 33 ! 36 $ 39 ’ 42 * 45 - 48 034 "
Page 201 and 202: R = 1; N = 200;[x,y] = meshgrid(lin
Page 203 and 204: (The plot might take a few seconds
Page 205: % Reproduce the layers at the diffe
show all

Basics of MATLAB and Beyond

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?