source: altifloat/matlab_toolbox/meanof4.m @ 199

Last change on this file since 199 was 129, checked in by jbrlod, 10 years ago

last version of Varanth

  • Property svn:executable set to *
File size: 1.9 KB
Line 
1ls;
2load obs_float_test.dat;
3whos;
4obs_float_test (:,:)
5x=obs_float_test (1:2:24,4)
6y=obs_float_test (2:2:24,4)
7x=x+1
8y=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
20load r_float_optim25float-filtre15-inc_16.dat;
21whos;
22r_float_optim25float_filtre15_inc_16(:,:)
23x1=r_float_optim25float_filtre15_inc_16(1:8:24,4)
24y1=r_float_optim25float_filtre15_inc_16(5:8:24,4)
25x1=x1+1
26y1=y1+1
27plot(x1(1),y1(1),'+')
28hold on;
29plot(x1(3),y1(3),'>')
30hold on;
31plot(x1(1:3),y1(1:3),'+')
32plot(x1(1:3),y1(1:3),'g')
33hold on;
34ls;
35r_float_optim25float_filtre15_inc_16(:,:)
36x2=r_float_optim25float_filtre15_inc_16(2:8:24,4)
37y2=r_float_optim25float_filtre15_inc_16(6:8:24,4)
38x2=x2+1
39y2=y2+1
40plot(x2(1),y2(1),'+')
41hold on;
42plot(x2(3),y2(3),'>')
43hold on;
44plot(x2(1:3),y2(1:3),'+')
45plot(x2(1:3),y2(1:3),'g')
46hold on;
47ls;
48r_float_optim25float_filtre15_inc_16(:,:)
49x3=r_float_optim25float_filtre15_inc_16(3:8:24,4)
50y3=r_float_optim25float_filtre15_inc_16(7:8:24,4)
51x3=x3+1
52y3=y3+1
53plot(x3(1),y3(1),'+')
54hold on;
55plot(x3(3),y3(3),'>')
56hold on;
57plot(x3(1:3),y3(1:3),'+')
58plot(x3(1:3),y3(1:3),'g')
59hold on;
60ls;
61r_float_optim25float_filtre15_inc_16(:,:)
62x4=r_float_optim25float_filtre15_inc_16(4:8:24,4)
63y4=r_float_optim25float_filtre15_inc_16(8:8:24,4)
64x4=x4+1
65y4=y4+1
66plot(x4(1),y4(1),'+')
67hold on;
68plot(x4(3),y4(3),'>')
69hold on;
70plot(x4(1:3),y4(1:3),'+')
71plot(x4(1:3),y4(1:3),'g')
72hold on;
73ls;
74
75
76load uv_total.dat
77   
78X=uv_total(1:5046,end-3);
79Y=uv_total(1:5046,end-2);
80U1=uv_total(1:5046,end-1);
81V1=uv_total(1:5046,end);
82U2=uv_total(5047:10092,end-1);
83V2=uv_total(5047:10092,end);
84U3=uv_total(10093:15138,end-1);
85V3=uv_total(10093:15138,end);
86   
87U=(U1+U2+U3)/3;
88V=(V1+V2+V3)/3;
89
90quiver(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
Note: See TracBrowser for help on using the repository browser.