MATLAB: Résolution numérique d’un système d’équations intégrales et exponentielles

J’essaie de trouver la solution numérique à un système, que j’ai besoin de donner de nombreuses variables en entrée. Comment faire? J’ai joint le problème ci-dessous. J’ai essayé juste en utilisant la méthode normale mais cela ne fonctionne pas pour ce type d’équation car j’ai besoin du codage pour résoudre directement l’équation lorsque je donne des valeurs d’entrée car je suppose que c’est trop compliqué. Je ne suis pas trop grand sur MATLAB, donc toute aide serait appréciée! Merci

Meilleure réponse

  •  

    E1 = 3300;
    E2 = 29700;
    v1 = 0.41;
    v2 = 0.28;
    E=((1-v1^2)/E1)+((1-v2^2)/E2);
    H = 7.2087e-07;
    beta = 6.5;
    sigma = 0.00331;
    psi=((H*beta)/(E*sigma));
    row = 0.1;
    z1 = 0.00168;
    z2 = 0.00387;
    z3 = 0.00265;
    c=2*z2-z1-z3;
    C=((-beta*log(row))^2)/(c*sigma);
    deltac=C/2;
    gamma1 = 22.14;
    gamma2 = 22.14;
    gamma12 = 0;
    gamma=gamma1+gamma2-gamma12;
    phi=((gamma*beta)/(E*(sigma^2)));
    deltac1 = 12.25207;
    %Z = 0.00423;
    h = 0.00289;
    syms Z
    P=0.532*(1/psi)*beta*((1/((C^1/2)*(-log(row))))*(int(exp(-Z^2/2)*(Z - h)^(3/2), Z, h, (h - deltac1)))-(3.76*phi^(1/2)*(1/((C^1/2)*(-log(row))^(1/2))))*(int(exp(-Z^2/2)*(Z - h)^(3/4), Z, h, (h - deltac1))))+(2.5)*((1/C)*((int(exp(-Z^2/2)*(Z - h), Z, (h+deltac), inf)))-(phi/psi)*(1/C)*((int(exp(-Z^2/2), Z, (h+deltac), inf))));
    Pval = vpa(P);
    disp(Pval)