source: MML/trunk/machine/SOLEIL/doc_html/machine/Soleil/StorageRing/modelnonlindisp.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: 5.3 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 modelnonlindisp</title>
6  <meta name="keywords" content="modelnonlindisp">
7  <meta name="description" content="">
8  <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
9  <meta name="generator" content="m2html v1.5 &copy; 2003-2005 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="#">machine</a> &gt; <a href="#">Soleil</a> &gt; <a href="index.html">StorageRing</a> &gt; modelnonlindisp.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 machine/Soleil/StorageRing&nbsp;<img alt=">" border="0" src="../../../right.png"></a></td></tr></table>-->
19
20<h1>modelnonlindisp
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></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 modelnonlindisp </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">
31% Written by Laurent S. Nadolski</pre></div>
32
33<!-- crossreference -->
34<h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
35This function calls:
36<ul style="list-style-image:url(../../../matlabicon.gif)">
37</ul>
38This function is called by:
39<ul style="list-style-image:url(../../../matlabicon.gif)">
40</ul>
41<!-- crossreference -->
42
43
44
45<h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../../../up.png"></a></h2>
46<div class="fragment"><pre>0001 <a name="_sub0" href="#_subfunctions" class="code">function modelnonlindisp</a>
470002
480003 <span class="comment">%</span>
490004 <span class="comment">%% Written by Laurent S. Nadolski</span>
500005
510006 kmax = 50;
520007 deltaRF = 50e-6; <span class="comment">% Hz</span>
530008 alpha1 = modelmcf;
540009
550010 hMat(:,:) = zeros((2*kmax+1),length(getx));
560011 vMat(:,:) = hMat(:,:);
570012 rfVect = zeros((2*kmax+1),1);
580013
590014
600015 steprf(-deltaRF*(kmax+1));
610016     
620017 <span class="keyword">for</span> k=1:(2*kmax+1),
630018     steprf(deltaRF)
640019     hMat(k,:) = getx;
650020     vMat(k,:) = getz;
660021     rfVect(k) = getrf;
670022 <span class="keyword">end</span>
680023
690024 steprf(-deltaRF*kmax);
700025
710026 <span class="comment">%%</span>
720027 plot((rfVect-rfVect(kmax+1))*1e6, hMat(:,1)*alpha1*getrf/(2*deltaRF),<span class="string">'.-'</span>);
730028 grid on
740029
750030 <span class="comment">%%</span>
760031
770032 refreshthering;
780033 <span class="keyword">global</span> THERING;
790034
800035 [CavityState, PassMethod, iCavity] = getcavity;
810036
820037 kmax = 50;
830038 deltaRF = 40; <span class="comment">% Hz</span>
840039 clear hMat rfVect vMat
850040 hMat(:,:) = zeros((2*kmax+1),length(THERING)+1);
860041 vMat(:,:) = hMat(:,:);
870042 rfVect = zeros((2*kmax+1),1);
880043
890044 CavityFrequency = THERING{iCavity(1)}.Frequency;
900045
910046 <span class="keyword">for</span> k=1:(2*kmax+1),
920047     <span class="keyword">for</span> kk = 1:length(iCavity)
930048         THERING{iCavity(kk)}.Frequency = CavityFrequency + deltaRF*(k - (kmax +1));
940049     <span class="keyword">end</span>
950050     rep  = findorbit6(THERING, 1:length(THERING)+1);
960051     hMat(k,:) = rep(1,:);
970052     vMat(k,:) = rep(3,:);
980053     rfVect(k) = deltaRF*(k - (kmax +1));
990054 <span class="keyword">end</span>
1000055
1010056 <span class="keyword">for</span> kk = 1:length(iCavity)
1020057     THERING{iCavity(kk)}.Frequency = CavityFrequency;
1030058 <span class="keyword">end</span>
1040059
1050060 <span class="comment">% for k=1:kmax,</span>
1060061 <span class="comment">%   deltahOrbit(k,:) = hMat(end-(k+1),:) - hMat(k,:);</span>
1070062 <span class="comment">%   deltaRFVect(k) =   rfVect(end-(k+1)) - rfVect(k);</span>
1080063 <span class="comment">% end</span>
1090064     
1100065 <span class="comment">%%</span>
1110066 order = 4;
1120067 order4 = order-3;
1130068 order3 = order-2;
1140069 order2 = order-1;
1150070 order1 = order;
1160071
1170072 spos = findspos(THERING,1:length(THERING)+1);
1180073 figure;
1190074 plot(rfVect, hMat(:,1),<span class="string">'.'</span>);
1200075
1210076 clear p;
1220077 <span class="keyword">for</span> k = 1:length(THERING)+1,
1230078     p(k,:) = polyfit(rfVect, hMat(:,k), order);
1240079 <span class="keyword">end</span>
1250080 <span class="comment">%%</span>
1260081
1270082 pdelta = -p*alpha1*CavityFrequency;
1280083
1290084 figure;
1300085 subplot(4,1,1)
1310086 plot(spos, pdelta(:,order1))
1320087 subplot(4,1,2)
1330088 plot(spos, pdelta(:,order2))
1340089 subplot(4,1,3)
1350090 plot(spos, pdelta(:,order3))
1360091 <span class="comment">% subplot(4,1,4)</span>
1370092 <span class="comment">% plot(spos, pdelta(:,order4))</span></pre></div>
138<hr><address>Generated on Fri 23-Jul-2010 00:42:07 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/" title="Matlab Documentation in HTML">m2html</a></strong> &copy; 2005</address>
139</body>
140</html>
Note: See TracBrowser for help on using the repository browser.