daily pastebin goal
39%
SHARE
TWEET

Untitled

a guest Dec 16th, 2018 69 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. close all
  2. clear all
  3. load acq.mat
  4.  
  5. %Initialisation variable
  6. w0 = 2*pi*(fin/fs);
  7. w=w0*linspace(1-0.0002,1+0.0002,200);
  8. C=zeros(1,200);
  9. A=zeros(1,200);
  10. phi=zeros(1,200);
  11. variance=zeros(1,200);
  12.  
  13. %Création de la matrice C,A,phi pour chaque valeur de w
  14. for index=1:length(w)
  15.  
  16.     Cj=Erreur(w(index),n,data); %ne pas écraser chacune des valeurs
  17.     % C,A,phi
  18.     a0=Cj(1);
  19.     a1=Cj(2);
  20.     a2=Cj(3);
  21.     C(index)=a0;
  22.     A(index)=sqrt((a1^2)+(a2^2));
  23.     phi(index)=atan(a2/a1);
  24.  
  25.     Vn = Cj(2,:)*sin(w(index)*n)+Cj(3,:)*cos(w(index)*n)+Cj(1,:); % Vn = a1*sin(w*n) + a2*cos(w*n) + a0
  26.     E=Vn-data;
  27.     variance(:,index)=var(E)
  28. end
  29.  
  30.  
  31. %tracés
  32. figure(1)
  33. plot(n,data,'*',n,Vn,'red') %valeurs expérimentales pour chaques échantillons
  34. legend('data','Vn');
  35. title('Variation de erreur en fonction de la pulsation')
  36.  
  37. figure(2)
  38. plot(w,variance) % signal d'erreur
  39. title('signal erreur')
  40. figure(3)
  41.  
  42. subplot(2,2,1)
  43. plot(w,variance) % signal d'erreur
  44. title('variance en fonction de la pulsation')
  45.  
  46. subplot(2,2,2)
  47. plot(w,A)
  48. xlabel('w');
  49. ylabel('A');
  50. title('Tracé de amplitude A en fonction de la pulsation');
  51.  
  52. subplot(2,2,3)
  53. plot(w,C)
  54. xlabel('w');
  55. ylabel('C');
  56. title('Tracé de offset C en fonction de la pulsation');
  57.  
  58. subplot(2,2,4)
  59. plot(w,phi)
  60. xlabel('w');
  61. ylabel('phi');
  62. title('Tracé de la phase phi en fonction de la pulsation');
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top