IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

MATLAB Discussion :

Aide pour compiler un code


Sujet :

MATLAB

Vue hybride

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Membre éprouvé
    Inscrit en
    Novembre 2006
    Messages
    1 073
    Détails du profil
    Informations forums :
    Inscription : Novembre 2006
    Messages : 1 073
    Par défaut Aide pour compiler un code
    Je n'ai pas Matlab. Or j'ai besoin d'utiliser un programme en Matlab, qui est donné ci dessous.
    est ce que quelq'un peut le compiler, avec les données que j'envoie ci dessous, et me renvoyer la réponse. (en message privé ou sur ce thread. La réponse est une matrice).
    merci

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
     
    function cwt = CWT_Wavelab(x,nvoice,wavelet,oct,scale)
    % CWT -- Continuous Wavelet Transform
    %  Usage
    %    cwt = CWT_Wavelab(x,nvoice,wavelet,oct,scale)
    %  Inputs
    %    x        signal, dyadic length n=2^J, real-valued
    %    nvoice   number of voices/octave
    %    wavelet  string 'Gauss', 'DerGauss','Sombrero', 'Morlet'
    %    octave   Default=2
    %    scale    Default=4
    %  Outputs
    %    cwt      matrix n by nscale where
    %             nscale = nvoice .* noctave
    %
    %  Description
    %    
    %
    	if nargin<4,
    		oct = 2;
    		scale = 4;
    	end	
    % preparation
    	x = ShapeAsRow(x);
    	n = length(x);
    	xhat = fft(x);
    	xi   = [ (0: (n/2)) (((-n/2)+1):-1) ] .* (2*pi/n);
     
    % root
    	omega0 = 5;
     
    %	noctave = floor(log2(n))-2;
    %	noctave = floor(log2(n))-1;
    	noctave = floor(log2(n))-oct;
    	nscale  = nvoice .* noctave;
     
    	cwt = zeros(n,nscale);
     
    	kscale  = 1;
    %	scale   = 4;
    %	scale = 16;
     
    	for jo = 1:noctave,
    	    for jv = 1:nvoice,
    		   qscale = scale .* (2^(jv/nvoice));
    		   omega =  n .* xi ./ qscale ;
    		   if strcmp(wavelet,'Gauss'),
    				window = exp(-omega.^2 ./2);
    		elseif strcmp(wavelet,'DerGauss'),
                                    window = i.*omega.*exp(-omega.^2 ./2);
    		   elseif strcmp(wavelet,'Sombrero'),
    				window = (omega.^2) .* exp(-omega.^2 ./2);
    		   elseif strcmp(wavelet,'Morlet'),
    				window = exp(-(omega - omega0).^2 ./2) - exp(-(omega.^2 + omega0.^2)/2);
    		   end
    		   % Renormalization
    		   window = window ./ sqrt(qscale);
    		   what = window .* xhat;
    		   w    = ifft(what);
    		   cwt(1:n,kscale) = real(w)';
    		   kscale = kscale+1;
    		end
    		scale  = scale .*2;
        end
    x=[1,1,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]
    nvoice=4
    oct=2
    wavelet=Sombrero
    scale=20

  2. #2
    Membre émérite
    Avatar de rostomus
    Homme Profil pro
    Doctorant électronique et traitement du signal
    Inscrit en
    Décembre 2006
    Messages
    791
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 42
    Localisation : France

    Informations professionnelles :
    Activité : Doctorant électronique et traitement du signal

    Informations forums :
    Inscription : Décembre 2006
    Messages : 791
    Par défaut
    Vous devez definir la fonction " ShapeAsRow"

  3. #3
    Membre éprouvé
    Inscrit en
    Novembre 2006
    Messages
    1 073
    Détails du profil
    Informations forums :
    Inscription : Novembre 2006
    Messages : 1 073
    Par défaut
    ok, je vais voir ce que c'est.

  4. #4
    Membre éprouvé
    Inscrit en
    Novembre 2006
    Messages
    1 073
    Détails du profil
    Informations forums :
    Inscription : Novembre 2006
    Messages : 1 073
    Par défaut
    apparamment, ca transforme toute colonne en ligne.
    et laisse une ligne invariante.

+ Répondre à la discussion
Cette discussion est résolue.

Discussions similaires

  1. Aide pour compiler du code pour DOS/XP
    Par jkepler dans le forum Fortran
    Réponses: 1
    Dernier message: 22/11/2007, 10h42
  2. [VB.NET] besoin d'aide pour déchiffrer un code
    Par pcdj dans le forum Windows Forms
    Réponses: 10
    Dernier message: 27/06/2006, 11h32
  3. Aide pour comprendre un code
    Par Spacy_green dans le forum Assembleur
    Réponses: 2
    Dernier message: 13/02/2006, 13h22
  4. Aide pour comprendre le code
    Par jfreuff dans le forum Assembleur
    Réponses: 2
    Dernier message: 31/01/2006, 17h54
  5. Je besoin d'aide pour terminer mon code
    Par Paulinho dans le forum C++
    Réponses: 7
    Dernier message: 06/11/2005, 23h30

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo