JoelSjogren

Untitled

Oct 2nd, 2020
216
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.19 KB | None | 0 0
  1. box% octave
  2. GNU Octave, version 5.2.0
  3. Copyright (C) 2020 John W. Eaton and others.
  4. This is free software; see the source code for copying conditions.
  5. There is ABSOLUTELY NO WARRANTY; not even for MERCHANTABILITY or
  6. FITNESS FOR A PARTICULAR PURPOSE. For details, type 'warranty'.
  7.  
  8. Octave was configured for "x86_64-pc-linux-gnu".
  9.  
  10. Additional information about Octave is available at https://www.octave.org.
  11.  
  12. Please contribute if you find this software useful.
  13. For more information, visit https://www.octave.org/get-involved.html
  14.  
  15. Read https://www.octave.org/bugs.html to learn how to submit bug reports.
  16. For information about changes from previous versions, type 'news'.
  17.  
  18. octave:1> a=exp(1j*2*pi*[0:0.01:1]);
  19. octave:2> b=1./a;
  20. octave:3> svd([real(a);imag(a);real(b);imag(b)])
  21. ans =
  22.  
  23. 1.0100e+01
  24. 1.0000e+01
  25. 2.6515e-15
  26. 7.6618e-16
  27.  
  28. octave:4> [a1,a2]=meshgrid(a,a);
  29. octave:5> b = 1./(a1.*a2);
  30. octave:6> a1_=reshape(a1,1,101^2);
  31. octave:7> a2_=reshape(a2,1,101^2);
  32. octave:8> b_=reshape(b,1,101^2);
  33. octave:9> svd([real(a1_);imag(a1_);real(a2_);imag(a2_);real(b_);imag(b_)])
  34. ans =
  35.  
  36. 72.132
  37. 71.764
  38. 71.764
  39. 71.063
  40. 71.063
  41. 70.711
  42.  
  43. octave:10>
Add Comment
Please, Sign In to add comment