How to extract data from txt file and plot spectogram?
19 次查看(过去 30 天)
显示 更早的评论
Hello,
I have a large sample of data (100000 samples) in txt file consisting of time and amplitude of the signal.
I want to plot spectogram. The code is not at all able to extract the .txt file at all along the path.
Kindly advice me as how to extract data from txt file when data is very large as it is in my case and hence to plot the spectogram.
The code is as below:
close all;
freq_sam = 100000; % samples in tt1 data
period_tt1 = 0.005; % time period
freq_tt1 = 1/period_tt1; % frequency
data = importdata('D:\Data_tokamak\data\obp.txt');
x_data = data(:,1);
y_data = data(:,2);
disp(x_data);
disp(y_data);
signal_tt1 = data;
fft_tt1 = fft(signal_tt1); % fourier transform
fft_tt1 = fftshift(fft_tt1); % performing fft shift
f = freq_sam/2*linspace(-1,1,freq_sam);
figure;
plot(f, fft_tt1);
采纳的回答
Star Strider
2023-6-3
I prefer the pspectrum function for these analyses. See specifically Spectrogram and Reassigned Spectrogram of Chirp for an example of returning the data and plotting it as well.
data = readmatrix('D:\Data_tokamak\data\obp.txt');
x_data = data(:,1);
y_data = data(:,2);
Ts = mean(diff(x_data))
Fs = 1/Ts;
figure
pspectrum(y_data, Fs, 'spectrogram') % Plot 'pspectrum' 'spectrogram'
[sp,fp,tp] = pspectrum(y_data,fs,"spectrogram"); % Return Values, Then Plot
figure
waterfall(fp,tp,sp')
set(gca,XDir="reverse",View=[60 60])
ylabel("Time (s)")
xlabel("Frequency (Hz)")
Be certain that the ‘x_data’ are sampled uniformly. The easiest way to do that is to calculate the standard deviation (std) of ‘diff(x_data)’ and compare it to ‘Ts’. It should be vanishingly small, and preferably less than 1E-6 of the ‘Ts’ value.
It would help to have ‘obp.txt’ to work with.
.
37 个评论
Rahul
2023-6-3
Hi
Star Strider,
The program gives an error
Error using pspectrum
Expected X to be finite
Is it due to too large sample of data or due to some other reasons?
rgds
Star Strider
2023-6-3
I doubt if it is too large, because 1E+5 samples is reasonable. Without having the data to work with, I suspect the problem is the presence of Inf or NaN values in ‘y_data’.
The fillmissing function searches only for NaN (or NaT — ‘Not a Time’) values, so the first thing to do is to search for those values in ‘y_data’ and first replace them with NaN values then use fillmissing to interpolate them —
data = [(0:10).' [rand(4,1); Inf; rand(3,1); NaN; rand(2,1)]]
data = 11×2
0 0.1625
1.0000 0.7347
2.0000 0.3144
3.0000 0.1415
4.0000 Inf
5.0000 0.6968
6.0000 0.6211
7.0000 0.9422
8.0000 NaN
9.0000 0.5554
data(~isfinite(data)) = NaN
data = 11×2
0 0.1625
1.0000 0.7347
2.0000 0.3144
3.0000 0.1415
4.0000 NaN
5.0000 0.6968
6.0000 0.6211
7.0000 0.9422
8.0000 NaN
9.0000 0.5554
data = fillmissing(data, 'linear')
data = 11×2
0 0.1625
1.0000 0.7347
2.0000 0.3144
3.0000 0.1415
4.0000 0.4191
5.0000 0.6968
6.0000 0.6211
7.0000 0.9422
8.0000 0.7488
9.0000 0.5554
I would not remove the missing value rows because spectral analysis functions (with the exception only of the nufft function that I mention here only to note its presence) require (and assume) regular sampling, and removing the rows with missing values violates that assumption. (The interpolation method can be whatever you want, I chose 'linear' here for convenience.)
The slightly modified version of ‘data’ should then work.
.
Rahul
2023-6-4
编辑:Rahul
2023-6-4
Hi Star Strider,
The data(~isfinite(data)) = NaN when implemented with if statement, I got it implemented with isnan(data(~isfinite(amplitude))), but it lead to converting all elements of amplitude into NaN.
The reason for using isnan is because inside if statement the syntax for isfinite gives error and couldn't be implemented.
plz advice.
if isnan(data(~isfinite(amplitude)))
amplitude = fillmissing(amplitude, 'linear');
end
%Plot the spectrogram
figure
pspectrum(amplitude, Fs, 'spectrogram') % Plot 'pspectrum' 'spectrogram'
[sp,fp,tp] = pspectrum(amplitude,Fs,"spectrogram");
Star Strider
2023-6-4
编辑:Star Strider
2023-6-4
The line:
data(~isfinite(data)) = NaN;
should not convert all the values to NaN. I demonstrated how it works earlier. In my code it does not need an if block to work correctly. It already detects all of them and then simply converts all the Inf values to NaN so that fillmissing can interpolate them.
Another option is:
data(isinf(data)) = NaN;
data = fillmissing(data, 'linear');
and then the rest of the code in my original Answer should work without problems.
Please share your data. Attach it using the ‘paperclip’ icon in the top toolbar. That is the easiest way to get this to work.
.
Rahul
2023-6-4
Hi Star Strider,
Its still not working. Gives error.
Error using pspectrum
Expected X to be finite
I have attached the data. plz find herewith the data as attached.
Star Strider
2023-6-4
Thank you!
Try this —
filecontents = fileread('obp3n.txt') % Check File Format & Contents
filecontents =
'DataBase = ***
ShotNo = ****
SignalName = OBP3N
SignalUnit = V
TriggerTime = 0.000000 ms
Period = 0.005000 ms
CreateTime = 2022/08/17 16:27:47
Samples = 100000
0.000000 -0.000408
0.005000 -0.001049
0.010000 -0.000088
0.015000 -0.001049
0.020000 0.003435
0.025000 -0.000408
0.030000 -0.000729
0.035000 -0.003932
0.040000 -0.000408
0.045000 0.002474
0.050000 -0.000088
0.054999 -0.001369
0.059999 -0.002651
0.064999 -0.001369
0.069999 0.000873
0.074999 0.001513
0.079999 0.000552
0.084999 -0.001049
0.089999 -0.001690
0.094999 0.002474
0.099999 0.004076
0.104999 0.000552
0.109999 -0.000408
0.114999 -0.001369
0.119999 0.000873
0.124999 0.000552
0.129999 -0.001369
0.134999 -0.000729
0.139999 0.001193
0.144999 0.001193
0.149999 0.000552
0.154998 0.001834
0.159998 -0.001369
0.164998 -0.001690
0.169998 -0.000408
0.174998 -0.000729
0.179998 0.001193
0.184998 0.000873
0.189998 -0.001369
0.194998 0.000873
0.199998 0.003435
0.204998 0.000232
0.209998 0.001193
0.214998 -0.000408
0.219998 0.000552
0.224998 0.002474
0.229998 0.001513
0.234998 -0.000408
0.239998 -0.000088
0.244998 0.000552
0.249998 0.000552
0.254997 0.001834
0.259997 0.001193
0.264997 -0.000729
0.269997 -0.000408
0.274997 0.000552
0.279997 0.001834
0.284997 -0.001690
0.289997 -0.001369
0.294997 -0.001690
0.299997 -0.001369
0.304997 -0.001369
0.309997 -0.001369
0.314997 0.000873
0.319997 -0.006174
0.324997 -0.001049
0.329997 -0.001049
0.334997 -0.001049
0.339997 -0.001369
0.344997 0.000552
0.349997 -0.003612
0.354996 0.001834
0.359996 0.002795
0.364996 -0.001369
0.369996 -0.001369
0.374996 0.000232
0.379996 -0.000088
0.384996 0.000552
0.389996 -0.002010
0.394996 -0.002010
0.399996 -0.000729
0.404996 0.000873
0.409996 0.002154
0.414996 -0.000729
0.419996 -0.000088
0.424996 -0.001049
0.429996 -0.000729
0.434996 -0.000729
0.439996 0.000552
0.444996 -0.001690
0.449996 -0.002971
0.454995 0.000232
0.459995 0.002154
0.464995 0.000552
0.469995 -0.000088
0.474995 -0.002651
0.479995 -0.001049
0.484995 0.001193
0.489995 0.000232
0.494995 -0.000729
0.499995 -0.000408
0.504995 0.000873
0.509995 0.000873
0.514995 0.002474
0.519995 0.002154
0.524995 0.002474
0.529995 0.000873
0.534995 0.000232
0.539995 0.002474
0.544995 0.000232
0.549995 0.000552
0.554994 -0.001049
0.559994 -0.000729
0.564994 0.001513
0.569994 -0.003612
0.574994 0.004717
0.579994 -0.000088
0.584994 -0.001369
0.589994 0.000873
0.594994 -0.000088
0.599994 0.001193
0.604994 -0.001049
0.609994 -0.000729
0.614994 0.001513
0.619994 -0.001369
0.624994 -0.001049
0.629994 -0.002651
0.634994 0.000873
0.639994 0.001513
0.644994 -0.000729
0.649994 0.000873
0.654993 0.001193
0.659993 -0.001690
0.664993 0.000232
0.669993 0.001834
0.674993 0.004396
0.679993 -0.001690
0.684993 0.000552
0.689993 -0.001369
0.694993 -0.000729
0.699993 0.000873
0.704993 -0.001369
0.709993 -0.000408
0.714993 0.000873
0.719993 0.000232
0.724993 0.000552
0.729993 -0.001369
0.734993 0.000552
0.739993 -0.000729
0.744993 0.001513
0.749993 -0.000408
0.754992 -0.000408
0.759992 0.000232
0.764992 -0.002010
0.769992 -0.001369
0.774992 0.001513
0.779992 -0.000408
0.784992 -0.001690
0.789992 0.000552
0.794992 -0.000408
0.799992 0.001513
0.804992 -0.000088
0.809992 -0.000729
0.814992 -0.003291
0.819992 -0.002971
0.824992 0.003756
0.829992 -0.002651
0.834992 -0.001369
0.839992 -0.000408
0.844992 -0.001049
0.849992 0.001834
0.854991 0.002154
0.859991 0.000873
0.864991 -0.001049
0.869991 0.000232
0.874991 0.001193
0.879991 -0.001049
0.884991 -0.001690
0.889991 0.000873
0.894991 -0.004252
0.899991 0.000552
0.904991 -0.001049
0.909991 0.000873
0.914991 0.000552
0.919991 -0.000729
0.924991 -0.001369
0.929991 0.000232
0.934991 0.001513
0.939991 0.000873
0.944991 -0.000729
0.949991 -0.000088
0.954990 0.000232
0.959990 0.000232
0.964990 -0.001690
0.969990 -0.001690
0.974990 0.001193
0.979990 0.001834
0.984990 -0.001690
0.989990 -0.000729
0.994990 0.000552
0.999990 -0.000408
1.004990 0.001513
1.009990 -0.000408
1.014990 0.001193
1.019990 -0.000729
1.024990 -0.000729
1.029990 0.000552
1.034990 0.002795
1.039990 0.000873
1.044990 -0.000408
1.049990 0.001193
1.054989 -0.000408
1.059989 0.003115
1.064989 -0.001369
1.069989 -0.001690
1.074989 0.000873
1.079989 0.001513
1.084989 0.000232
1.089989 -0.000729
1.094989 0.000552
1.099989 -0.000408
1.104989 0.000232
1.109989 0.000232
1.114989 -0.000088
1.119989 -0.000729
1.124989 -0.000088
1.129989 0.000232
1.134989 0.002795
1.139989 0.001193
1.144989 0.001193
1.149989 0.002154
1.154988 -0.000408
1.159988 -0.000729
1.164988 -0.000729
1.169988 -0.000408
1.174988 0.001193
1.179988 -0.001369
1.184988 0.000873
1.189988 0.000873
1.194988 0.000873
1.199988 -0.000408
1.204988 0.000232
1.209988 -0.000729
1.214988 0.001513
1.219988 0.000873
1.224988 -0.001690
1.229988 0.001193
1.234988 -0.000729
1.239988 -0.000408
1.244988 0.000873
1.249988 -0.000408
1.254987 0.000552
1.259987 -0.001690
1.264987 -0.000408
1.269987 -0.000729
1.274987 -0.000088
1.279987 -0.001369
1.284987 -0.001690
1.289987 0.000873
1.294987 0.001834
1.299987 -0.001049
1.304987 -0.000088
1.309987 -0.001049
1.314987 0.000232
1.319987 -0.001369
1.324987 0.001193
1.329987 -0.001369
1.334987 -0.000729
1.339987 -0.000088
1.344987 -0.000088
1.349987 -0.001049
1.354986 -0.000088
1.359986 -0.000088
1.364986 -0.000729
1.369986 0.000552
1.374986 -0.001049
1.379986 -0.001369
1.384986 -0.001049
1.389986 -0.000088
1.394986 0.000552
1.399986 -0.001690
1.404986 0.001834
1.409986 -0.002651
1.414986 -0.000408
1.419986 -0.000408
1.424986 0.000552
1.429986 -0.002010
1.434986 0.000873
1.439986 0.000552
1.444986 0.000873
1.449986 0.001513
1.454985 -0.001049
1.459985 0.001193
1.464985 0.000232
1.469985 0.001193
1.474985 -0.000088
1.479985 -0.000729
1.484985 0.001513
1.489985 0.001193
1.494985 0.000873
1.499985 0.001193
1.504985 0.000232
1.509985 0.000552
1.514985 -0.000088
1.519985 0.001193
1.524985 -0.001369
1.529985 -0.000729
1.534985 -0.001369
1.539985 0.000873
1.544985 0.000232
1.549985 0.002154
1.554984 0.001193
1.559984 0.000873
1.564984 -0.000408
1.569984 -0.000408
1.574984 0.000232
1.579984 0.000232
1.584984 0.001513
1.589984 0.000552
1.594984 0.003435
1.599984 0.001513
1.604984 0.000552
1.609984 -0.001369
1.614984 -0.001369
1.619984 -0.000408
1.624984 0.000232
1.629984 -0.001690
1.634984 0.001513
1.639984 0.000552
1.644984 -0.000729
1.649984 0.000552
1.654983 -0.000729
1.659983 -0.002010
1.664983 0.000232
1.669983 0.002154
1.674983 -0.000729
1.679983 -0.001049
1.684983 -0.001369
1.689983 0.000232
1.694983 -0.001369
1.699983 -0.000408
1.704983 -0.000408
1.709983 -0.000729
1.714983 -0.001049
1.719983 -0.001369
1.724983 0.001193
1.729983 0.001834
1.734983 -0.000729
1.739983 -0.001369
1.744983 -0.000088
1.749983 0.001193
1.754982 -0.001369
1.759982 -0.001369
1.764982 -0.001369
1.769982 0.001513
1.774982 -0.000088
1.779982 -0.000408
1.784982 -0.000729
1.789982 -0.000408
1.794982 -0.001049
1.799982 -0.000729
1.804982 0.001513
1.809982 -0.001369
1.814982 -0.000088
1.819982 0.001834
1.824982 0.000552
1.829982 0.000232
1.834982 -0.001369
1.839982 -0.000729
1.844982 -0.000729
1.849982 -0.001369
1.854981 -0.001369
1.859981 -0.001690
1.864981 0.000552
1.869981 -0.000088
1.874981 -0.001690
1.879981 -0.000729
1.884981 -0.000408
1.889981 -0.000088
1.894981 0.000873
1.899981 -0.001049
1.904981 0.001834
1.909981 -0.000408
1.914981 0.001193
1.919981 0.001513
1.924981 -0.001690
1.929981 -0.000088
1.934981 0.001513
1.939981 0.000232
1.944981 0.000873
1.949981 0.001834
1.954980 -0.000408
1.959980 -0.001049
1.964980 -0.000729
1.969980 0.000552
1.974980 0.000232
1.979980 -0.002330
1.984980 -0.000408
1.989980 -0.001690
1.994980 -0.001049
1.999980 -0.000088
2.004980 0.000552
2.009980 -0.001369
2.014980 0.002154
2.019980 -0.000408
2.024980 -0.000088
2.029980 0.000873
2.034980 -0.001690
2.039980 0.001193
2.044980 -0.001049
2.049980 0.000873
2.054979 -0.001369
2.059979 0.000552
2.064979 -0.000088
2.069979 -0.000088
2.074979 0.001193
2.079979 0.000552
2.084979 -0.001369
2.089979 -0.001369
2.094979 0.000552
2.099979 0.001193
2.104979 0.000873
2.109979 0.000873
2.114979 -0.002971
2.119979 -0.000088
2.124979 -0.000408
2.129979 0.000873
2.134979 -0.000408
2.139979 0.000232
2.144979 -0.000408
2.149979 0.000232
2.154978 -0.000088
2.159978 -0.001369
2.164978 0.000873
2.169978 -0.001049
2.174978 0.001513
2.179978 -0.000088
2.184978 0.000873
2.189978 0.000873
2.194978 0.000552
2.199978 -0.001049
2.204978 -0.000408
2.209978 0.000552
2.214978 -0.000408
2.219978 0.000552
2.224978 0.000232
2.229978 0.001193
2.234978 -0.001049
2.239978 0.000873
2.244978 0.000232
2.249978 -0.002010
2.254977 -0.000408
2.259977 0.000552
2.264977 -0.001690
2.269977 -0.000729
2.274977 -0.000729
2.279977 -0.000729
2.284977 0.001513
2.289977 -0.003612
2.294977 -0.000408
2.299977 -0.002010
2.304977 0.001834
2.309977 -0.000729
2.314977 -0.000729
2.319977 0.002474
2.324977 -0.000408
2.329977 0.002154
2.334977 -0.000088
2.339977 -0.001369
2.344977 0.000873
2.349977 -0.001369
2.354976 -0.000408
2.359976 -0.001049
2.364976 0.001193
2.369976 0.000232
2.374976 0.000873
2.379976 0.000552
2.384976 0.001193
2.389976 -0.001690
2.394976 -0.001049
2.399976 -0.000729
2.404976 0.001834
2.409976 0.000232
2.414976 0.001513
2.419976 -0.000408
2.424976 -0.000088
2.429976 0.001513
2.434976 0.001193
2.439976 -0.000729
2.444976 0.000552
2.449976 -0.002010
2.454975 0.000873
2.459975 0.001834
2.464975 -0.001049
2.469975 -0.000088
2.474975 -0.000729
2.479975 -0.000729
2.484975 -0.000408
2.489975 -0.000729
2.494975 0.001834
2.499975 -0.000088
2.504975 0.000552
2.509975 -0.000088
2.514975 0.001834
2.519975 -0.000088
2.524975 -0.000408
2.529975 -0.000408
2.534975 -0.002010
2.539975 0.000873
2.544975 -0.000729
2.549975 -0.000408
2.554974 -0.000408
2.559974 -0.000088
2.564974 -0.000408
2.569974 -0.001049
2.574974 0.000873
2.579974 -0.000088
2.584974 0.000873
2.589974 -0.002010
2.594974 0.000873
2.599974 -0.001690
2.604974 -0.000408
2.609974 0.001193
2.614974 -0.001049
2.619974 -0.000088
2.624974 -0.002010
2.629974 -0.001369
2.634974 0.001193
2.639974 0.001513
2.644974 -0.002010
2.649974 -0.000088
2.654973 0.000232
2.659973 -0.001369
2.664973 0.000873
2.669973 -0.000408
2.674973 -0.004252
2.679973 -0.002330
2.684973 0.000232
2.689973 -0.002330
2.694973 -0.000088
2.699973 -0.000408
2.704973 -0.000408
2.709973 0.000552
2.714973 -0.001369
2.719973 0.001513
2.724973 0.000232
2.729973 -0.001049
2.734973 0.001513
2.739973 0.000873
2.744973 -0.000088
2.749973 -0.000408
2.754972 0.001193
2.759972 -0.000408
2.764972 0.000232
2.769972 0.001193
2.774972 -0.000729
2.779972 0.000232
2.784972 -0.000729
2.789972 0.001193
2.794972 0.000873
2.799972 -0.001049
2.804972 -0.003932
2.809972 -0.000088
2.814972 -0.000088
2.819972 0.000873
2.824972 -0.000408
2.829972 -0.002010
2.834972 -0.000408
2.839972 -0.000088
2.844972 -0.001369
2.849972 0.001834
2.854971 -0.000088
2.859971 -0.000408
2.864971 -0.000408
2.869971 0.000552
2.874971 0.001513
2.879971 -0.000088
2.884971 -0.000408
2.889971 -0.000729
2.894971 0.000552
2.899971 -0.000729
2.904971 0.001193
2.909971 -0.001690
2.914971 0.000232
2.919971 -0.000408
2.924971 -0.000408
2.929971 -0.000408
2.934971 -0.000088
2.939971 0.002154
2.944971 0.001513
2.949971 -0.000729
2.954970 0.001834
2.959970 0.001193
2.964970 0.000873
2.969970 0.002154
2.974970 -0.000408
2.979970 -0.001690
2.984970 0.001193
2.989970 -0.000088
2.994970 -0.000088
2.999970 -0.001369
3.004970 0.000552
3.009970 -0.001369
3.014970 0.000873
3.019970 -0.000408
3.024970 0.000552
3.029970 -0.001049
3.034970 -0.000088
3.039970 -0.000729
3.044970 -0.001369
3.049970 0.000552
3.054969 -0.000088
3.059969 -0.000729
3.064969 -0.000088
3.069969 -0.000729
3.074969 -0.001369
3.079969 -0.000088
3.084969 -0.000088
3.089969 -0.000408
3.094969 0.000552
3.099969 -0.000088
3.104969 -0.001049
3.109969 -0.000088
3.114969 -0.000088
3.119969 0.002154
3.124969 -0.001369
3.129969 -0.000408
3.134969 -0.001690
3.139969 -0.000729
3.144969 -0.000729
3.149969 0.000552
3.154968 -0.000408
3.159968 -0.000729
3.164968 0.000552
3.169968 -0.001369
3.174968 -0.001369
3.179968 -0.001049
3.184968 -0.000729
3.189968 -0.000088
3.194968 0.000232
3.199968 0.000873
3.204968 0.000873
3.209968 -0.000088
3.214968 -0.001049
3.219968 -0.000408
3.224968 -0.000088
3.229968 -0.004252
3.234968 0.001513
3.239968 -0.002010
3.244968 -0.002651
3.249968 -0.001049
3.254967 0.000873
3.259967 0.000232
3.264967 0.000873
3.269967 -0.000088
3.274967 0.000873
3.279967 -0.000408
3.284967 -0.000408
3.289967 0.002474
3.294967 -0.001369
3.299967 0.001193
3.304967 0.000873
3.309967 -0.001690
3.314967 -0.000088
3.319967 -0.000729
3.324967 -0.000088
3.329967 -0.000088
3.334967 0.000873
3.339967 0.000552
3.344967 -0.000088
3.349967 -0.001049
3.354966 -0.001369
3.359966 -0.001049
3.364966 0.000873
3.369966 0.004396
3.374966 -0.000408
3.379966 -0.001049
3.384966 -0.000088
3.389966 -0.000088
3.394966 0.002795
3.399966 0.001193
3.404966 -0.001369
3.409966 -0.002651
3.414966 0.000232
3.419966 0.000873
3.424966 -0.000088
3.429966 -0.000729
3.434966 0.000552
3.439966 0.001193
3.444966 0.000873
3.449966 0.001193
3.454965 -0.000088
3.459965 -0.000729
3.464965 -0.000088
3.469965 0.001513
3.474965 -0.000088
3.479965 -0.000729
3.484965 -0.000729
3.489965 -0.000408
3.494965 0.002474
3.499965 -0.000729
3.504965 -0.000729
3.509965 -0.000088
3.514965 -0.000088
3.519965 -0.000088
3.524965 -0.000408
3.529965 -0.000088
3.534965 -0.000088
3.539965 0.000232
3.544965 -0.000408
3.549965 -0.000408
3.554964 -0.000408
3.559964 -0.000088
3.564964 0.000552
3.569964 0.000552
3.574964 0.001513
3.579964 0.000232
3.584964 -0.001049
3.589964 0.000873
3.594964 0.001193
3.599964 -0.001690
3.604964 0.000232
3.609964 -0.001690
3.614964 0.001513
3.619964 -0.002010
3.624964 -0.000088
3.629964 -0.000408
3.634964 -0.002330
3.639964 0.002154
3.644964 -0.000088
3.649964 0.002154
3.654963 -0.000408
3.659963 0.001834
3.664963 -0.001049
3.669963 -0.000729
3.674963 -0.000408
3.679963 0.000552
3.684963 -0.000729
3.689963 -0.001369
3.694963 -0.001369
3.699963 -0.002010
3.704963 0.000552
3.709963 0.000552
3.714963 -0.000408
3.719963 0.000232
3.724963 -0.000729
3.729963 -0.003612
3.734963 0.001513
3.739963 0.000232
3.744963 -0.000408
3.749963 -0.000729
3.754962 0.000232
3.759962 0.001193
3.764962 -0.000088
3.769962 -0.001049
3.774962 0.000232
3.779962 0.001513
3.784962 -0.000729
3.789962 -0.000729
3.794962 0.000873
3.799962 0.000552
3.804962 -0.001049
3.809962 -0.001049
3.814962 -0.000088
3.819962 -0.000088
3.824962 -0.002010
3.829962 -0.000088
3.834962 0.000873
3.839962 0.002474
3.844962 0.000552
3.849962 -0.000729
3.854961 -0.001369
3.859961 0.000873
3.864961 -0.000408
3.869961 0.000552
3.874961 0.001193
3.879961 0.000232
3.884961 0.000232
3.889961 -0.000088
3.894961 0.000232
3.899961 -0.000729
3.904961 -0.001690
3.909961 -0.000088
3.914961 -0.000088
3.919961 0.000552
3.924961 -0.000408
3.929961 0.000552
3.934961 -0.001690
3.939961 0.001834
3.944961 0.000873
3.949961 -0.000729
3.954960 -0.000408
3.959960 0.000552
3.964960 -0.001369
3.969960 0.000873
3.974960 -0.000088
3.979960 -0.001049
3.984960 -0.000088
3.989960 -0.000088
3.994960 -0.000729
3.999960 -0.000729
4.004960 0.000232
4.009960 -0.001690
4.014960 0.000873
4.019960 0.001834
4.024960 -0.002330
4.029960 0.000552
4.034960 -0.001049
4.039960 0.000873
4.044960 0.000552
4.049960 -0.000088
4.054959 -0.000729
4.059959 -0.001369
4.064959 0.000232
4.069959 -0.000088
4.074959 0.000552
4.079959 -0.000408
4.084959 0.002154
4.089959 -0.001049
4.094959 0.000552
4.099959 0.001834
4.104959 -0.000408
4.109959 -0.000408
4.114959 -0.000088
4.119959 -0.001369
4.124959 0.000552
4.129959 -0.000408
4.134959 0.000232
4.139959 -0.000729
4.144959 -0.002330
4.149959 -0.001049
4.154958 -0.000729
4.159958 -0.000729
4.164958 -0.000729
4.169958 0.001513
4.174958 -0.001049
4.179958 -0.000729
4.184958 -0.000088
4.189958 -0.001369
4.194958 0.000552
4.199958 -0.001690
4.204958 -0.001049
4.209958 -0.001049
4.214958 -0.000729
4.219958 -0.000408
4.224958 0.000552
4.229958 0.003435
4.234958 -0.001049
4.239958 0.000552
4.244958 0.001513
4.249958 0.000873
4.254957 0.001834
4.259957 0.000232
4.264957 -0.000088
4.269957 0.000873
4.274957 -0.002651
4.279957 0.002154
4.284957 0.000873
4.289957 0.000232
4.294957 0.000552
4.299957 -0.001369
4.304957 -0.001049
4.309957 0.000232
4.314957 0.000873
4.319957 -0.001690
4.324957 -0.000408
4.329957 -0.001049
4.334957 -0.001049
4.339957 0.001193
4.344957 0.001513
4.349957 -0.000729
4.354956 0.000873
4.359956 0.000232
4.364956 0.000552
4.369956 0.000552
4.374956 -0.000088
4.379956 0.000873
4.384956 -0.000729
4.389956 0.000552
4.394956 -0.000729
4.399956 -0.000408
4.404956 0.000873
4.409956 0.001834
4.414956 -0.000408
4.419956 -0.002010
4.424956 -0.001369
4.429956 -0.000729
4.434956 0.002795
4.439956 -0.000088
4.444956 -0.000408
4.449956 0.000232
4.454955 -0.001049
4.459955 0.001513
4.464955 -0.000088
4.469955 0.000873
4.474955 -0.000088
4.479955 0.000552
4.484955 0.000873
4.489955 -0.000729
4.494955 -0.000729
4.499955 -0.001690
4.504955 -0.000088
4.509955 0.000232
4.514955 -0.001049
4.519955 -0.000729
4.524955 0.000873
4.529955 -0.001369
4.534955 0.000552
4.539955 -0.001690
4.544955 0.000232
4.549955 0.001513
4.554954 0.000232
4.559954 0.000873
4.564954 0.001834
4.569954 -0.001369
4.574954 -0.001049
4.579954 0.000873
4.584954 -0.000408
4.589954 -0.001049
4.594954 0.001193
4.599954 -0.001049
4.604954 -0.001049
4.609954 0.000552
4.614954 0.000873
4.619954 -0.000729
4.624954 0.000552
4.629954 0.001513
4.634954 -0.000408
4.639954 -0.001049
4.644954 -0.002651
4.649954 -0.001690
4.654953 -0.000088
4.659953 -0.001049
4.664953 0.000232
4.669953 -0.000408
4.674953 0.001193
4.679953 -0.000408
4.684953 -0.000088
4.689953 -0.001049
4.694953 -0.001369
4.699953 0.000873
4.704953 0.002795
4.709953 -0.000729
4.714953 0.001513
4.719953 -0.000088
4.724953 0.000232
4.729953 -0.000088
4.734953 0.001193
4.739953 -0.000408
4.744953 0.002154
4.749953 -0.000408
4.754952 0.000873
4.759952 0.000552
4.764952 0.000873
4.769952 -0.000088
4.774952 -0.000729
4.779952 0.000873
4.784952 -0.001049
4.789952 -0.001049
4.794952 0.000552
4.799952 0.000232
4.804952 0.001513
4.809952 -0.000729
4.814952 0.000552
4.819952 -0.001049
4.824952 -0.000088
4.829952 0.000232
4.834952 -0.001690
4.839952 0.000232
4.844952 0.000873
4.849952 -0.001369
4.854951 0.000232
4.859951 -0.001049
4.864951 -0.000408
4.869951 0.003115
4.874951 0.000552
4.879951 -0.001369
4.884951 0.000873
4.889951 0.002474
4.894951 -0.001049
4.899951 0.001193
4.904951 -0.000408
4.909951 -0.000408
4.914951 0.001513
4.919951 -0.000729
4.924951 0.000232
4.929951 -0.000088
4.934951 -0.000088
4.939951 0.001193
4.944951 -0.000088
4.949951 -0.000088
4.954950 0.001193
4.959950 -0.000088
4.964950 -0.001049
4.969950 -0.000088
4.974950 0.000232
4.979950 0.000552
4.984950 -0.000408
4.989950 -0.000088
4.994950 -0.002010
4.999950 0.000552
5.004950 0.000873
5.009950 0.000873
5.014950 -0.000729
5.019950 0.000552
5.024950 -0.000408
5.029950 -0.000088
5.034950 0.000552
5.039950 0.000232
5.044950 0.000232
5.049950 0.000552
5.054949 0.000873
5.059949 -0.000088
5.064949 0.000873
5.069949 -0.000408
5.074949 0.000873
5.079949 0.000232
5.084949 0.000873
5.089949 0.000552
5.094949 -0.000408
5.099949 -0.000088
5.104949 -0.001049
5.109949 -0.001049
5.114949 0.001193
5.119949 -0.000088
5.124949 -0.000088
5.129949 -0.000408
5.134949 -0.000408
5.139949 0.000552
5.144949 -0.000408
5.149949 -0.000729
5.154948 0.000552
5.159948 -0.000408
5.164948 -0.001049
5.169948 0.000873
5.174948 0.000232
5.179948 0.001193
5.184948 -0.000088
5.189948 -0.000088
5.194948 -0.000729
5.199948 0.001193
5.204948 0.000552
5.209948 0.000552
5.214948 0.000552
5.219948 0.001193
5.224948 0.000232
5.229948 0.000552
5.234948 -0.000088
5.239948 0.000232
5.244948 0.000232
5.249948 0.000552
5.254947 0.001513
5.259947 -0.003612
5.264947 0.001513
5.269947 0.000552
5.274947 0.000552
5.279947 -0.001690
5.284947 -0.000088
5.289947 -0.000088
5.294947 -0.001369
5.299947 -0.001369
5.304947 -0.000088
5.309947 -0.000408
5.314947 0.001193
5.319947 -0.001049
5.324947 0.001513
5.329947 -0.002010
5.334947 -0.000729
5.339947 -0.000088
5.344947 0.000552
5.349947 0.000552
5.354946 0.000873
5.359946 0.000873
5.364946 0.001193
5.369946 -0.000729
5.374946 0.001513
5.379946 0.001193
5.384946 -0.000088
5.389946 0.004076
5.394946 -0.001690
5.399946 -0.000408
5.404946 -0.002651
5.409946 0.001193
5.414946 -0.000088
5.419946 0.003115
5.424946 0.001193
5.429946 -0.000408
5.434946 -0.000729
5.439946 0.000552
5.444946 0.000232
5.449946 -0.000729
5.454945 -0.000408
5.459945 -0.000408
5.464945 -0.000729
5.469945 -0.000088
5.474945 -0.000729
5.479945 -0.000088
5.484945 -0.002010
5.489945 0.001193
5.494945 0.000552
5.499945 -0.000088
5.504945 0.000552
5.509945 0.000232
5.514945 -0.000088
5.519945 -0.000408
5.524945 0.000552
5.529945 0.001513
5.534945 0.000552
5.539945 -0.001369
5.544945 -0.000729
5.549945 -0.001049
5.554944 -0.000408
5.559944 0.000552
5.564944 0.003115
5.569944 -0.000408
5.574944 0.000232
5.579944 -0.001369
5.584944 0.000552
5.589944 0.000873
5.594944 -0.000408
5.599944 0.000552
5.604944 -0.001049
5.609944 -0.000088
5.614944 0.000232
5.619944 -0.001049
5.624944 -0.000408
5.629944 -0.002010
5.634944 -0.000088
5.639944 -0.000729
5.644944 -0.000408
5.649944 0.000552
5.654943 0.001193
5.659943 0.002795
5.664943 -0.001690
5.669943 0.000232
5.674943 -0.001049
5.679943 -0.000729
5.684943 -0.001369
5.689943 0.000552
5.694943 -0.000729
5.699943 -0.001369
5.704943 -0.001049
5.709943 0.000873
5.714943 0.002154
5.719943 0.000552
5.724943 0.000873
5.729943 0.000873
5.734943 -0.000088
5.739943 -0.000408
5.744943 0.002154
5.749943 0.000873
5.754942 -0.000088
5.759942 0.002795
5.764942 -0.000088
5.769942 -0.000088
5.774942 0.002795
5.779942 0.001513
5.784942 -0.001049
5.789942 0.000873
5.794942 -0.000088
5.799942 0.000873
5.804942 0.002795
5.809942 -0.001690
5.814942 -0.000729
5.819942 -0.002330
5.824942 -0.000088
5.829942 -0.002651
5.834942 -0.001690
5.839942 -0.001049
5.844942 0.000552
5.849942 -0.000088
5.854941 -0.000408
5.859941 -0.000408
5.864941 0.000552
5.869941 0.000232
5.874941 0.000232
5.879941 -0.001690
5.884941 -0.000729
5.889941 0.001193
5.894941 0.000552
5.899941 0.000552
5.904941 0.000552
5.909941 0.000232
5.914941 0.001193
5.919941 -0.000408
5.924941 -0.001690
5.929941 0.000552
5.934941 -0.000088
5.939941 -0.001049
5.944941 -0.000408
5.949941 -0.000408
5.954940 0.001513
5.959940 0.000232
5.964940 0.000232
5.969940 -0.002010
5.974940 -0.001049
5.979940 0.001193
5.984940 0.000873
5.989940 0.001193
5.994940 -0.001690
5.999940 -0.000088
6.004940 0.000873
6.009940 -0.001369
6.014940 0.001193
6.019940 -0.000408
6.024940 -0.000729
6.029940 -0.001369
6.034940 0.001513
6.039940 -0.000088
6.044940 0.000552
6.049940 0.000552
6.054939 0.001834
6.059939 0.000232
6.064939 -0.001369
6.069939 0.000552
6.074939 -0.000408
6.079939 -0.000408
6.084939 0.000552
6.089939 -0.001369
6.094939 -0.000729
6.099939 0.001834
6.104939 -0.001690
6.109939 0.001513
6.114939 0.000873
6.119939 -0.000408
6.124939 0.001513
6.129939 -0.002010
6.134939 0.000232
6.139939 -0.002010
6.144939 -0.000408
6.149939 -0.001049
6.154938 -0.000729
6.159938 -0.000088
6.164938 0.000552
6.169938 0.000232
6.174938 0.000232
6.179938 0.001193
6.184938 -0.000729
6.189938 -0.001049
6.194938 -0.001690
6.199938 0.000232
6.204938 0.000232
6.209938 0.001513
6.214938 -0.000088
6.219938 -0.000729
6.224938 -0.000729
6.229938 -0.000088
6.234938 0.000232
6.239938 0.001193
6.244938 -0.000088
6.249938 -0.000088
6.254937 -0.000729
6.259937 0.001193
6.264937 0.001513
6.269937 -0.000088
6.274937 -0.000088
6.279937 -0.000729
6.284937 0.004076
6.289937 0.001834
6.294937 0.000232
6.299937 -0.003932
6.304937 -0.000088
6.309937 0.000873
6.314937 0.002795
6.319937 0.001834
6.324937 -0.001049
6.329937 -0.001690
6.334937 -0.000408
6.339937 0.000873
6.344937 0.002154
6.349937 -0.000729
6.354936 -0.000729
6.359936 0.000232
6.364936 -0.000408
6.369936 -0.001049
6.374936 -0.002010
6.379936 -0.001049
6.384936 -0.000729
6.389936 0.000873
6.394936 0.001834
6.399936 -0.000088
6.404936 -0.002010
6.409936 0.001513
6.414936 -0.002010
6.419936 0.000873
6.424936 -0.001369
6.429936 -0.000088
6.434936 -0.000088
6.439936 0.000552
6.444936 0.000552
6.449936 0.000552
6.454935 -0.003932
6.459935 0.000232
6.464935 -0.001049
6.469935 -0.002971
6.474935 0.001513
6.479935 -0.001049
6.484935 -0.000729
6.489935 -0.000408
6.494935 0.001193
6.499935 -0.002010
6.504935 -0.001369
6.509935 -0.000408
6.514935 0.001513
6.519935 -0.000088
6.524935 0.000873
6.529935 0.001513
6.534935 -0.001049
6.539935 -0.001049
6.544935 0.001513
6.549935 -0.001049
6.554934 -0.000088
6.559934 -0.000729
6.564934 -0.002010
6.569934 0.000873
6.574934 0.000232
6.579934 0.002154
6.584934 -0.000729
6.589934 -0.003612
6.594934 -0.000088
6.599934 0.003115
6.604934 0.002154
6.609934 -0.001369
6.614934 -0.000729
6.619934 -0.000729
6.624934 -0.000729
6.629934 0.002474
6.634934 -0.000729
6.639934 -0.001690
6.644934 -0.002010
6.649934 -0.000088
6.654933 0.000232
6.659933 -0.000729
6.664933 -0.000408
6.669933 -0.002651
6.674933 -0.000408
6.679933 -0.000408
6.684933 -0.000088
6.689933 -0.001690
6.694933 -0.001690
6.699933 -0.000088
6.704933 0.001513
6.709933 0.002474
6.714933 0.002474
6.719933 0.000873
6.724933 0.000552
6.729933 -0.000088
6.734933 -0.000408
6.739933 -0.000729
6.744933 -0.001690
6.749933 -0.000729
6.754932 0.001193
6.759932 0.000552
6.764932 -0.000729
6.769932 -0.002330
6.774932 -0.000408
6.779932 -0.000729
6.784932 0.001513
6.789932 0.001513
6.794932 -0.000408
6.799932 -0.002971
6.804932 0.000552
6.809932 -0.001369
6.814932 0.002795
6.819932 0.001513
6.824932 -0.002010
6.829932 -0.000408
6.834932 0.001193
6.839932 0.000552
6.844932 -0.002330
6.849932 -0.000729
6.854931 0.001193
6.859931 0.000873
6.864931 -0.000088
6.869931 0.000552
6.874931 -0.001690
6.879931 -0.000088
6.884931 0.000232
6.889931 -0.001690
6.894931 -0.000088
6.899931 -0.000408
6.904931 -0.000088
6.909931 0.000232
6.914931 -0.000088
6.919931 -0.000088
6.924931 -0.000408
6.929931 -0.001690
6.934931 0.000552
6.939931 -0.000729
6.944931 0.000552
6.949931 -0.001369
6.954930 -0.000408
6.959930 0.000552
6.964930 0.002795
6.969930 -0.002010
6.974930 -0.001690
6.979930 -0.002010
6.984930 0.000232
6.989930 0.001834
6.994930 0.000232
6.999930 -0.001049
7.004930 0.001834
7.009930 0.000552
7.014930 -0.000088
7.019930 -0.001049
7.024930 -0.000408
7.029930 -0.001049
7.034930 -0.001369
7.039930 -0.000088
7.044930 0.000552
7.049930 -0.001369
7.054929 -0.002971
7.059929 -0.000408
7.064929 -0.000408
7.069929 0.002154
7.074929 -0.000088
7.079929 -0.000088
7.084929 -0.000088
7.089929 0.001193
7.094929 -0.000729
7.099929 0.001834
7.104929 0.000873
7.109929 -0.000729
7.114929 0.000873
7.119929 -0.000088
7.124929 0.002154
7.129929 -0.000088
7.134929 -0.000729
7.139929 0.001193
7.144929 0.000232
7.149929 0.002154
7.154928 -0.001369
7.159928 0.000552
7.164928 0.001834
7.169928 0.000552
7.174928 -0.000088
7.179928 0.002474
7.184928 -0.003291
7.189928 -0.000729
7.194928 0.002474
7.199928 0.001834
7.204928 -0.000729
7.209928 0.003115
7.214928 -0.000729
7.219928 -0.000408
7.224928 0.000552
7.229928 -0.000408
7.234928 0.000552
7.239928 -0.001049
7.244928 0.000873
7.249928 -0.000729
7.254927 -0.001049
7.259927 0.000873
7.264927 -0.000088
7.269927 -0.000408
7.274927 0.001193
7.279927 -0.000729
7.284927 -0.000088
7.289927 0.001193
7.294927 0.000873
7.299927 -0.000408
7.304927 0.001193
7.309927 -0.001369
7.314927 -0.000088
7.319927 -0.002010
7.324927 -0.000088
7.329927 -0.000088
7.334927 -0.001369
7.339927 0.001193
7.344927 -0.001049
7.349927 0.000873
7.354926 -0.002010
7.359926 -0.003612
7.364926 0.003115
7.369926 0.000552
7.374926 -0.000088
7.379926 -0.001690
7.384926 0.000232
7.389926 0.002154
7.394926 0.000552
7.399926 -0.000088
7.404926 0.001193
7.409926 0.000232
7.414926 -0.000088
7.419926 0.002154
7.424926 -0.002010
7.429926 -0.000729
7.434926 0.001193
7.439926 -0.000729
7.444926 0.000873
7.449926 -0.001690
7.454925 0.001834
7.459925 -0.001369
7.464925 -0.000729
7.469925 -0.002010
7.474925 -0.000729
7.479925 -0.002971
7.484925 -0.000088
7.489925 0.001513
7.494925 -0.000408
7.499925 0.000232
7.504925 -0.000408
7.509925 -0.000408
7.514925 0.000873
7.519925 -0.001049
7.524925 0.001834
7.529925 -0.002010
7.534925 -0.000729
7.539925 -0.000088
7.544925 0.000232
7.549925 -0.000088
7.554924 0.001193
7.559924 0.000552
7.564924 0.000873
7.569924 0.000552
7.574924 -0.000408
7.579924 0.000873
7.584924 0.001193
7.589924 -0.000408
7.594924 0.002474
7.599924 -0.001049
7.604924 -0.000408
7.609924 -0.000408
7.614924 0.000552
7.619924 0.001834
7.624924 -0.000729
7.629924 0.001193
7.634924 -0.002330
7.639924 -0.001690
7.644924 -0.000408
7.649924 -0.000088
7.654923 -0.001369
7.659923 0.000232
7.664923 -0.001369
7.669923 -0.001369
7.674923 0.003435
7.679923 -0.000729
7.684923 -0.001369
7.689923 -0.000088
7.694923 -0.001369
7.699923 -0.001049
7.704923 0.000552
7.709923 -0.000729
7.714923 -0.000408
7.719923 0.001513
7.724923 0.001193
7.729923 0.001193
7.734923 0.001193
7.739923 -0.001690
7.744923 0.001834
7.749923 -0.000729
7.754922 -0.001049
7.759922 -0.001369
7.764922 -0.001690
7.769922 0.000232
7.774922 -0.000088
7.779922 -0.000408
7.784922 0.001834
7.789922 0.002474
7.794922 -0.001049
7.799922 -0.000408
7.804922 -0.001690
7.809922 -0.000088
7.814922 0.000552
7.819922 -0.000729
7.824922 0.000552
7.829922 -0.000729
7.834922 -0.001369
7.839922 -0.000088
7.844922 0.000232
7.849922 -0.000088
7.854921 -0.000729
7.859921 -0.001369
7.864921 -0.002330
7.869921 -0.001049
7.874921 0.000552
7.879921 0.000232
7.884921 -0.000408
7.889921 -0.000729
7.894921 0.003756
7.899921 -0.000729
7.904921 0.000232
7.909921 -0.001690
7.914921 0.001193
7.919921 0.002154
7.924921 -0.001049
7.929921 -0.002010
7.934921 0.001193
7.939921 0.000552
7.944921 -0.000729
7.949921 0.000552
7.954920 -0.000729
7.959920 0.001513
7.964920 -0.000729
7.969920 0.000873
7.974920 -0.000088
7.979920 0.000232
7.984920 0.000873
7.989920 -0.000088
7.994920 -0.000408
7.999920 0.000552
8.004920 0.001193
8.009920 0.002154
8.014920 -0.000088
8.019920 0.000873
8.024920 0.002474
8.029920 0.001193
8.034920 -0.000729
8.039920 0.000552
8.044920 0.002154
8.049920 0.000232
8.054919 -0.000408
8.059919 0.000873
8.064919 0.001513
8.069919 0.000873
8.074919 0.000873
8.079919 0.000232
8.084919 -0.000408
8.089919 -0.000729
8.094919 0.002795
8.099919 0.000873
8.104919 -0.001369
8.109919 -0.002010
8.114919 -0.001049
8.119919 0.001513
8.124919 0.003435
8.129919 -0.001690
8.134919 0.001193
8.139919 -0.002971
8.144919 -0.001049
8.149919 -0.000729
8.154918 0.002474
8.159918 0.000873
8.164918 -0.000729
8.169918 -0.002010
8.174918 0.000232
8.179918 0.002154
8.184918 -0.000408
8.189918 0.000552
8.194918 -0.001049
8.199918 -0.001369
8.204918 -0.000088
8.209918 -0.001049
8.214918 0.001513
8.219918 -0.000408
8.224918 -0.000088
8.229918 -0.000408
8.234918 0.000552
8.239918 -0.001049
8.244918 -0.001049
8.249918 -0.000088
8.254917 0.001193
8.259917 0.000552
8.264917 -0.000729
8.269917 0.000232
8.274917 0.000552
8.279917 0.000873
8.284917 -0.001369
8.289917 -0.001049
8.294917 -0.001690
8.299917 -0.001369
8.304917 -0.000729
8.309917 -0.000408
8.314917 -0.001690
8.319917 0.000552
8.324917 -0.000088
8.329917 0.000232
8.334917 -0.000408
8.339917 -0.000088
8.344917 0.002154
8.349917 -0.000088
8.354916 -0.000408
8.359916 0.002474
8.364916 0.000873
8.369916 -0.000729
8.374916 0.000873
8.379916 0.000873
8.384916 -0.000729
8.389916 0.000873
8.394916 0.000873
8.399916 -0.000408
8.404916 -0.002330
8.409916 0.002474
8.414916 -0.000408
8.419916 -0.001690
8.424916 0.001193
8.429916 -0.002651
8.434916 -0.000088
8.439916 -0.000088
8.444916 -0.002010
8.449916 -0.000088
8.454915 -0.001690
8.459915 -0.000408
8.464915 -0.000408
8.469915 -0.001690
8.474915 0.000552
8.479915 -0.001049
8.484915 -0.000408
8.489915 -0.000729
8.494915 0.000552
8.499915 -0.000408
8.504915 -0.000408
8.509915 0.000552
8.514915 -0.001369
8.519915 -0.000088
8.524915 -0.000088
8.529915 -0.001369
8.534915 -0.001049
8.539915 -0.000729
8.544915 0.000552
8.549915 -0.000729
8.554914 0.001193
8.559914 -0.000408
8.564914 0.000232
8.569914 -0.001690
8.574914 0.003115
8.579914 -0.001049
8.584914 0.000873
8.589914 -0.000729
8.594914 -0.000729
8.599914 -0.000408
8.604914 -0.002330
8.609914 -0.000408
8.614914 -0.000729
8.619914 -0.000729
8.624914 0.001834
8.629914 -0.002010
8.634914 0.000552
8.639914 -0.001690
8.644914 -0.000088
8.649914 -0.000408
8.654913 0.000552
8.659913 -0.000088
8.664913 -0.001690
8.669913 0.000552
8.674913 0.001193
8.679913 0.001513
8.684913 -0.001049
8.689913 -0.000729
8.694913 -0.000408
8.699913 0.000232
8.704913 -0.000088
8.709913 -0.000729
8.714913 -0.001690
8.719913 -0.000088
8.724913 0.000552
8.729913 -0.002010
8.734913 0.001834
8.739913 0.002154
8.744913 -0.002010
8.749913 0.000232
8.754912 -0.001690
8.759912 -0.000088
8.764912 -0.001049
8.769912 0.000873
8.774912 -0.001369
8.779912 0.000873
8.784912 -0.001369
8.789912 -0.001690
8.794912 0.000873
8.799912 -0.002010
8.804912 -0.000088
8.809912 -0.002010
8.814912 -0.000088
8.819912 -0.004893
8.824912 -0.001369
8.829912 -0.000729
8.834912 0.002154
8.839912 0.002795
8.844912 -0.001049
8.849912 -0.000088
8.854911 -0.000408
8.859911 -0.000088
8.864911 0.000873
8.869911 -0.000729
8.874911 -0.000729
8.879911 -0.000408
8.884911 0.001193
8.889911 0.002154
8.894911 -0.001369
8.899911 -0.000408
8.904911 0.000873
8.909911 -0.000088
8.914911 -0.000088
8.919911 -0.000408
8.924911 -0.001690
8.929911 -0.002971
8.934911 -0.002010
8.939911 0.002795
8.944911 0.001193
8.949911 -0.001690
8.954910 -0.002971
8.959910 -0.000729
8.964910 0.002795
8.969910 0.001513
8.974910 -0.000408
8.979910 -0.000408
8.984910 -0.001690
8.989910 -0.000088
8.994910 0.000552
8.999910 0.000552
9.004910 0.001513
9.009910 -0.002010
9.014910 -0.000729
9.019910 0.002474
9.024910 -0.002010
9.029910 -0.000729
9.034910 -0.001049
9.039910 -0.002651
9.044910 0.001834
9.049910 -0.001690
9.054909 0.001193
9.059909 -0.001049
9.064909 -0.000408
9.069909 0.001513
9.074909 0.001834
9.079909 0.000873
9.084909 -0.001690
9.089909 0.000232
9.094909 -0.000088
9.099909 -0.001369
9.104909 -0.002330
9.109909 -0.003291
9.114909 -0.001690
9.119909 -0.002010
9.124909 0.000873
9.129909 -0.001049
9.134909 0.000552
9.139909 -0.001049
9.144909 -0.000729
9.149909 -0.000729
9.154908 0.001834
9.159908 0.000873
9.164908 0.000552
9.169908 -0.002330
9.174908 0.001834
9.179908 0.001193
9.184908 -0.000408
9.189908 -0.000729
9.194908 -0.000408
9.199908 -0.001369
9.204908 0.001513
9.209908 0.000552
9.214908 0.000873
9.219908 0.000552
9.224908 0.000873
9.229908 0.001834
9.234908 -0.000408
9.239908 -0.001369
9.244908 0.001193
9.249908 -0.000408
9.254907 0.003435
9.259907 -0.001369
9.264907 0.000232
9.269907 -0.001690
9.274907 0.000873
9.279907 -0.000088
9.284907 0.000552
9.289907 -0.001049
9.294907 -0.001690
9.299907 0.004717
9.304907 0.000552
9.309907 0.001834
9.314907 0.000873
9.319907 -0.003291
9.324907 0.000552
9.329907 0.000873
9.334907 -0.001049
9.339907 -0.001690
9.344907 -0.000408
9.349907 -0.001690
9.354906 -0.000088
9.359906 0.001513
9.364906 0.001193
9.369906 0.001193
9.374906 -0.002010
9.379906 0.000552
9.384906 0.001834
9.389906 -0.000088
9.394906 -0.000088
9.399906 -0.002010
9.404906 0.001513
9.409906 0.002795
9.414906 0.001193
9.419906 -0.002010
9.424906 0.000232
9.429906 -0.002651
9.434906 0.001834
9.439906 0.000232
9.444906 0.000552
9.449906 -0.001049
9.454905 -0.001049
9.459905 0.000552
9.464905 0.003115
9.469905 0.000552
9.474905 -0.001690
9.479905 -0.003612
9.484905 -0.000729
9.489905 0.000232
9.494905 -0.001690
9.499905 0.001513
9.504905 0.000873
9.509905 0.001513
9.514905 0.001193
9.519905 0.000873
9.524905 0.000873
9.529905 -0.002651
9.534905 -0.001369
9.539905 0.001193
9.544905 0.000873
9.549905 -0.001049
9.554904 -0.001690
9.559904 -0.000088
9.564904 0.000873
9.569904 -0.000408
9.574904 -0.001049
9.579904 -0.001049
9.584904 -0.001049
9.589904 -0.001049
9.594904 0.003115
9.599904 -0.002010
9.604904 -0.001690
9.609904 -0.000088
9.614904 0.000232
9.619904 -0.000088
9.624904 -0.000729
9.629904 -0.000408
9.634904 -0.000088
9.639904 0.000552
9.644904 0.003435
9.649904 0.001513
9.654903 -0.000729
9.659903 -0.002651
9.664903 -0.000088
9.669903 0.000873
9.674903 0.000232
9.679903 0.000552
9.684903 -0.001690
9.689903 -0.001690
9.694903 0.001193
9.699903 -0.001049
9.704903 0.001834
9.709903 -0.002010
9.714903 -0.001049
9.719903 -0.002010
9.724903 0.000552
9.729903 0.000552
9.734903 0.000873
9.739903 -0.002010
9.744903 -0.000408
9.749903 -0.001690
9.754902 0.001193
9.759902 -0.000408
9.764902 -0.000408
9.769902 0.000232
9.774902 0.000873
9.779902 0.000232
9.784902 -0.001369
9.789902 0.000552
9.794902 -0.000088
9.799902 0.000232
9.804902 -0.000408
9.809902 -0.001049
9.814902 -0.001049
9.819902 -0.001369
9.824902 -0.002651
9.829902 0.000232
9.834902 0.003115
9.839902 0.000873
9.844902 -0.001049
9.849902 0.003115
9.854901 0.003115
9.859901 -0.000088
9.864901 0.002154
9.869901 -0.001049
9.874901 0.000552
9.879901 0.002474
9.884901 0.000873
9.889901 0.001513
9.894901 -0.000729
9.899901 -0.001049
9.904901 0.002154
9.909901 0.000552
9.914901 0.000873
9.919901 -0.000088
9.924901 0.000232
9.929901 -0.001369
9.934901 -0.002010
9.939901 0.000873
9.944901 0.000552
9.949901 -0.000729
9.954900 0.001193
9.959900 0.000232
9.964900 -0.001690
9.969900 0.000873
9.974900 -0.000408
9.979900 -0.003612
9.984900 -0.000408
9.989900 0.000232
9.994900 -0.000088
9.999900 -0.001690
10.004900 0.000873
10.009900 -0.000088
10.014900 -0.001049
10.019900 0.000873
10.024900 -0.000729
10.029900 0.000232
10.034900 0.001513
10.039900 0.000552
10.044900 0.000552
10.049900 -0.000408
10.054899 -0.000408
10.059899 0.000232
10.064899 0.000873
10.069899 -0.000408
10.074899 -0.000729
10.079899 -0.002010
10.084899 -0.000088
10.089899 -0.000729
10.094899 0.001193
10.099899 -0.001049
10.104899 0.000552
10.109899 0.000232
10.114899 -0.001369
10.119899 0.000552
10.124899 0.000552
10.129899 0.000232
10.134899 -0.001049
10.139899 -0.000729
10.144899 0.000552
10.149899 0.001193
10.154898 -0.000088
10.159898 -0.001049
10.164898 -0.002010
10.169898 0.000552
10.174898 0.000232
10.179898 0.000873
10.184898 -0.001369
10.189898 0.001193
10.194898 0.000232
10.199898 -0.000729
10.204898 0.003756
10.209898 -0.001690
10.214898 0.000873
10.219898 0.000552
10.224898 -0.000729
10.229898 -0.000729
10.234898 -0.001049
10.239898 0.002474
10.244898 0.000232
10.249898 0.001513
10.254897 0.000873
10.259897 -0.001369
10.264897 0.000552
10.269897 0.001834
10.274897 0.000873
10.279897 0.000873
10.284897 0.000873
10.289897 0.000552
10.294897 -0.000729
10.299897 -0.000088
10.304897 -0.000088
10.309897 0.000552
10.314897 0.000232
10.319897 0.000873
10.324897 -0.000408
10.329897 0.002154
10.334897 -0.001369
10.339897 -0.001369
10.344897 0.002474
10.349897 -0.001049
10.354896 -0.000088
10.359896 0.000552
10.364896 -0.001369
10.369896 -0.001369
10.374896 0.001513
10.379896 0.001193
10.384896 -0.001690
10.389896 -0.000088
10.394896 0.001513
10.399896 0.000232
10.404896 0.000232
10.409896 -0.000088
10.414896 0.001193
10.419896 -0.000729
10.424896 0.003756
10.429896 -0.000088
10.434896 0.001513
10.439896 -0.002651
10.444896 -0.001369
10.449896 -0.000729
10.454895 0.000552
10.459895 -0.000088
10.464895 0.001513
10.469895 0.000552
10.474895 -0.001369
10.479895 -0.002651
10.484895 0.000552
10.489895 -0.001369
10.494895 0.001193
10.499895 -0.000729
10.504895 0.001834
10.509895 -0.000408
10.514895 -0.000088
10.519895 -0.001690
10.524895 0.000552
10.529895 0.000873
10.534895 -0.000088
10.539895 -0.000408
10.544895 -0.000088
10.549895 0.001513
10.554894 -0.002010
10.559894 0.000873
10.564894 0.000873
10.569894 -0.001369
10.574894 -0.000729
10.579894 -0.000088
10.584894 0.000873
10.589894 -0.002330
10.594894 -0.000408
10.599894 -0.000088
10.604894 0.001193
10.609894 -0.000088
10.614894 0.000552
10.619894 0.003756
10.624894 0.002474
10.629894 0.002154
10.634894 -0.003291
10.639894 0.002795
10.644894 0.000232
10.649894 -0.000408
10.654893 -0.000729
10.659893 0.000873
10.664893 -0.000408
10.669893 -0.000729
10.674893 -0.001369
10.679893 0.000552
10.684893 -0.000408
10.689893 -0.000729
10.694893 0.000552
10.699893 -0.001049
10.704893 0.001513
10.709893 0.000873
10.714893 -0.002971
10.719893 -0.000088
10.724893 -0.000408
10.729893 0.000232
10.734893 -0.000729
10.739893 -0.002010
10.744893 0.001193
10.749893 0.000552
10.754892 0.001193
10.759892 -0.000408
10.764892 0.000232
10.769892 0.000552
10.774892 0.000232
10.779892 0.000873
10.784892 -0.001690
10.789892 0.001834
10.794892 -0.000088
10.799892 -0.000088
10.804892 0.001834
10.809892 0.002154
10.814892 0.002154
10.819892 0.002154
10.824892 0.001193
10.829892 -0.000408
10.834892 -0.000408
10.839892 -0.001049
10.844892 0.003756
10.849892 0.000232
10.854891 -0.001369
10.859891 -0.000729
10.864891 -0.000729
10.869891 -0.000088
10.874891 -0.001049
10.879891 -0.000088
10.884891 -0.002010
10.889891 -0.002651
10.894891 -0.000408
10.899891 0.001513
10.904891 -0.000408
10.909891 -0.001690
10.914891 -0.001690
10.919891 -0.000729
10.924891 0.003115
10.929891 0.001834
10.934891 -0.002010
10.939891 -0.001690
10.944891 0.001193
10.949891 0.001513
10.954890 -0.000408
10.959890 -0.003291
10.964890 -0.000408
10.969890 0.001193
10.974890 -0.001690
10.979890 -0.000729
10.984890 0.001193
10.989890 -0.000088
10.994890 -0.000408
10.999890 -0.002010
11.004890 0.001193
11.009890 -0.001690
11.014890 0.000873
11.019890 -0.001049
11.024890 -0.000088
11.029890 -0.000408
11.034890 0.001834
11.039890 -0.000408
11.044890 -0.000088
11.049890 -0.005534
11.054889 -0.000088
11.059889 0.000873
11.064889 0.000873
11.069889 -0.000729
11.074889 -0.002010
11.079889 -0.001049
11.084889 0.000232
11.089889 0.000873
11.094889 -0.000088
11.099889 -0.002330
11.104889 0.000552
11.109889 -0.001690
11.114889 0.000873
11.119889 0.001834
11.124889 -0.000408
11.129889 -0.000088
11.134889 -0.001369
11.139889 -0.000088
11.144889 -0.000408
11.149889 0.002154
11.154888 0.000552
11.159888 0.000552
11.164888 -0.002010
11.169888 0.000873
11.174888 -0.000408
11.179888 -0.001049
11.184888 -0.000088
11.189888 -0.002971
11.194888 -0.000408
11.199888 -0.000088
11.204888 0.001193
11.209888 -0.000729
11.214888 -0.000088
11.219888 0.001513
11.224888 -0.002651
11.229888 -0.000408
11.234888 0.004076
11.239888 0.000232
11.244888 -0.001369
11.249888 -0.000088
11.254887 0.001513
11.259887 0.000873
11.264887 -0.000408
11.269887 0.002154
11.274887 -0.001049
11.279887 -0.002330
11.284887 -0.000729
11.289887 -0.000408
11.294887 0.001193
11.299887 -0.002010
11.304887 -0.000088
11.309887 0.000232
11.314887 0.000552
11.319887 -0.002651
11.324887 -0.000729
11.329887 -0.001369
11.334887 -0.000088
11.339887 -0.001690
11.344887 0.001513
11.349887 -0.000088
11.354886 -0.001690
11.359886 -0.000408
11.364886 -0.001049
11.369886 -0.000408
11.374886 0.002474
11.379886 -0.000408
11.384886 0.001834
11.389886 -0.000408
11.394886 -0.000088
11.399886 -0.000408
11.404886 -0.002010
11.409886 -0.001690
11.414886 0.000552
11.419886 0.000552
11.424886 -0.000088
11.429886 -0.001369
11.434886 0.002154
11.439886 0.001193
11.444886 0.001193
11.449886 -0.002010
11.454885 0.000552
11.459885 -0.001369
11.464885 -0.000088
11.469885 -0.000729
11.474885 -0.002010
11.479885 -0.001690
11.484885 -0.002651
11.489885 0.000873
11.494885 -0.002971
11.499885 -0.000729
11.504885 -0.002010
11.509885 -0.002971
11.514885 -0.001369
11.519885 -0.000088
11.524885 0.000552
11.529885 0.001193
11.534885 -0.001369
11.539885 0.000873
11.544885 0.002154
11.549885 0.000232
11.554884 0.003435
11.559884 0.001513
11.564884 -0.002010
11.569884 0.000552
11.574884 0.001834
11.579884 0.001834
11.584884 -0.002651
11.589884 0.001193
11.594884 -0.002010
11.599884 0.002795
11.604884 0.000873
11.609884 0.001513
11.614884 0.002154
11.619884 -0.000088
11.624884 -0.000088
11.629884 -0.003612
11.634884 0.000552
11.639884 0.000873
11.644884 0.000552
11.649884 -0.000088
11.654883 0.000552
11.659883 -0.000408
11.664883 0.001193
11.669883 0.001193
11.674883 0.001193
11.679883 -0.000088
11.684883 0.000232
11.689883 -0.001690
11.694883 -0.001049
11.699883 0.001193
11.704883 0.001193
11.709883 0.000552
11.714883 0.001513
11.719883 -0.001690
11.724883 -0.000088
11.729883 -0.001690
11.734883 -0.000729
11.739883 -0.001049
11.744883 0.001513
11.749883 -0.000729
11.754882 0.001193
11.759882 -0.000088
11.764882 0.002154
11.769882 -0.000088
11.774882 0.000873
11.779882 0.001193
11.784882 0.000552
11.789882 -0.001049
11.794882 -0.000729
11.799882 0.000232
11.804882 -0.000729
11.809882 -0.001369
11.814882 0.002795
11.819882 -0.002971
11.824882 -0.000088
11.829882 -0.000088
11.834882 -0.002330
11.839882 0.000873
11.844882 0.000552
11.849882 0.001513
11.854881 -0.000408
11.859881 0.000552
11.864881 0.001193
11.869881 0.000552
11.874881 0.002154
11.879881 -0.000088
11.884881 0.002474
11.889881 0.000873
11.894881 -0.000088
11.899881 0.000552
11.904881 -0.002010
11.909881 -0.000088
11.914881 -0.000088
11.919881 0.003435
11.924881 -0.000408
11.929881 -0.000729
11.934881 -0.002010
11.939881 0.000552
11.944881 -0.000088
11.949881 -0.000729
11.954880 -0.000408
11.959880 0.001193
11.964880 -0.001049
11.969880 0.001193
11.974880 -0.003291
11.979880 -0.000729
11.984880 -0.000408
11.989880 0.001513
11.994880 -0.000408
11.999880 -0.000088
12.004880 0.001193
12.009880 -0.001049
12.014880 0.000552
12.019880 -0.000408
12.024880 -0.001369
12.029880 0.000873
12.034880 0.000552
12.039880 -0.000729
12.044880 0.000232
12.049880 -0.000408
12.054879 0.001513
12.059879 0.000232
12.064879 0.000873
12.069879 0.001513
12.074879 0.000552
12.079879 -0.000088
data = readmatrix('obp3n.txt', 'HeaderLines',8) % Skip First 8 Header Lines
data = 100000×2
0 -0.0004
0.0050 -0.0010
0.0100 -0.0001
0.0150 -0.0010
0.0200 0.0034
0.0250 -0.0004
0.0300 -0.0007
0.0350 -0.0039
0.0400 -0.0004
0.0450 0.0025
nonfinite_rows = nnz(~any(isfinite(data),2)) % Check To See If Any Rows Have Non-Finite Values
nonfinite_rows = 0
if nonfinite_rows > 0
data(~isfinite(data)) = NaN;
data = fillmissing(data, 'linear');
end
x_data = data(:,1);
y_data = data(:,2);
Ts = mean(diff(x_data))
Ts = 0.0050
Tsd = std(diff(x_data))
Tsd = 2.1795e-07
Fs = 1/Ts
Fs = 200.0020
[y_data,x_data] = resample(y_data, x_data, round(Fs)); % Resample To Constant Sampling Frequency
Fs = round(Fs);
figure
pspectrum(y_data, Fs, 'spectrogram') % Plot 'pspectrum' 'spectrogram'
colormap(turbo) % Introduced: R2020b
[sp,fp,tp] = pspectrum(y_data,Fs,"spectrogram"); % Return Values, Then Plot
figure
waterfall(fp,tp,sp')
set(gca,XDir="reverse",View=[60 60])
colormap(turbo) % Introduced: R2020b
ylabel("Time (s)")
xlabel("Frequency (Hz)")
figure
waterfall(fp,tp,sp')
set(gca,XDir="reverse",View=[60 60])
Ax = gca;
Ax.ZScale = 'log';
colormap(turbo) % Introduced: R2020b
ylabel("Time (s)")
xlabel("Frequency (Hz)")
It was necessary to change the readmatrix call arguments to import the file correctly. It now works, and gives reasonable results.
There were enough variation in the sampling intervals thaat I considered resampling the data to be worthwhile. That is not necessary, however I believe it improved the result.
.
Rahul
2023-6-4
Thanks a lot Star Strider.
It is working.
But when I change the filename by another filename, the code cannot detect it. It only detects the file OBP3N.txt
Error using fileread
Could not open file OBP4N.txt. No such file or directory.
Error in spectrum_tt1_obp4n (line 6)
filecontents = fileread('OBP4N.txt'); % Check File Format & Contents
Can you please tell me how to generalize this code which can detect all .txt files?
rgds
Star Strider
2023-6-4
My pleasure!
You will probably need to use the fullfile function to create a file name that includes the relevant path information to your other files. Once you have the directory information, simply change the file name argument in fullfile to create the complete path, and then use that as the file name.
For this file here, that might go something like this:
f = @(filename) fullfile('https://www.mathworks.com', 'matlabcentral/answers/uploaded_files/1402564', filename);
file = f('obp3n.txt')
data = readmatrix(file, 'HeaderLines',8) % Skip First 8 Header Lines
That creates:
file = 'https://www.mathworks.com/matlabcentral/answers/uploaded_files/1402564/obp3n.txt'
You will have to experiment with it on your computer, although it would be easier if all your files were added to your MATLAB search path, then the fullfile function (and the anonymous funciton I created here to make its use more straightforward) would not be necessary. See: What Is the MATLAB Search Path? for those details.
If my Answer helped you solve your problem, please Accept it!
.
Rahul
2023-6-5
Hello Star Strider,
The plots are generated and also i can access all .txt files. But I still have some queries regarding the flow of code.
- How are NaN/Inf values generated when actually there are no NaN/Inf values in the real-time machine generated data, which needsto be linearized?
Regards,
Star Strider
2023-6-5
No NaN values are generated.
These lines:
nonfinite_rows = nnz(~any(isfinite(data),2)) % Check To See If Any Rows Have Non-Finite Values
if nonfinite_rows > 0
data(~isfinite(data)) = NaN;
data = fillmissing(data, 'linear');
end
check to see if there are any non-finite values in the file, and if there are, it fills them so that the structure of the data remains intact. It does not linearize anything, it uses linear interpolation (other methods are possible) to fill any non-finite values. That avoids problems.
If you are certain that the data will never have any non-finite values, you can eliminate those lines. The code will run correctly without them.
.
Rahul
2023-6-6
编辑:Rahul
2023-6-6
Hello Star Strider,
Thank you. I got it now. Had some confusion earlier which got cleared now.
Regarding the spectrogram plot, my professor told me that it is wrong. This is because the Time period is in minutes, should have to be in milli-seconds, since the device ran only for few milli-seconds and thus the data generated is also for few milli-seconds.
But, since axes levels are generated by pspectrum itself, how can we change it?
Kindly advice.
Star Strider
2023-6-6
Multiply ‘data(:,1)’ by 1E-3 and it works —
filecontents = fileread('obp3n.txt') % Check File Format & Contents
filecontents =
'DataBase = ***
ShotNo = ****
SignalName = OBP3N
SignalUnit = V
TriggerTime = 0.000000 ms
Period = 0.005000 ms
CreateTime = 2022/08/17 16:27:47
Samples = 100000
0.000000 -0.000408
0.005000 -0.001049
0.010000 -0.000088
0.015000 -0.001049
0.020000 0.003435
0.025000 -0.000408
0.030000 -0.000729
0.035000 -0.003932
0.040000 -0.000408
0.045000 0.002474
0.050000 -0.000088
0.054999 -0.001369
0.059999 -0.002651
0.064999 -0.001369
0.069999 0.000873
0.074999 0.001513
0.079999 0.000552
0.084999 -0.001049
0.089999 -0.001690
0.094999 0.002474
0.099999 0.004076
0.104999 0.000552
0.109999 -0.000408
0.114999 -0.001369
0.119999 0.000873
0.124999 0.000552
0.129999 -0.001369
0.134999 -0.000729
0.139999 0.001193
0.144999 0.001193
0.149999 0.000552
0.154998 0.001834
0.159998 -0.001369
0.164998 -0.001690
0.169998 -0.000408
0.174998 -0.000729
0.179998 0.001193
0.184998 0.000873
0.189998 -0.001369
0.194998 0.000873
0.199998 0.003435
0.204998 0.000232
0.209998 0.001193
0.214998 -0.000408
0.219998 0.000552
0.224998 0.002474
0.229998 0.001513
0.234998 -0.000408
0.239998 -0.000088
0.244998 0.000552
0.249998 0.000552
0.254997 0.001834
0.259997 0.001193
0.264997 -0.000729
0.269997 -0.000408
0.274997 0.000552
0.279997 0.001834
0.284997 -0.001690
0.289997 -0.001369
0.294997 -0.001690
0.299997 -0.001369
0.304997 -0.001369
0.309997 -0.001369
0.314997 0.000873
0.319997 -0.006174
0.324997 -0.001049
0.329997 -0.001049
0.334997 -0.001049
0.339997 -0.001369
0.344997 0.000552
0.349997 -0.003612
0.354996 0.001834
0.359996 0.002795
0.364996 -0.001369
0.369996 -0.001369
0.374996 0.000232
0.379996 -0.000088
0.384996 0.000552
0.389996 -0.002010
0.394996 -0.002010
0.399996 -0.000729
0.404996 0.000873
0.409996 0.002154
0.414996 -0.000729
0.419996 -0.000088
0.424996 -0.001049
0.429996 -0.000729
0.434996 -0.000729
0.439996 0.000552
0.444996 -0.001690
0.449996 -0.002971
0.454995 0.000232
0.459995 0.002154
0.464995 0.000552
0.469995 -0.000088
0.474995 -0.002651
0.479995 -0.001049
0.484995 0.001193
0.489995 0.000232
0.494995 -0.000729
0.499995 -0.000408
0.504995 0.000873
0.509995 0.000873
0.514995 0.002474
0.519995 0.002154
0.524995 0.002474
0.529995 0.000873
0.534995 0.000232
0.539995 0.002474
0.544995 0.000232
0.549995 0.000552
0.554994 -0.001049
0.559994 -0.000729
0.564994 0.001513
0.569994 -0.003612
0.574994 0.004717
0.579994 -0.000088
0.584994 -0.001369
0.589994 0.000873
0.594994 -0.000088
0.599994 0.001193
0.604994 -0.001049
0.609994 -0.000729
0.614994 0.001513
0.619994 -0.001369
0.624994 -0.001049
0.629994 -0.002651
0.634994 0.000873
0.639994 0.001513
0.644994 -0.000729
0.649994 0.000873
0.654993 0.001193
0.659993 -0.001690
0.664993 0.000232
0.669993 0.001834
0.674993 0.004396
0.679993 -0.001690
0.684993 0.000552
0.689993 -0.001369
0.694993 -0.000729
0.699993 0.000873
0.704993 -0.001369
0.709993 -0.000408
0.714993 0.000873
0.719993 0.000232
0.724993 0.000552
0.729993 -0.001369
0.734993 0.000552
0.739993 -0.000729
0.744993 0.001513
0.749993 -0.000408
0.754992 -0.000408
0.759992 0.000232
0.764992 -0.002010
0.769992 -0.001369
0.774992 0.001513
0.779992 -0.000408
0.784992 -0.001690
0.789992 0.000552
0.794992 -0.000408
0.799992 0.001513
0.804992 -0.000088
0.809992 -0.000729
0.814992 -0.003291
0.819992 -0.002971
0.824992 0.003756
0.829992 -0.002651
0.834992 -0.001369
0.839992 -0.000408
0.844992 -0.001049
0.849992 0.001834
0.854991 0.002154
0.859991 0.000873
0.864991 -0.001049
0.869991 0.000232
0.874991 0.001193
0.879991 -0.001049
0.884991 -0.001690
0.889991 0.000873
0.894991 -0.004252
0.899991 0.000552
0.904991 -0.001049
0.909991 0.000873
0.914991 0.000552
0.919991 -0.000729
0.924991 -0.001369
0.929991 0.000232
0.934991 0.001513
0.939991 0.000873
0.944991 -0.000729
0.949991 -0.000088
0.954990 0.000232
0.959990 0.000232
0.964990 -0.001690
0.969990 -0.001690
0.974990 0.001193
0.979990 0.001834
0.984990 -0.001690
0.989990 -0.000729
0.994990 0.000552
0.999990 -0.000408
1.004990 0.001513
1.009990 -0.000408
1.014990 0.001193
1.019990 -0.000729
1.024990 -0.000729
1.029990 0.000552
1.034990 0.002795
1.039990 0.000873
1.044990 -0.000408
1.049990 0.001193
1.054989 -0.000408
1.059989 0.003115
1.064989 -0.001369
1.069989 -0.001690
1.074989 0.000873
1.079989 0.001513
1.084989 0.000232
1.089989 -0.000729
1.094989 0.000552
1.099989 -0.000408
1.104989 0.000232
1.109989 0.000232
1.114989 -0.000088
1.119989 -0.000729
1.124989 -0.000088
1.129989 0.000232
1.134989 0.002795
1.139989 0.001193
1.144989 0.001193
1.149989 0.002154
1.154988 -0.000408
1.159988 -0.000729
1.164988 -0.000729
1.169988 -0.000408
1.174988 0.001193
1.179988 -0.001369
1.184988 0.000873
1.189988 0.000873
1.194988 0.000873
1.199988 -0.000408
1.204988 0.000232
1.209988 -0.000729
1.214988 0.001513
1.219988 0.000873
1.224988 -0.001690
1.229988 0.001193
1.234988 -0.000729
1.239988 -0.000408
1.244988 0.000873
1.249988 -0.000408
1.254987 0.000552
1.259987 -0.001690
1.264987 -0.000408
1.269987 -0.000729
1.274987 -0.000088
1.279987 -0.001369
1.284987 -0.001690
1.289987 0.000873
1.294987 0.001834
1.299987 -0.001049
1.304987 -0.000088
1.309987 -0.001049
1.314987 0.000232
1.319987 -0.001369
1.324987 0.001193
1.329987 -0.001369
1.334987 -0.000729
1.339987 -0.000088
1.344987 -0.000088
1.349987 -0.001049
1.354986 -0.000088
1.359986 -0.000088
1.364986 -0.000729
1.369986 0.000552
1.374986 -0.001049
1.379986 -0.001369
1.384986 -0.001049
1.389986 -0.000088
1.394986 0.000552
1.399986 -0.001690
1.404986 0.001834
1.409986 -0.002651
1.414986 -0.000408
1.419986 -0.000408
1.424986 0.000552
1.429986 -0.002010
1.434986 0.000873
1.439986 0.000552
1.444986 0.000873
1.449986 0.001513
1.454985 -0.001049
1.459985 0.001193
1.464985 0.000232
1.469985 0.001193
1.474985 -0.000088
1.479985 -0.000729
1.484985 0.001513
1.489985 0.001193
1.494985 0.000873
1.499985 0.001193
1.504985 0.000232
1.509985 0.000552
1.514985 -0.000088
1.519985 0.001193
1.524985 -0.001369
1.529985 -0.000729
1.534985 -0.001369
1.539985 0.000873
1.544985 0.000232
1.549985 0.002154
1.554984 0.001193
1.559984 0.000873
1.564984 -0.000408
1.569984 -0.000408
1.574984 0.000232
1.579984 0.000232
1.584984 0.001513
1.589984 0.000552
1.594984 0.003435
1.599984 0.001513
1.604984 0.000552
1.609984 -0.001369
1.614984 -0.001369
1.619984 -0.000408
1.624984 0.000232
1.629984 -0.001690
1.634984 0.001513
1.639984 0.000552
1.644984 -0.000729
1.649984 0.000552
1.654983 -0.000729
1.659983 -0.002010
1.664983 0.000232
1.669983 0.002154
1.674983 -0.000729
1.679983 -0.001049
1.684983 -0.001369
1.689983 0.000232
1.694983 -0.001369
1.699983 -0.000408
1.704983 -0.000408
1.709983 -0.000729
1.714983 -0.001049
1.719983 -0.001369
1.724983 0.001193
1.729983 0.001834
1.734983 -0.000729
1.739983 -0.001369
1.744983 -0.000088
1.749983 0.001193
1.754982 -0.001369
1.759982 -0.001369
1.764982 -0.001369
1.769982 0.001513
1.774982 -0.000088
1.779982 -0.000408
1.784982 -0.000729
1.789982 -0.000408
1.794982 -0.001049
1.799982 -0.000729
1.804982 0.001513
1.809982 -0.001369
1.814982 -0.000088
1.819982 0.001834
1.824982 0.000552
1.829982 0.000232
1.834982 -0.001369
1.839982 -0.000729
1.844982 -0.000729
1.849982 -0.001369
1.854981 -0.001369
1.859981 -0.001690
1.864981 0.000552
1.869981 -0.000088
1.874981 -0.001690
1.879981 -0.000729
1.884981 -0.000408
1.889981 -0.000088
1.894981 0.000873
1.899981 -0.001049
1.904981 0.001834
1.909981 -0.000408
1.914981 0.001193
1.919981 0.001513
1.924981 -0.001690
1.929981 -0.000088
1.934981 0.001513
1.939981 0.000232
1.944981 0.000873
1.949981 0.001834
1.954980 -0.000408
1.959980 -0.001049
1.964980 -0.000729
1.969980 0.000552
1.974980 0.000232
1.979980 -0.002330
1.984980 -0.000408
1.989980 -0.001690
1.994980 -0.001049
1.999980 -0.000088
2.004980 0.000552
2.009980 -0.001369
2.014980 0.002154
2.019980 -0.000408
2.024980 -0.000088
2.029980 0.000873
2.034980 -0.001690
2.039980 0.001193
2.044980 -0.001049
2.049980 0.000873
2.054979 -0.001369
2.059979 0.000552
2.064979 -0.000088
2.069979 -0.000088
2.074979 0.001193
2.079979 0.000552
2.084979 -0.001369
2.089979 -0.001369
2.094979 0.000552
2.099979 0.001193
2.104979 0.000873
2.109979 0.000873
2.114979 -0.002971
2.119979 -0.000088
2.124979 -0.000408
2.129979 0.000873
2.134979 -0.000408
2.139979 0.000232
2.144979 -0.000408
2.149979 0.000232
2.154978 -0.000088
2.159978 -0.001369
2.164978 0.000873
2.169978 -0.001049
2.174978 0.001513
2.179978 -0.000088
2.184978 0.000873
2.189978 0.000873
2.194978 0.000552
2.199978 -0.001049
2.204978 -0.000408
2.209978 0.000552
2.214978 -0.000408
2.219978 0.000552
2.224978 0.000232
2.229978 0.001193
2.234978 -0.001049
2.239978 0.000873
2.244978 0.000232
2.249978 -0.002010
2.254977 -0.000408
2.259977 0.000552
2.264977 -0.001690
2.269977 -0.000729
2.274977 -0.000729
2.279977 -0.000729
2.284977 0.001513
2.289977 -0.003612
2.294977 -0.000408
2.299977 -0.002010
2.304977 0.001834
2.309977 -0.000729
2.314977 -0.000729
2.319977 0.002474
2.324977 -0.000408
2.329977 0.002154
2.334977 -0.000088
2.339977 -0.001369
2.344977 0.000873
2.349977 -0.001369
2.354976 -0.000408
2.359976 -0.001049
2.364976 0.001193
2.369976 0.000232
2.374976 0.000873
2.379976 0.000552
2.384976 0.001193
2.389976 -0.001690
2.394976 -0.001049
2.399976 -0.000729
2.404976 0.001834
2.409976 0.000232
2.414976 0.001513
2.419976 -0.000408
2.424976 -0.000088
2.429976 0.001513
2.434976 0.001193
2.439976 -0.000729
2.444976 0.000552
2.449976 -0.002010
2.454975 0.000873
2.459975 0.001834
2.464975 -0.001049
2.469975 -0.000088
2.474975 -0.000729
2.479975 -0.000729
2.484975 -0.000408
2.489975 -0.000729
2.494975 0.001834
2.499975 -0.000088
2.504975 0.000552
2.509975 -0.000088
2.514975 0.001834
2.519975 -0.000088
2.524975 -0.000408
2.529975 -0.000408
2.534975 -0.002010
2.539975 0.000873
2.544975 -0.000729
2.549975 -0.000408
2.554974 -0.000408
2.559974 -0.000088
2.564974 -0.000408
2.569974 -0.001049
2.574974 0.000873
2.579974 -0.000088
2.584974 0.000873
2.589974 -0.002010
2.594974 0.000873
2.599974 -0.001690
2.604974 -0.000408
2.609974 0.001193
2.614974 -0.001049
2.619974 -0.000088
2.624974 -0.002010
2.629974 -0.001369
2.634974 0.001193
2.639974 0.001513
2.644974 -0.002010
2.649974 -0.000088
2.654973 0.000232
2.659973 -0.001369
2.664973 0.000873
2.669973 -0.000408
2.674973 -0.004252
2.679973 -0.002330
2.684973 0.000232
2.689973 -0.002330
2.694973 -0.000088
2.699973 -0.000408
2.704973 -0.000408
2.709973 0.000552
2.714973 -0.001369
2.719973 0.001513
2.724973 0.000232
2.729973 -0.001049
2.734973 0.001513
2.739973 0.000873
2.744973 -0.000088
2.749973 -0.000408
2.754972 0.001193
2.759972 -0.000408
2.764972 0.000232
2.769972 0.001193
2.774972 -0.000729
2.779972 0.000232
2.784972 -0.000729
2.789972 0.001193
2.794972 0.000873
2.799972 -0.001049
2.804972 -0.003932
2.809972 -0.000088
2.814972 -0.000088
2.819972 0.000873
2.824972 -0.000408
2.829972 -0.002010
2.834972 -0.000408
2.839972 -0.000088
2.844972 -0.001369
2.849972 0.001834
2.854971 -0.000088
2.859971 -0.000408
2.864971 -0.000408
2.869971 0.000552
2.874971 0.001513
2.879971 -0.000088
2.884971 -0.000408
2.889971 -0.000729
2.894971 0.000552
2.899971 -0.000729
2.904971 0.001193
2.909971 -0.001690
2.914971 0.000232
2.919971 -0.000408
2.924971 -0.000408
2.929971 -0.000408
2.934971 -0.000088
2.939971 0.002154
2.944971 0.001513
2.949971 -0.000729
2.954970 0.001834
2.959970 0.001193
2.964970 0.000873
2.969970 0.002154
2.974970 -0.000408
2.979970 -0.001690
2.984970 0.001193
2.989970 -0.000088
2.994970 -0.000088
2.999970 -0.001369
3.004970 0.000552
3.009970 -0.001369
3.014970 0.000873
3.019970 -0.000408
3.024970 0.000552
3.029970 -0.001049
3.034970 -0.000088
3.039970 -0.000729
3.044970 -0.001369
3.049970 0.000552
3.054969 -0.000088
3.059969 -0.000729
3.064969 -0.000088
3.069969 -0.000729
3.074969 -0.001369
3.079969 -0.000088
3.084969 -0.000088
3.089969 -0.000408
3.094969 0.000552
3.099969 -0.000088
3.104969 -0.001049
3.109969 -0.000088
3.114969 -0.000088
3.119969 0.002154
3.124969 -0.001369
3.129969 -0.000408
3.134969 -0.001690
3.139969 -0.000729
3.144969 -0.000729
3.149969 0.000552
3.154968 -0.000408
3.159968 -0.000729
3.164968 0.000552
3.169968 -0.001369
3.174968 -0.001369
3.179968 -0.001049
3.184968 -0.000729
3.189968 -0.000088
3.194968 0.000232
3.199968 0.000873
3.204968 0.000873
3.209968 -0.000088
3.214968 -0.001049
3.219968 -0.000408
3.224968 -0.000088
3.229968 -0.004252
3.234968 0.001513
3.239968 -0.002010
3.244968 -0.002651
3.249968 -0.001049
3.254967 0.000873
3.259967 0.000232
3.264967 0.000873
3.269967 -0.000088
3.274967 0.000873
3.279967 -0.000408
3.284967 -0.000408
3.289967 0.002474
3.294967 -0.001369
3.299967 0.001193
3.304967 0.000873
3.309967 -0.001690
3.314967 -0.000088
3.319967 -0.000729
3.324967 -0.000088
3.329967 -0.000088
3.334967 0.000873
3.339967 0.000552
3.344967 -0.000088
3.349967 -0.001049
3.354966 -0.001369
3.359966 -0.001049
3.364966 0.000873
3.369966 0.004396
3.374966 -0.000408
3.379966 -0.001049
3.384966 -0.000088
3.389966 -0.000088
3.394966 0.002795
3.399966 0.001193
3.404966 -0.001369
3.409966 -0.002651
3.414966 0.000232
3.419966 0.000873
3.424966 -0.000088
3.429966 -0.000729
3.434966 0.000552
3.439966 0.001193
3.444966 0.000873
3.449966 0.001193
3.454965 -0.000088
3.459965 -0.000729
3.464965 -0.000088
3.469965 0.001513
3.474965 -0.000088
3.479965 -0.000729
3.484965 -0.000729
3.489965 -0.000408
3.494965 0.002474
3.499965 -0.000729
3.504965 -0.000729
3.509965 -0.000088
3.514965 -0.000088
3.519965 -0.000088
3.524965 -0.000408
3.529965 -0.000088
3.534965 -0.000088
3.539965 0.000232
3.544965 -0.000408
3.549965 -0.000408
3.554964 -0.000408
3.559964 -0.000088
3.564964 0.000552
3.569964 0.000552
3.574964 0.001513
3.579964 0.000232
3.584964 -0.001049
3.589964 0.000873
3.594964 0.001193
3.599964 -0.001690
3.604964 0.000232
3.609964 -0.001690
3.614964 0.001513
3.619964 -0.002010
3.624964 -0.000088
3.629964 -0.000408
3.634964 -0.002330
3.639964 0.002154
3.644964 -0.000088
3.649964 0.002154
3.654963 -0.000408
3.659963 0.001834
3.664963 -0.001049
3.669963 -0.000729
3.674963 -0.000408
3.679963 0.000552
3.684963 -0.000729
3.689963 -0.001369
3.694963 -0.001369
3.699963 -0.002010
3.704963 0.000552
3.709963 0.000552
3.714963 -0.000408
3.719963 0.000232
3.724963 -0.000729
3.729963 -0.003612
3.734963 0.001513
3.739963 0.000232
3.744963 -0.000408
3.749963 -0.000729
3.754962 0.000232
3.759962 0.001193
3.764962 -0.000088
3.769962 -0.001049
3.774962 0.000232
3.779962 0.001513
3.784962 -0.000729
3.789962 -0.000729
3.794962 0.000873
3.799962 0.000552
3.804962 -0.001049
3.809962 -0.001049
3.814962 -0.000088
3.819962 -0.000088
3.824962 -0.002010
3.829962 -0.000088
3.834962 0.000873
3.839962 0.002474
3.844962 0.000552
3.849962 -0.000729
3.854961 -0.001369
3.859961 0.000873
3.864961 -0.000408
3.869961 0.000552
3.874961 0.001193
3.879961 0.000232
3.884961 0.000232
3.889961 -0.000088
3.894961 0.000232
3.899961 -0.000729
3.904961 -0.001690
3.909961 -0.000088
3.914961 -0.000088
3.919961 0.000552
3.924961 -0.000408
3.929961 0.000552
3.934961 -0.001690
3.939961 0.001834
3.944961 0.000873
3.949961 -0.000729
3.954960 -0.000408
3.959960 0.000552
3.964960 -0.001369
3.969960 0.000873
3.974960 -0.000088
3.979960 -0.001049
3.984960 -0.000088
3.989960 -0.000088
3.994960 -0.000729
3.999960 -0.000729
4.004960 0.000232
4.009960 -0.001690
4.014960 0.000873
4.019960 0.001834
4.024960 -0.002330
4.029960 0.000552
4.034960 -0.001049
4.039960 0.000873
4.044960 0.000552
4.049960 -0.000088
4.054959 -0.000729
4.059959 -0.001369
4.064959 0.000232
4.069959 -0.000088
4.074959 0.000552
4.079959 -0.000408
4.084959 0.002154
4.089959 -0.001049
4.094959 0.000552
4.099959 0.001834
4.104959 -0.000408
4.109959 -0.000408
4.114959 -0.000088
4.119959 -0.001369
4.124959 0.000552
4.129959 -0.000408
4.134959 0.000232
4.139959 -0.000729
4.144959 -0.002330
4.149959 -0.001049
4.154958 -0.000729
4.159958 -0.000729
4.164958 -0.000729
4.169958 0.001513
4.174958 -0.001049
4.179958 -0.000729
4.184958 -0.000088
4.189958 -0.001369
4.194958 0.000552
4.199958 -0.001690
4.204958 -0.001049
4.209958 -0.001049
4.214958 -0.000729
4.219958 -0.000408
4.224958 0.000552
4.229958 0.003435
4.234958 -0.001049
4.239958 0.000552
4.244958 0.001513
4.249958 0.000873
4.254957 0.001834
4.259957 0.000232
4.264957 -0.000088
4.269957 0.000873
4.274957 -0.002651
4.279957 0.002154
4.284957 0.000873
4.289957 0.000232
4.294957 0.000552
4.299957 -0.001369
4.304957 -0.001049
4.309957 0.000232
4.314957 0.000873
4.319957 -0.001690
4.324957 -0.000408
4.329957 -0.001049
4.334957 -0.001049
4.339957 0.001193
4.344957 0.001513
4.349957 -0.000729
4.354956 0.000873
4.359956 0.000232
4.364956 0.000552
4.369956 0.000552
4.374956 -0.000088
4.379956 0.000873
4.384956 -0.000729
4.389956 0.000552
4.394956 -0.000729
4.399956 -0.000408
4.404956 0.000873
4.409956 0.001834
4.414956 -0.000408
4.419956 -0.002010
4.424956 -0.001369
4.429956 -0.000729
4.434956 0.002795
4.439956 -0.000088
4.444956 -0.000408
4.449956 0.000232
4.454955 -0.001049
4.459955 0.001513
4.464955 -0.000088
4.469955 0.000873
4.474955 -0.000088
4.479955 0.000552
4.484955 0.000873
4.489955 -0.000729
4.494955 -0.000729
4.499955 -0.001690
4.504955 -0.000088
4.509955 0.000232
4.514955 -0.001049
4.519955 -0.000729
4.524955 0.000873
4.529955 -0.001369
4.534955 0.000552
4.539955 -0.001690
4.544955 0.000232
4.549955 0.001513
4.554954 0.000232
4.559954 0.000873
4.564954 0.001834
4.569954 -0.001369
4.574954 -0.001049
4.579954 0.000873
4.584954 -0.000408
4.589954 -0.001049
4.594954 0.001193
4.599954 -0.001049
4.604954 -0.001049
4.609954 0.000552
4.614954 0.000873
4.619954 -0.000729
4.624954 0.000552
4.629954 0.001513
4.634954 -0.000408
4.639954 -0.001049
4.644954 -0.002651
4.649954 -0.001690
4.654953 -0.000088
4.659953 -0.001049
4.664953 0.000232
4.669953 -0.000408
4.674953 0.001193
4.679953 -0.000408
4.684953 -0.000088
4.689953 -0.001049
4.694953 -0.001369
4.699953 0.000873
4.704953 0.002795
4.709953 -0.000729
4.714953 0.001513
4.719953 -0.000088
4.724953 0.000232
4.729953 -0.000088
4.734953 0.001193
4.739953 -0.000408
4.744953 0.002154
4.749953 -0.000408
4.754952 0.000873
4.759952 0.000552
4.764952 0.000873
4.769952 -0.000088
4.774952 -0.000729
4.779952 0.000873
4.784952 -0.001049
4.789952 -0.001049
4.794952 0.000552
4.799952 0.000232
4.804952 0.001513
4.809952 -0.000729
4.814952 0.000552
4.819952 -0.001049
4.824952 -0.000088
4.829952 0.000232
4.834952 -0.001690
4.839952 0.000232
4.844952 0.000873
4.849952 -0.001369
4.854951 0.000232
4.859951 -0.001049
4.864951 -0.000408
4.869951 0.003115
4.874951 0.000552
4.879951 -0.001369
4.884951 0.000873
4.889951 0.002474
4.894951 -0.001049
4.899951 0.001193
4.904951 -0.000408
4.909951 -0.000408
4.914951 0.001513
4.919951 -0.000729
4.924951 0.000232
4.929951 -0.000088
4.934951 -0.000088
4.939951 0.001193
4.944951 -0.000088
4.949951 -0.000088
4.954950 0.001193
4.959950 -0.000088
4.964950 -0.001049
4.969950 -0.000088
4.974950 0.000232
4.979950 0.000552
4.984950 -0.000408
4.989950 -0.000088
4.994950 -0.002010
4.999950 0.000552
5.004950 0.000873
5.009950 0.000873
5.014950 -0.000729
5.019950 0.000552
5.024950 -0.000408
5.029950 -0.000088
5.034950 0.000552
5.039950 0.000232
5.044950 0.000232
5.049950 0.000552
5.054949 0.000873
5.059949 -0.000088
5.064949 0.000873
5.069949 -0.000408
5.074949 0.000873
5.079949 0.000232
5.084949 0.000873
5.089949 0.000552
5.094949 -0.000408
5.099949 -0.000088
5.104949 -0.001049
5.109949 -0.001049
5.114949 0.001193
5.119949 -0.000088
5.124949 -0.000088
5.129949 -0.000408
5.134949 -0.000408
5.139949 0.000552
5.144949 -0.000408
5.149949 -0.000729
5.154948 0.000552
5.159948 -0.000408
5.164948 -0.001049
5.169948 0.000873
5.174948 0.000232
5.179948 0.001193
5.184948 -0.000088
5.189948 -0.000088
5.194948 -0.000729
5.199948 0.001193
5.204948 0.000552
5.209948 0.000552
5.214948 0.000552
5.219948 0.001193
5.224948 0.000232
5.229948 0.000552
5.234948 -0.000088
5.239948 0.000232
5.244948 0.000232
5.249948 0.000552
5.254947 0.001513
5.259947 -0.003612
5.264947 0.001513
5.269947 0.000552
5.274947 0.000552
5.279947 -0.001690
5.284947 -0.000088
5.289947 -0.000088
5.294947 -0.001369
5.299947 -0.001369
5.304947 -0.000088
5.309947 -0.000408
5.314947 0.001193
5.319947 -0.001049
5.324947 0.001513
5.329947 -0.002010
5.334947 -0.000729
5.339947 -0.000088
5.344947 0.000552
5.349947 0.000552
5.354946 0.000873
5.359946 0.000873
5.364946 0.001193
5.369946 -0.000729
5.374946 0.001513
5.379946 0.001193
5.384946 -0.000088
5.389946 0.004076
5.394946 -0.001690
5.399946 -0.000408
5.404946 -0.002651
5.409946 0.001193
5.414946 -0.000088
5.419946 0.003115
5.424946 0.001193
5.429946 -0.000408
5.434946 -0.000729
5.439946 0.000552
5.444946 0.000232
5.449946 -0.000729
5.454945 -0.000408
5.459945 -0.000408
5.464945 -0.000729
5.469945 -0.000088
5.474945 -0.000729
5.479945 -0.000088
5.484945 -0.002010
5.489945 0.001193
5.494945 0.000552
5.499945 -0.000088
5.504945 0.000552
5.509945 0.000232
5.514945 -0.000088
5.519945 -0.000408
5.524945 0.000552
5.529945 0.001513
5.534945 0.000552
5.539945 -0.001369
5.544945 -0.000729
5.549945 -0.001049
5.554944 -0.000408
5.559944 0.000552
5.564944 0.003115
5.569944 -0.000408
5.574944 0.000232
5.579944 -0.001369
5.584944 0.000552
5.589944 0.000873
5.594944 -0.000408
5.599944 0.000552
5.604944 -0.001049
5.609944 -0.000088
5.614944 0.000232
5.619944 -0.001049
5.624944 -0.000408
5.629944 -0.002010
5.634944 -0.000088
5.639944 -0.000729
5.644944 -0.000408
5.649944 0.000552
5.654943 0.001193
5.659943 0.002795
5.664943 -0.001690
5.669943 0.000232
5.674943 -0.001049
5.679943 -0.000729
5.684943 -0.001369
5.689943 0.000552
5.694943 -0.000729
5.699943 -0.001369
5.704943 -0.001049
5.709943 0.000873
5.714943 0.002154
5.719943 0.000552
5.724943 0.000873
5.729943 0.000873
5.734943 -0.000088
5.739943 -0.000408
5.744943 0.002154
5.749943 0.000873
5.754942 -0.000088
5.759942 0.002795
5.764942 -0.000088
5.769942 -0.000088
5.774942 0.002795
5.779942 0.001513
5.784942 -0.001049
5.789942 0.000873
5.794942 -0.000088
5.799942 0.000873
5.804942 0.002795
5.809942 -0.001690
5.814942 -0.000729
5.819942 -0.002330
5.824942 -0.000088
5.829942 -0.002651
5.834942 -0.001690
5.839942 -0.001049
5.844942 0.000552
5.849942 -0.000088
5.854941 -0.000408
5.859941 -0.000408
5.864941 0.000552
5.869941 0.000232
5.874941 0.000232
5.879941 -0.001690
5.884941 -0.000729
5.889941 0.001193
5.894941 0.000552
5.899941 0.000552
5.904941 0.000552
5.909941 0.000232
5.914941 0.001193
5.919941 -0.000408
5.924941 -0.001690
5.929941 0.000552
5.934941 -0.000088
5.939941 -0.001049
5.944941 -0.000408
5.949941 -0.000408
5.954940 0.001513
5.959940 0.000232
5.964940 0.000232
5.969940 -0.002010
5.974940 -0.001049
5.979940 0.001193
5.984940 0.000873
5.989940 0.001193
5.994940 -0.001690
5.999940 -0.000088
6.004940 0.000873
6.009940 -0.001369
6.014940 0.001193
6.019940 -0.000408
6.024940 -0.000729
6.029940 -0.001369
6.034940 0.001513
6.039940 -0.000088
6.044940 0.000552
6.049940 0.000552
6.054939 0.001834
6.059939 0.000232
6.064939 -0.001369
6.069939 0.000552
6.074939 -0.000408
6.079939 -0.000408
6.084939 0.000552
6.089939 -0.001369
6.094939 -0.000729
6.099939 0.001834
6.104939 -0.001690
6.109939 0.001513
6.114939 0.000873
6.119939 -0.000408
6.124939 0.001513
6.129939 -0.002010
6.134939 0.000232
6.139939 -0.002010
6.144939 -0.000408
6.149939 -0.001049
6.154938 -0.000729
6.159938 -0.000088
6.164938 0.000552
6.169938 0.000232
6.174938 0.000232
6.179938 0.001193
6.184938 -0.000729
6.189938 -0.001049
6.194938 -0.001690
6.199938 0.000232
6.204938 0.000232
6.209938 0.001513
6.214938 -0.000088
6.219938 -0.000729
6.224938 -0.000729
6.229938 -0.000088
6.234938 0.000232
6.239938 0.001193
6.244938 -0.000088
6.249938 -0.000088
6.254937 -0.000729
6.259937 0.001193
6.264937 0.001513
6.269937 -0.000088
6.274937 -0.000088
6.279937 -0.000729
6.284937 0.004076
6.289937 0.001834
6.294937 0.000232
6.299937 -0.003932
6.304937 -0.000088
6.309937 0.000873
6.314937 0.002795
6.319937 0.001834
6.324937 -0.001049
6.329937 -0.001690
6.334937 -0.000408
6.339937 0.000873
6.344937 0.002154
6.349937 -0.000729
6.354936 -0.000729
6.359936 0.000232
6.364936 -0.000408
6.369936 -0.001049
6.374936 -0.002010
6.379936 -0.001049
6.384936 -0.000729
6.389936 0.000873
6.394936 0.001834
6.399936 -0.000088
6.404936 -0.002010
6.409936 0.001513
6.414936 -0.002010
6.419936 0.000873
6.424936 -0.001369
6.429936 -0.000088
6.434936 -0.000088
6.439936 0.000552
6.444936 0.000552
6.449936 0.000552
6.454935 -0.003932
6.459935 0.000232
6.464935 -0.001049
6.469935 -0.002971
6.474935 0.001513
6.479935 -0.001049
6.484935 -0.000729
6.489935 -0.000408
6.494935 0.001193
6.499935 -0.002010
6.504935 -0.001369
6.509935 -0.000408
6.514935 0.001513
6.519935 -0.000088
6.524935 0.000873
6.529935 0.001513
6.534935 -0.001049
6.539935 -0.001049
6.544935 0.001513
6.549935 -0.001049
6.554934 -0.000088
6.559934 -0.000729
6.564934 -0.002010
6.569934 0.000873
6.574934 0.000232
6.579934 0.002154
6.584934 -0.000729
6.589934 -0.003612
6.594934 -0.000088
6.599934 0.003115
6.604934 0.002154
6.609934 -0.001369
6.614934 -0.000729
6.619934 -0.000729
6.624934 -0.000729
6.629934 0.002474
6.634934 -0.000729
6.639934 -0.001690
6.644934 -0.002010
6.649934 -0.000088
6.654933 0.000232
6.659933 -0.000729
6.664933 -0.000408
6.669933 -0.002651
6.674933 -0.000408
6.679933 -0.000408
6.684933 -0.000088
6.689933 -0.001690
6.694933 -0.001690
6.699933 -0.000088
6.704933 0.001513
6.709933 0.002474
6.714933 0.002474
6.719933 0.000873
6.724933 0.000552
6.729933 -0.000088
6.734933 -0.000408
6.739933 -0.000729
6.744933 -0.001690
6.749933 -0.000729
6.754932 0.001193
6.759932 0.000552
6.764932 -0.000729
6.769932 -0.002330
6.774932 -0.000408
6.779932 -0.000729
6.784932 0.001513
6.789932 0.001513
6.794932 -0.000408
6.799932 -0.002971
6.804932 0.000552
6.809932 -0.001369
6.814932 0.002795
6.819932 0.001513
6.824932 -0.002010
6.829932 -0.000408
6.834932 0.001193
6.839932 0.000552
6.844932 -0.002330
6.849932 -0.000729
6.854931 0.001193
6.859931 0.000873
6.864931 -0.000088
6.869931 0.000552
6.874931 -0.001690
6.879931 -0.000088
6.884931 0.000232
6.889931 -0.001690
6.894931 -0.000088
6.899931 -0.000408
6.904931 -0.000088
6.909931 0.000232
6.914931 -0.000088
6.919931 -0.000088
6.924931 -0.000408
6.929931 -0.001690
6.934931 0.000552
6.939931 -0.000729
6.944931 0.000552
6.949931 -0.001369
6.954930 -0.000408
6.959930 0.000552
6.964930 0.002795
6.969930 -0.002010
6.974930 -0.001690
6.979930 -0.002010
6.984930 0.000232
6.989930 0.001834
6.994930 0.000232
6.999930 -0.001049
7.004930 0.001834
7.009930 0.000552
7.014930 -0.000088
7.019930 -0.001049
7.024930 -0.000408
7.029930 -0.001049
7.034930 -0.001369
7.039930 -0.000088
7.044930 0.000552
7.049930 -0.001369
7.054929 -0.002971
7.059929 -0.000408
7.064929 -0.000408
7.069929 0.002154
7.074929 -0.000088
7.079929 -0.000088
7.084929 -0.000088
7.089929 0.001193
7.094929 -0.000729
7.099929 0.001834
7.104929 0.000873
7.109929 -0.000729
7.114929 0.000873
7.119929 -0.000088
7.124929 0.002154
7.129929 -0.000088
7.134929 -0.000729
7.139929 0.001193
7.144929 0.000232
7.149929 0.002154
7.154928 -0.001369
7.159928 0.000552
7.164928 0.001834
7.169928 0.000552
7.174928 -0.000088
7.179928 0.002474
7.184928 -0.003291
7.189928 -0.000729
7.194928 0.002474
7.199928 0.001834
7.204928 -0.000729
7.209928 0.003115
7.214928 -0.000729
7.219928 -0.000408
7.224928 0.000552
7.229928 -0.000408
7.234928 0.000552
7.239928 -0.001049
7.244928 0.000873
7.249928 -0.000729
7.254927 -0.001049
7.259927 0.000873
7.264927 -0.000088
7.269927 -0.000408
7.274927 0.001193
7.279927 -0.000729
7.284927 -0.000088
7.289927 0.001193
7.294927 0.000873
7.299927 -0.000408
7.304927 0.001193
7.309927 -0.001369
7.314927 -0.000088
7.319927 -0.002010
7.324927 -0.000088
7.329927 -0.000088
7.334927 -0.001369
7.339927 0.001193
7.344927 -0.001049
7.349927 0.000873
7.354926 -0.002010
7.359926 -0.003612
7.364926 0.003115
7.369926 0.000552
7.374926 -0.000088
7.379926 -0.001690
7.384926 0.000232
7.389926 0.002154
7.394926 0.000552
7.399926 -0.000088
7.404926 0.001193
7.409926 0.000232
7.414926 -0.000088
7.419926 0.002154
7.424926 -0.002010
7.429926 -0.000729
7.434926 0.001193
7.439926 -0.000729
7.444926 0.000873
7.449926 -0.001690
7.454925 0.001834
7.459925 -0.001369
7.464925 -0.000729
7.469925 -0.002010
7.474925 -0.000729
7.479925 -0.002971
7.484925 -0.000088
7.489925 0.001513
7.494925 -0.000408
7.499925 0.000232
7.504925 -0.000408
7.509925 -0.000408
7.514925 0.000873
7.519925 -0.001049
7.524925 0.001834
7.529925 -0.002010
7.534925 -0.000729
7.539925 -0.000088
7.544925 0.000232
7.549925 -0.000088
7.554924 0.001193
7.559924 0.000552
7.564924 0.000873
7.569924 0.000552
7.574924 -0.000408
7.579924 0.000873
7.584924 0.001193
7.589924 -0.000408
7.594924 0.002474
7.599924 -0.001049
7.604924 -0.000408
7.609924 -0.000408
7.614924 0.000552
7.619924 0.001834
7.624924 -0.000729
7.629924 0.001193
7.634924 -0.002330
7.639924 -0.001690
7.644924 -0.000408
7.649924 -0.000088
7.654923 -0.001369
7.659923 0.000232
7.664923 -0.001369
7.669923 -0.001369
7.674923 0.003435
7.679923 -0.000729
7.684923 -0.001369
7.689923 -0.000088
7.694923 -0.001369
7.699923 -0.001049
7.704923 0.000552
7.709923 -0.000729
7.714923 -0.000408
7.719923 0.001513
7.724923 0.001193
7.729923 0.001193
7.734923 0.001193
7.739923 -0.001690
7.744923 0.001834
7.749923 -0.000729
7.754922 -0.001049
7.759922 -0.001369
7.764922 -0.001690
7.769922 0.000232
7.774922 -0.000088
7.779922 -0.000408
7.784922 0.001834
7.789922 0.002474
7.794922 -0.001049
7.799922 -0.000408
7.804922 -0.001690
7.809922 -0.000088
7.814922 0.000552
7.819922 -0.000729
7.824922 0.000552
7.829922 -0.000729
7.834922 -0.001369
7.839922 -0.000088
7.844922 0.000232
7.849922 -0.000088
7.854921 -0.000729
7.859921 -0.001369
7.864921 -0.002330
7.869921 -0.001049
7.874921 0.000552
7.879921 0.000232
7.884921 -0.000408
7.889921 -0.000729
7.894921 0.003756
7.899921 -0.000729
7.904921 0.000232
7.909921 -0.001690
7.914921 0.001193
7.919921 0.002154
7.924921 -0.001049
7.929921 -0.002010
7.934921 0.001193
7.939921 0.000552
7.944921 -0.000729
7.949921 0.000552
7.954920 -0.000729
7.959920 0.001513
7.964920 -0.000729
7.969920 0.000873
7.974920 -0.000088
7.979920 0.000232
7.984920 0.000873
7.989920 -0.000088
7.994920 -0.000408
7.999920 0.000552
8.004920 0.001193
8.009920 0.002154
8.014920 -0.000088
8.019920 0.000873
8.024920 0.002474
8.029920 0.001193
8.034920 -0.000729
8.039920 0.000552
8.044920 0.002154
8.049920 0.000232
8.054919 -0.000408
8.059919 0.000873
8.064919 0.001513
8.069919 0.000873
8.074919 0.000873
8.079919 0.000232
8.084919 -0.000408
8.089919 -0.000729
8.094919 0.002795
8.099919 0.000873
8.104919 -0.001369
8.109919 -0.002010
8.114919 -0.001049
8.119919 0.001513
8.124919 0.003435
8.129919 -0.001690
8.134919 0.001193
8.139919 -0.002971
8.144919 -0.001049
8.149919 -0.000729
8.154918 0.002474
8.159918 0.000873
8.164918 -0.000729
8.169918 -0.002010
8.174918 0.000232
8.179918 0.002154
8.184918 -0.000408
8.189918 0.000552
8.194918 -0.001049
8.199918 -0.001369
8.204918 -0.000088
8.209918 -0.001049
8.214918 0.001513
8.219918 -0.000408
8.224918 -0.000088
8.229918 -0.000408
8.234918 0.000552
8.239918 -0.001049
8.244918 -0.001049
8.249918 -0.000088
8.254917 0.001193
8.259917 0.000552
8.264917 -0.000729
8.269917 0.000232
8.274917 0.000552
8.279917 0.000873
8.284917 -0.001369
8.289917 -0.001049
8.294917 -0.001690
8.299917 -0.001369
8.304917 -0.000729
8.309917 -0.000408
8.314917 -0.001690
8.319917 0.000552
8.324917 -0.000088
8.329917 0.000232
8.334917 -0.000408
8.339917 -0.000088
8.344917 0.002154
8.349917 -0.000088
8.354916 -0.000408
8.359916 0.002474
8.364916 0.000873
8.369916 -0.000729
8.374916 0.000873
8.379916 0.000873
8.384916 -0.000729
8.389916 0.000873
8.394916 0.000873
8.399916 -0.000408
8.404916 -0.002330
8.409916 0.002474
8.414916 -0.000408
8.419916 -0.001690
8.424916 0.001193
8.429916 -0.002651
8.434916 -0.000088
8.439916 -0.000088
8.444916 -0.002010
8.449916 -0.000088
8.454915 -0.001690
8.459915 -0.000408
8.464915 -0.000408
8.469915 -0.001690
8.474915 0.000552
8.479915 -0.001049
8.484915 -0.000408
8.489915 -0.000729
8.494915 0.000552
8.499915 -0.000408
8.504915 -0.000408
8.509915 0.000552
8.514915 -0.001369
8.519915 -0.000088
8.524915 -0.000088
8.529915 -0.001369
8.534915 -0.001049
8.539915 -0.000729
8.544915 0.000552
8.549915 -0.000729
8.554914 0.001193
8.559914 -0.000408
8.564914 0.000232
8.569914 -0.001690
8.574914 0.003115
8.579914 -0.001049
8.584914 0.000873
8.589914 -0.000729
8.594914 -0.000729
8.599914 -0.000408
8.604914 -0.002330
8.609914 -0.000408
8.614914 -0.000729
8.619914 -0.000729
8.624914 0.001834
8.629914 -0.002010
8.634914 0.000552
8.639914 -0.001690
8.644914 -0.000088
8.649914 -0.000408
8.654913 0.000552
8.659913 -0.000088
8.664913 -0.001690
8.669913 0.000552
8.674913 0.001193
8.679913 0.001513
8.684913 -0.001049
8.689913 -0.000729
8.694913 -0.000408
8.699913 0.000232
8.704913 -0.000088
8.709913 -0.000729
8.714913 -0.001690
8.719913 -0.000088
8.724913 0.000552
8.729913 -0.002010
8.734913 0.001834
8.739913 0.002154
8.744913 -0.002010
8.749913 0.000232
8.754912 -0.001690
8.759912 -0.000088
8.764912 -0.001049
8.769912 0.000873
8.774912 -0.001369
8.779912 0.000873
8.784912 -0.001369
8.789912 -0.001690
8.794912 0.000873
8.799912 -0.002010
8.804912 -0.000088
8.809912 -0.002010
8.814912 -0.000088
8.819912 -0.004893
8.824912 -0.001369
8.829912 -0.000729
8.834912 0.002154
8.839912 0.002795
8.844912 -0.001049
8.849912 -0.000088
8.854911 -0.000408
8.859911 -0.000088
8.864911 0.000873
8.869911 -0.000729
8.874911 -0.000729
8.879911 -0.000408
8.884911 0.001193
8.889911 0.002154
8.894911 -0.001369
8.899911 -0.000408
8.904911 0.000873
8.909911 -0.000088
8.914911 -0.000088
8.919911 -0.000408
8.924911 -0.001690
8.929911 -0.002971
8.934911 -0.002010
8.939911 0.002795
8.944911 0.001193
8.949911 -0.001690
8.954910 -0.002971
8.959910 -0.000729
8.964910 0.002795
8.969910 0.001513
8.974910 -0.000408
8.979910 -0.000408
8.984910 -0.001690
8.989910 -0.000088
8.994910 0.000552
8.999910 0.000552
9.004910 0.001513
9.009910 -0.002010
9.014910 -0.000729
9.019910 0.002474
9.024910 -0.002010
9.029910 -0.000729
9.034910 -0.001049
9.039910 -0.002651
9.044910 0.001834
9.049910 -0.001690
9.054909 0.001193
9.059909 -0.001049
9.064909 -0.000408
9.069909 0.001513
9.074909 0.001834
9.079909 0.000873
9.084909 -0.001690
9.089909 0.000232
9.094909 -0.000088
9.099909 -0.001369
9.104909 -0.002330
9.109909 -0.003291
9.114909 -0.001690
9.119909 -0.002010
9.124909 0.000873
9.129909 -0.001049
9.134909 0.000552
9.139909 -0.001049
9.144909 -0.000729
9.149909 -0.000729
9.154908 0.001834
9.159908 0.000873
9.164908 0.000552
9.169908 -0.002330
9.174908 0.001834
9.179908 0.001193
9.184908 -0.000408
9.189908 -0.000729
9.194908 -0.000408
9.199908 -0.001369
9.204908 0.001513
9.209908 0.000552
9.214908 0.000873
9.219908 0.000552
9.224908 0.000873
9.229908 0.001834
9.234908 -0.000408
9.239908 -0.001369
9.244908 0.001193
9.249908 -0.000408
9.254907 0.003435
9.259907 -0.001369
9.264907 0.000232
9.269907 -0.001690
9.274907 0.000873
9.279907 -0.000088
9.284907 0.000552
9.289907 -0.001049
9.294907 -0.001690
9.299907 0.004717
9.304907 0.000552
9.309907 0.001834
9.314907 0.000873
9.319907 -0.003291
9.324907 0.000552
9.329907 0.000873
9.334907 -0.001049
9.339907 -0.001690
9.344907 -0.000408
9.349907 -0.001690
9.354906 -0.000088
9.359906 0.001513
9.364906 0.001193
9.369906 0.001193
9.374906 -0.002010
9.379906 0.000552
9.384906 0.001834
9.389906 -0.000088
9.394906 -0.000088
9.399906 -0.002010
9.404906 0.001513
9.409906 0.002795
9.414906 0.001193
9.419906 -0.002010
9.424906 0.000232
9.429906 -0.002651
9.434906 0.001834
9.439906 0.000232
9.444906 0.000552
9.449906 -0.001049
9.454905 -0.001049
9.459905 0.000552
9.464905 0.003115
9.469905 0.000552
9.474905 -0.001690
9.479905 -0.003612
9.484905 -0.000729
9.489905 0.000232
9.494905 -0.001690
9.499905 0.001513
9.504905 0.000873
9.509905 0.001513
9.514905 0.001193
9.519905 0.000873
9.524905 0.000873
9.529905 -0.002651
9.534905 -0.001369
9.539905 0.001193
9.544905 0.000873
9.549905 -0.001049
9.554904 -0.001690
9.559904 -0.000088
9.564904 0.000873
9.569904 -0.000408
9.574904 -0.001049
9.579904 -0.001049
9.584904 -0.001049
9.589904 -0.001049
9.594904 0.003115
9.599904 -0.002010
9.604904 -0.001690
9.609904 -0.000088
9.614904 0.000232
9.619904 -0.000088
9.624904 -0.000729
9.629904 -0.000408
9.634904 -0.000088
9.639904 0.000552
9.644904 0.003435
9.649904 0.001513
9.654903 -0.000729
9.659903 -0.002651
9.664903 -0.000088
9.669903 0.000873
9.674903 0.000232
9.679903 0.000552
9.684903 -0.001690
9.689903 -0.001690
9.694903 0.001193
9.699903 -0.001049
9.704903 0.001834
9.709903 -0.002010
9.714903 -0.001049
9.719903 -0.002010
9.724903 0.000552
9.729903 0.000552
9.734903 0.000873
9.739903 -0.002010
9.744903 -0.000408
9.749903 -0.001690
9.754902 0.001193
9.759902 -0.000408
9.764902 -0.000408
9.769902 0.000232
9.774902 0.000873
9.779902 0.000232
9.784902 -0.001369
9.789902 0.000552
9.794902 -0.000088
9.799902 0.000232
9.804902 -0.000408
9.809902 -0.001049
9.814902 -0.001049
9.819902 -0.001369
9.824902 -0.002651
9.829902 0.000232
9.834902 0.003115
9.839902 0.000873
9.844902 -0.001049
9.849902 0.003115
9.854901 0.003115
9.859901 -0.000088
9.864901 0.002154
9.869901 -0.001049
9.874901 0.000552
9.879901 0.002474
9.884901 0.000873
9.889901 0.001513
9.894901 -0.000729
9.899901 -0.001049
9.904901 0.002154
9.909901 0.000552
9.914901 0.000873
9.919901 -0.000088
9.924901 0.000232
9.929901 -0.001369
9.934901 -0.002010
9.939901 0.000873
9.944901 0.000552
9.949901 -0.000729
9.954900 0.001193
9.959900 0.000232
9.964900 -0.001690
9.969900 0.000873
9.974900 -0.000408
9.979900 -0.003612
9.984900 -0.000408
9.989900 0.000232
9.994900 -0.000088
9.999900 -0.001690
10.004900 0.000873
10.009900 -0.000088
10.014900 -0.001049
10.019900 0.000873
10.024900 -0.000729
10.029900 0.000232
10.034900 0.001513
10.039900 0.000552
10.044900 0.000552
10.049900 -0.000408
10.054899 -0.000408
10.059899 0.000232
10.064899 0.000873
10.069899 -0.000408
10.074899 -0.000729
10.079899 -0.002010
10.084899 -0.000088
10.089899 -0.000729
10.094899 0.001193
10.099899 -0.001049
10.104899 0.000552
10.109899 0.000232
10.114899 -0.001369
10.119899 0.000552
10.124899 0.000552
10.129899 0.000232
10.134899 -0.001049
10.139899 -0.000729
10.144899 0.000552
10.149899 0.001193
10.154898 -0.000088
10.159898 -0.001049
10.164898 -0.002010
10.169898 0.000552
10.174898 0.000232
10.179898 0.000873
10.184898 -0.001369
10.189898 0.001193
10.194898 0.000232
10.199898 -0.000729
10.204898 0.003756
10.209898 -0.001690
10.214898 0.000873
10.219898 0.000552
10.224898 -0.000729
10.229898 -0.000729
10.234898 -0.001049
10.239898 0.002474
10.244898 0.000232
10.249898 0.001513
10.254897 0.000873
10.259897 -0.001369
10.264897 0.000552
10.269897 0.001834
10.274897 0.000873
10.279897 0.000873
10.284897 0.000873
10.289897 0.000552
10.294897 -0.000729
10.299897 -0.000088
10.304897 -0.000088
10.309897 0.000552
10.314897 0.000232
10.319897 0.000873
10.324897 -0.000408
10.329897 0.002154
10.334897 -0.001369
10.339897 -0.001369
10.344897 0.002474
10.349897 -0.001049
10.354896 -0.000088
10.359896 0.000552
10.364896 -0.001369
10.369896 -0.001369
10.374896 0.001513
10.379896 0.001193
10.384896 -0.001690
10.389896 -0.000088
10.394896 0.001513
10.399896 0.000232
10.404896 0.000232
10.409896 -0.000088
10.414896 0.001193
10.419896 -0.000729
10.424896 0.003756
10.429896 -0.000088
10.434896 0.001513
10.439896 -0.002651
10.444896 -0.001369
10.449896 -0.000729
10.454895 0.000552
10.459895 -0.000088
10.464895 0.001513
10.469895 0.000552
10.474895 -0.001369
10.479895 -0.002651
10.484895 0.000552
10.489895 -0.001369
10.494895 0.001193
10.499895 -0.000729
10.504895 0.001834
10.509895 -0.000408
10.514895 -0.000088
10.519895 -0.001690
10.524895 0.000552
10.529895 0.000873
10.534895 -0.000088
10.539895 -0.000408
10.544895 -0.000088
10.549895 0.001513
10.554894 -0.002010
10.559894 0.000873
10.564894 0.000873
10.569894 -0.001369
10.574894 -0.000729
10.579894 -0.000088
10.584894 0.000873
10.589894 -0.002330
10.594894 -0.000408
10.599894 -0.000088
10.604894 0.001193
10.609894 -0.000088
10.614894 0.000552
10.619894 0.003756
10.624894 0.002474
10.629894 0.002154
10.634894 -0.003291
10.639894 0.002795
10.644894 0.000232
10.649894 -0.000408
10.654893 -0.000729
10.659893 0.000873
10.664893 -0.000408
10.669893 -0.000729
10.674893 -0.001369
10.679893 0.000552
10.684893 -0.000408
10.689893 -0.000729
10.694893 0.000552
10.699893 -0.001049
10.704893 0.001513
10.709893 0.000873
10.714893 -0.002971
10.719893 -0.000088
10.724893 -0.000408
10.729893 0.000232
10.734893 -0.000729
10.739893 -0.002010
10.744893 0.001193
10.749893 0.000552
10.754892 0.001193
10.759892 -0.000408
10.764892 0.000232
10.769892 0.000552
10.774892 0.000232
10.779892 0.000873
10.784892 -0.001690
10.789892 0.001834
10.794892 -0.000088
10.799892 -0.000088
10.804892 0.001834
10.809892 0.002154
10.814892 0.002154
10.819892 0.002154
10.824892 0.001193
10.829892 -0.000408
10.834892 -0.000408
10.839892 -0.001049
10.844892 0.003756
10.849892 0.000232
10.854891 -0.001369
10.859891 -0.000729
10.864891 -0.000729
10.869891 -0.000088
10.874891 -0.001049
10.879891 -0.000088
10.884891 -0.002010
10.889891 -0.002651
10.894891 -0.000408
10.899891 0.001513
10.904891 -0.000408
10.909891 -0.001690
10.914891 -0.001690
10.919891 -0.000729
10.924891 0.003115
10.929891 0.001834
10.934891 -0.002010
10.939891 -0.001690
10.944891 0.001193
10.949891 0.001513
10.954890 -0.000408
10.959890 -0.003291
10.964890 -0.000408
10.969890 0.001193
10.974890 -0.001690
10.979890 -0.000729
10.984890 0.001193
10.989890 -0.000088
10.994890 -0.000408
10.999890 -0.002010
11.004890 0.001193
11.009890 -0.001690
11.014890 0.000873
11.019890 -0.001049
11.024890 -0.000088
11.029890 -0.000408
11.034890 0.001834
11.039890 -0.000408
11.044890 -0.000088
11.049890 -0.005534
11.054889 -0.000088
11.059889 0.000873
11.064889 0.000873
11.069889 -0.000729
11.074889 -0.002010
11.079889 -0.001049
11.084889 0.000232
11.089889 0.000873
11.094889 -0.000088
11.099889 -0.002330
11.104889 0.000552
11.109889 -0.001690
11.114889 0.000873
11.119889 0.001834
11.124889 -0.000408
11.129889 -0.000088
11.134889 -0.001369
11.139889 -0.000088
11.144889 -0.000408
11.149889 0.002154
11.154888 0.000552
11.159888 0.000552
11.164888 -0.002010
11.169888 0.000873
11.174888 -0.000408
11.179888 -0.001049
11.184888 -0.000088
11.189888 -0.002971
11.194888 -0.000408
11.199888 -0.000088
11.204888 0.001193
11.209888 -0.000729
11.214888 -0.000088
11.219888 0.001513
11.224888 -0.002651
11.229888 -0.000408
11.234888 0.004076
11.239888 0.000232
11.244888 -0.001369
11.249888 -0.000088
11.254887 0.001513
11.259887 0.000873
11.264887 -0.000408
11.269887 0.002154
11.274887 -0.001049
11.279887 -0.002330
11.284887 -0.000729
11.289887 -0.000408
11.294887 0.001193
11.299887 -0.002010
11.304887 -0.000088
11.309887 0.000232
11.314887 0.000552
11.319887 -0.002651
11.324887 -0.000729
11.329887 -0.001369
11.334887 -0.000088
11.339887 -0.001690
11.344887 0.001513
11.349887 -0.000088
11.354886 -0.001690
11.359886 -0.000408
11.364886 -0.001049
11.369886 -0.000408
11.374886 0.002474
11.379886 -0.000408
11.384886 0.001834
11.389886 -0.000408
11.394886 -0.000088
11.399886 -0.000408
11.404886 -0.002010
11.409886 -0.001690
11.414886 0.000552
11.419886 0.000552
11.424886 -0.000088
11.429886 -0.001369
11.434886 0.002154
11.439886 0.001193
11.444886 0.001193
11.449886 -0.002010
11.454885 0.000552
11.459885 -0.001369
11.464885 -0.000088
11.469885 -0.000729
11.474885 -0.002010
11.479885 -0.001690
11.484885 -0.002651
11.489885 0.000873
11.494885 -0.002971
11.499885 -0.000729
11.504885 -0.002010
11.509885 -0.002971
11.514885 -0.001369
11.519885 -0.000088
11.524885 0.000552
11.529885 0.001193
11.534885 -0.001369
11.539885 0.000873
11.544885 0.002154
11.549885 0.000232
11.554884 0.003435
11.559884 0.001513
11.564884 -0.002010
11.569884 0.000552
11.574884 0.001834
11.579884 0.001834
11.584884 -0.002651
11.589884 0.001193
11.594884 -0.002010
11.599884 0.002795
11.604884 0.000873
11.609884 0.001513
11.614884 0.002154
11.619884 -0.000088
11.624884 -0.000088
11.629884 -0.003612
11.634884 0.000552
11.639884 0.000873
11.644884 0.000552
11.649884 -0.000088
11.654883 0.000552
11.659883 -0.000408
11.664883 0.001193
11.669883 0.001193
11.674883 0.001193
11.679883 -0.000088
11.684883 0.000232
11.689883 -0.001690
11.694883 -0.001049
11.699883 0.001193
11.704883 0.001193
11.709883 0.000552
11.714883 0.001513
11.719883 -0.001690
11.724883 -0.000088
11.729883 -0.001690
11.734883 -0.000729
11.739883 -0.001049
11.744883 0.001513
11.749883 -0.000729
11.754882 0.001193
11.759882 -0.000088
11.764882 0.002154
11.769882 -0.000088
11.774882 0.000873
11.779882 0.001193
11.784882 0.000552
11.789882 -0.001049
11.794882 -0.000729
11.799882 0.000232
11.804882 -0.000729
11.809882 -0.001369
11.814882 0.002795
11.819882 -0.002971
11.824882 -0.000088
11.829882 -0.000088
11.834882 -0.002330
11.839882 0.000873
11.844882 0.000552
11.849882 0.001513
11.854881 -0.000408
11.859881 0.000552
11.864881 0.001193
11.869881 0.000552
11.874881 0.002154
11.879881 -0.000088
11.884881 0.002474
11.889881 0.000873
11.894881 -0.000088
11.899881 0.000552
11.904881 -0.002010
11.909881 -0.000088
11.914881 -0.000088
11.919881 0.003435
11.924881 -0.000408
11.929881 -0.000729
11.934881 -0.002010
11.939881 0.000552
11.944881 -0.000088
11.949881 -0.000729
11.954880 -0.000408
11.959880 0.001193
11.964880 -0.001049
11.969880 0.001193
11.974880 -0.003291
11.979880 -0.000729
11.984880 -0.000408
11.989880 0.001513
11.994880 -0.000408
11.999880 -0.000088
12.004880 0.001193
12.009880 -0.001049
12.014880 0.000552
12.019880 -0.000408
12.024880 -0.001369
12.029880 0.000873
12.034880 0.000552
12.039880 -0.000729
12.044880 0.000232
12.049880 -0.000408
12.054879 0.001513
12.059879 0.000232
12.064879 0.000873
12.069879 0.001513
12.074879 0.000552
12.079879 -0.000088
data = readmatrix('obp3n.txt', 'HeaderLines',8) % Skip First 8 Header Lines
data = 100000×2
0 -0.0004
0.0050 -0.0010
0.0100 -0.0001
0.0150 -0.0010
0.0200 0.0034
0.0250 -0.0004
0.0300 -0.0007
0.0350 -0.0039
0.0400 -0.0004
0.0450 0.0025
data(:,1) = data(:,1) * 1E-3;
nonfinite_rows = nnz(~any(isfinite(data),2)) % Check To See If Any Rows Have Non-Finite Values
nonfinite_rows = 0
if nonfinite_rows > 0
data(~isfinite(data)) = NaN;
data = fillmissing(data, 'linear');
end
x_data = data(:,1);
y_data = data(:,2);
Ts = mean(diff(x_data));
Tsd = std(diff(x_data));
Fs = 1/Ts
Fs = 2.0000e+05
[y_data,x_data] = resample(y_data, x_data, round(Fs)); % Resample To Constant Sampling Frequency
Fs = round(Fs);
figure
pspectrum(y_data, Fs, 'spectrogram') % Plot 'pspectrum' 'spectrogram'
colormap(turbo) % Introduced: R2020b
[sp,fp,tp] = pspectrum(y_data,Fs,"spectrogram"); % Return Values, Then Plot
figure
waterfall(fp,tp,sp')
set(gca,XDir="reverse",View=[60 60])
colormap(turbo) % Introduced: R2020b
ylabel("Time (s)")
xlabel("Frequency (Hz)")
figure
waterfall(fp,tp,sp')
set(gca,XDir="reverse",View=[60 60])
Ax = gca;
Ax.ZScale = 'log';
colormap(turbo) % Introduced: R2020b
ylabel("Time (s)")
xlabel("Frequency (Hz)")
Fixed!
.
Rahul
2023-6-7
Hi Star Strider,
I want to study the activity at a particular frequency of the spectrogram. Hence, I'm trying to design a band pass filter. I want to study say for example at frequency band 35-50 KHz.
I'm unable to implement filtfilt.
Can you plz advice on the code below:
% Design filter for studying activities at particular frequency
w_p = 50/Fs/2;
w_s = 150/Fs/2;
[n,w_n]=buttord(w_p,w_s,3,60);
[z,p,k]=butter(n,w_n);
sos=zp2sos(z,p,k);
figure
freqz(sos,512,'whole')
Star Strider
2023-6-7
I generally prefer to use elliptic filters. If you have R2018a or later, you can just do this:
y_data_filt = bandpass(y_data, [35 50], Fs, 'ImpulseResponse','iir');
Otherwise, this works —
data = readmatrix('obp3n.txt', 'HeaderLines',8) % Skip First 8 Header Lines
data = 100000×2
0 -0.0004
0.0050 -0.0010
0.0100 -0.0001
0.0150 -0.0010
0.0200 0.0034
0.0250 -0.0004
0.0300 -0.0007
0.0350 -0.0039
0.0400 -0.0004
0.0450 0.0025
data(:,1) = data(:,1) * 1E-3; % Time In Milliseconds
x_data = data(:,1);
y_data = data(:,2);
Ts = mean(diff(x_data));
Tsd = std(diff(x_data));
Fs = 1/Ts
Fs = 2.0000e+05
[y_data,x_data] = resample(y_data, x_data, round(Fs)); % Resample To Constant Sampling Frequency
Fs = round(Fs);
Fn = Fs/2;
Wp = [35 50]/Fn; % Passband Frequency (Normalised)
Ws = [0.98 1.02].*Wp; % Stopband Frequency (Normalised)
Rp = 1; % Passband Ripple
Rs = 60; % Passband Ripple (Attenuation)
[n,Wp] = ellipord(Wp,Ws,Rp,Rs); % Elliptic Order Calculation
[z,p,k] = ellip(n,Rp,Rs,Wp); % Elliptic Filter Design: Zero-Pole-Gain
[sos,g] = zp2sos(z,p,k); % Second-Order Section For Stability
figure
freqz(sos, 2^18, Fs) % Filter Bode Plot
set(subplot(2,1,1), 'XLim',[0 100]) % OPTTIONAL
set(subplot(2,1,2), 'XLim',[0 100]) % OPTTIONAL
y_data_filt = filtfilt(sos, g, y_data); % Filter Signal
figure
plot(x_data, y_data_filt)
grid
xlabel('Time')
ylabel('Amplitude')
figure
pspectrum(y_data_filt, Fs, 'spectrogram') % Plot 'pspectrum' 'spectrogram'
colormap(turbo) % Introduced: R2020b
[sp,fp,tp] = pspectrum(y_data,Fs,"spectrogram"); % Return Values, Then Plot
figure
waterfall(fp,tp,sp')
set(gca,XDir="reverse",View=[60 60])
colormap(turbo) % Introduced: R2020b
ylabel("Time (s)")
xlabel("Frequency (Hz)")
figure
waterfall(fp,tp,sp')
set(gca,XDir="reverse",View=[60 60])
Ax = gca;
Ax.ZScale = 'log';
colormap(turbo) % Introduced: R2020b
ylabel("Time (s)")
xlabel("Frequency (Hz)")
.
Rahul
2023-6-8
Hi Star Strider,
The profiles are a bit confusing. The spectrogram shows frequencies in KHz (before implementing filter).
I wished to study the frequency 35-50 KHz range and after implementing the filter the profiles shows frequency in Hz and not in KHz.
Can you kindly elaborate tis to me plz.
Star Strider
2023-6-8
That relates to the way the axes are scaled. They all go from 0 Hz (D-C) to Hz, equivalent to kHz. The first plot is scaled in kHz and goes from 0 to kHz, the others are scaled in Hz and go from 0 Hz to (equivalent to ) Hz.
So, they are all the same.
Rahul
2023-6-11
Hi Star Strider,
I tried to do some modification in the code as per my requirement, but got stuck.
I request you to kindly advice me how to get the two subplots executed along with figure (5)
figure(4)
[h,w]=freqz(sos, 10,Fs) ; % Filter Bode Plot
set(subplot(2,1,1), 'XLim',[0 100]) % OPTTIONAL
set(subplot(2,1,2), 'XLim',[0 100]) % OPTTIONAL
figure(5)
ax = gca;
ax.YLim=[0,100];
ax.XTick = 0:10:100;
plot(w, pi/180*unwrap(angle(h)));
xlabel('Frequency (Hz)');
ylabel('Phase (rad)');
Star Strider
2023-6-11
The freqz function with outputs will not also plot the Bode plot, so two separate calls to it are necessary here, however the second argument to it in this code limits the fft to only length 10 (2^3.3), while I usually use something on the order of 2^18.
The result is that will not be anywhere close to being an accurate Bode plot for this filter —
Fs = 2.0000e+05
Fs = 200000
Fn = Fs/2;
Wp = [35 50]/Fn; % Passband Frequency (Normalised)
Ws = [0.98 1.02].*Wp; % Stopband Frequency (Normalised)
Rp = 1; % Passband Ripple
Rs = 60; % Passband Ripple (Attenuation)
[n,Wp] = ellipord(Wp,Ws,Rp,Rs); % Elliptic Order Calculation
[z,p,k] = ellip(n,Rp,Rs,Wp); % Elliptic Filter Design: Zero-Pole-Gain
[sos,g] = zp2sos(z,p,k); % Second-Order Section For Stability
figure(4)
freqz(sos, 10,Fs) ; % Filter Bode Plot
% set(subplot(2,1,1), 'XLim',[0 100]) % OPTTIONAL
% set(subplot(2,1,2), 'XLim',[0 100]) % OPTTIONAL
[h,w]=freqz(sos, 10,Fs) ;
figure(5)
ax = gca;
ax.YLim=[0,100];
ax.XTick = 0:10:100;
plot(w, pi/180*unwrap(angle(h)));
xlabel('Frequency (Hz)');
ylabel('Phase (rad)');
I am not certain what you want to do here.
.
Rahul
2023-6-11
Hi,
I want to plot the phase (in radians) versus frequency plot (previously phase was in degree).
So just conversion of phase from degree to radians.
regards
Star Strider
2023-6-11
Rahul
2023-6-13
移动:Star Strider
2023-6-13
Hi,
I need suggestion again.
I want to access 12 no. of txt files one after another and plot them as a single plot.
I wrote the code but fopen is not working and fid = -1.
Kindly advice.
figure(5)
for i=1:12
myfilename = sprintf('OBP%i.txt',i);
f = @(filename) fullfile('D:\Data_tokamak\TT1 data',filename );
fid = fopen(filename,'r');
Mydata = readmatrix(file, 'HeaderLines',8);
fid = fclose(fid);
ax = gca;
ax.YLim=[0,100];
ax.XTick = [0:10:100];
plot(w, pi/180*(angle(h)));
grid on;
xlabel('Frequency (Hz)');
ylabel('Phase (rad)');
hold on
end
with regards
Star Strider
2023-6-13
I have no idea how you generate ‘h’ (that previously was created by freqz) or how it relates to ‘Mydata’ that your code never uses or refers to after creating it. If you are supposed to filter it after reading it, you need to add a filtfilt call using the previously-created filter, and then do sometihing with the fitlered data.
There is no freqz call anywhere. If it is called before the loop, the code plots the same information 12 times.
I also changed the file name format descriptor to match the previous file names. MATLAB is case-sensitive, and the file names must match exactly.
I made some minor changes, however I cannot fix it because I do not know what you want to do with ‘Mydata’ or even what it is —
figure(5)
hold on
for i=1:12
myfilename = sprintf('obp%in.txt',i);
file = fullfile('D:\Data_tokamak\TT1 data',myfilename );
Mydata = readmatrix(file, 'HeaderLines',8);
ax = gca;
ax.YLim=[0,100];
ax.XTick = [0:10:100];
plot(w, pi/180*(angle(h)));
grid on;
xlabel('Frequency (Hz)');
ylabel('Phase (rad)');
end
hold off
This code honestly does not make sense to me.
.
Rahul
2023-6-13
Hi Star Strider,
Yes, I understood my mistake now.
Moreover I sent you only part of the code.
I will have to make a separate function may be. I will do it.
Thanks a lot for your support.
rgds
Rahul
2023-6-13
Hi,
I can't open the txt file, (fopen not working), since fid is showing -1.
This is probably because I'm not using the syntax
file = f('OBP1N.txt');
that was used in the previous program.
Since I want to access multiple files how can I do this,I want to know from you please, but only one phase versus frequency plot (phase in radians).
% Phase (in radians) versus Frequency (in Hz) plot for
% multiple input data from txt files
numfiles = 12;
mydata = cell(1, numfiles);
for k = 1:numfiles
myfilename = sprintf('OBP%d.txt', k);
f = @(filename)fullfile('D:\Data_tokamak\TT1 data',filename);
fid = fopen(myfilename,'r');
%mydata{k} = importdata(myfilename);
mydata(k) = readmatrix(myfilename, 'HeaderLines',8);
mydata(:,1) = mydata(:,1) * 1E-3;
fid = fclose(fid);
end
for k=1:numfiles
nonfinite_rows(k) = nnz(~any(isfinite(mydata),2));
if nonfinite_rows(k) > 0
data(~isfinite(mydata)) = NaN;
data = fillmissing(mydata, 'linear');
end
x_data(k) = mydata(:,1);
y_data(k) = mydata(:,2);
Ts = mean(diff(x_data));
Tsd = std(diff(x_data));
Fs = 1/Ts;
end
pspectrum(y_data, Fs, 'spectrogram');
colormap(turbo);
[sp,fp,tp] = pspectrum(y_data,Fs,"spectrogram");
wp = [39.8 40.2]/Fn;
ws = [0.98 1.02].*wp;
Rp = 1;
Rs = 60;
[n,wp] = ellipord(wp,ws,Rp,Rs);
[z,p,k] = ellip(n,Rp,Rs,wp);
[sos,g] = zp2sos(z,p,k);
freqz(sos, 2^18,Fs) ;
set(subplot(2,1,1), 'XLim',[0 100])
set(subplot(2,1,2), 'XLim',[0 100])
[h w] = freqz(sos, 10,Fs) ;
%Implementation of filtfilt function without plotting
y_data_filt = filtfilt(sos, g, y_data);
% Plot for phase (in radians) versus frequency (in hz)
for i=1:numfiles
ax = gca;
ax.YLim=[0,100];
ax.XTick = [0:10:100];
plot(w, pi/180*(angle(h)));
grid on;
xlabel('Frequency (Hz)');
ylabel('Phase (rad)');
hold on
end
Sincerely
Star Strider
2023-6-13
编辑:Star Strider
2023-6-13
There are several problems with the first loop.
The first loop should be something like this —
f = @(filename)fullfile('D:\Data_tokamak\TT1 data',filename);
for k = 1:numfiles
myfilename = sprintf('obp%dn.txt', k);
mydata{k} = readmatrix(f(myfilename), 'HeaderLines',8);
mydata{k}(:,1) = mydata(:,1) * 1E-3;
end
That looks like it should work.
Using this construction, this needs to be changed to:
x_data(:,k) = mydata{k}(:,1);
y_data(:,k) = mydata{k}(:,2);
The rest of the code seems to be what has already been demonstrated to work, at least with the provided file.
I cannot test it beyond wiat I already wrote.
.
Rahul
2023-6-14
Hi Star Strider,
Apology for disturbing you time and again, but it still does not run but gives error:
Too many input arguments
I thought may be because mydata contains only one output variable (k) whereas the input contains two columns of data (time and amplitude). I hence introduced two variables in output k1 and k2. But still it didn't solve the problem.
Please advice.
Star Strider
2023-6-14
‘Too many input arguments’
What line is throwing that error? (Please copy and paste all the red text from the Command Window.)
What arguments are you passing to that function call?
It is UCT-6 here (22:25) so I will deal with this in the morning.
.
Rahul
2023-6-14
Here it is as below:
Error using phase_plot_rad_tt1>@(filename)fullfile('D:\Data_tokamak\TT1 data',filename)
Too many input arguments.
Error in phase_plot_rad_tt1 (line 10)
mydata{k} = readmatrix(f(myfilename, 'HeaderLines',8));
with regards
Star Strider
2023-6-14
O.K.
mydata{k} = readmatrix(f(myfilename), 'HeaderLines',8);
That should work.
Rahul
2023-6-14
Hi,
Error using readmatrix
Unable to find or open 'D:\Data_tokamak\TT1 data\OBP1.txt'. Check the path and filename or file permissions.
Error in phase_plot_rad_tt1 (line 10)
mydata{k} = readmatrix(f(myfilename), 'HeaderLines',8);
regards
Star Strider
2023-6-14
Rahul
2023-6-14
Hi,
How to solve this error.
Unable to perform assignment because the size of the left side is 100000-by-1 and the size of the right side is
99999-by-1.
Error in phase_plot_rad_tt1 (line 30)
[y_data(:,k),x_data(:,k)] = resample(y_data(:,k), x_data(:,k), round(Fs));
plz advice.
with rgds
Star Strider
2023-6-14
[y_datac{:,k},x_datac{:,k}] = resample(y_data(:,k), x_data(:,k), round(Fs));
That should work. The only complication is that you will have to use cell array indexing on ‘x_datac’ and ‘y_datac’ after that. It is not possible to store them in ‘x_data’ and ‘y_data’ because those were not originally cell arrays, and the sizes do not match, so they cannot be stored as numeric arrays.
.
Star Strider
2023-6-14
The problem seems to be whatever ‘k’ is, since it works here for different values of ‘k’ —
Fs = 200000;
Fn = Fs/2;
Wp = [35 50]/Fn; % Passband Frequency (Normalised)
Ws = [0.98 1.02].*Wp; % Stopband Frequency (Normalised)
Rp = 1; % Passband Ripple
Rs = 60; % Passband Ripple (Attenuation)
[n,Wp] = ellipord(Wp,Ws,Rp,Rs); % Elliptic Order Calculation
[z,p,k] = ellip(n,Rp,Rs,Wp); % Elliptic Filter Design: Zero-Pole-Gain
[sos,g] = zp2sos(z,p,k); % Second-Order Section For Stability
k = 5;
y_datac{:,k} = randn(50,1);
y_data_filt = filtfilt(sos, g, y_datac{:,k});
CheckSize = size(y_data_filt)
CheckSize = 1×2
50 1
So ‘k’ must be an integer greater than zero. (The code does not use logical indexing, so we can rule that out.)
.
Rahul
2023-6-15
Hi,
The value of k in my code is never initialized to zero, its going from 1 to numfiles(=12).
I did as you adviced. But none of the plots are getting executed.
Even the spectrogram is also not getting executed.
kindly advice.
with regards.
% Phase (in radians) versus Frequency (in Hz) plot for
% multiple input data from txt files
numfiles = 12;
x_datac{:,1}=[];x_datac{:,2}=[];x_datac{:,3}=[];x_datac{:,4}=[];x_datac{:,5}=[];x_datac{:,6}=[];x_datac{:,7}=[];
x_datac{:,8}=[];x_datac{:,9}=[];x_datac{:,10}=[];x_datac{:,11}=[];x_datac{:,12}=[];
y_datac{:,1}=[];y_datac{:,2}=[];y_datac{:,3}=[];y_datac{:,4}=[];y_datac{:,5}=[];y_datac{:,6}=[];y_datac{:,7}=[];
y_datac{:,8}=[];y_datac{:,9}=[];y_datac{:,10}=[];y_datac{:,11}=[];y_datac{:,12}=[];
mydata = cell(1, numfiles);
f = @(filename)fullfile('D:\Data_tokamak\TT1 data',filename);
for k = 1:numfiles
myfilename = sprintf('OBP%dN.txt', k);
mydata{k} = readmatrix(f(myfilename), 'HeaderLines',8);
mydata{k}(:,1) = mydata{k}(:,1) * 1E-3;
end
for k=1:numfiles
nonfinite_rows{k} = nnz(~any(isfinite(mydata{k}),2)); % Check To See If Any Rows Have Non-Finite Values
if nonfinite_rows{k} > 0
data(~isfinite(mydata{k})) = NaN;
data = fillmissing(mydata{k}, 'linear');
end
x_data(:,k) = mydata{k}(:,1);
y_data(:,k) = mydata{k}(:,2);
Ts = mean(diff(x_data(:,k)));
Tsd = std(diff(x_data(:,k)));
Fs = 1/Ts;
[y_datac{:,k},x_datac{:,k}] = resample(y_data(:,k), x_data(:,k), round(Fs)); % Resample To Constant Sampling Frequency and using cell arrays
end
Fs = round(Fs);
Fn = Fs/2;
figure(1)
for k=1:numfiles
pspectrum(y_datac{:,k}, Fs, 'spectrogram'); % Plot 'pspectrum' 'spectrogram'
colormap(turbo); % Introduced: R2020b
[sp,fp,tp] = pspectrum(y_datac{:,k},Fs,"spectrogram"); % Return Values, Then Plot
end
wp = [39.8 40.2]/Fn; % Passband frequency Normalized
ws = [0.98 1.02].*wp; % Stopband frequency Normalized
Rp = 1; % Passband Ripple
Rs = 60; % Passband Ripple (Attenuation)
[n,wp] = ellipord(wp,ws,Rp,Rs); % Elliptic order calculation
[z,p,k] = ellip(n,Rp,Rs,wp); % Elliptic filter design: Zero-Pole-Gain
[sos,g] = zp2sos(z,p,k); % Second order selection for stability
figure(2)
freqz(sos, 2^18,Fs) ; % Filter Bode Plot
set(subplot(2,1,1), 'XLim',[0 100]) % OPTTIONAL
set(subplot(2,1,2), 'XLim',[0 100]) % OPTTIONAL
[h w] = freqz(sos, 10,Fs) ;
%Implementation of filtfilt function without plotting
k=12;
y_data_filt = filtfilt(sos, g, y_datac{:,k}); % Filter Signal
% Plot for phase (in radians) versus frequency (in hz)
for k=1:numfiles
ax = gca;
ax.YLim=[0,100];
ax.XTick = [0:10:100];
plot(w, pi/180*(angle(h)));
grid on;
xlabel('Frequency (Hz)');
ylabel('Phase (rad)');
hold on
end
Star Strider
2023-6-15
Are you certain this is correct:
myfilename = sprintf('OBP%dN.txt', k);
MATLAB is case-sensitive, and the file names must match exactly, even with respect to case, and the other file names were all lower-case. If that is the problem, the code should throw an error about not being able to find the file. Are there any errors?
So the first thing to do is to check to see what ‘mydata’ contains. Be certain that it is not empty.
更多回答(1 个)
Diwakar Diwakar
2023-6-3
Try this code. May be this code will help you.
% Specify the path to your text file
file_path = 'path/to/your/file.txt';
% Read the data from the text file
data = dlmread(file_path);
% Extract the time and amplitude columns
time = data(:, 1);
amplitude = data(:, 2);
% Set the parameters for the spectrogram
window_size = 256; % Size of the window for the spectrogram
overlap = 128; % Overlap between consecutive windows
% Compute the spectrogram
[~, frequency, time_spectrogram, power] = spectrogram(amplitude, window_size, overlap, [], 'onesided', 'yaxis');
% Plot the spectrogram
figure;
imagesc(time_spectrogram, frequency, 10*log10(power));
axis xy; % Flip the y-axis direction
colormap jet; % Choose a colormap
colorbar; % Add a colorbar for reference
xlabel('Time');
ylabel('Frequency');
title('Spectrogram');
4 个评论
Rahul
2023-6-3
Hi Diwakar,
Thank you for your help. I ran the code. But it gives error for line 2 for dlmread.
rgds,
rc
Diwakar Diwakar
2023-6-4
data = readmatrix(file_path);
If you are using MATLAB R2019b or later, you can use the readmatrix function instead of dlmread
Rahul
2023-6-4
Its not working
it gives error:
Error using pwelch
Expected x to be finite.
Plz advice.
rgds
Star Strider
2023-6-4
@Rahul — We already covered this and I told you how to fix it.
If we are going to solve it, you need to attach your data.
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!发生错误
由于页面发生更改,无法完成操作。请重新加载页面以查看其更新后的状态。
您也可以从以下列表中选择网站:
如何获得最佳网站性能
选择中国网站(中文或英文)以获得最佳网站性能。其他 MathWorks 国家/地区网站并未针对您所在位置的访问进行优化。
美洲
- América Latina (Español)
- Canada (English)
- United States (English)
欧洲
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
亚太
- Australia (English)
- India (English)
- New Zealand (English)
- 中国
- 日本Japanese (日本語)
- 한국Korean (한국어)