estimation de paramètres

**sibiath10** · 26/08/2013, 16h06

Bonjour à tous
j'ai un petit souci avec un programme d'estimation et j'aimerais avoir votre aide
je veux faire tourner ce programme mais en vain

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
ts=[];
Ei=40;S=0.0104;mc=2;Rel=41.3;U=23;e=0.03;lam=3;Rc1=1e-8;E=2258;es=3e-4;lami=0.033;ei=0.05;
tf=20;
a=(lam/E)^2;Rc1=Rc1/S;fi=U^2/Rel;roc=lam/a;
mc=1.2e6*S*es;
 
 
pas=1;
j=1;
for i=1:pas:tf
    ts(j)=i;
    j=j+1;
 end
 
 
B0=[E;lam;Rc1;mc]
 
Tc22=PCInvlap2(ts,B0);
plot(ts,Tc22,'k')

avec les sous programmes suivants:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
function [Tc22] = PCInvlap2(ts,B0)
 
E=B0(1,1);
lam=B0(2,1);
Rc1=B0(3,1);
mc=B0(4,1);
 
Tc22=Invlap('TetaPC2',ts,0,1e-9,E,lam,Rc1,mc);

et

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
function fp = TetaPC2(p,E,lam,Rc1,mc)
 
global e fi lam a ai lami ei
 
flag=0;
taille=size(p);
S=0.0105;
rce=2/S;
a=(lam/E)^2;
ai=(lami/Ei)^2;
Rc2=1e-8;Rc3=1e-8;
 
if taille(1)>taille(2),
	flag=1;
	p=p';
end;
 
for i=1:taille(1)
     q = (p(i) / a) ^ 0.5;
     qi = (p(i) / ai) ^ 0.5;
     mcp=rce*p(i); 
     A=cosh(q*e);B=sinh(q*e)/(lam*q);C=lam*q*sinh(q*e);D=cosh(q*e);
 
      Ai=cosh(qi*ei);Bi=sinh(qi*ei)/lami/qi;Ci=lami*qi*sinh(qi*ei);Di=cosh(qi*ei);
      M10=[1,Rc3;0,1];M20=[Ai,Bi;Ci,Di];
      M0=M10*M20;
      A1=M0(1,1);B1=M0(1,2);C1=M0(2,1);D1=M0(2,2);
 
      M12=[1,0;mcp,1];M22=[1,Rc1;0,1];M32=[A,B;C,D];
      M2=M12*M22*M32;
      B02=M2(1,2);D02=M2(2,2);
 
      M3=[1,Rc2;0,1];
      M00=M12*M22*M32*M3*M20;
      B0=M00(1,2);D0=M00(2,2);
 
     Y0=(fi/p(i))/((D1/B1)+(D0/B0));
 
 
     Y2=Y0*(D02-((D0/B0)*B02));
     THETA2(i)=Y2;
 
end
 
if flag==1,
   fp=THETA2';
end;