Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- box% octave
- GNU Octave, version 5.2.0
- Copyright (C) 2020 John W. Eaton and others.
- This is free software; see the source code for copying conditions.
- There is ABSOLUTELY NO WARRANTY; not even for MERCHANTABILITY or
- FITNESS FOR A PARTICULAR PURPOSE. For details, type 'warranty'.
- Octave was configured for "x86_64-pc-linux-gnu".
- Additional information about Octave is available at https://www.octave.org.
- Please contribute if you find this software useful.
- For more information, visit https://www.octave.org/get-involved.html
- Read https://www.octave.org/bugs.html to learn how to submit bug reports.
- For information about changes from previous versions, type 'news'.
- octave:1> a=exp(1j*2*pi*[0:0.01:1]);
- octave:2> b=1./a;
- octave:3> svd([real(a);imag(a);real(b);imag(b)])
- ans =
- 1.0100e+01
- 1.0000e+01
- 2.6515e-15
- 7.6618e-16
- octave:4> [a1,a2]=meshgrid(a,a);
- octave:5> b = 1./(a1.*a2);
- octave:6> a1_=reshape(a1,1,101^2);
- octave:7> a2_=reshape(a2,1,101^2);
- octave:8> b_=reshape(b,1,101^2);
- octave:9> svd([real(a1_);imag(a1_);real(a2_);imag(a2_);real(b_);imag(b_)])
- ans =
- 72.132
- 71.764
- 71.764
- 71.063
- 71.063
- 70.711
- octave:10>
Add Comment
Please, Sign In to add comment