이 질문을 팔로우합니다.
- 팔로우하는 게시물 피드에서 업데이트를 확인할 수 있습니다.
- 정보 수신 기본 설정에 따라 이메일을 받을 수 있습니다.
Error with data processing
조회 수: 2 (최근 30일)
이전 댓글 표시
I get an error if I run this matlab code of processing data for Discrete Wavelet Transformation analysis
댓글 수: 4
Jonas
2022년 7월 19일
please post your error message and supply all required data files used in your code
채택된 답변
Walter Roberson
2022년 7월 19일
You need getmswtfeat from https://www.mathworks.com/matlabcentral/fileexchange/33146-feature-extraction-using-multisignal-wavelet-packet-decomposition
댓글 수: 31
Walter Roberson
2022년 7월 20일
you appear to load() a file, and use parts of some variables from the file. However we are not told what datatype is in the file.
Based on the error messages we can tell that you are expecting that the mat file contains double precision data, but that instead it has a table() object.
You need to extract data from the table object, or else you need to change the mat file to have double precision data instead of table.
It is as-if you used readtable() to read the data when readmatrix() might possibly have been a better choice. But that depends on whether the file is all numeric or not
John Amoo-Otoo
2022년 7월 21일
편집: Walter Roberson
2022년 7월 22일
Please attached the data in table format
temp = load('John_FaultData0km0ohm.mat');
fieldnames(temp)
ans = 2×1 cell array
{'DataPG0km0hmsFaultData' }
{'SteadyStateNoneFaultState'}
temp.DataPG0km0hmsFaultData(1:5,:)
ans = 5×2 table
Time DataPG0km0hmsFaultData
______ ______________________
19.8 0.00030982
19.801 0.00030954
19.802 0.00030947
19.803 0.00031002
19.804 0.00030967
temp.SteadyStateNoneFaultState(1:5,:)
ans = 5×2 table
Time SteadyStateNoneFaultState
_______ _________________________
0 0
0.00099 6.4183e-09
0.00198 2.0619e-07
0.00297 1.0102e-06
0.00396 2.7701e-06
Walter Roberson
2022년 7월 22일
I do not know what you want to do with the features extracted.
Note that there is a theoretical fault in the code. You remove data that has nan, but you do not take into account that that leaves a hole in the timing. You never use the timing information.
In practice the only nan is at the very last entry in the SteadyState, so removing it makes no timing difference for this data.
%% load the data first
load('John_FaultData0km0ohm.mat')
%% does this signal has any NaNs, if so remove
SteadyStateNoneFaultState = rmmissing(SteadyStateNoneFaultState);
DataPG0km0hmsFaultData = rmmissing(DataPG0km0hmsFaultData);
SSTime = SteadyStateNoneFaultState.Time;
SSNFS = SteadyStateNoneFaultState.SteadyStateNoneFaultState;
DPTime = DataPG0km0hmsFaultData.Time;
DPFD = DataPG0km0hmsFaultData.DataPG0km0hmsFaultData;
%% Define the filtering and inspect
n = 8; %% defines window length
w = [-ones(n,1); ones(n,1)];
SSNFS = filter(w, n, SSNFS);
DPFD = filter(w, n, DPFD);
%% Normalize signals
med_training = prctile(SSNFS,50);
iqr_training = iqr(SSNFS);
SSNFS = (SSNFS-med_training)./iqr_training;
med_fault = prctile(DPFD,50);
iqr_fault = iqr(DPFD);
DPFD = (DPFD-med_fault)./iqr_fault;
% Plot & Observe the data
subplot(2,1,1)
plot(DPTime, DPFD)
title('filtered DataPG0km0hmsFaultData')
subplot(2,1,2)
plot(SSTime, SSNFS)
title('filtered SteadyStateNoneFaultState')
%% Let's observe the FFT power spectrum for differences
feat_fault = getmswtfeat(DPFD,32,16,100000);
feat_Good = getmswtfeat(SSNFS,32,16,100000);
John Amoo-Otoo
2022년 7월 24일
Walter you are a Guru and super consultant in MATLAB. the script in data proicessing runs perfectly. The only script now that requires debugging is getmsfeat.m. After that I will be squared away
John Amoo-Otoo
2022년 7월 24일
The getsmsfeat.m is the Discrete Wavelet Transform Decomposition script using different mother wavelets( Daubechies, Symlets, Coiflets,etc)
Walter Roberson
2022년 7월 24일
What error are you seeing? Earlier when I ran your code after fixing the other problems, the only problem that I had with the feature extraction that you posted earlier, was that the code needed Wavelet Toolbox. When I installed that, it ran to completion. If it runs without error but does not produce the output you want then you will need to be specific about what I should be looking for.
Walter Roberson
2022년 7월 26일
MATLAB would have stopped after a single error, not continued to tell you about all the different problems encountered in getmswtfeat .
I do not have any problems on my system when I run the processdata that I posted.
Which MATLAB release are you using?
That error you are getting about nargin would only apply if somehow your getmswtfeat is a script instead of a function. What you attached is a function. You should
which -all getmswtfeat
to check in case somehow you have an additional getmswtfeat that is interfering.
John Amoo-Otoo
2022년 7월 30일
Walter, Good day, from the attached, for an 8 level deocmposition for daubechies 2, I am expecting 9 features from column 1 through column 9. In the columns are these the detailed coefficients from column 1 through 8 and the column 9 should be the approximate. Are these entropies or detailed coefficients?
Walter Roberson
2022년 7월 31일
Those are entropies.
feature_out = zeros(numwin,(J+1)*Nsignals);
tab_entropy(:,k) = -sum(prob.*log(prob),2);%./size(percentENER(:,st:en),2);
feature_out(:,(1:((J+1)))+(dims-1)*(J+1)) =tab_entropy;
feature_out is not written to anywhere else, and tab_entropy is not written to anywhere else (other than pre-allocation with zeros)
John Amoo-Otoo
2022년 8월 2일
Walter, it looks like the column one is the detailed coefficients and the plots is the entropy. Just checking if that is true. Also if I am looking for the highest detailed coefficients I can see negative and positive number the column, which one is the highest.
Walter Roberson
2022년 8월 2일
Line 51: feature_out = zeros(numwin,(J+1)*Nsignals);
Line 97: tab_entropy = zeros(numOfSIGs,level+1);
Line 110: tab_entropy(:,k) = -sum(prob.*log(prob),2);%./size(percentENER(:,st:en),2);
Line 114: feature_out(:,(1:((J+1)))+(dims-1)*(J+1)) =tab_entropy;
Those are the only lines that write to feature_out or tab_entropy that are not commented out. You can see that column 1 of tab_entropy is not written to differently than the other columns, and you can see that nothing other than tab_entropy is written into feature_out (other than the initial zeros.)
Therefore, column 1 of the output is an entropy, not a detailed coefficient.
John Amoo-Otoo
2022년 8월 20일
Walt, please could you look into this attached error. My last simulation run perfectly and I changed the data to process another data and feature fault was given me an error.
Also for the feature out code if I want to add a script for tab detailed coefficients juts like tab entropy how do I do that
John Amoo-Otoo
2022년 8월 20일
Walter, I was able to resolve the error. could you send me a script line for tab calculating Detailed Coefficients
John Amoo-Otoo
2022년 8월 21일
Also if you have a script for tabing Detailed coefficient you can send it to me
Walter Roberson
2022년 8월 21일
As I do not have any script for that purpose, should I be arranging to send you the lack of script, or should I be arranging to not send you the script that I do not have?
John Amoo-Otoo
2022년 8월 22일
Walter, please could you look into this error. The script is for calculating the detailed coefficients using Matlab Discrete Wavelet Transform
Walter Roberson
2022년 8월 28일
As I reported to you in https://www.mathworks.com/matlabcentral/answers/1763300-error-with-data-processing#comment_2278845 your very last entry in your data is NaN. wavedec() cannot handle NaN, so you need to filter your input before calling wavedec() ... like I showed in that comment.
John Amoo-Otoo
2022년 8월 28일
Walter, Thank you it worked. I have been able to add a single line script that will remove the NANs
John Amoo-Otoo
2022년 9월 27일
Walt, a quick question. I am processing another script with a signal that I added Gaussian Noise to mimic the intrusion of lightning and switchig disturbance during the fault condition signal. and attached is the error I got. Please could you look into it at your convenience
Walter Roberson
2022년 9월 27일
The code before the error line does not define any variable named signal but it does define a variable named signals
추가 답변 (1개)
John Amoo-Otoo
2022년 7월 20일
Jonas, please attached is my error for data pocessing
댓글 수: 2
참고 항목
카테고리
Help Center 및 File Exchange에서 Get Started with Signal Processing Toolbox에 대해 자세히 알아보기
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 (한국어)