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

Signal Discussion :

probleme au niveau du bruit


Sujet :

Signal

Vue hybride

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Membre régulier
    Homme Profil pro
    Étudiant
    Inscrit en
    Avril 2013
    Messages
    9
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Étudiant
    Secteur : High Tech - Opérateur de télécommunications

    Informations forums :
    Inscription : Avril 2013
    Messages : 9
    Par défaut probleme au niveau du bruit
    salut j'essaye de compiler un programme sur matlab qui traite l’Égaliseur ZF et MMSE en transmission MIMO ,mais ça me donne erreur au niveau de cette ligne v = NoiseGen(nr,d,2/snr(z));voilà ci-joint le programme ;veuillez m'aider s'il vous plait j'ai besoin d'une solution c'est pour mon PFE 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
    clc; clear all;
    nt = 4; nr = 4; d = 1000; MinErrs = 100; MaxData = 1e6; chs = 10;
    snr_dB = 0:3:21; snr = 10.^(snr_dB/10); Lsnr = length(snr);
    berZF = zeros(1,Lsnr); berMMSE = zeros(1,Lsnr);
    for e = 1:chs
        disp(['ch ' num2str(e) '/' num2str(chs)]);
        H = (randn(nr,nt)+sqrt(-1)*randn(nr,nt))/sqrt(2);
        G_ZF = inv(H'*H)*H';
        nErZF = zeros(1,Lsnr); nErMMSE = zeros(1,Lsnr); data = zeros(1,Lsnr);
        for z = 1:Lsnr
            disp(['step ' num2str(z) '/' num2str(Lsnr) ': SNR(dB) = ' num2str(snr_dB(z))]);
            G_MMSE = inv(H'*H + (1/snr(z))*eye(nt))*H';
            while((nErZF(z)<MinErrs || nErMMSE(z)<MinErrs) && data(z)<MaxData)
                data(z) = data(z) + 2*nt*d;
                a = randsrc(nt,d,[-1 1]) + sqrt(-1)*randsrc(nt,d,[-1 1]);
                %v = NoiseGen(nr,d,2/snr(z));
                r=H*a + v; r_zf=r; r_mmse=r;
                aTempZF = G_ZF*r; aTempMMSE = G_MMSE*r;
                aHatZF = FnDecodeQAM(aTempZF,'QPS',1);
                aHatMMSE = FnDecodeQAM(aTempMMSE,'QPSK',1);
                [dumy dumy nBerZF dumy] = FnSerBerQAM(a,aHatZF,4,'gray');
                [dumy dumy nBerMMSE dumy] = FnSerBerQAM(a,aHatMMSE,4,'gray');
                nErZF(z) = nErZF(z) + nBerZF;
                nErMMSE(z) = nErMMSE(z) + nBerMMSE;
            end
            berZF(z) = berZF(z) + nErZF(z)/data(z);
            berMMSE(z) = berMMSE(z) + nErMMSE(z)/data(z);
        end
    end
    berZF = berZF/chs; berMMSE = berMMSE/chs;
    semilogy(snr_dB,berZF,'k-s',snr_dB,berMMSE,'m-o'); grid on
    legend('ZF','MMSE'); grid on; xlabel('SNR (dB)'); ylabel('BER');

  2. #2
    Invité
    Invité(e)
    Par défaut
    Bonjour,

    La fonction NoiseGen n'étant pas une fonction MATLAB, ce serait utile d'avoir le message d'erreur.

  3. #3
    Membre régulier
    Homme Profil pro
    Étudiant
    Inscrit en
    Avril 2013
    Messages
    9
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Étudiant
    Secteur : High Tech - Opérateur de télécommunications

    Informations forums :
    Inscription : Avril 2013
    Messages : 9
    Par défaut erreur afficher
    voilà ce qu'il m'affiche comme erreur :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    ch 1/10
    step 1/8: SNR(dB) = 0
    ??? Error using ==> mrdivide
    Matrix dimensions must agree.
    
    Error in ==> Untitled3 at 16
    v = NoiseGen(nr,d,2/(snr));

  4. #4
    Invité
    Invité(e)
    Par défaut
    Ton code initial était v = NoiseGen(nr,d,2/snr(z));.
    Code avec lequel tu ne devrais pas avoir cette erreur.

  5. #5
    Membre régulier
    Homme Profil pro
    Étudiant
    Inscrit en
    Avril 2013
    Messages
    9
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Étudiant
    Secteur : High Tech - Opérateur de télécommunications

    Informations forums :
    Inscription : Avril 2013
    Messages : 9
    Par défaut
    j'ai oublié d'enlevé le % dans le code que j'ai met sur ce forum , j'ai toujours la même erreur avec ce code je comprend pas pourquoi
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    v = NoiseGen(nr,d,2/snr(z));

  6. #6
    Invité
    Invité(e)
    Par défaut
    Exactement la même... j'en doute snz(z) n'est qu'un scalaire.
    Matrix dimensions must agree.

Discussions similaires

  1. Débutant probleme au niveau "tutorial java"
    Par carton dans le forum Java ME
    Réponses: 3
    Dernier message: 22/07/2006, 10h54
  2. Probleme au niveau de la connexion de tomcat avec oracle8i
    Par hamska2 dans le forum Tomcat et TomEE
    Réponses: 1
    Dernier message: 08/05/2006, 01h53
  3. [JMF] Problème au niveau de RTP
    Par Alienx dans le forum Multimédia
    Réponses: 2
    Dernier message: 05/03/2006, 16h50
  4. [Thread] petit problème au niveau du reveil d'un thread
    Par sagitarium dans le forum Concurrence et multi-thread
    Réponses: 3
    Dernier message: 13/01/2006, 12h35
  5. PROBLEME haut niveau BACKUP
    Par Dev_Michael dans le forum Administration
    Réponses: 7
    Dernier message: 10/03/2005, 16h32

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