質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

89.69%

matlabでの音場グラフの作成について(境界の条件について)

受付中

回答 0

投稿

  • 評価
  • クリップ 0
  • VIEW 985

sasataka

score 9

前提・実現したいこと

球面上のトランスデューサを4つのセクタに分割し異なる初期位相(具体的には0,,π/2,π,2π/3)を与えて下記の様なグラフを作りたいのですが、そこで境界部分をどのように設定していいのかわからず角度のずれたグラフになってしまいます。きれいに90度十字の様なグラフにするにはどのように直せばいいのでしょうか?
上が作りたいグラフ、下が私が作ったグラフです。また実際にはmatlabで角度を調整しています。

![イメージ説明](f2261591c0ca1b19ac069f863084276d.png)

イメージ説明

該当のソースコード

clear all;
close all;

h=0;
rho=1000;
Ae=0.001;
t=1;
v=1500;
omega=2*pi*1000000;
k=omega/v;
r=0.1;

x_start=-0.003;
x_end=0.003;
num_div=51;
delta_x=(x_end-x_start)/num_div;
x=x_start:delta_x:x_end;
y_start=-0.003;
y_end=0.003;
delta_y=(y_end-y_start)/num_div;
y=y_start:delta_y:y_end;
[X,Y]=meshgrid(x,y);

P1=0;
P2=0;
P3=0;
P4=0;
SdeltaS=0;
for i=0:pi/60:pi/3
for j=0:pi/10:(pi/2)-(pi/10)

a=i;
b=j;

x1=r*cos(b)*sin(a);
y1=r*sin(b)*sin(a);
z1=r*cos(a);
deltaS=r^2*sin(a)*(pi/60)*(pi/10);
R=sqrt((X-x1).^2+(Y-y1).^2+(h+z1).^2);
theta=0;
P=(rho*Ae.*exp(1i*omega*t).*deltaS.*exp(1i*theta))./(2*pi.*R.*exp(1i*k*R));

P1=P1+P;

end
end
for i=0:pi/60:pi/3
for j=pi/2:pi/10:pi-(pi/10)

a=i;
b=j;

x1=r*cos(b)*sin(a);
y1=r*sin(b)*sin(a);
z1=r*cos(a);
deltaS=r^2*sin(a)*(pi/60)*(pi/10);
R=sqrt((X-x1).^2+(Y-y1).^2+(h+z1).^2);
theta=pi/2;
P=(rho*Ae.*exp(1i*omega*t).*deltaS.*exp(1i*theta))./(2*pi.*R.*exp(1i*k*R));

P2=P2+P;

end
end
for i=0:pi/60:pi/3
for j=pi:pi/10:(3*pi/2)-(pi/10)

a=i;
b=j;

x1=r*cos(b)*sin(a);
y1=r*sin(b)*sin(a);
z1=r*cos(a);
deltaS=r^2*sin(a)*(pi/60)*(pi/10);
R=sqrt((X-x1).^2+(Y-y1).^2+(h+z1).^2);
theta=pi;
P=(rho*Ae.*exp(1i*omega*t).*deltaS.*exp(1i*theta))./(2*pi.*R.*exp(1i*k*R));

P3=P3+P;

end
end
for i=0:pi/60:pi/3
for j=3*pi/2:pi/10:(2*pi)-(pi/10)

a=i;
b=j;

x1=r*cos(b)*sin(a);
y1=r*sin(b)*sin(a);
z1=r*cos(a);
deltaS=r^2*sin(a)*(pi/60)*(pi/10);
R=sqrt((X-x1).^2+(Y-y1).^2+(h+z1).^2);
theta=3*pi/2;
P=(rho*Ae.*exp(1i*omega*t).*deltaS.*exp(1i*theta))./(2*pi.*R.*exp(1i*k*R));

P4=P4+P;

end
end
P=P1+P2+P3+P4;
P=abs(P);
mesh(X,Y,P)

grid on

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

まだ回答がついていません

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 89.69%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる