how can i plot in the while loop
显示 更早的评论
hey gays ive written a code @ it dosent have an error but it doesnt plot what i want what can i do with this ?
clear
t=16;
d=23;
M=5;
N1=floor(((275*(1/M))-30+d)-2);
s1=23.45*sind((360/365)*(284+(N1+2)));
l=35.44;
B=(N1-81)*360/364;
h=0.25*(t-12)*60;
z(1)=0;
while z<pi
z=0.1+z;
x=sin(l)*sin(s1)*cos(B)-cos(l)*sin(s1)*sin(B)*cos(z)+cos(l)*cos(s1)*cos(h)*cos(B)+sin(l)*cos(s1)*cos(h)*sin(B)*cos(z)+cos(s1)*sin(h)*sin(B)*sin(z)
theta=acos(x)
plot(theta,z)
end
6 个评论
Bob Thompson
2019-4-8
It looks like theta and x are both single values, so your results will just be a series of points. Matlab doesn't usually show points of lines in a very obvious manner. Are you sure it isn't plotting, and the problem is just not the visibility. You can test this by changing the format of the points of the plot.
Alex Mcaulley
2019-4-8
What do you want to plot? In this code you are plotting one point (the last value of theta, z)
Bob Thompson
2019-4-8
Also, as Alex mentioned, you are currently only plotting the last point of theta and x. Put 'hold on' before the while loop if you want to plot all points.
fatemeh jan
2019-4-8
fatemeh jan
2019-4-8
fatemeh jan
2019-4-8
采纳的回答
Alex Mcaulley
2019-4-8
Then, you simply need this:
t=16;
d=23;
M=5;
N1=floor(((275*(1/M))-30+d)-2);
s1=23.45*sind((360/365)*(284+(N1+2)));
l=35.44;
B=(N1-81)*360/364;
h=0.25*(t-12)*60;
z = 0:0.1:pi;
theta = acos(sin(l)*sin(s1)*cos(B)-cos(l)*sin(s1)*sin(B)*cos(z)+cos(l)*cos(s1)*cos(h)*cos(B)+sin(l)*cos(s1)*cos(h)*sin(B)*cos(z)+cos(s1)*sin(h)*sin(B)*sin(z));
plot(theta,z)
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 MATLAB 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
