implantation d'une fonction contenant une integrale sur matlab

[invite164df21a]

bonjour,
j'ai une fonction qui contient une intégrale que j'arrive pas à implémenter sur matlab, la voici

avec:
n et B deux inconnues à déterminer ultérieurement (c'est le but du code: optimisation de n et B)
t est le temps: entier positif
je suis débutante sur matlab et je confond tous, et toujours des messages d'erreur
-----

[fiatlux]

salut

tu sais à quoi tu veux arriver? ce que je veux savoir est-ce que par exemple tu sais que cette intégrale vaut (par ex.) 37 et tu veux connaître n et B en fonction de ça?
Et si j'ai bien compris, t est connu?

j'ai essayé d'avoir l'expression approchée de l'intégrale avec:

syms y n B t
f = y*((n-1)*B*1400*t+(-35.95*y)^(1-n))^(1/(1-n));
int(f,-1.8,1.8)

mais malheureusement Matlab ne trouve pas de solution d'expression explicite de l'intégrale. Message d'erreur: Warning: Explicit integral could not be found.

Par contre si tu sais à quoi tu veux arriver, tu peux résoudre numériquement avec une fonction comme lsqonlin ou lsqcurvefit

[invite164df21a]

bonsoir,

tu sais à quoi tu veux arriver? ce que je veux savoir est-ce que par exemple tu sais que cette intégrale vaut (par ex.) 37 et tu veux connaître n et B en fonction de ça?
Et si j'ai bien compris, t est connu?

en fait non je sais pas, ce que j'ai c'est la valeur de f en certain instant t, bien sur t est connu, n et b sont inconnus je cherche à les déterminer,
pour ça j'ai un code d'optimisation qui marche correctement avec des fonction simples, j'ai voulu l'appliquer à cette fonction, mais pour ça, il faut implémenter la fonction j'ai déjà essayé ta proposition ainsi que la méthode du trapèze mais je reçoit toujours des messages d'erreurs

[invite164df21a]

en fait il ne reconnait pas syms

??? Undefined function or method 'syms' for input arguments of type 'char'.

Error in ==> main>f at 10
syms v a(1) a(2) t

[A voir en vidéo sur Futura]

[fiatlux]

je reçoit toujours des messages d'erreurs

Mais si tu reçois des messages d'erreurs, ça veut sans doute dire qu'il y a une erreur, pas forcément que la méthode est fausse... A la limite montre-nous ton code, comme ça on pourra voir quel message d'erreur ça génère et voir où ça cloche.

A part ça, si tu sais par exemple que pour t=3 tu as f=34, tu peux sans problème déterminer n et B de façon itérative, genre avec lsqonlin ou lsqcurvefit comme j'ai dit.


EDIT:

il ne reconnait pas syms

je pense plutôt qu'il faut genre nommer c=a(1) et d=a(2) et faire syms v c d t

[invite164df21a]

bon voici le code

function main
global ti yi	 
ti = [10800 43200 86400 172800 604800];
yi = [0.6 0.7 1 1.2 1.4];
a0 = [0.1 5]; 
asol = fminsearch ( @F, a0)
plot (ti, yi, '*', ti, f(asol, ti))

function y = f(a,t) 
syms v a(1) a(2) t
fc = v.*((a(2)-1).*a(1).*1400.*t+(-35.95.*v).^(1-a(2))).^(1/(1-a(2)));
int(fc,-1.8,1.8);
y=0.035.*int(fc,-1.8,1.8);

function out = F(a)
global ti yi
out = sum ( (yi - f(a,ti)).^2 );