How to implement the math equation correctly to make a plot ?

Question

Vahram Voskerchyan 2021-6-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/862350-how-to-implement-the-math-equation-correctly-to-make-a-plot

回答： Walter Roberson 2021-6-22

采纳的回答： Walter Roberson

在 MATLAB Online 中打开

I have a equation for irradiance distribution that says.

it's an radiance distribution where the k is the wave vector , Λ is the period between the emitting elements. Implementing the equaiton can give the plot showing below. I am struggling to implement this equation when I plug everything in all I get is a single value. How do I loop this equation and create a matrix so it can give me the plot that is shown below ?

%rectangular aperture array factor. 
I_0 = 1;
lambda =1.55e-6;
 k=(2*pi)/lambda;
period = 2e-6;
theta = 20;
numel=33;
I = I_0*abs(((sind(k*period*sind(theta))/k*period*sind(theta))*(sind(numel/2*k*sind(theta))/sind(0.5*k*period*sind(theta)))))

For now I am trying to implement the equation to have the plot like below. But this will only give a signle value. How can I implement the right equation to have the plot shown above.

2 个评论
显示无隐藏无

Walter Roberson 2021-6-22

I suggest you read https://www.mathworks.com/matlabcentral/answers/395816-why-do-i-get-array-indices-must-be-positive-integers-or-logical-values-error-when-using-angle-fu#comment_882605

KSSV 2021-6-22

You need to show your code and thn specify your problem,

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2021-6-22

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/862350-how-to-implement-the-math-equation-correctly-to-make-a-plot#answer_730530

在 MATLAB Online 中打开

%rectangular aperture array factor.

N = 150;

I_0 = 1;

lambda =1.55e-6;

k=(2*pi)/lambda;

period = 2e-6;

theta = linspace(-25,25,N);

numel=33;

I = I_0 .* abs(((sind(k .* period .* sind(theta)) ./ k .* period .* sind(theta)) .* (sind(numel ./ 2 .* k .* sind(theta)) ./ sind(0.5 .* k .* period .* sind(theta)))));

plot(theta, I)

As you increase N, you will get a denser and denser plot that is more and more unreadable. That is because you are plotting a high frequency signal.