本文共 1345 字，大约阅读时间需要 4 分钟。
| 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 95 96 | %-------------- outer loop for x= 1 : 40 for y = 1 : 48 for z = 1 : 34 %----------inter loop x= 20 ; y= 30 ; z= 15 ; k = 1 ; for ztmp =- 1 : 1 z_ztmp = z+ztmp; if z_ztmp < 1 || z_ztmp > 34 continue ; end for ytmp =- 1 : 1 y_ytmp = y+ytmp; if y_ytmp < 1 || y_ytmp > 48 continue ; end for xtmp =- 1 : 1 if xtmp == 0 && ytmp == 0 && ztmp == 0 continue ; end x_xtmp = x+xtmp; if x_xtmp < 1 || x_xtmp > 48 continue ; end index = 40 * 48 *(z+ztmp- 1 )+ 40 *( y+ytmp- 1 )+(x+xtmp); if index < 1 || index > 27524 continue ; end data(k,:) = full_nt_matrix(index , :); k = k+ 1 ; end end end %% --------处理data l = size(data, 1 ); if l <= 1 continue ; end index = 40 * 48 *(z- 1 )+ 40 *(y- 1 )+x; time_series = full_nt_matrix(index,:); j = 1 ; for i= 1 :l r = corr2(time_series,data(i,:)); if r> 0.5 sim_data(j,:) = data(i,:); j = j+ 1 ; end end l = size(sim_data, 1 ); if l> 0 result_time_series = sum( sim_data, 1 )/ 9 ; end tmp_V(x,y,z,:) = result_time_series; %---------inter loop end end end |
本文转自二郎三郎博客园博客，原文链接：http://www.cnblogs.com/haore147/p/3795802.html，如需转载请自行联系原作者