这个公式怎样用MATLAB作图,三维,求程序

来源：百度知道 编辑：UC知道 时间：2024/06/10 20:01:25
d1=d2=7.5mm,N=2 ,M=2 ,k=2*PI/波长, 波长=7.5mm ,求程序,万分感谢!
感谢两位的帮助,第一个程序运行不对...第二个是极坐标系的,但是声场仿真如果能建立直角坐标系能更好地分析声场的特性,还请指导..

d1=7.5;
d2=d1;
N=2;
M=2;
lumda=7.5;
k=2*pi/lumda;

[alpha,seta]=meshgrid(0:0.01:pi);
D=sin(k*M*d1/2*cos(alpha).*sin(seta)).*sin(k*M*d2/2*sin(alpha).*sin(seta))/M./sin(k*d1/2*cos(alpha).*sin(seta))/N./sin(k*d2/2*sin(alpha).*sin(seta));
mesh(alpha,seta,D)
xlabel('\alpha');
ylabel('\theta');
zlabel('D')

clear;
d1=7.5;
d2=7.5;
n=2;
m=2;
b=7.5; % 波长
k=2*pi/b;
a=(0.1:0.1:0.9)*pi/2; % alpha 取值范围根据你的要求可自己改
c=a; % theta 取值范围根据你的要求可自己改
for i=1:9
for j=1:9
cs1=k*d1/2*cos(a(i))*sin(c(j));
cs2=k*d2/2*sin(a(i))*sin(c(j));
D(i,j)=sin(m*cs1)*sin(n*cs2)/(m*sin(cs1)*n*sin(cs2));
end
end
mesh(a,c,D)