I have the following code but i am not able to get the graph to join the points together:
a few of the files are attached.
for i = 1:62
T=(i-1)*30;
FTA=load(['XYpoints_timeFromDeparture_',num2str(T),'.txt']);
FTA = fopen(['XYpoints_timeFromDeparture_',num2str(T),'.txt']);
Xa=fscanf(FTA,'%f %*f %*f %*f',Inf);
fclose(FTA);
FTA = fopen(['XYpoints_timeFromDeparture_',num2str(T),'.txt']);
Ya=fscanf(FTA,'%*f %f %*f %*f',Inf);
fclose(FTA);
FTA = fopen(['XYpoints_timeFromDeparture_',num2str(T),'.txt']);
Xb=fscanf(FTA,'%*f %*f %f %*f',Inf);
fclose(FTA);
FTA = fopen(['XYpoints_timeFromDeparture_',num2str(T),'.txt']);
Yb=fscanf(FTA,'%*f %*f %*f %f',Inf);
fclose(FTA);
Op=([Xa]-[Xb]);
Ad=([Ya]-[Yb]);
Angle = atand(Op./Ad);
A = (i*30);
TTxt =sprintf('For Time At %.0f: ', A);
disp(TTxt)
AvgAngle = mean(Angle)*-1;
AvgAngleDISP = sprintf("The Average Angle is: %.3f°",AvgAngle);
disp(AvgAngleDISP)
V_MPH = ((tand(AvgAngle)*10*cosd(30)-10*sind(30))*2.237);
VT = sprintf("The Velocity is: %.3fmph",V_MPH);
disp(VT)
hold on
figure(3);
x = A/60;
y = V_MPH;
plot(x,y,'k-x')
xlabel('Time (Minutes)')
ylabel('Velocity (MPH)')
end

6 件のコメント

Ameer Hamza
Ameer Hamza 2020 年 4 月 19 日
We cannot run your code because we don't have the txt files. Also, you haven't attached any image. How can we suggest anything useful?
SkyRider44
SkyRider44 2020 年 4 月 19 日
i have now added an image and just a few of the files
Star Strider
Star Strider 2020 年 4 月 19 日
Please describe what you want to do with the data. (Your code is, at best, confusing.)
Ameer Hamza
Ameer Hamza 2020 年 4 月 19 日
SkyRider44, the txt files you shared are not the ones required by this portion of your code.
Can you add a breakpoint at this line
plot(x,y,'k-x')
and paste the output of
size(x)
size(y)
SkyRider44
SkyRider44 2020 年 4 月 19 日
Each file contains 200+ rows and 4 columns of coordinates. 2 of the coloumns represent a raindrop at point x1,y1 and the other two columns represent a point x2,y2, i used the coordinates to find the angle of the line a in repspect to a verticaal line, from this i calulate the velocity and put this onto a graph on the y-axis and the time in minutes on the x-axis, however the points on the graph do not connect with a line.
SkyRider44
SkyRider44 2020 年 4 月 19 日
Did you want the output like so? (attached image)

サインインしてコメントする。

 採用された回答

Ameer Hamza
Ameer Hamza 2020 年 4 月 19 日
編集済み: Ameer Hamza 2020 年 4 月 19 日

0 投票

As you can see, x and y are 1x1, which indicates that for each file, you are just plotting a single point. If you want to draw a line, then you will need to move the plot() statement out of the for-loop and store the values of x and y in an array. Try the following code.
for i = 1:62
T=(i-1)*30;
FTA=load(['XYpoints_timeFromDeparture_',num2str(T),'.txt']);
FTA = fopen(['XYpoints_timeFromDeparture_',num2str(T),'.txt']);
Xa=fscanf(FTA,'%f %*f %*f %*f',Inf);
fclose(FTA);
FTA = fopen(['XYpoints_timeFromDeparture_',num2str(T),'.txt']);
Ya=fscanf(FTA,'%*f %f %*f %*f',Inf);
fclose(FTA);
FTA = fopen(['XYpoints_timeFromDeparture_',num2str(T),'.txt']);
Xb=fscanf(FTA,'%*f %*f %f %*f',Inf);
fclose(FTA);
FTA = fopen(['XYpoints_timeFromDeparture_',num2str(T),'.txt']);
Yb=fscanf(FTA,'%*f %*f %*f %f',Inf);
fclose(FTA);
Op=([Xa]-[Xb]);
Ad=([Ya]-[Yb]);
Angle = atand(Op./Ad);
A = (i*30);
TTxt =sprintf('For Time At %.0f: ', A);
disp(TTxt)
AvgAngle = mean(Angle)*-1;
AvgAngleDISP = sprintf("The Average Angle is: %.3f°",AvgAngle);
disp(AvgAngleDISP)
V_MPH = ((tand(AvgAngle)*10*cosd(30)-10*sind(30))*2.237);
VT = sprintf("The Velocity is: %.3fmph",V_MPH);
disp(VT)
x(i) = A/60;
y(i) = V_MPH;
end
figure;
plot(x,y,'k-x')
xlabel('Time (Minutes)')
ylabel('Velocity (MPH)')

2 件のコメント

SkyRider44
SkyRider44 2020 年 4 月 19 日
Thank you very helpful!
Ameer Hamza
Ameer Hamza 2020 年 4 月 19 日
Glad to be of help.

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

ヘルプ センター および File Exchange2-D and 3-D Plots についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by