source: altifloat/matlab_toolbox/nbdepas7.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: 10.4 KB
Line 
1ls;
2load obs_float_test.dat;
3whos;
4obs_float_test (:,:)
5x=obs_float_test (1:2:16,4)
6y=obs_float_test (2:2:16,4)
7x=x+1
8y=y+1
9for p=1:16
10    figure (p)
11   plot(x(1),y(1),'+')
12   hold on;
13   plot(x(8),y(8),'>')
14   hold on;
15   plot(x(1:8),y(1:8),'+')
16   hold on;
17   plot(x(1:8),y(1:8),'-r')
18   hold on;
19end
20
21load r_float_optim25float-filtre15-inc_1.dat;
22whos;
23r_float_optim25float_filtre15_inc_1(:,:)
24x1=r_float_optim25float_filtre15_inc_1(1:2:16,4)
25y1=r_float_optim25float_filtre15_inc_1(2:2:16,4)
26x1=x1+1
27y1=y1+1
28figure(1)
29plot(x1(1),y1(1),'+')
30hold on;
31plot(x1(8),y1(8),'>')
32hold on;
33plot(x1(1:8),y1(1:8),'+')
34plot(x1(1:8),y1(1:8),'g')
35hold on;
36ls;
37
38load r_float_optim25float-filtre15-inc_2.dat;
39whos;
40r_float_optim25float_filtre15_inc_2(:,:)
41x2=r_float_optim25float_filtre15_inc_2(1:2:16,4)
42y2=r_float_optim25float_filtre15_inc_2(2:2:16,4)
43x2=x2+1
44y2=y2+1
45figure(2)
46plot(x2(1),y2(1),'+')
47hold on;
48plot(x2(8),y2(8),'>')
49hold on;
50plot(x2(1:8),y2(1:8),'+')
51plot(x2(1:8),y2(1:8),'g')
52hold on;
53ls;
54
55load r_float_optim25float-filtre15-inc_3.dat;
56whos;
57r_float_optim25float_filtre15_inc_3(:,:)
58x3=r_float_optim25float_filtre15_inc_3(1:2:16,4)
59y3=r_float_optim25float_filtre15_inc_3(2:2:16,4)
60x3=x3+1
61y3=y3+1
62figure(3)
63plot(x3(1),y3(1),'+')
64hold on;
65plot(x3(8),y3(8),'>')
66hold on;
67plot(x3(1:8),y3(1:8),'+')
68plot(x3(1:8),y3(1:8),'g')
69hold on;
70ls;
71
72load r_float_optim25float-filtre15-inc_4.dat;
73whos;
74r_float_optim25float_filtre15_inc_4(:,:)
75x4=r_float_optim25float_filtre15_inc_4(1:2:16,4)
76y4=r_float_optim25float_filtre15_inc_4(2:2:16,4)
77x4=x4+1
78y4=y4+1
79figure(4)
80plot(x4(1),y4(1),'+')
81hold on;
82plot(x4(8),y4(8),'>')
83hold on;
84plot(x4(1:8),y4(1:8),'+')
85plot(x4(1:8),y4(1:8),'g')
86hold on;
87ls;
88
89load r_float_optim25float-filtre15-inc_5.dat;
90whos;
91r_float_optim25float_filtre15_inc_5(:,:)
92x5=r_float_optim25float_filtre15_inc_5(1:2:16,4)
93y5=r_float_optim25float_filtre15_inc_5(2:2:16,4)
94x5=x5+1
95y5=y5+1
96figure(5)
97plot(x5(1),y5(1),'+')
98hold on;
99plot(x5(8),y5(8),'>')
100hold on;
101plot(x5(1:8),y5(1:8),'+')
102plot(x5(1:8),y5(1:8),'g')
103hold on;
104ls;
105
106load r_float_optim25float-filtre15-inc_6.dat;
107whos;
108r_float_optim25float_filtre15_inc_6(:,:)
109x6=r_float_optim25float_filtre15_inc_6(1:2:16,4)
110y6=r_float_optim25float_filtre15_inc_6(2:2:16,4)
111x6=x6+1
112y6=y6+1
113figure(6)
114plot(x6(1),y6(1),'+')
115hold on;
116plot(x6(8),y6(8),'>')
117hold on;
118plot(x6(1:8),y6(1:8),'+')
119plot(x6(1:8),y6(1:8),'g')
120hold on;
121ls;
122
123load r_float_optim25float-filtre15-inc_7.dat;
124whos;
125r_float_optim25float_filtre15_inc_7(:,:)
126x7=r_float_optim25float_filtre15_inc_7(1:2:16,4)
127y7=r_float_optim25float_filtre15_inc_7(2:2:16,4)
128x7=x7+1
129y7=y7+1
130figure(7)
131plot(x7(1),y7(1),'+')
132hold on;
133plot(x7(8),y7(8),'>')
134hold on;
135plot(x7(1:8),y7(1:8),'+')
136plot(x7(1:8),y7(1:8),'g')
137hold on;
138ls;
139
140load r_float_optim25float-filtre15-inc_8.dat;
141whos;
142r_float_optim25float_filtre15_inc_8(:,:)
143x8=r_float_optim25float_filtre15_inc_8(1:2:16,4)
144y8=r_float_optim25float_filtre15_inc_8(2:2:16,4)
145x8=x8+1
146y8=y8+1
147figure(8)
148plot(x8(1),y8(1),'+')
149hold on;
150plot(x8(8),y8(8),'>')
151hold on;
152plot(x8(1:8),y8(1:8),'+')
153plot(x8(1:8),y8(1:8),'g')
154hold on;
155ls;
156
157load r_float_optim25float-filtre15-inc_9.dat;
158whos;
159r_float_optim25float_filtre15_inc_9(:,:)
160x9=r_float_optim25float_filtre15_inc_9(1:2:16,4)
161y9=r_float_optim25float_filtre15_inc_9(2:2:16,4)
162x9=x9+1
163y9=y9+1
164figure(9)
165plot(x9(1),y9(1),'+')
166hold on;
167plot(x9(8),y9(8),'>')
168hold on;
169plot(x9(1:8),y9(1:8),'+')
170plot(x9(1:8),y9(1:8),'g')
171hold on;
172ls;
173
174load r_float_optim25float-filtre15-inc_10.dat;
175whos;
176r_float_optim25float_filtre15_inc_10(:,:)
177x10=r_float_optim25float_filtre15_inc_10(1:2:16,4)
178y10=r_float_optim25float_filtre15_inc_10(2:2:16,4)
179x10=x10+1
180y10=y10+1
181figure(10)
182plot(x10(1),y10(1),'+')
183hold on;
184plot(x10(8),y10(8),'>')
185hold on;
186plot(x10(1:8),y10(1:8),'+')
187plot(x10(1:8),y10(1:8),'g')
188hold on;
189ls;
190
191load r_float_optim25float-filtre15-inc_11.dat;
192whos;
193r_float_optim25float_filtre15_inc_11(:,:)
194x11=r_float_optim25float_filtre15_inc_11(1:2:16,4)
195y11=r_float_optim25float_filtre15_inc_11(2:2:16,4)
196x11=x11+1
197y11=y11+1
198figure(11)
199plot(x11(1),y11(1),'+')
200hold on;
201plot(x11(8),y11(8),'>')
202hold on;
203plot(x11(1:8),y11(1:8),'+')
204plot(x11(1:8),y11(1:8),'g')
205hold on;
206ls;
207
208load r_float_optim25float-filtre15-inc_12.dat;
209whos;
210r_float_optim25float_filtre15_inc_12(:,:)
211x12=r_float_optim25float_filtre15_inc_12(1:2:16,4)
212y12=r_float_optim25float_filtre15_inc_12(2:2:16,4)
213x12=x12+1
214y12=y12+1
215figure(12)
216plot(x12(1),y12(1),'+')
217hold on;
218plot(x12(8),y12(8),'>')
219hold on;
220plot(x12(1:8),y12(1:8),'+')
221plot(x12(1:8),y12(1:8),'g')
222hold on;
223ls;
224
225load r_float_optim25float-filtre15-inc_13.dat;
226whos;
227r_float_optim25float_filtre15_inc_13(:,:)
228x13=r_float_optim25float_filtre15_inc_13(1:2:16,4)
229y13=r_float_optim25float_filtre15_inc_13(2:2:16,4)
230x13=x13+1
231y13=y13+1
232figure(13)
233plot(x13(1),y13(1),'+')
234hold on;
235plot(x13(8),y13(8),'>')
236hold on;
237plot(x13(1:8),y13(1:8),'+')
238plot(x13(1:8),y13(1:8),'g')
239hold on;
240ls;
241
242load r_float_optim25float-filtre15-inc_14.dat;
243whos;
244r_float_optim25float_filtre15_inc_14(:,:)
245x14=r_float_optim25float_filtre15_inc_14(1:2:16,4)
246y14=r_float_optim25float_filtre15_inc_14(2:2:16,4)
247x14=x14+1
248y14=y14+1
249figure(14)
250plot(x14(1),y14(1),'+')
251hold on;
252plot(x14(8),y14(8),'>')
253hold on;
254plot(x14(1:8),y14(1:8),'+')
255plot(x14(1:8),y14(1:8),'g')
256hold on;
257ls;
258
259load r_float_optim25float-filtre15-inc_15.dat;
260whos;
261r_float_optim25float_filtre15_inc_15(:,:)
262x15=r_float_optim25float_filtre15_inc_15(1:2:16,4)
263y15=r_float_optim25float_filtre15_inc_15(2:2:16,4)
264x15=x15+1
265y15=y15+1
266figure(15)
267plot(x15(1),y15(1),'+')
268hold on;
269plot(x15(8),y15(8),'>')
270hold on;
271plot(x15(1:8),y15(1:8),'+')
272plot(x15(1:8),y15(1:8),'g')
273hold on;
274ls;
275
276load r_float_optim25float-filtre15-inc_16.dat;
277whos;
278r_float_optim25float_filtre15_inc_16(:,:)
279x16=r_float_optim25float_filtre15_inc_16(1:2:16,4)
280y16=r_float_optim25float_filtre15_inc_16(2:2:16,4)
281x16=x16+1
282y16=y16+1
283figure(16)
284plot(x16(1),y16(1),'+')
285hold on;
286plot(x16(8),y16(8),'>')
287hold on;
288plot(x16(1:8),y16(1:8),'+')
289plot(x16(1:8),y16(1:8),'g')
290hold on;
291ls;
292
293load uoptim25float-filtre15-inc1.dat;
294whos;
295X=uoptim25float_filtre15_inc1(:,end-2);
296Y=uoptim25float_filtre15_inc1(:,end-1);
297U=uoptim25float_filtre15_inc1(:,end);
298ls;
299load voptim25float-filtre15-inc1.dat;
300whos;
301V=voptim25float_filtre15_inc1(:,end);
302figure(1)
303quiver(X,Y,U,V);
304
305load uoptim25float-filtre15-inc2.dat;
306whos;
307X=uoptim25float_filtre15_inc2(:,end-2);
308Y=uoptim25float_filtre15_inc2(:,end-1);
309U=uoptim25float_filtre15_inc2(:,end);
310ls;
311load voptim25float-filtre15-inc2.dat;
312whos;
313V=voptim25float_filtre15_inc2(:,end);
314figure(2)
315quiver(X,Y,U,V);
316
317load uoptim25float-filtre15-inc3.dat;
318whos;
319X=uoptim25float_filtre15_inc3(:,end-2);
320Y=uoptim25float_filtre15_inc3(:,end-1);
321U=uoptim25float_filtre15_inc3(:,end);
322ls;
323load voptim25float-filtre15-inc3.dat;
324whos;
325V=voptim25float_filtre15_inc3(:,end);
326figure(3)
327quiver(X,Y,U,V);
328
329load uoptim25float-filtre15-inc4.dat;
330whos;
331X=uoptim25float_filtre15_inc4(:,end-2);
332Y=uoptim25float_filtre15_inc4(:,end-1);
333U=uoptim25float_filtre15_inc4(:,end);
334ls;
335load voptim25float-filtre15-inc4.dat;
336whos;
337V=voptim25float_filtre15_inc4(:,end);
338figure(4)
339quiver(X,Y,U,V);
340
341load uoptim25float-filtre15-inc5.dat;
342whos;
343X=uoptim25float_filtre15_inc5(:,end-2);
344Y=uoptim25float_filtre15_inc5(:,end-1);
345U=uoptim25float_filtre15_inc5(:,end);
346ls;
347load voptim25float-filtre15-inc5.dat;
348whos;
349V=voptim25float_filtre15_inc5(:,end);
350figure(5)
351quiver(X,Y,U,V);
352
353load uoptim25float-filtre15-inc6.dat;
354whos;
355X=uoptim25float_filtre15_inc6(:,end-2);
356Y=uoptim25float_filtre15_inc6(:,end-1);
357U=uoptim25float_filtre15_inc6(:,end);
358ls;
359load voptim25float-filtre15-inc6.dat;
360whos;
361V=voptim25float_filtre15_inc6(:,end);
362figure(6)
363quiver(X,Y,U,V);
364
365load uoptim25float-filtre15-inc7.dat;
366whos;
367X=uoptim25float_filtre15_inc7(:,end-2);
368Y=uoptim25float_filtre15_inc7(:,end-1);
369U=uoptim25float_filtre15_inc7(:,end);
370ls;
371load voptim25float-filtre15-inc7.dat;
372whos;
373V=voptim25float_filtre15_inc7(:,end);
374figure(7)
375quiver(X,Y,U,V);
376
377load uoptim25float-filtre15-inc8.dat;
378whos;
379X=uoptim25float_filtre15_inc8(:,end-2);
380Y=uoptim25float_filtre15_inc8(:,end-1);
381U=uoptim25float_filtre15_inc8(:,end);
382ls;
383load voptim25float-filtre15-inc8.dat;
384whos;
385V=voptim25float_filtre15_inc8(:,end);
386figure(8)
387quiver(X,Y,U,V);
388
389load uoptim25float-filtre15-inc9.dat;
390whos;
391X=uoptim25float_filtre15_inc9(:,end-2);
392Y=uoptim25float_filtre15_inc9(:,end-1);
393U=uoptim25float_filtre15_inc9(:,end);
394ls;
395load voptim25float-filtre15-inc9.dat;
396whos;
397V=voptim25float_filtre15_inc9(:,end);
398figure(9)
399quiver(X,Y,U,V);
400
401load uoptim25float-filtre15-inc10.dat;
402whos;
403X=uoptim25float_filtre15_inc10(:,end-2);
404Y=uoptim25float_filtre15_inc10(:,end-1);
405U=uoptim25float_filtre15_inc10(:,end);
406ls;
407load voptim25float-filtre15-inc10.dat;
408whos;
409V=voptim25float_filtre15_inc10(:,end);
410figure(10)
411quiver(X,Y,U,V);
412
413load uoptim25float-filtre15-inc11.dat;
414whos;
415X=uoptim25float_filtre15_inc11(:,end-2);
416Y=uoptim25float_filtre15_inc11(:,end-1);
417U=uoptim25float_filtre15_inc11(:,end);
418ls;
419load voptim25float-filtre15-inc11.dat;
420whos;
421V=voptim25float_filtre15_inc11(:,end);
422figure(11)
423quiver(X,Y,U,V);
424
425load uoptim25float-filtre15-inc12.dat;
426whos;
427X=uoptim25float_filtre15_inc12(:,end-2);
428Y=uoptim25float_filtre15_inc12(:,end-1);
429U=uoptim25float_filtre15_inc12(:,end);
430ls;
431load voptim25float-filtre15-inc12.dat;
432whos;
433V=voptim25float_filtre15_inc12(:,end);
434figure(12)
435quiver(X,Y,U,V);
436
437load uoptim25float-filtre15-inc13.dat;
438whos;
439X=uoptim25float_filtre15_inc13(:,end-2);
440Y=uoptim25float_filtre15_inc13(:,end-1);
441U=uoptim25float_filtre15_inc13(:,end);
442ls;
443load voptim25float-filtre15-inc13.dat;
444whos;
445V=voptim25float_filtre15_inc13(:,end);
446figure(13)
447quiver(X,Y,U,V);
448
449load uoptim25float-filtre15-inc14.dat;
450whos;
451X=uoptim25float_filtre15_inc14(:,end-2);
452Y=uoptim25float_filtre15_inc14(:,end-1);
453U=uoptim25float_filtre15_inc14(:,end);
454ls;
455load voptim25float-filtre15-inc14.dat;
456whos;
457V=voptim25float_filtre15_inc14(:,end);
458figure(14)
459quiver(X,Y,U,V);
460
461load uoptim25float-filtre15-inc15.dat;
462whos;
463X=uoptim25float_filtre15_inc15(:,end-2);
464Y=uoptim25float_filtre15_inc15(:,end-1);
465U=uoptim25float_filtre15_inc15(:,end);
466ls;
467load voptim25float-filtre15-inc15.dat;
468whos;
469V=voptim25float_filtre15_inc15(:,end);
470figure(15)
471quiver(X,Y,U,V);
472
473load uoptim25float-filtre15-inc16.dat;
474whos;
475X1=uoptim25float_filtre15_inc16(:,end-2);
476Y1=uoptim25float_filtre15_inc16(:,end-1);
477U2=uoptim25float_filtre15_inc16(:,end);
478ls;
479load voptim25float-filtre15-inc16.dat;
480whos;
481V2=voptim25float_filtre15_inc16(:,end);
482figure(16)
483quiver(X1,Y1,U2,V2,'-r');
484
485% for k=1:16
486% load uzero.dat;
487% whos;
488% X=uzero(:,end-2);
489% Y=uzero(:,end-1);
490% U=uzero(:,end);
491% ls;
492% load vzero.dat;
493% whos;
494% V=vzero(:,end);
495% figure(k);
496% quiver(X,Y,U,V);
497% end
Note: See TracBrowser for help on using the repository browser.