source: MML/trunk/mml/doc_html/mml/plotbpmpsd.html @ 4

Last change on this file since 4 was 4, checked in by zhangj, 10 years ago

Initial import--MML version from SOLEIL@2013
File size: 8.6 KB
Line 
1<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"
2                "http://www.w3.org/TR/REC-html40/loose.dtd">
3<html>
4<head>
5  <title>Description of plotbpmpsd</title>
6  <meta name="keywords" content="plotbpmpsd">
7  <meta name="description" content="PLOTBPMPSD - Plots data taken will getbpmpsd">
8  <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
9  <meta name="generator" content="m2html &copy; 2003 Guillaume Flandin">
10  <meta name="robots" content="index, follow">
11  <link type="text/css" rel="stylesheet" href="../m2html.css">
12</head>
13<body>
14<a name="_top"></a>
15<div><a href="../index.html">Home</a> &gt;  <a href="index.html">mml</a> &gt; plotbpmpsd.m</div>
16
17<!--<table width="100%"><tr><td align="left"><a href="../index.html"><img alt="<" border="0" src="../left.png">&nbsp;Master index</a></td>
18<td align="right"><a href="index.html">Index for mml&nbsp;<img alt=">" border="0" src="../right.png"></a></td></tr></table>-->
19
20<h1>plotbpmpsd
21</h1>
22
23<h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
24<div class="box"><strong>PLOTBPMPSD - Plots data taken will getbpmpsd</strong></div>
25
26<h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
27<div class="box"><strong>function [Pyy, fy, Pyy_Int, FileName] = plotbpmpsd(FileName, LineColorInput, HoldFlag1) </strong></div>
28
29<h2><a name="_description"></a>DESCRIPTION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
30<div class="fragment"><pre class="comment">PLOTBPMPSD - Plots data taken will getbpmpsd
31  [Pxx, f, PxxInt, FileName] = plotbpmpsd(FileName, LineType, HoldFlag)
32
33  INPUTS
34  1. FileName - Filename or '' to browse
35  2. LineType - Line type, like '-r'
36  3. HoldFlag - 0   -&gt; new plot
37               else -&gt; hold last plot
38
39  OUTPUTS
40  1. Pxx - Power spectrum
41  2. f - Frequency vector
42  3. PxxInt - Cumulated integrated power spectrum
43  4. FileName - Filename used
44
45   Written by Greg Portmann</pre></div>
46
47<!-- crossreference -->
48<h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
49This function calls:
50<ul style="list-style-image:url(../matlabicon.gif)">
51</ul>
52This function is called by:
53<ul style="list-style-image:url(../matlabicon.gif)">
54</ul>
55<!-- crossreference -->
56
57
58<h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>
59<div class="fragment"><pre>0001 <a name="_sub0" href="#_subfunctions" class="code">function [Pyy, fy, Pyy_Int, FileName] = plotbpmpsd(FileName, LineColorInput, HoldFlag1)</a>
600002 <span class="comment">%PLOTBPMPSD - Plots data taken will getbpmpsd</span>
610003 <span class="comment">%  [Pxx, f, PxxInt, FileName] = plotbpmpsd(FileName, LineType, HoldFlag)</span>
620004 <span class="comment">%</span>
630005 <span class="comment">%  INPUTS</span>
640006 <span class="comment">%  1. FileName - Filename or '' to browse</span>
650007 <span class="comment">%  2. LineType - Line type, like '-r'</span>
660008 <span class="comment">%  3. HoldFlag - 0   -&gt; new plot</span>
670009 <span class="comment">%               else -&gt; hold last plot</span>
680010 <span class="comment">%</span>
690011 <span class="comment">%  OUTPUTS</span>
700012 <span class="comment">%  1. Pxx - Power spectrum</span>
710013 <span class="comment">%  2. f - Frequency vector</span>
720014 <span class="comment">%  3. PxxInt - Cumulated integrated power spectrum</span>
730015 <span class="comment">%  4. FileName - Filename used</span>
740016 <span class="comment">%</span>
750017 <span class="comment">%   Written by Greg Portmann</span>
760018
770019
780020 <span class="keyword">if</span> nargin &lt; 1
790021     FileName = [];
800022 <span class="keyword">end</span>
810023 <span class="keyword">if</span> nargin &lt; 2
820024     LineColorInput = <span class="string">'b'</span>;
830025 <span class="keyword">end</span>
840026 <span class="keyword">if</span> nargin &lt; 3
850027     HoldFlag1 = 0;
860028 <span class="keyword">end</span>
870029
880030 <span class="keyword">if</span> isempty(FileName)
890031     [FileName, DirectoryName] = uigetfile(<span class="string">'*.mat'</span>, <span class="string">'Select a configuration file to load'</span>);
900032     <span class="keyword">if</span> FileName == 0
910033         <span class="keyword">return</span>
920034     <span class="keyword">else</span>
930035         load([DirectoryName FileName]);
940036     <span class="keyword">end</span>
950037 <span class="keyword">else</span>
960038     load(FileName);
970039 <span class="keyword">end</span>
980040
990041 <span class="keyword">if</span> exist(<span class="string">'PxxAvg'</span>, <span class="string">'var'</span>)
1000042     Pxx = PxxAvg;
1010043     Pyy = PyyAvg;
1020044     Xrms = XAvgRMS;
1030045     Yrms = YAvgRMS;
1040046 <span class="keyword">end</span>
1050047
1060048 newList = editlist(DeviceList, <span class="string">'BPM'</span>, zeros(size(DeviceList,1),1));
1070049 <span class="keyword">if</span> isempty(newList)
1080050     <span class="keyword">return</span>
1090051 <span class="keyword">end</span>
1100052 j = findrowindex(newList(1,:), DeviceList);
1110053
1120054 <span class="keyword">if</span> ~HoldFlag1
1130055     clf reset
1140056 <span class="keyword">end</span>
1150057
1160058 subplot(2,2,1);
1170059 <span class="keyword">if</span> HoldFlag1
1180060     hold on;
1190061 <span class="keyword">end</span>
1200062 NumFreqRemove = length(find(Pxx(j,:)==0));
1210063 loglog(fx(NumFreqRemove+1:end), 1e12*Pxx(j,NumFreqRemove+1:end), LineColorInput);
1220064 <span class="keyword">if</span> HoldFlag1
1230065     title(sprintf(<span class="string">'BPMx POWER SPECTRUM (%d points)'</span>,length(fx)));
1240066 <span class="keyword">else</span>
1250067     title(sprintf(<span class="string">'BPMx(%d,%d) POWER SPECTRUM (%d points)'</span>,DeviceList(j,:),length(fx)));
1260068 <span class="keyword">end</span>
1270069 xlabel(<span class="string">'Frequency [Hz]'</span>);
1280070 ylabel(<span class="string">'BPMx PSD [\mum^2/Hz]'</span>);
1290071 grid on;
1300072 legend off;
1310073 aa=axis;
1320074 axis([1 fx(end) aa(3) aa(4)]);
1330075 <span class="comment">%axis([1 fx(end) 1e-3 10]);</span>
1340076
1350077 <span class="comment">% Position spectrum</span>
1360078 subplot(2,2,3);
1370079 <span class="keyword">if</span> HoldFlag1
1380080     hold on;
1390081 <span class="keyword">end</span>
1400082 semilogx(fx(NumFreqRemove:end), 1e12*Pxx_Int(j,NumFreqRemove:end), LineColorInput);
1410083 <span class="keyword">if</span> HoldFlag1
1420084     title(sprintf(<span class="string">'BPMx Integrated PSD'</span>));
1430085 <span class="keyword">else</span>
1440086     title(sprintf(<span class="string">'BPMx(%d,%d) Integrated PSD (RMS=%.1f \\mum)'</span>,DeviceList(j,:),1e6*Xrms(j)));
1450087 <span class="keyword">end</span>
1460088 xlabel(<span class="string">'Frequency [Hz]'</span>);
1470089 ylabel(<span class="string">'Mean Square Displacement [\mum^2]'</span>);
1480090 grid on;
1490091 legend off;
1500092 aa=axis;
1510093 axis([1 fx(end) aa(3) aa(4)]);
1520094
1530095 subplot(2,2,2);
1540096 <span class="keyword">if</span> HoldFlag1
1550097     hold on;
1560098 <span class="keyword">end</span>
1570099 loglog(fy(NumFreqRemove+1:end), 1e12*Pyy(j,NumFreqRemove+1:end), LineColorInput);
1580100 <span class="keyword">if</span> HoldFlag1
1590101     title(sprintf(<span class="string">'BPMy POWER SPECTRUM (%d points)'</span>,length(fy)));
1600102 <span class="keyword">else</span>
1610103     title(sprintf(<span class="string">'BPMy(%d,%d) POWER SPECTRUM (%d points)'</span>,DeviceList(j,:),length(fy)));
1620104 <span class="keyword">end</span>
1630105 xlabel(<span class="string">'Frequency [Hz]'</span>);
1640106 ylabel(<span class="string">'BPMy PSD [\mum^2/Hz]'</span>);
1650107 grid on;
1660108 legend off;
1670109 aa=axis;
1680110 axis([1 fx(end) aa(3) aa(4)]);
1690111 <span class="comment">%axis([1 fx(end) 1e-3 10]);</span>
1700112
1710113 <span class="comment">% Position spectrum</span>
1720114 subplot(2,2,4);
1730115 <span class="keyword">if</span> HoldFlag1
1740116     hold on;
1750117 <span class="keyword">end</span>
1760118 semilogx(fy(NumFreqRemove:end), 1e12*Pyy_Int(j,NumFreqRemove:end), LineColorInput);
1770119 <span class="keyword">if</span> HoldFlag1
1780120     title(sprintf(<span class="string">'BPMy Integrated PSD'</span>));
1790121 <span class="keyword">else</span>
1800122     title(sprintf(<span class="string">'BPMy(%d,%d) Integrated PSD (RMS=%.1f \\mum)'</span>,DeviceList(j,:),1e6*Yrms(j)));
1810123     addlabel(1,0,datestr(TimeStart));
1820124 <span class="keyword">end</span>
1830125 xlabel(<span class="string">'Frequency [Hz]'</span>);
1840126 ylabel(<span class="string">'Mean Square Displacement [\mum^2]'</span>);
1850127 grid on;
1860128 legend off;
1870129 aa=axis;
1880130 axis([1 fx(end) aa(3) aa(4)]);
1890131</pre></div>
190<hr><address>Generated on Mon 21-May-2007 15:29:18 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> &copy; 2003</address>
191</body>
192</html>
Note: See TracBrowser for help on using the repository browser.