###前提・実現したいこと
matlabでの球面状トランスデューサが作る音場についてのシミュレーションについての課題が出されましたがうまくいきません
(x,y)平面に音圧Pの音場を下記の3次元のグラフのようにしなければならないのですが私の作ったものではトゲトゲのものになってしまいます。どこを直せばいいか教えてください。
式: P=(ρAee^jωt/2π)Σ(e^-jkri/ri)dSi
それぞれの値はコードに記入しているものと同様です
ここでΣ以降のri,dSiはΣで和を取ります
###発生している問題・エラーメッセージ
エラーメッセージ
###該当のソースコード
clear all;
close all;
h=0.2;
rho=1000;
Ae=1;
t=1;
v=1500;
omega=2pi1000000;
k=omega/v;
r=1;
x_start=-3;
x_end=3;
num_div=50;
delta_x=(x_end-x_start)/num_div;
x=x_start:delta_x:x_end;
y_start=-3;
y_end=3;
delta_y=(y_end-y_start)/num_div;
y=y_start:delta_y:y_end;
[X,Y]=meshgrid(x,y);
P1=0;
for i=0:pi/120:pi/6
for j=0:pi/10:2*pi
a=i;
b=j;
x1=rcos(b)sin(a);
y1=rsin(b)sin(a);
z1=rcos(a);
deltaS=r^2sin(a)(pi/100)(pi/10);
R=sqrt((X-x1).^2+(Y-y1).^2+(h+z1).^2);
P=(rhoAe.exp(1iomegat).deltaS)./(2pi.R.exp(1ikR));
P1=P1+P;
end
end
P1=abs(P1);
mesh(X,Y,P1)
grid on
###試したこと
課題に対してアプローチしたことを記載してください
###補足情報(言語/FW/ツール等のバージョンなど)
より詳細な情報