Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
function radrms(fconf)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Function:
% Package: HoP
% Author: Franco Buffa
% Purpose:
% Date:
%
% GNU Octave Version: 4.0.2
%
% revised by
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% is a function that
%
% Output:
%
%
% Input:
%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if nargin>0
run(fconf);
else
run('conf.m');
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
load zern.dat
% i0=find(abs(u3)>tol4);
% u3(i0)=sign(u3(i0))*tol4;
x0=((1:mp)-mp/2);
[xm ym]=meshgrid(x0,x0);
R=sqrt(xm.^2+ym.^2);
d=5;
x0=y=[];
for r=r0:r1-d
z=(R>=r);
z=z.*(R<r+d);
i0=find(z>0);
% imagesc(z)
y=[y ;std(u3(i0))];
x0=[x0 ; r];
end
dx=D/(mp-1)/k0;
x=x0*dx;
y(end-d:end)=y(end-d);
figure(9);
plot(x,y,'*-');
grid minor on;
axis([4.5 32]);
end