Bonjour, cliquez-ici pour vous inscrire et participer au forum.
  • Login:



+ Répondre à la discussion
Affichage des résultats 1 à 4 sur 4

Svvp Aide code Matlab

  1. bochra017

    Date d'inscription
    juin 2017
    Messages
    3

    Svvp Aide code Matlab

    Bonjour,

    j'ai un code Matlab pour 1 dimension, j'ai besoin de votre aide pour écrit même code pour 2 dimension.

    J’ai besoin de montrer que l'erreur Monte-Carlo est fonction de 1/sqrt(N) . Mon prof me dit que : Pour obtenir les erreurs MC, pour chaque nombre de points N (10^2, 10^3, 10^4...10^7) vous refaites le calcul (30 fois par exemple) pour calculer sigma, et l'erreur sera sigma/sqrt(N). vous aurez l'erreur pour chaque N, vous tracer la courbe pour obtenir la courbe en bleu. Et tracer 1/sqrt(N) qui est la courbe verte. Il me donne un exemple d’intégrale pour 1D ( code+commentaires explicatif ) et il me demande de faire la même chose pour 2D, pour montré que l’erreur converge vers 1/sqrt(N) n'importent qu’il dimension. Mais je n’arrive pas affaire , j’ai besoin l'aide svvvp

    code 1D :
    Code:
    close all
    clc
    format long
    
    N = [10^2 10^3 10^4 10^5 10^6 10^7]';
    C = 0.25; %fmax
    a = 0;
    b = 1;
    
    %Nombre des essais
    MCTries = 7;
    
    %On aura besoin de ces 2 tableaux pour garder enregister lerreur 
    %correspondante à chaque essaye (pour cet exemp: pour chaque N on refera le calcul 7 fois.
    %on prendera à la fin la moyenne des erreurs.
    array_sigmas = zeros(MCTries,1);
    array_errors = zeros(length(N),1);
    
    for i=1:length(N) % boucle sur les N
        for j=1:MCTries % boucle sur les essais
            I = 0;
            f = zeros(N(i,1),1);
            for k = 1:N(i,1) % boucle sur les xi
                x = rand();
                fx = x*(1-x)*sin(200*x*(1-x))^2;
                I = I + fx;
                f(i,1) = fx;%stockage des fx
            end
            % I = (b-a)*I/N(i,1) % valeur de l'integrale
        array_sigmas(j,1) = sqrt(var(f)); % calcul de l'erreur (la racine de la variance)  
        end
        % on prends la moyenne des essais
        array_errors(i,1) = mean(array_sigmas);
    end
    
    
    %data to plot
     x = N(:,1);
    y0 = 1./sqrt(N(:,1));
    y1 = array_errors(:,1)
    
    %plotting
    figure1 = figure;
    axes1 = axes('Parent',figure1);
    hold(axes1,'on');
    box(axes1,'on');
    grid(axes1,'on');
    plot(x,y0,'-r',x,y1,'-bo');
    set(axes1,'XMinorTick','on','XScale','log','YMinorTick','on','YScale','log');
    xlabel('N') % x-axis label
    ylabel('errors') % y-axis label
    legend('1/sqrt(N)','error in MC integration')
    
    
    
    <
    THISONE.png

    pour 2D : par exemple cette fonction f(x,y)=xye^((-x^2)y) ou bien une fonction de votre choix.pic1.png

    Et merci

    -----

     


    • Publicité



  2.  

  3. bochra017

    Date d'inscription
    juin 2017
    Messages
    3

    Re : Svvp Aide code Matlab

    J'ai vraiment besoin de votre aide svp! ,tout ce que je veux c'est :même code pour 2D.
     

  4. Antoane

    Date d'inscription
    août 2007
    Localisation
    IDF
    Âge
    26
    Messages
    12 259

    Re : Svvp Aide code Matlab

    Bonjour,

    EXERCICES et FORUM
    Pas grand monde n'a l'air motivé pour faire ton travail ; tu auras peut-être davantage de réponses une fois que tu auras proposé un code que les participants pourront t'aider à améliorer/finir.
    Deux pattes c'est une diode, trois pattes c'est un transistor, quatre pattes c'est une vache.
     


    • Publicité







Sur le même thème :


    Sur le même sujet

    Définitions



 

Discussions similaires

  1. [matlab] programmation : peut-on générer un arbre des fonctions d'un code Matlab ?
    Par bratisla dans le forum Logiciel - Software - Open Source
    Réponses: 3
    Dernier message: 09/07/2015, 11h12
  2. [Microbiologie] interprétation svvp
    Par alyal dans le forum Biologie
    Réponses: 1
    Dernier message: 07/01/2013, 13h11
  3. Aide code matlab
    Par dalida_73 dans le forum Programmation et langages, Algorithmique
    Réponses: 3
    Dernier message: 10/04/2012, 23h35
  4. code matlab
    Par boch dans le forum Logiciel - Software - Open Source
    Réponses: 0
    Dernier message: 29/11/2010, 10h34
  5. est ce du code matlab?
    Par ABN84 dans le forum Logiciel - Software - Open Source
    Réponses: 2
    Dernier message: 17/11/2008, 10h01