-
Notifications
You must be signed in to change notification settings - Fork 0
/
creator_of_vectors.m
94 lines (77 loc) · 3.44 KB
/
creator_of_vectors.m
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
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
%%
%Normal vectors for fit_v4
matrix_I_low_ln=log(matrix_I_low/255);
matrix_I_high_ln=log(matrix_I_high/255);
N=13;
vector_al=[];
vector_al_I_low_ln=[];
vector_al_I_high_ln=[];
for n=1:N
vector_al=[vector_al matrix_al(n,:)];
vector_al_I_low_ln=[vector_al_I_low_ln matrix_I_low_ln(n,:)];
vector_al_I_high_ln=[vector_al_I_high_ln matrix_I_high_ln(n,:)];
end
plot(vector_al,vector_al_I_low_ln)
%vector_al=[vector_al(1:11) 1.92 1.92 1.92 1.92 1.92 1.92 1.92 1.92 1.92 1.92 vector_al(12:143)];
N=11;
vector_lu=[];
vector_lu_I_low_ln=[];
vector_lu_I_high_ln=[];
for n=1:N
vector_lu=[vector_lu matrix_lu(:,n)'];
vector_lu_I_low_ln=[vector_lu_I_low_ln matrix_I_low_ln(:,n)'];
vector_lu_I_high_ln=[vector_lu_I_high_ln matrix_I_high_ln(:,n)'];
end
%%
%Extra data for fit_v5
matrix_al_extra_I_high_ln=log(matrix_al_extra_I_high/255);
matrix_al_extra_I_low_ln=log(matrix_al_extra_I_low/255);
for n=1:10
vector_al=[vector_al(1:18+3*n) matrix_al_extra(n,:) vector_al(18+3*n+1:end)];
vector_al_I_low_ln=[vector_al_I_low_ln(1:18+3*n) matrix_al_extra_I_low_ln(n,:) vector_al_I_low_ln(18+3*n+1:end)];
vector_al_I_high_ln=[vector_al_I_high_ln(1:18+3*n) matrix_al_extra_I_high_ln(n,:) vector_al_I_high_ln(18+3*n+1:end)];
end
matrix_lu_extra_I_high_ln=log(matrix_lu_extra_I_high/255);
matrix_lu_extra_I_low_ln=log(matrix_lu_extra_I_low/255);
for n=1:8
vector_lu=[vector_lu(1:5+8*n) matrix_lu_extra(:,n)' vector_lu(5+8*n+1:end)];
vector_lu_I_low_ln=[vector_lu_I_low_ln(1:5+8*n) matrix_lu_extra_I_low_ln(:,n)' vector_lu_I_low_ln(5+8*n+1:end)];
vector_lu_I_high_ln=[vector_lu_I_high_ln(1:5+8*n) matrix_lu_extra_I_high_ln(:,n)' vector_lu_I_high_ln(5+8*n+1:end)];
end
plot(vector_al,vector_al_I_low_ln)
%%
%al_1
for n=1:4
pre_x=find(vector_al==1.92);
x=max(pre_x);
vector_al=[vector_al(1:x) matrix_al_extra_extra_1(n,:) vector_al(x+1:end)];
vector_al_I_high_ln=[vector_al_I_high_ln(1:x) log(matrix_I_high_extra_extra_1(n,:)/255) vector_al_I_high_ln(x+1:end)];
vector_al_I_low_ln=[vector_al_I_low_ln(1:x) log(matrix_I_low_extra_extra_1(n,:)/255) vector_al_I_low_ln(x+1:end)];
end
%lu_1
vector_auxiliar=[8 8 7.2 6.4 5.6 4.8 4 3.2 2.4 1.6 0.800000000000001 0 0];
for n=1:13
pre_x=find(vector_lu==aux);
x=min(pre_x);
vector_lu=[vector_lu(1:x) matrix_lu_extra_extra_1(:,n)' vector_lu(x+1:end)];
vector_lu_I_high_ln=[vector_lu_I_high_ln(1:x) log(matrix_I_high_extra_extra_1(:,n)'/255) vector_lu_I_high_ln(x+1:end)];
vector_lu_I_low_ln=[vector_lu_I_low_ln(1:x) log(matrix_I_low_extra_extra_1(:,n)'/255) vector_lu_I_low_ln(x+1:end)];
end
%al_2
vector_auxiliar=[1.76 1.6 1.44 1.28 1.12 0.96 0.8 0.64 0.48 0.32 0.16];
for n=1:11
pre_x=find(vector_al==vector_auxiliar(n));
x=min(pre_x);
vector_al=[vector_al(1:x) matrix_al_extra_extra_2(n,:) vector_al(x+1:end)];
vector_al_I_high_ln=[vector_al_I_high_ln(1:x) log(matrix_I_high_extra_extra_2(n,:)/255) vector_al_I_high_ln(x+1:end)];
vector_al_I_low_ln=[vector_al_I_low_ln(1:x) log(matrix_I_low_extra_extra_2(n,:)/255) vector_al_I_low_ln(x+1:end)];
end
%lu_2
for n=1:5
pre_x=find(vector_lu==8);
x=min(pre_x);
vector_lu=[vector_lu(1:x) matrix_lu_extra_extra_2(:,n)' vector_lu(x+1:end)];
vector_lu_I_high_ln=[vector_lu_I_high_ln(1:x) log(matrix_I_high_extra_extra_2(:,n)'/255) vector_lu_I_high_ln(x+1:end)];
vector_lu_I_low_ln=[vector_lu_I_low_ln(1:x) log(matrix_I_low_extra_extra_2(:,n)'/255) vector_lu_I_low_ln(x+1:end)];
end
plot(vector_al,vector_al_I_low_ln)