source: altifloat/matlab_toolbox/drawing.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: 2.0 KB
Line 
1close all
2clear all
3
4load obs_float/nemed_observations_2.dat
5
6figure (1)
7    I1=find( nemed_observations_2(:,3)==1);
8    II1=nemed_observations_2(I1,:)
9    Iu1=find(II1(:,1)==1);
10    Iv1=find(II1(:,1)==2);
11    size1u=size(Iu1)
12    size1v=size(Iv1)   
13    plot(II1(Iu1,4),II1(Iv1,4),'--r');
14    hold on;
15
16    I2=find( nemed_observations_2(:,3)==2);
17    II2=nemed_observations_2(I2,:)
18    Iu2=find(II2(:,1)==1);
19    Iv2=find(II2(:,1)==2);
20    size2u=size(Iu2)
21    size2v=size(Iv2)   
22    plot(II2(Iu2,4),II2(Iv2,4),'r');
23    hold on;
24   
25    I3=find( nemed_observations_2(:,3)==3);
26    II3=nemed_observations_2(I3,:)
27    Iu3=find(II3(:,1)==1);
28    Iv3=find(II3(:,1)==2);
29    size3u=size(Iu3)
30    size3v=size(Iv3)   
31    plot(II3(Iu3,4),II3(Iv3,4),'b');
32    hold on;
33   
34    I4=find( nemed_observations_2(:,3)==4);
35    II4=nemed_observations_2(I4,:)
36    Iu4=find(II4(:,1)==1);
37    Iv4=find(II4(:,1)==2);
38    size4u=size(Iu4)
39    size4v=size(Iv4)   
40    plot(II4(Iu4,4),II4(Iv4,4),'m');
41    hold on;
42   
43    I5=find( nemed_observations_2(:,3)==5);
44    II5=nemed_observations_2(I5,:)
45    Iu5=find(II5(:,1)==1);
46    Iv5=find(II5(:,1)==2);
47    size5u=size(Iu5)
48    size5v=size(Iv5)   
49    plot(II5(Iu5,4),II5(Iv5,4),'g');
50    hold on;
51   
52    I6=find( nemed_observations_2(:,3)==6);
53    II6=nemed_observations_2(I6,:)
54    Iu6=find(II6(:,1)==1);
55    Iv6=find(II6(:,1)==2);
56    size6u=size(Iu6)
57    size6v=size(Iv6)   
58    plot(II6(Iu6,4),II6(Iv6,4),'c');
59    hold on;
60   
61    I7=find( nemed_observations_2(:,3)==7);
62    II7=nemed_observations_2(I7,:)
63    Iu7=find(II7(:,1)==1);
64    Iv7=find(II7(:,1)==2);
65    size7u=size(Iu7)
66    size7v=size(Iv7)   
67    plot(II7(Iu7,4),II7(Iv7,4),'k');
68    hold on;
69
70% % % % % %     I8=find( nemed_observations_2(:,3)==8);
71% % % % % %     II8=nemed_observations_2(I8,:)
72% % % % % %     Iu8=find(II8(:,1)==1);
73% % % % % %     Iv8=find(II8(:,1)==2);
74% % % % % %     size8u=size(Iu8)
75% % % % % %     size8v=size(Iv8)   
76% % % % % %     plot(II8(Iu8,4),II8(Iv8,4),'--b');
77% % % % % %     hold on;
Note: See TracBrowser for help on using the repository browser.