hi very valuable group members. I wish you all a good day. I ask you to help Hilbert Huan transformation point. I have two different nonlinear signals. I obtained imfs for these signals. I need the matlab code for the hilbert huang transform of these

EXAMPLE SİGNAL 1 CODE

>> nsamples=10000;
fs=1e3;                                                                                                                                  
 
time=(1:nsamples-1)/fs;
 
signal=0.5*time+sin(pi*time)+sin(2*pi*time)+sin(6*pi*time);%
>> plot(signal)
>> 

EXAMPLE SİGNAL 1 RESULT

EXAMPLE SİGNAL 1 İMF CODE

>> clear all
close all
 
nsamples=10000;
fs=1e3;
 
time=(1:nsamples-1)/fs;
 
signal=0.5*time+sin(pi*time)+sin(2*pi*time)+sin(6*pi*time);%
threshouldSifting=0.1;
monotone=0;
indexIMF=0;
signalmem=signal;
while (not(monotone))
    h=signal;
    SD=1;
    while(SD>threshouldSifting)
        %% upper
        [pksMax,locsMax] = findpeaks(h);
        pksMax=[h(1), pksMax, h(end)];
        locsMax=[1, locsMax, length(h)];
        upperSpline=spline(time(locsMax),pksMax);
        spmax=ppval(upperSpline,time);
        
        %% lower
        [pksMin,locsMin] = findpeaks(-1*h);
        pksMin=-1*pksMin;
        pksMin=[h(1), pksMin, h(end)];
        locsMin=[1, locsMin, length(h)];
        lowerSpline=spline(time(locsMin),pksMin);
        spmin=ppval(lowerSpline,time);
        
        %% evaluate the mean
        
        mspline=(spmax+spmin)/2;
 
        hprec=h;
        h=h-mspline;
        SD=sum((hprec-h).^2)/sum(hprec.^2);  %criterio di stop proposto da by Huang et al. (1998)
    end
 
    indexIMF=indexIMF+1;
    IMF(indexIMF,:)=h;
    signal=signal-h;
    monotone=all(diff(signal)>0)||all(diff(signal)<0);
end
subplot(indexIMF+2,1,1), plot(time,signalmem),title('Segnale di ingresso'),xlabel('Tempo [s]');
 
for index=1:indexIMF
    subplot(indexIMF+2,1,index+1), plot(time,IMF(index,:)),title(['IMF_',num2str(index)]),xlabel('Tempo [s]');
end
 
subplot(indexIMF+2,1,indexIMF+2), plot(time,signal),title('Residuo'),xlabel('Tempo [s]');
>> 
>> 

EXAMPLE SİGNAL 1 IMF RESULT

HELP TOPİC

WHAT I WANT : HİLBERT HUANG TRANSFORM OF THİS SİGNAL AND RESULT GRAPHS. SO; FREQUENCY , TİME AND ENERGY RESULTS.

MATLAM CODE AND GRAPHS ????

………………………………………………………………………………………………………..

EXAMPLE SİGNAL 2 CODE

>> Fs= 10000;
Ts=2/Fs;
L=150;
t=(0:L-1);
x=36*sin(2*pi*50*t)+sin(2*pi*120*t);
plot(t,x);
>> 

EXAMPLE 2 SİGNAL RESULT

EXAMPLE 2 SİGNAL İMF CODES

>> clear all
close all
 
nsamples=10000;
fs=1e3;
 
time=(1:nsamples-1)/fs;
 
signal=36*sin(2*pi*50*time)+sin(2*pi*120*time);
%
threshouldSifting=0.1;
monotone=0;
indexIMF=0;
signalmem=signal;
while (not(monotone))
    h=signal;
    SD=1;
    while(SD>threshouldSifting)
        %% upper
        [pksMax,locsMax] = findpeaks(h);
        pksMax=[h(1), pksMax, h(end)];
        locsMax=[1, locsMax, length(h)];
        upperSpline=spline(time(locsMax),pksMax);
        spmax=ppval(upperSpline,time);
        
        %% lower
        [pksMin,locsMin] = findpeaks(-1*h);
        pksMin=-1*pksMin;
        pksMin=[h(1), pksMin, h(end)];
        locsMin=[1, locsMin, length(h)];
        lowerSpline=spline(time(locsMin),pksMin);
        spmin=ppval(lowerSpline,time);
        
        %% evaluate the mean
        
        mspline=(spmax+spmin)/2;
 
        hprec=h;
        h=h-mspline;
        SD=sum((hprec-h).^2)/sum(hprec.^2);  %criterio di stop proposto da by Huang et al. (1998)
    end
 
    indexIMF=indexIMF+1;
    IMF(indexIMF,:)=h;
    signal=signal-h;
    monotone=all(diff(signal)>0)||all(diff(signal)<0);
end
subplot(indexIMF+2,1,1), plot(time,signalmem),title('Segnale di ingresso'),xlabel('Tempo [s]');
 
for index=1:indexIMF
    subplot(indexIMF+2,1,index+1), plot(time,IMF(index,:)),title(['IMF_',num2str(index)]),xlabel('Tempo [s]');
end
 
subplot(indexIMF+2,1,indexIMF+2), plot(time,signal),title('Residuo'),xlabel('Tempo [s]');
>> 
>> 

EXAMPLE 2 İMF RESULTS

HELP TOPİC

WHAT I WANT : HİLBERT HUANG TRANSFORM OF THİS SİGNAL AND RESULT GRAPHS. SO; FREQUENCY , TİME AND ENERGY RESULTS.

MATLAM CODE AND GRAPHS ????

………………………………………………………………………..

SAME THİS ARTİCLE APPLİCATİON

MAİL : [e-mail]

[Walter Roberson]

sukru ozdemir you flagged your own Question as unclear. If you feel your question is unclear, then you should revise your question to clarify it.

[sukru ozdemir]

Thank you very much for your historic reply. My question is; one of my nonlinear signals is necessary for the hilbert huang transformation. I need the matlab code for that

[sukru ozdemir]

I need the matlab code for a work just like in this article. codes; should be written as a continuation of the code I have written

Demagnetization Fault Detection by means of Hilbert

Huang Transform of the stator current Decomposition

in PMSM

[sukru ozdemir]

https://www.researchgate.net/publication/224110913_Fault_Detection_by_Means_of_Hilbert-Huang_Transform_of_the_Stator_Current_in_a_PMSM_With_Demagnetization

the article can be accessed from this link. As in this article; hilbert huang transformation of any of the two nonlinear signals that I have written. finally time, frequency and energy graph

[sukru ozdemir]

At this point, if you help me and send the necessary code and files to my mail I would be very happy. Thank you very much again.

[Walter Roberson]

It is not clear why using hht() would not work for you?

[sukru ozdemir]

for signals; i don't know the matlab code required to do hilbert huang transform.

[sukru ozdemir]

so; For the above two different signals, EMD must first be applied and the IMFs must be found. I found the IMFs. but I don't know the mablab code needed for HHT and I couldn't write it.

[sukru ozdemir]

or my request, if it is possible for you, can you write a matlab code for me to be a hilbert huang conversion for any nonlinear signal

[Walter Roberson]

? You yourself linked to hht() that appears to be the transform you are looking for. That page gives examples of the use. It is not clear what else is needed?

[sukru ozdemir]

what i want is this: i need matlab code for hilbert huang transformation analysis of this signal

x=36*sin(2*pi*50*t)+sin(2*pi*120*t);

[sukru ozdemir]

my assignment is that gentleman; my teacher said; write one nonlinear signal and find me his hilbert huang transformation

[Walter Roberson]

It is not clear what the difficulty is in calling emd() on your signal and passing the result to hht()?

[sukru ozdemir]

I need code I can't

[sukru ozdemir]

>> 
clear variables
close all
fs=100;
t=0:1/fs:10;
 
 
fresolution=0.1;
signal=cos(2*pi*t+0.5*sin(4*pi*t))+sin(pi/10*t);
 
 
spec=zeros(length(signal),length(signal));
 
 
for indexIMF=1:min(size(IMF))
    Z=hilbert(IMF(indexIMF,:));
    modZ=abs(Z);
    phaseZ=unwrap(2*atan(imag(Z)./real(Z)));
    phaseZ=phaseZ/2;
    freq=diff(phaseZ)/2/pi*fs/fresolution;
    
    indexFreqVal=find(freq>0);
    
    for index=1:length(indexFreqVal)
   spec(indexFreqVal(index),ceil(freq(indexFreqVal(index))))=spec(indexFreqVal(index),ceil(freq(indexFreqVal(index))))+modZ(indexFreqVal(index));
        
    end
    
end
 
mesh((0:length(t)-1)*fresolution, t,spec)
xlim([0 10])
 
xlabel('Freq (rad/s)')
ylabel('time (s)')

....................................................................................

I want to write code like the code above. for my signal. ie; for the hilbert transformation of my signal

[sukru ozdemir]

when I add this to the code that finds the imf, it gives an error so I need help

[Walter Roberson]

Given a variable named signal how would you call empirical mode decomposition on the variable?

[sukru ozdemir]

So, sir, I mean, can you write me this code?

[sukru ozdemir]

I'm asking because I can't do that anyway because I said this, sir

[sukru ozdemir]

I'll need the following help: I ​​have a nonlinear signal. the imfs of this signal will be found by decomposition methods (decomposiition). and then the hilbert transformation. time frequency and energy dimension. That 's all my problem. can't write the required codes

[Walter Roberson]

So, sir, I mean, can you write me this code?

Will I get the academic credit for doing the assignment, or will you get the academic credit?

I'm asking because I can't do that

In MATLAB, what is the syntax for passing a parameter whose variable name you know ("signal") to a function whose name you know ("edf") because you read it out of the documentation at https://www.mathworks.com/help/signal/ref/hht.html ?

[sukru ozdemir]

it's just an assignment, but it's very important to me. I read that page, but I don't know the code that will apply to my signal

[Walter Roberson]

In MATLAB, what is the syntax for passing a parameter whose variable name you know ("signal") to a function whose name you know ("edf") because you read it out of the documentation at https://www.mathworks.com/help/signal/ref/hht.html ?

[sukru ozdemir]

https://www.researchgate.net/publication/224110913_Fault_Detection_by_Means_of_Hilbert-Huang_Transform_of_the_Stator_Current_in_a_PMSM_With_Demagnetization

So I want to do the same work on my link on this link. that is all. first: nonlinear signal generated.

second: I found the imfs of this signal. Third: I need to find the hilbert transformation of this signal, as in the link. I stayed in the third step.

[sukru ozdemir]

first signal, then imf to be found, (I found), then hilbert transformation. I stayed in the Hilbert transformation. that is all. I still do not understand what you do not understand, my precious brother

[Walter Roberson]

The imf (intrinsic mode functions) is what you would get as output of HHT, so you would appear to have already done the work.

Is there are reason you did not use edm and hht ?

[sukru ozdemir]

don't understand what you mean, please excuse me. I need the hht code to see my signal in the frequency, time and energy domain

Answer 1

[sukru ozdemir]

https://www.mathworks.com/help/signal/ref/hht.html

same application