plotting Intensity function which consists of an exponential sum

Question

Elinor Ginzburg 2020-3-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/512948-plotting-intensity-function-which-consists-of-an-exponential-sum

评论： Sriram Tadavarty 2020-3-26

在 MATLAB Online 中打开

hello,

Im trying to plot the followin function in matlab:

where k=10,d=1,A=1,0.1,0.01, epsilon is a random number between [0,1], the summation n is from 1 to 20 and theta goes from [-pi/2,pi/2]

I've written this code:

% param defenition %
L=100; N=20; d=1; k=10; Ar=1;
A=0.1;
theta=linspace(-pi/2,pi/2,200);
phase=(d+A*epsilon)*sin(theta);
epsilon=rand(0,1);
function psi = Psi(x)
    psi=1;
    for n=1:N
        psi=psi+(exp(i*k*phase).^n);
    end
end
I=([abs(psi)].^2);
fplot(I,[-pi/2,pi/2])

but when I try to run it, I get nothing printed in the figure. what did I do wrong?

Thank you very much for your time and attention.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Elinor Ginzburg 2020-3-26

在 MATLAB Online 中打开

thank you for your reply, but it still doesn't work. perhaps you have another idea? now it looks like this:

%% Question 03 %%
% param defenition %
L=100; N=20; d=1; k=10; Ar=1;
A=0.1;
theta=linspace(-pi/2,pi/2,200);
epsilon=rand(0,1);
phase=(d+A*epsilon)*sin(theta);
psi=1;
for n=1:N
    psi=psi+(exp(i*k*phase).^n);
end
I=((abs(psi)).^2);
plot(I)

Sriram Tadavarty 2020-3-26

rand(0,1) returns empty vector. So instead provide rand(1,1)

请先登录，再进行评论。

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

Answer 1

Birdman 2020-3-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/512948-plotting-intensity-function-which-consists-of-an-exponential-sum#answer_422066

在 MATLAB Online 中打开

Try this:

L=100; N=20; d=1; k=10; Ar=1;
A=0.1;
theta=linspace(-pi/2,pi/2,200);
epsilon=rand;
phase=(d+A*epsilon)*sin(theta);
psi=1;
for n=1:N
    psi=psi+(exp(1i*k*phase).^n);
end
I=((abs(psi)).^2);
plot(I)

2 个评论
显示无隐藏无

Elinor Ginzburg 2020-3-26

it worked! thank you so much! I had no idea I need to use rand instead of rand(0,1). how can you keep track between all the changes in different versions? thanks once again!