# Plot data from txt file

5 ビュー (過去 30 日間)
Angel Lozada 2023 年 6 月 9 日
コメント済み: Star Strider 2023 年 6 月 13 日
Hello.
I plotted a set of data (860000 samples) getting the followed plots:
Plot 1: Emitted Signal from NAA Station (Phase)
Plot 2: Emitted Signal from NAA Station (Amplitude)
The data were taking as 10 samples per second (see attached file named Figure 1), i.e., in second 0 the are 10 numbers, in second 1 another 10 numbers, and so on.
However, when I analyze the plots, is possible to see that there are more than 10 data (>>10 numbers) in 0s, the same for second 1 and so on.
Regards.
clear;
T = data(:,1);
T = seconds(data(:,1));
T.Format = 'hh:mm';
T2 = data(:,2);
figure;
plot (T,T2,'x');
title 'Emitted Signal from NAA Station';
ylabel ('Phase');
xlabel ('Samples'); T = data(:,1);
T = seconds(data(:,1));
T.Format = 'hh:mm';
T3 = data(:,3);
figure;
plot (T,T3,'x');
title 'Emitted Signal from NAA Station';
ylabel ('Apmlitude');
xlabel ('Samples'); ##### 2 件のコメント表示 1 件の古いコメント非表示 1 件の古いコメント
Angel Lozada 2023 年 6 月 10 日
Dear Star Strider.
I will try to explain in a better way my inquiry.
As you can see in previous attached .txt file there are 3 columns indicating a giving information for a signal.
Column 1: seconds and sample in each second; Column: 2 Phase and Column 3: Amplitude.
Let say we take a look to Phase data. In column 1 the collected data begin at 0s and during this time we have 10 samples (from 0.0 to 0.9), then at 1s we have other ten samples (1.0 to 1.9), and so on.
Attached file Figure 3 shows the expected plot (plotting 10 samples in each second) and also shows what I am getting in current Phase plot (a lot of samples between second 0 and 1). You can see in this Figure 3 that is a kind of discrete plot, because there is information (samples) only in each second and not in between.
Again, thanks for your support and patience.
Regards.

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

### 採用された回答

Star Strider 2023 年 6 月 10 日
If you simply want to aggregate the data in each second and plot them that way, try something like this —
Uzp = unzip('Data 2.zip');
A1 = 864000×3
0 -101.0300 23.5600 0.1000 -100.2100 23.8600 0.2000 -101.7000 23.5600 0.3000 -107.6600 23.4700 0.4000 -110.5400 22.7800 0.5000 -112.2900 21.6600 0.6000 -107.4600 21.6700 0.7000 -115.9300 21.2700 0.8000 -114.0500 22.2300 0.9000 -110.9900 22.2400
t = seconds(A1(:,1));
ixv = ceil(A1(:,1) + 1E-12); % Index Vector For 'accumarray'
secv = accumarray(ixv, (1:size(A1,1)).', [], @(x){A1(x,[2 3])}) % Aggregatte Points By Second
secv = 86400×1 cell array
{10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double} {10×2 double}
A2 = reshape(cell2mat(secv).', 2,10,[]); % Convert To Numeric MAtrix From Cell Array & Reshape ARray
tv = unique(ixv); % Corresponding Time Vector
figure
yyaxis left
plot(tv, squeeze(A2(1,:,tv)), '.')
ylabel('Phase (°)')
grid
yyaxis right
plot(tv, squeeze(A2(2,:,tv)), '.')
ylabel('Amplitude')
xlabel('Time')
% xlim([tv(1) tv(end)])
xlim([0 5.5]) % cs = A2(2,:,:) .* exp(1j*deg2rad(A2(1,:,:))); % Complex Signal
%
% figure
% plot(tv, squeeze(abs(cs)))
% grid
% xlabel('Time')
% ylabel('Magnitude')
% % xlim([0 5.5])
% xlim([tv(1) tv(end)])
This plots them as dots rather than asterisks. Choose whatever marker you want.
.
##### 2 件のコメント表示 1 件の古いコメント非表示 1 件の古いコメント
Star Strider 2023 年 6 月 13 日
As always, my pleasure!

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

### カテゴリ

Help Center および File Exchange2-D and 3-D Plots についてさらに検索

R2022a

### Community Treasure Hunt

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

Start Hunting!