别介意,我设法找到一个解决方案。它不是很优雅,但至少工作原理:
v =兰德(1000 1)* 15;
φ=兰德(1000 1)* 360;
Data = [v,φ);
θ= 30;
theta_bounds = (theta-45,θ+ 45);
如果theta_bounds (1) < 0;%为θ略高于0
L_bound = rem (360 + theta_bounds (1), 360);
U_bound = rem (360 + theta_bounds (2), 360);
%数据可以基于一个或声明
idx =找到(数据(:,2)> = L_bound |数据(:,2)< = U_bound);
Data_cut =数据(idx:);
elseiftheta_bounds (2) > 360;%为θ略低于360
L_bound =快速眼动(theta_bounds (1), 360);
U_bound =快速眼动(theta_bounds (2), 360);
%数据可以基于一个或声明
idx =找到(数据(:,2)> = L_bound |数据(:,2)< = U_bound);
Data_cut =数据(idx:);
其他的
% 330θ之间
L_bound =快速眼动(theta_bounds (1), 360);
U_bound =快速眼动(theta_bounds (2), 360);
%可以减少数据基于和声明
idx =找到(数据(:,2)> = L_bound &数据(:,2)< = U_bound);
Data_cut =数据(idx:);
结束
图()
散射(数据(:,2),数据(:1),10日“填充”)
持有在
散射(Data_cut (:, 2), Data_cut(: 1), 10日“填充”)
xlim ([0360])