edo.m

3种伪随机数列发生器的matlab编码以及其各种特性的分析

% Premi鑢e partie du TD4

clear all

disp('Entrez le pas de temps h (valeur initiale sugg閞閑 0.025 : ');
H = input(' ');

Horizon = 4;    %  Dur閑 de la simulation
Vs_0 = 1;       % Condition initiale
Lambda = 10;     % 1/(RC), inverse de la constante de temps du syst鑝e
N = Horizon/H;
T = 0:H:Horizon;
Vs=ones(size(T));
Sigma=H*Lambda;

%-----------------------------------------------------------------------
  
% m閠hode d'Adams-Bashforth (naturellement explicite)
      % initialisation en supposant la sortie nulle pour t<0
      F_n = Lambda*(-Vs(1) + sin(T(1))) + cos(T(1));
      F_nmoinsun = %A completer
      Vs(2) = Vs(1) + (H/2)*(3*F_n - F_nmoinsun);   
      

  for n=2:N
      
      % A COMPLETER
  end

  figure; plot(T,Vs); 
  hold on; plot(T,Vs,'*');
  Tfin = 0:0.001:Horizon;
  plot(Tfin, sin(Tfin) + exp(-Lambda*Tfin),'r');
  title(['r閜onse r閑lle (rouge), Adams-Bashforth (bleu), Sigma=',num2str(Sigma)]);
  hold off;

%-----------------------------------------------------------------------
  
% m閠hode d'Adams-Moulton (mise sous forme explicite)
 
  for n=1:N
     
     % A COMPLETER
   
  end
  
  figure; plot(T,Vs); 
  hold on; plot(T,Vs,'*')
  Tfin = 0:0.001:Horizon;
  hold on; plot(Tfin, sin(Tfin) + exp(-Lambda*Tfin),'r');
  title(['r閜onse r閑lle (rouge), Adams-Moulton (bleu), Sigma=',num2str(Sigma)]);
  hold off;