how to plot a vector

17 次查看（过去 30 天）
Elinor Ginzburg
Elinor Ginzburg 2020-5-24
回答： Elinor Ginzburg 2020-5-24
hello,
I'm trying to plot the points of a vector to see if it looks like I want it to, but I get a blank figure. can anyone please advice me on what's wrong with my code?
% units %
lambda = 1; L = 10; N = 500; a = 1; V0 = 5;
% Create lattice points %
dx = L/N;
x = zeros([1 N]);
for n = 1:N
x(n) = -L/2+(n-0.5)*dx;
end
% define potential %
Vho = zeros([1 N]);
for i = 1:N
Vsw = -V0*(heaviside(x(i)+a/2)-heaviside(x(i)-a/2));
end
figure
plot(Vsw)
Thank you very much for your help!
  1 个评论
Stephen23
Stephen23 2020-5-24
Given that heaviside belongs to the Symbolic Toolbox, you will probably need to consider it a bit more carefully.

请先登录，再进行评论。

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

采纳的回答

Stephen23
Stephen23 2020-5-24
编辑：Stephen23 2020-5-24
"what's wrong with my code?"
  • you preallocate a variable named Vho which you never use again.
  • inside the loop you allocate to a variable named Vso , but you do not use any indexing, so on each loop iteration you simply overwrite Vso, thus leaving only the (scalar value) from the final iteration.
So even though you write that you are plotting a vector, if you actually take a look at Vso it is actually scalar (and plotting a scalar gives the appearance of an empty plot, although it isn't actually, just the default linespec does not include any point marker).
I fixed those mistakes (and replaced heaviside with a simple logical comparison):
Vsw = zeros(1,N);
for k = 1:N
Vsw(k) = -V0*((x(k)+a/2)>0 - (x(k)-a/2)>0);
end
plot(Vsw,'-*')

更多回答（2 个）

Tommy
Tommy 2020-5-24
Vsw is a scalar, so nothing will show when you plot it. You are overwriting the value of Vsw in each loop of your for loop. Do you mean to instead store the values within Vho and plot that?
% define potential %
Vho = zeros([1 N]);
for i = 1:N
Vho(i) = -V0*(heaviside(x(i)+a/2)-heaviside(x(i)-a/2));
end
figure
plot(Vho)
Elinor Ginzburg
Elinor Ginzburg 2020-5-24
Thank you very much guys!
I mistyped and didn't index Vsw as you all mentioned.

类别

在 Help Center 和 File Exchange 中查找有关 Annotations 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by