1 | ls; |
---|
2 | load obs_float_test.dat; |
---|
3 | whos; |
---|
4 | obs_float_test (:,:) |
---|
5 | x=obs_float_test (1:2:24,4) |
---|
6 | y=obs_float_test (2:2:24,4) |
---|
7 | x=x+1 |
---|
8 | y=y+1 |
---|
9 | |
---|
10 | plot(x(1),y(1),'+') |
---|
11 | hold on; |
---|
12 | plot(x(12),y(12),'>') |
---|
13 | hold on; |
---|
14 | plot(x(1:12),y(1:12),'+') |
---|
15 | hold on; |
---|
16 | plot(x(1:12),y(1:12),'-b') |
---|
17 | hold on;1 |
---|
18 | |
---|
19 | |
---|
20 | load r_float_optim25float-filtre15-inc_16.dat; |
---|
21 | whos; |
---|
22 | r_float_optim25float_filtre15_inc_16(:,:) |
---|
23 | x1=r_float_optim25float_filtre15_inc_16(1:8:24,4) |
---|
24 | y1=r_float_optim25float_filtre15_inc_16(5:8:24,4) |
---|
25 | x1=x1+1 |
---|
26 | y1=y1+1 |
---|
27 | plot(x1(1),y1(1),'+') |
---|
28 | hold on; |
---|
29 | plot(x1(3),y1(3),'>') |
---|
30 | hold on; |
---|
31 | plot(x1(1:3),y1(1:3),'+') |
---|
32 | plot(x1(1:3),y1(1:3),'g') |
---|
33 | hold on; |
---|
34 | ls; |
---|
35 | r_float_optim25float_filtre15_inc_16(:,:) |
---|
36 | x2=r_float_optim25float_filtre15_inc_16(2:8:24,4) |
---|
37 | y2=r_float_optim25float_filtre15_inc_16(6:8:24,4) |
---|
38 | x2=x2+1 |
---|
39 | y2=y2+1 |
---|
40 | plot(x2(1),y2(1),'+') |
---|
41 | hold on; |
---|
42 | plot(x2(3),y2(3),'>') |
---|
43 | hold on; |
---|
44 | plot(x2(1:3),y2(1:3),'+') |
---|
45 | plot(x2(1:3),y2(1:3),'g') |
---|
46 | hold on; |
---|
47 | ls; |
---|
48 | r_float_optim25float_filtre15_inc_16(:,:) |
---|
49 | x3=r_float_optim25float_filtre15_inc_16(3:8:24,4) |
---|
50 | y3=r_float_optim25float_filtre15_inc_16(7:8:24,4) |
---|
51 | x3=x3+1 |
---|
52 | y3=y3+1 |
---|
53 | plot(x3(1),y3(1),'+') |
---|
54 | hold on; |
---|
55 | plot(x3(3),y3(3),'>') |
---|
56 | hold on; |
---|
57 | plot(x3(1:3),y3(1:3),'+') |
---|
58 | plot(x3(1:3),y3(1:3),'g') |
---|
59 | hold on; |
---|
60 | ls; |
---|
61 | r_float_optim25float_filtre15_inc_16(:,:) |
---|
62 | x4=r_float_optim25float_filtre15_inc_16(4:8:24,4) |
---|
63 | y4=r_float_optim25float_filtre15_inc_16(8:8:24,4) |
---|
64 | x4=x4+1 |
---|
65 | y4=y4+1 |
---|
66 | plot(x4(1),y4(1),'+') |
---|
67 | hold on; |
---|
68 | plot(x4(3),y4(3),'>') |
---|
69 | hold on; |
---|
70 | plot(x4(1:3),y4(1:3),'+') |
---|
71 | plot(x4(1:3),y4(1:3),'g') |
---|
72 | hold on; |
---|
73 | ls; |
---|
74 | |
---|
75 | |
---|
76 | load uv_total.dat |
---|
77 | |
---|
78 | X=uv_total(1:5046,end-3); |
---|
79 | Y=uv_total(1:5046,end-2); |
---|
80 | U1=uv_total(1:5046,end-1); |
---|
81 | V1=uv_total(1:5046,end); |
---|
82 | U2=uv_total(5047:10092,end-1); |
---|
83 | V2=uv_total(5047:10092,end); |
---|
84 | U3=uv_total(10093:15138,end-1); |
---|
85 | V3=uv_total(10093:15138,end); |
---|
86 | |
---|
87 | U=(U1+U2+U3)/3; |
---|
88 | V=(V1+V2+V3)/3; |
---|
89 | |
---|
90 | quiver(X,Y,U,V) |
---|
91 | |
---|
92 | |
---|
93 | |
---|
94 | |
---|
95 | % for k=1:16 |
---|
96 | % load uzero.dat; |
---|
97 | % whos; |
---|
98 | % X=uzero(:,end-2); |
---|
99 | % Y=uzero(:,end-1); |
---|
100 | % U=uzero(:,end); |
---|
101 | % ls; |
---|
102 | % load vzero.dat; |
---|
103 | % whos; |
---|
104 | % V=vzero(:,end); |
---|
105 | % figure(k); |
---|
106 | % quiver(X,Y,U,V); |
---|
107 | % end |
---|