Advertisement
Guest User

Untitled

a guest
Mar 22nd, 2018
72
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.44 KB | None | 0 0
  1. clear all
  2. close all
  3. t=0:0.02:6*pi;
  4. w=1;
  5. A=5;
  6. i_r=[];
  7. i_m=[];
  8. B=3;
  9. y=A*((sin(w*t).^3).*sin(2*w*t).^2).*exp(-B);
  10. yw=y;
  11. for i=2:length(t)-1;
  12. if yw(i)>y(i+1)&yw(i)>y(i-1);
  13. i_r=[i_r i];
  14. elseif yw(i)<y(i+1)&yw(i)<y(i-1);
  15. i_m=[i_m i];
  16. else
  17. end
  18. end
  19. plot(t(i_r),yw(i_r),'om')
  20. ylabel('Amplitude')
  21. xlabel('t')
  22. grid on
  23. hold on
  24. plot(t(i_m),yw(i_m),'<g')
  25. plot(t,yw,'r')
  26. legend({'y_m_a_x','y_m_i_n','y'})
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement