MATLAB：1つのループのサブプロット

Question

軌道速度Vs Depthをプロットしたいと思います。これは次のスクリプトで成功します：

clear all; clc; close all; 

u = [ 
    0.1841 0.6168 0.6889 0.3531 -0.2019 -0.6549 -0.7366 -0.3882 0.1988; 
    0.1627 0.5404 0.6010 0.3045 -0.1819 -0.5763 -0.6439 -0.3353 0.1799; 
    0.1436 0.4732 0.5240 0.2623 -0.1639 -0.5071 -0.5628 -0.2895 0.1629; 
    0.1266 0.4141 0.4566 0.2258 -0.1475 -0.4463 -0.4919 -0.2498 0.1475; 
    0.1113 0.3622 0.3977 0.1941 -0.1328 -0.3927 -0.4299 -0.2154 0.1336; 
    0.0977 0.3165 0.3461 0.1666 -0.1196 -0.3455 -0.3756 -0.1855 0.1211]; 

w = [ 
    0.3543 0.0801 -0.2553 -0.4510 -0.3879 -0.1006 0.2431 0.4419 0.3824; 
    0.3055 0.0660 -0.2262 -0.3956 -0.3391 -0.0874 0.2129 0.3866 0.3351; 
    0.2630 0.0538 -0.2005 -0.3472 -0.2966 -0.0760 0.1866 0.3385 0.2939; 
    0.2262 0.0434 -0.1780 -0.3050 -0.2596 -0.0661 0.1635 0.2965 0.2579; 
    0.1942 0.0345 -0.1583 -0.2680 -0.2273 -0.0576 0.1434 0.2598 0.2264; 
    0.1666 0.0269 -0.1409 -0.2358 -0.1992 -0.0502 0.1258 0.2278 0.1990]; 


xx = 0; yy = 0; 
scale = 0; p = 1; 
dz = 2; 
figure(1) 
for d = 1:6; 
%d 
    dz = dz - 2; 
    for i =1:9; 
    quiver(xx, yy+dz, p*u(d,i),p*w(d,i),scale) 
    hold on 
    end 
end  
grid on  
axis equal;  
axis([-1 1 -10.5 0.5])  
%subplot(1,2,1)  
%axis([-1 1 -4.5 0.5])  
%subplot(1,2,2)  
%axis([-1 1 -10.5 -5.5]) 

私の質問：このプロットを2つのサブプロットに分割するには？

subplot(1,2,1)  
axis([-1 1 -4.5 0.5])  
subplot(1,2,2)  
axis([-1 1 -10.5 -5.5]) 

Answer 1

あなたは、単に1つのサブプロットおよびその他のサブプロット内の値の他の半分の値の半分をプロットすることにより、これを行うことができます。間隔がd=1:6であるため、ifの条件を使用して、1つのサブプロット内にd=1:3の値をプロットし、2番目のサブプロットにd=4:6のプロット値をプロットすることができます。すなわち

if d<=3 subplot(1,2,1); axis([-1 1 -4.5 0.5]); grid on; end 
if d>3 subplot(1,2,2); axis([-1 1 -10.5 -5.5]); grid on; end 

ので、全体的にあなたのコードは次のようになります。

clear all; clc; close all; 

u = [ 
    0.1841 0.6168 0.6889 0.3531 -0.2019 -0.6549 -0.7366 -0.3882 0.1988; 
    0.1627 0.5404 0.6010 0.3045 -0.1819 -0.5763 -0.6439 -0.3353 0.1799; 
    0.1436 0.4732 0.5240 0.2623 -0.1639 -0.5071 -0.5628 -0.2895 0.1629; 
    0.1266 0.4141 0.4566 0.2258 -0.1475 -0.4463 -0.4919 -0.2498 0.1475; 
    0.1113 0.3622 0.3977 0.1941 -0.1328 -0.3927 -0.4299 -0.2154 0.1336; 
    0.0977 0.3165 0.3461 0.1666 -0.1196 -0.3455 -0.3756 -0.1855 0.1211]; 

w = [ 
    0.3543 0.0801 -0.2553 -0.4510 -0.3879 -0.1006 0.2431 0.4419 0.3824; 
    0.3055 0.0660 -0.2262 -0.3956 -0.3391 -0.0874 0.2129 0.3866 0.3351; 
    0.2630 0.0538 -0.2005 -0.3472 -0.2966 -0.0760 0.1866 0.3385 0.2939; 
    0.2262 0.0434 -0.1780 -0.3050 -0.2596 -0.0661 0.1635 0.2965 0.2579; 
    0.1942 0.0345 -0.1583 -0.2680 -0.2273 -0.0576 0.1434 0.2598 0.2264; 
    0.1666 0.0269 -0.1409 -0.2358 -0.1992 -0.0502 0.1258 0.2278 0.1990]; 


xx = 0; yy = 0; 
scale = 0; p = 1; 
dz = 2; 
figure(1) 

for d = 1:6; 
    dz = dz - 2; 
    for i =1:9; 
    if d<=3 subplot(1,2,1); axis([-1 1 -4.5 0.5]); grid on; end 
    if d>3 subplot(1,2,2); axis([-1 1 -10.5 -5.5]); grid on; end 
    quiver(xx, yy+dz, p*u(d,i),p*w(d,i),scale) 
    hold on; 
    end 
end