Discussion :

[paulino7526]

Bonjour la famille,

Svp, j'aimerai realiser un fonction qui permet de tracer un signal rampe(en dans de scie) voir figure ci-dessous.

sous Matlab il y a la fonction sawtooth qui permet de le fait. Donc j'aimerai savoir s'il existe une telle fonction sous fortran.

Exemple sous Matlab :

Code MATLAB :

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
function [t,x1,x2]=ramp
 
format long 
clc;
clear all;
%parametre d'integration
a=0;
b=1 ;
n=1e7;
h=(b-a)/n;
hh=h/2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fonction dent de scie
u=0:h:1;
Vramp=6.0+2.2*sawtooth(320*50*u);
figure(1)
plot(u,Vramp)
xlabel('t ','fontsize',16')
ylabel(' v(t)','fontsize',16')
grid

courbe:

[paulino7526]

Bonjour la famille,

suite au problème ci-dessus, j'ai réalisé un code sous matlab qui donne la fonction ramp.

Code matlab :

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
function [Vramp,u]=ramp_fortran
 
format long 
clc
clear all;
%parametre d'integration
a=0;
b=1 ;
n=1e7;
h=(b-a)/n;
hh=h/2;
T = 400e-6; 
g=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fonction dent de scie
v=1;
 for u=0:h:1
 mm(v)=u;
Vramp(v)=0.011e6*(u-g*T)+3.8;
 if Vramp(v)>= 8.2
     g=g+1;
Vramp(v)=0.011e6*(u-g*T)+3.8;
 else
Vramp(v)=0.011e6*(u-g*T)+3.8;
 end
v=v+1;
 end
 figure
plot(mm,Vramp)
xlabel('E','fontsize',16')
ylabel('v','fontsize',16')
grid
 end

j'aimerai donc obtenir la même fonction sous fortran, et voici mon code ci-dessous:

Code fortran :

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
48
49
50
51
52
53
program ramp
 
	real v,u,h,hh,g,mm,ramp,T
 
 
 
!	interval de variation du parametre de controle
 
 
	integer i,n												   
	print*,'													  '
	print*,'													  '
	print*,'													  '
	print*,'													  '
	print*,'GOD IS LOVE, BE PATIENT ......... FER ..........'
 
        open (50,file='C:\test\b50.dat')
        open (51,file='C:\test\b51.dat')
 
!	==== this bifurcation is done with the parameter delta =====
 
 
            a=0
            b=1
            n=1e7
            h=(b-a)/n
            hh=h/2
            T = 400e-6
            g=0
 
            v=1
 
             u=0
              do while(u<=1)
 
             mm(v)=u
              ramp(v)=0.011e6*(u-g*T)+3.8
              if (ramp(v).ge.8.2) then
               g=g+1
               ramp(v)=0.011e6*(u-g*T)+3.8
               else
               ramp(v)=0.011e6*(u-g*T)+3.8
                endif
               v=v+1;
 
               write(51,*)ramp,u
               print*,ramp,u
               u=u+h
             end do
       	close(50)
	close(51)
 
        end

les erreurs générées sont les suivantes:

E:\SAVE PAULAIN\DATA\00 PHD AND MASTER DSCHANG\PHD\DATA ROGER\ramp.f90:3.24:

 real v,u,h,hh,g,mm,ramp,T
                        1
Error: Symbol 'ramp' at (1) cannot have a type
E:\SAVE PAULAIN\DATA\00 PHD AND MASTER DSCHANG\PHD\DATA ROGER\ramp.f90:36.20:

             mm(v)=u
                    1
Error: Unexpected STATEMENT FUNCTION statement at (1)
E:\SAVE PAULAIN\DATA\00 PHD AND MASTER DSCHANG\PHD\DATA ROGER\ramp.f90:37.18:

              ramp(v)=0.011e6*(u-g*T)+3.8
                  1
Error: 'ramp' at (1) is not a variable
E:\SAVE PAULAIN\DATA\00 PHD AND MASTER DSCHANG\PHD\DATA ROGER\ramp.f90:38.22:

              if (ramp(v).ge.8.2) then
                      1
Error: Symbol at (1) is not appropriate for an expression
E:\SAVE PAULAIN\DATA\00 PHD AND MASTER DSCHANG\PHD\DATA ROGER\ramp.f90:40.19:

               ramp(v)=0.011e6*(u-g*T)+3.8
                   1
Error: 'ramp' at (1) is not a variable
E:\SAVE PAULAIN\DATA\00 PHD AND MASTER DSCHANG\PHD\DATA ROGER\ramp.f90:41.19:

               else
                   1
Error: Unexpected ELSE statement at (1)
E:\SAVE PAULAIN\DATA\00 PHD AND MASTER DSCHANG\PHD\DATA ROGER\ramp.f90:42.19:

               ramp(v)=0.011e6*(u-g*T)+3.8
                   1
Error: 'ramp' at (1) is not a variable
E:\SAVE PAULAIN\DATA\00 PHD AND MASTER DSCHANG\PHD\DATA ROGER\ramp.f90:43.19:

                endif
                   1
Error: Expecting END DO statement at (1)
E:\SAVE PAULAIN\DATA\00 PHD AND MASTER DSCHANG\PHD\DATA ROGER\ramp.f90:46.30:

               write(51,*)ramp,u
                              1
Error: Symbol at (1) is not appropriate for an expression
E:\SAVE PAULAIN\DATA\00 PHD AND MASTER DSCHANG\PHD\DATA ROGER\ramp.f90:47.26:

               print*,ramp,u
                          1
Error: Symbol at (1) is not appropriate for an expression

SVP quelqu'un peut-il m'aidé à corriger ses problèmes? Merci.

[Sylvain Bergeron]

La principale erreur, récurrente, est que tu déclare à la fois un programme ramp, une variable scalaire ramp, puis tu utilises ensuite ramp comme un vecteur...

Tu devrais possiblement appeler ton programme PrgRamp puis déclarer ramp comme un vecteur de dimension suffisante.

Une fois ce problème réglé, le nombre d'erreur devrait être nettement plus faible.

Ah oui, après la ligne program, tu dois insérer la ligne implicite none.

[paulino7526]

Bjr à toutes la famille,

Merci infiniment à toi Sylvain Bergeron pour ton aide.

J'ai corrigé mon code selon tes instructions et j'ai pu avoir une courbe.

cependant elle ne se répète pas comme dans Matlab, voir ci-dessous:




ci-dessous le code qui m'a permit de tracer cette courbe:

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
48
49
50
51
52
53
54
	program PrgRamp
        implicit none
	real v,u,h,hh,g,t,a,b
	real, DIMENSION(10000)::mm
        real, DIMENSION(10000)::r
 
 
 
!	interval de variation du parametre de controle
 
 
	integer i,n												   
	print*,'													  '
	print*,'													  '
	print*,'													  '
	print*,'													  '
	print*,'GOD IS LOVE, BE PATIENT ......... FER ..........'
 
        open (50,file='C:\test\b50.dat')
        open (51,file='C:\test\b51.dat')
 
!	==== this bifurcation is done with the parameter delta =====
 
 
            a=0
            b=1
            n=1e7
            h=(b-a)/n
            hh=h/2
            t = 400e-6
            g=0
 
            v=0
 
             u=0
              do while(u<=1)
 
             mm(v)=u
              r(v)=0.011e6*(u-g*t)+3.8
              if (r(v).ge.8.2) then
               g=g+1
               r(v)=0.011e6*(u-g*t)+3.8
               else
               r(v)=0.011e6*(u-g*t)+3.8
                endif
               write(51,*)mm(v),r(v)
               print*,r,u
               v=v+1
               u=u+h
             end do
       	close(50)
	close(51)
 
        end

Après plusieurs simulations, je n'arrive pas à produire le signal en dent de scie. SVP quelqu'un peut-il m'aider dans ce sens?

Merci d'avant.

[Sylvain Bergeron]

La seule chose que je vois, c'est que dans ton graphique, tu n'affiches pas toutes les données disponibles. L'axe des X arrête à 3.5e-4 alors que la première dent pointe à 4e-4.

Petit commentaire moins important : la variable v devrait être entière et non réelle puisqu'elle sert d'indice...

[paulino7526]

Bonjours tous,

Une fois de plus merci à toi Sylvain Bergeron pour ta disponibilité, j'ai effectué des modifications suite à tes remarques.

Mais toujours est-il que je n'obtient pas l'affichage de tout mon signal sous Matlab.

Je rappel que je compile mon programme sous Fortran et que j'importe les données sous Matlab pour le représentation graphique.

Ci-dessous mon code modifier et la représentation graphique:

code:

Code fortran :

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
48
49
50
51
52
53
54
55
56
	program PrgRamp
        implicit none
	real h,hh,g,t,a,b,n
	real, DIMENSION(10000000)::mm
        real, DIMENSION(10000000)::r
        real, DIMENSION(10000000)::u
 
 
 
!	interval de variation du parametre de controle
 
 
	integer i,v
	print*,'													  '
	print*,'													  '
	print*,'													  '
	print*,'													  '
	print*,'GOD IS LOVE, BE PATIENT ......... FER ..........'
 
        open (50,file='C:\test\b50.dat')
        open (51,file='C:\test\b51.dat')
 
!	==== this bifurcation is done with the parameter delta =====
 
 
            a=0
            b=1
            n=1e7
            h=(b-a)/n
            hh=h/2
            t = 400e-6
            g=0
 
            v=0
             i=0
             u=0
              do while(u(v).lt.1)
 
             mm(v)=u(v)
              r(v)=0.011e6*(u(v)-g*t)+3.8
              if (r(v).ge.8.2) then
               g=g+1
               r(v)=0.011e6*(u(v)-g*t)+3.8
               else
               r(v)=0.011e6*(u(v)-g*t)+3.8
                endif
               write(51,*)mm(v),r(v),i
               print*,r,u
               v=v+1
               u(v)=u(v)+h
               i=i+1
             end do
       	close(50)
	close(51)
 
        end

graphe:




Comme on peut le constater sur le graphe, puis j'augmente la dimension des variables mm, r, et u, puis le nombre d'itération diminuer. pour la figure ci-dessus nous avons pour une dimension de 10 000 000 juste deux itérations.

Tant dit que, en effectuant le même programme sous Matlab avec le même part de 1e-7 j'ai eu 10 000 000 itération dans la boucle for, vois ci-dessous:

Code matlab :

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
function [Vramp,u]=ramp_fortran
 
format long 
clc
clear all;
%parametre d'integration
a=0;
b=1 ;
n=1e7;
h=(b-a)/n;
hh=h/2;
T = 400e-6; 
i=0;
g=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fonction dent de scie
v=1;
 for u=0:h:1
 mm(v)=u;
Vramp(v)=0.011e6*(u-g*T)+3.8;
 if Vramp(v)>= 8.2
     g=g+1;
Vramp(v)=0.011e6*(u-g*T)+3.8;
 else
Vramp(v)=0.011e6*(u-g*T)+3.8;
 end
v=v+1;
i=i+1
 end
 figure(1)
plot(mm,Vramp)
xlabel('E','fontsize',16')
ylabel('v','fontsize',16')
grid
z=length(i)
 end

Ainsi donc, je ne sais pas sur quel paramètre je peux jouer pour obtenir toute la représentation de la fonction ramp. Merci our votre aide.