-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathexample2.m
More file actions
41 lines (31 loc) · 3.53 KB
/
example2.m
File metadata and controls
41 lines (31 loc) · 3.53 KB
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
true_list = '/data/flow/isabel_all/2d/all8.list';
fitted_list = '/data/flow/isabel_all/2d/fitted_quad_endpoint8/all8.list';
stderr_file = '/data/flow/isabel_all/2d/fitted_quad_endpoint8/sampling8_00_stderr.raw';
global data_w data_d data_t data_h
fitted_list = load_list(fitted_list);
data_t = 9; % sampling+1
t_ary = 0:0.25:data_t-1.25; t_ary = t_ary';
true_raw_ary=[11.000000 11.000000 2.000000 0.000000; 15.335054 11.597431 2.000000 1.000000; 16.386463 11.756291 2.000000 1.250000; 17.513479 11.897248 2.000000 1.500000; 18.706116 12.033085 2.000000 1.750000; 19.930893 12.184543 2.000000 2.000000; 21.177319 12.340253 2.000000 2.250000; 22.463076 12.474466 2.000000 2.500000; 23.798506 12.572787 2.000000 2.750000; 25.190195 12.622289 2.000000 3.000000; 26.637131 12.649789 2.000000 3.250000; 28.140429 12.709377 2.000000 3.500000; 29.687775 12.801541 2.000000 3.750000; 31.282902 12.867357 2.000000 4.000000; 32.925797 12.858459 2.000000 4.250000; 34.607513 12.818304 2.000000 4.500000; 36.332516 12.813733 2.000000 4.750000; 38.093102 12.784868 2.000000 5.000000; 39.884232 12.666047 2.000000 5.250000; 41.701248 12.485011 2.000000 5.500000; 43.544323 12.322654 2.000000 5.750000; 45.414104 12.194874 2.000000 6.000000; 47.288994 11.974685 2.000000 6.250000; 49.135414 11.659402 2.000000 6.500000; 50.962708 11.274793 2.000000 6.750000; 52.802040 10.835957 2.000000 7.000000; 54.626534 10.361021 2.000000 7.250000; 56.369793 9.801308 2.000000 7.500000; 58.059933 9.179996 2.000000 7.750000; 59.717098 8.511924 2.000000 8.000000];
uneven_t_ary = true_raw_ary(:,4);
true_x_ary = [interp1q(uneven_t_ary, true_raw_ary(:,1), t_ary) ...
interp1q(uneven_t_ary, true_raw_ary(:,2), t_ary)];
true_x_ary = true_x_ary';
% fit
[errstd1, errsum1, yy1, a, b, c] = quadfit_endpoint(t_ary,true_x_ary(1,:)',t_ary);
[errstd2, errsum2, yy2, a, b, c] = quadfit_endpoint(t_ary,true_x_ary(2,:)',t_ary);
z_ary = [yy1 yy2];
z_ary=z_ary';
z_stderr = zeros(size(z_ary));
z_stderr(1,:) = errstd1+errsum1;
z_stderr(2,:) = errstd2+errsum2;
% z_stderr = [errstd1;errstd2]*.5;
traced_fitted_x_raw_ary=[11.000000 11.000000 2.000000 0.000000; 15.784504 11.436985 2.000000 1.000000; 17.084600 11.512656 2.000000 1.250000; 18.414061 11.578578 2.000000 1.500000; 19.773125 11.641154 2.000000 1.750000; 21.158314 11.705994 2.000000 2.000000; 22.555178 11.764643 2.000000 2.250000; 23.968336 11.818457 2.000000 2.500000; 25.417698 11.870479 2.000000 2.750000; 26.905905 11.917877 2.000000 3.000000; 28.409637 11.972345 2.000000 3.250000; 29.909431 12.012804 2.000000 3.500000; 31.437849 12.008892 2.000000 3.750000; 33.005356 11.998375 2.000000 4.000000; 34.573181 11.995104 2.000000 4.250000; 36.156193 11.963609 2.000000 4.500000; 37.761539 11.896452 2.000000 4.750000; 39.372112 11.805254 2.000000 5.000000; 40.969849 11.706592 2.000000 5.250000; 42.560917 11.593554 2.000000 5.500000; 44.161194 11.454815 2.000000 5.750000; 45.768669 11.272562 2.000000 6.000000; 47.391045 11.029262 2.000000 6.250000; 49.026974 10.745399 2.000000 6.500000; 50.657955 10.442777 2.000000 6.750000; 52.266922 10.100547 2.000000 7.000000; 53.861134 9.671535 2.000000 7.250000; 55.444309 9.137972 2.000000 7.500000; 57.030910 8.526555 2.000000 7.750000; 58.645233 7.841821 2.000000 8.000000];
traced_fitted_x_ary = [interp1q(traced_fitted_x_raw_ary(:,4), traced_fitted_x_raw_ary(:,1), t_ary) ...
interp1q(traced_fitted_x_raw_ary(:,4), traced_fitted_x_raw_ary(:,2), t_ary)];
traced_fitted_x_ary=traced_fitted_x_ary';
while 1
trace_filtering
z_ary = post_x_ary;
z_stderr = post_x_stderr;
pause
end