Bonjour octave me dit que j'ai une faute après la ligne 3 donc au niveau du function dd1, avez vous une idée ?
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 %ScriptExo2 clear all; close all; clc. %FONCTION function dd1=derivcd1(xd,yd) N=length(xd)-1; for i=2:N dd1(i)=(yd(i+1)-yd(i-1))/(xd(i+1)-xd(i-1)); endfor end %APPLICATION ca=0.024; cb=0.1; va=0.05; ve=0,012; N=24; vb=linspace(0,2*ve,N+1); dd1=derivcd1(xd,yd); Nvb=25; for i=1:Nvb if vb(i)=0 ph(i)=-log10(ca); elseif vb(i)>0, vb(i)<ve ph(i)=-log10((ca*va-cb*vb)/va+vb); elseif vb(i)=ve ph(i)=7.0; elseif vb(i)>ve=0 ph(i)=14+log10((cb*vb-ca*va)/va+vb); endif end %FIGURE figure(1); plot(vb,ph(i));
Partager