Last change
on this file since 17 was
17,
checked in by zhangj, 11 years ago
|
To have a stable version on the server.
|
-
Property svn:executable set to
*
|
File size:
978 bytes
|
Line | |
---|
1 | %pertes sur chambre |
---|
2 | |
---|
3 | dx=[-12 -8 -3 -2 0 2 3 5 ]; |
---|
4 | dr=[0.12 8.2 8.5 8.5 8 3.4 1 0 ]; |
---|
5 | dr=dr/8.5; |
---|
6 | |
---|
7 | |
---|
8 | a=10; |
---|
9 | sig=2.8; |
---|
10 | i=0; |
---|
11 | |
---|
12 | for x=-a:0.1:a |
---|
13 | i=i+1; |
---|
14 | p(i)=x-4; |
---|
15 | ax=a-abs(x); |
---|
16 | r(i)=(1-exp(-0.5*(ax/sig)*(ax/sig))); |
---|
17 | end |
---|
18 | |
---|
19 | plot(dx,dr,'-ob' ,p,r,'-k') |
---|
20 | |
---|
21 | dn= [ 3 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21]; |
---|
22 | dx1=[5 5 5 5 5 3 5 5 5 5 5 5 5 5 4 5 5 5]; |
---|
23 | dr1=[5.8 3.3 1.9 3.45 3 0.95 5. 3.8 7.3 2.5 7.3 1.4 6.4 0.2 1.7 5.6 0 4]; |
---|
24 | |
---|
25 | x0=9-sig*sqrt(-2*log(1-dr1/8.5)); |
---|
26 | if (abs(x0-dx1)<abs(-x0-dx1)) |
---|
27 | dx0=x0-dx1; |
---|
28 | else |
---|
29 | dx0=-x0-dx1; |
---|
30 | end |
---|
31 | |
---|
32 | %mesure BPM |
---|
33 | s_bpmx=[14.2382 28.4764 50.2337 64.4719 92.5483 106.7865 128.5438 142.7820]; |
---|
34 | Zm=[-0.5356 0.2893 -0.3448 0.1526 -0.096 -0.2161 0.0730 -0.0045]; |
---|
35 | |
---|
36 | %ds1=[1:22]*3.55*2; |
---|
37 | dx=(dn-1)*3.55*2; |
---|
38 | s=[err_beta(:,1) ; err_beta(:,1)+78.31]; |
---|
39 | sigz1=[err_beta(:,2) ; err_beta(:,2)]*1000*1.2; |
---|
40 | sigz2=-sigz1; |
---|
41 | figure(2) |
---|
42 | plot(dx,dx0,'-ob',s_bpmx,Zm,'-or',s,sigz1,'-k',s,sigz2,'-k') |
---|
43 | |
---|
Note: See
TracBrowser
for help on using the repository browser.