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

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

新規登録して質問してみよう
ただいま回答率
85.50%
MATLAB

MATLABはMathWorksで開発された数値計算や数値の視覚化のための高水準の対話型プログラミング環境です。

Q&A

0回答

2151閲覧

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

sasataka

総合スコア11

MATLAB

MATLABはMathWorksで開発された数値計算や数値の視覚化のための高水準の対話型プログラミング環境です。

0グッド

0クリップ

投稿2017/09/26 07:47

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

イメージ説明

イメージ説明

###該当のソースコード
clear all;
close all;

h=0;
rho=1000;
Ae=0.001;
t=1;
v=1500;
omega=2pi1000000;
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=rcos(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=(rhoAe.exp(1iomegat).deltaS.exp(1itheta))./(2pi.R.exp(1ikR));

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=rcos(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=(rhoAe.exp(1iomegat).deltaS.exp(1itheta))./(2pi.R.exp(1ikR));

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=rcos(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=(rhoAe.exp(1iomegat).deltaS.exp(1itheta))./(2pi.R.exp(1ikR));

P3=P3+P;

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

a=i;
b=j;

x1=rcos(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=3pi/2;
P=(rho
Ae.exp(1iomega*t).deltaS.exp(1itheta))./(2pi.R.exp(1ikR));

P4=P4+P;

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

grid on

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

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

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

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

会員登録して回答してみよう

アカウントをお持ちの方は

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

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問