Articles
CWT matlab-funktion [stängd]
On februari 17, 2021 by admin <åt sidan class = "s-notice s-notice__info js-post-notice mb16" role = "status">
Stängd. Denna fråga är utanför ämnet . För närvarande accepteras inte svar.
Kommentarer
- Spaghetti händer på grund av den komplexa attityden hos vissa CWT, som du redan lyckades med Fourier.
- betyder detta att tillämpa cwt till den här signalen är inte möjlig?
- Det är möjligt, och möjligen mer informativt än med Fourier, eftersom signalen inte är stillastående.
- din spaghetti är förmodligen en artefakt av ditt tomtsamtal. , inte cwt
- I ' jag röstar för att stänga den här frågan som off-topic eftersom problemet har att göra med att anropa matlabs plotfunktion med komplexa värden och inte ett problem med wavelets
Svar
Försöker gissa vilken signal du analyserar och syftet, här är en demo, på en riktig signal, med hälften av Fourier-spektrumet och motsvarande kontinuerliga wavelet-transform-skalogram.
Här misstänker jag att signalen är för kort (utan ytterligare mål) för att FFT och CWT ska ge tolkbara resultat. Matlab-koden är:
nsample = 64; % An odd number timeSampling = 1/nsample; time = (0:nsample-1)*timeSampling; ratioSecondHalf = 20; data = zeros(nsample,1); data(1:nsample/2,1) = rand(nsample/2,1)-0.5; data = medfilt1(data,5); data(nsample/2+1:end,1) = rand(nsample/2,1)/ratioSecondHalf; [fftR,fftAxe] = FFTR(data,timeSampling); [cwtCoeff,cwfFreq] = cwt(data,1:64,"morl",timeSampling); figure(1);clf subplot(3,1,1) plot(time,data,"x-");;axis tight xlabel("Time (a. u.)") ylabel("Amplitude (a. u.)") subplot(3,1,2) plot(fftAxe,fftR,"x-");axis tight xlabel("Frequency (a. u.)") ylabel("Amplitude (a. u.)") subplot(3,1,3) h=imagesc("XData",time,"YData",2*cwfFreq/pi,"CData",abs(cwtCoeff));axis tight xlabel("Time (a. u.)") ylabel("Frequency (a. u.)")
FFTR.m erhålls från här .
Kommentarer
- i ' Vi har också provat detta på en längre signal men det gav mig ändå spagetti tror jag det ' s eftersom det ' är icke-stationärt som du sa.
- Abslute-värde saknas någonstans
- OCH längre längd handlar om tolkbarhet, inte spagetti
- jag gjorde det med absolut värde och jag fick bättre resultat, inte mer spagetti.
- Bra, spagetti uppstår från plottningskomplex värden som 2D-poäng (real + imaginär)
Lämna ett svar