データが正方行列ではない場合の3Dグラフの作成について
显示 更早的评论
卒業研究で急遽matlabを使用することとなったため、matlabの知識に関してとても疎いです。
76×3行列の3Dグラフ(surfなど)の作成方法について教えていただきたいです。
1列ずつ3つのデータに分けて、それをそれぞれ正方行列に変換してからsurf関数を使うのかと考えているのですがどうも上手くいきません。
4 个评论
michio
2019-10-29
方向性は良い気がしますが、データやコードがないとコメントは難しいです。できれば、
- これが試したコードです
- ここでこんなエラーがでます
- 実行はできるがこの表示が意図と違う
などの追加情報を纏めてくれると助かります。偽データでもよいので、データまでつけてすぐに作業できると特に嬉しいです。。
rei shimizu
2019-10-29
Kazuya
2019-10-29
できればデータはファイル(csv, Excel, mat など) の添付で、コードとエラーメッセージはコメント蘭にコピペで投稿頂けると助かります。
実際にMATLABで試してみた方が速いのですが、画像から情報を吸い出すのは多少手間でして。
rei shimizu
2019-10-30
采纳的回答
michio
2019-10-29
matlab質問2.PNG のコードの中で、
x = A(:,1);
y = A(:,1);
z = A(:,1);
となっているところ、
x = A(:,1);
y = A(:,2);
z = A(:,3);
に変えるとよいように思います。
適当なダミーデータで実行してみると、
x = rand(100,1);
y = rand(100,1);
z = sin(2*pi*x)+cos(2*pi*y);
xlin = linspace(min(x),max(x),76);
ylin = linspace(min(y),max(y),76);
[X,Y] = meshgrid(xlin,ylin);
Z = griddata(x,y,z,X,Y,'cubic');
surf(X,Y,Z)
こんな感じです。griddata に 'cubic' と入れているのは、データの補間方法を指定しています。何も指定しなければ 'linear' ということで、線形補間を行います。
7 个评论
rei shimizu
2019-10-30
michio
2019-10-30
その為の MATLAB Answers ですので、もちろんです。もし周りもに困っている人がいたら利用されるように勧めてくださいね。質問される際にはこちらも参考にして頂けると幸いです。
rei shimizu
2019-10-30
ありがとうございます。
試しにやってみたところ3Dグラフの作成に成功しました。
x6 = A6(:,1);
y6 = A6(:,2);
z6 = A6(:,3);
xlin6 = linspace(min(x6),max(x6),76);
ylin6 = linspace(min(y6),max(y6),76);
[X6,Y6] = meshgrid(xlin6,ylin6);
Z6 = griddata(x6,y6,z6,X6,Y6,'cubic');
surf(X6,Y6,Z6)
しかし、griddataについてMathWorksさんのドキュメンテーションを読みましたが、理解が足りなかったためわからない点が1つできました。
以下の画像のようにZ6の中身にNaNの文字がいくつか出てしまうのですがこれは正常なのでしょうか？
michio
2019-10-30
もともとデータ点が存在しない範囲（＝外挿部分）は NaN を返す仕様になっています。
x6 = A6(:,1);
y6 = A6(:,2);
z6 = A6(:,3);
xlin6 = linspace(min(x6),max(x6),76);
ylin6 = linspace(min(y6),max(y6),76);
[X6,Y6] = meshgrid(xlin6,ylin6);
Z6 = griddata(x6,y6,z6,X6,Y6,'v4');
surf(X6,Y6,Z6)
と 'v4' を使うと外挿部分もデータを返すので、こちらを使っても良いかもしれません。
rei shimizu
2019-10-30
rei shimizu
2019-10-30
度々申し訳ございません。
以下の図のように、座標のデータを表示させてみると、元データには存在しないx座標、y座標、ｚ座標の数値の組み合わせが表示されました。
これはもともとこのような仕様なのでしょうか？
念のため元データも添付しておきます。
michio
2019-10-31
そうですね。。
xlin6 = linspace(min(x6),max(x6),76);
ylin6 = linspace(min(y6),max(y6),76);
[X6,Y6] = meshgrid(xlin6,ylin6);
の部分、xlin6, ylin6 それぞれ個別にみるとデータに存在するデータ点ですが、2次元で見るとデータに含まれない部分も出てくると思います。'v4' 以外の手法だと、その存在しない部分が NaN で返されているはずです。
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 2 次元および 3 次元プロット 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
