用matlab中的simulink画出这个公式的立体方向图！

simulink中有个f(U)模块，多个输入加上一个竖线输出上连着这个模块就可以编辑公式了 可以在里面放几个子模块简化公式。

clear
clc
close
thetaB=1; %初值设置
N=10;
lamda=10^(-3);
d=0.1;
theta=0:2*pi/100:2*pi; %变量范围设置
Ftheta=sin((N*pi)/lamda*d*(sin(theta)-sin(thetaB)))./sin((N*pi)/lamda*d*(sin(theta)-sin(thetaB)))...
*exp(1j*(N-1)*pi/lamda).*(sin(theta)-sin(thetaB));
figure(1)
plot(theta,abs(Ftheta)); %幅度与theta的关系图
figure(2)
plot(theta,real(Ftheta)); %实部与theta的关系图
figure(3)
plot(theta,imag(Ftheta)); %虚部与theta的关系图
希望对你有用咯。

这是个复函数，用real和imag分别取实部和虚部，看看matlab入门的书吧，p.s.这个公式是多缝衍射公式吗？