no plot appearing for loop

3 views (last 30 days)
Seb apple
Seb apple on 13 Aug 2020
Edited: Stephen23 on 13 Aug 2020
when i run this code i have the plot screen appear but no data on it. can anbody help me please?
f = 1:0.1:1;
h = 0:0.1:4;
for i=1:length(f)
for j=1:length(h)
t=0;
L = 1;
m = 7+L;
u = 0.042;
Um = h(i) * (f(i));
r = Um * -cos(f(i) * t);
dy = Um*r(i)*sin(f(i)*t);
df = 1+L;
p = ((df * u)/m) * Um;
B = (df * (u^2) * f(i))/(2 * pi * m);
a = p/B;
g = 1 + ((pi * B)^(-0.5));
c = df*u^2;
Cd = (pi*B)^(-0.5);
n = r*abs(r)*u*Cd;
e=((c + n));
plot(f,e(i),'Linewidth',1.5);
hold on
xlabel('x');
ylabel('y');
end
end
  1 Comment
Stephen23
Stephen23 on 13 Aug 2020
Edited: Stephen23 on 13 Aug 2020
1- Your code only runs without error because this line happens to define a scalar value:
f = 1:0.1:1;
As soon as f has two or more elements (as that colon operator and the loop itself implies, otherwise they would be both completely superfluous) then your code will try to access the second element of many scalar variables and you will get a stream of errors (which is what happened when I tried it just now).
2- Probably most/all of the matrix operations should be replaced with array operations: https://www.mathworks.com/help/matlab/matlab_prog/array-vs-matrix-operations.html
3- I suspect that most/all of this code can be replaced with some simpler vectorized code: https://www.mathworks.com/help/matlab/matlab_prog/vectorization.html
4- Probably you should replace (...)^(-0.5) with 1./sqrt(...)
5- Plotting lots of separate points like that is not a good use of MATLAB: at best it can give you a set of disconnected points. But because the default linestyle does not have a point marker, by default you will not get anything displayed in the axes. You could fix this by changing the linestyle to include a marker, but you still won't get a line (which you seem to expect, due to the use of the 'Linewidth' option).
The MATLAB approach is to store the data in an array and then plot it all at once after the loop.
6- There might be other bugs...

Sign in to comment.

Answers (2)

Alan Stevens
Alan Stevens on 13 Aug 2020
You have t = 0 every time you enter the inner loop, and you don't update it. This gives rise to zeros in subsequent lines.

hosein Javan
hosein Javan on 13 Aug 2020
"f" was scalar rather than vector. some constants needed to be defined outside of the loop. and the iteration-depandant variables are fixed. I just don't know what is the inner loop for. you need to clarify what you want.
f = 0:0.1:1;
h = 0:0.1:4;
t=0;
L = 1;
m = 7+L;
u = 0.042;
df = 1+L;
c = df*u^2;
for i=1:length(f)
for j=1:length(h)
Um(i) = h(i) * (f(i));
r(i) = Um(i) * -cos(f(i) * t);
dy(i) = Um(i)*r(i)*sin(f(i)*t);
p(i) = ((df * u)/m) * Um(i);
B(i) = (df * (u^2) * f(i))/(2 * pi * m);
a(i) = p(i)/B(i);
g(i) = 1 + ((pi * B(i))^(-0.5));
Cd(i) = (pi*B(i))^(-0.5);
n(i) = r(i)*abs(r(i))*u*Cd(i);
e(i) = ((c + n(i)));
end
end
plot(f,e,'Linewidth',1.5);
xlabel('x');
ylabel('y');

Community Treasure Hunt

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

Start Hunting!