Facebook
From lud, 3 Years ago, written in Matlab M-file.
Embed
Download Paste or View Raw
Hits: 45
  1.  
  2. plot(sin(x))
  3. t=linspace(0,1e-5,10000);
  4.  
  5. Harm = 50;
  6.  
  7. A=0.3;
  8. T=1/10^6;
  9. tr=10*10^-9;
  10. tau=T-2*tr;
  11. k=1/(3*10^7);
  12. m=5*10^-7;
  13. wo=2*pi/T;
  14.  
  15. a=(k*tr^2)/(2);
  16. b=A*((tr+tau)-tr);
  17. c=(((tr+tau)-(T/2))*(k*((tr-tau)+(T/2)-2*m)))/2;
  18.  
  19. a0=(1/T)*(a+b+c);
  20.  
  21. f = a0/2;
  22.  
  23. for n = 1:2:Harm
  24.     a1 = (k*sin(n*wo*tr))/(n*wo);
  25.     b1 = 0;
  26.     c1 = 0;
  27.     f = f + 2/T*(a1+b1+c1);
  28. end
  29.  
  30. plot(t,f)
  31.  
  32.  
  33.  
  34.