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
| %ScriptProjet2
clear all; close all; clc.
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.0024
cb=0.1
va=0.05
ve=(ca*va)/cb;
for vb(i)=[0:2(ve)]
vb(i)=vb(i)/10; %pdiscretisation
vb(i)=zeros(N+1,1); %tableau commande zeros
Nvb=vb(i);
%2eme partie
for i=1:vb(i)
if vb(i)=0
ph=-log10(ca);
elseif vb(i)=0<vb<ve
ph=-log10((ca*va-cb*vb)/va+vb);
elseif vb(i)=ve
ph=7.0;
else vb(i)>ve=0
ph=14+log10((cb*vb-ca*va)/va+vb);
endif ph=14.0
end |
Partager