表面としてプロットすると周期的なデータが周期的に見えない

Question

波が多層球面に入射する問題をシミュレートしています。

forceZ1という名前の変数は、phaseとfrequencyの2つの変数の関数です。私は方程式を解いて、forceZ1の行列を見つけました。また、私はfrequencyが1 * 3000マトリックスであり、phaseが同じサイズであることを知っています。だから明らかにforceZ1は3000 * 3000になります。

この場合、位相は0と2*piの間で変化し、データに基づいて、力は周期2πで周期的であるはずです。これは、角度に関して任意の列ベクトルをforceZ1にプロットすることによって検証することができる。

figure(); plot(ang1, forceZ1(:,567)); 

います：

を私は表面forceZ1対frequencyとphaseを描くしようとしたが、図では、定期的な挙動を示していません。私は間違って何をしていますか？

ここに私のプロットのコードは次のとおりです。

Z=forceZ1; 
X=ang1; 
Y=frequency; 

figure 
surf(X,Y,Z,'edgecolor','none') 
shading interp 
view(3) 
xlim([0,2*pi]); 
xlabel('X=phase') 
ylabel('Y=nondimensional frequency') 
zlabel('Z=force in Z direction') 
title('Force-Frequency-Phase') 
colormap jet 

私はあなたの参考のために.mat file of my dataをアップロードしました。

Answer 1

あなたのforceZ1行列が途中で転置されたようです。データは確かに周期的ですが、周期性の方向はY軸です。

surf(X,Y,Z.','edgecolor','none'); 

をにつながる...：あなたはこのようなあなたのデータをプロットした場合

それは意味を作り始めるだろう、あなたのデータは、位相によって異なる高調波を含んでいるので、

しかし、大きさとオフセットを調整するには、ribbon：

など、高調波の形を強調する別のタイプのチャートを使用してみてください。

figure(); ribbon(Z(1:40:end,1:100:end)); 

...または示された周波数のためのスライダーを使用して、アニメーション、またはインタラクティブなチャートを作成してみてください。

また、データが転置されただけでなく反転された可能性もありますので、正しい結果を得るにはflipを使用する必要があります。