|
|
Wavelettransformation komplexer Signale |
|
Jonas45 |
Forum-Newbie
|
|
Beiträge: 1
|
|
|
|
Anmeldedatum: 29.10.12
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 29.10.2012, 01:31
Titel: Wavelettransformation komplexer Signale
|
|
|
|
|
Hallo Leute,
ich habe da leider ein Problem mit der Wavelet Transformation komplexer Signale. Wenn ich ein komplexes Testsignal mit "wavedec" einer Wavelet-Transformation unterziehe und dann mittels "waverec" zurücktransformiere stelle ich fest, dass der Fehler zwischen Rekonstruktion und Originalsignal ziemlich groß ist.
Weiß jemand, was man da machen kann????????????
Habe außer dem Haar-Wavelet auch alle anderen mir bekannten Wavelets ausprobiert, der Fehler tritt trotzdem noch auch.
Wäre für Hilfe sehr dankbar.
Der Code ist im Folgenden beigefügt:
clear;
close all;
%% generate chirp signal
f0 = 4000;%Hz
Tp = 0.005; %seconds
fs = 40000;%Hz
fc = 8000;%Hz
%chirp rate alpha
alpha = pi*f0/Tp;
% % % chirp duration in samples
Nc = floor(Tp*fs);
%generate time vector
t = (0:Nc-1)/fs;
wcm = 2*pi*fc-alpha*Tp;
s = exp( 1i*(wcm*t + alpha*t.^2) );% .* (td1>=0 & td1<=Tp) ;
% baseband reference signal
s = s .* exp(-1i*2*pi*fc*t);
%% generate complex baseband signal representing SAS echo assuming one
%% ideal point target
%fast time length
N_ft = 2048;
%round- trip delay in samples
N_d = 500;
%shift chirp signal to round-trip delay of the point target;
%round-trip delay set to 500 samples
signal = zeros(1,N_ft);
signal(N_d:N_d+length(s)-1) = s * exp(-j* 2*pi*fc*N_d /fs);
%% perform wavelet decomposition and reconstruction
wname = 'db5';
N_dec = 20;
[signal_wavelet, L] = wavedec(signal,N_dec,wname);
signal_rec = waverec(signal_wavelet,L,wname);
%% plot reconstruction error
figure;plot(abs(signal -signal_rec));
|
|
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
| RSS
Hosted by:
Copyright © 2007 - 2024
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|