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>
Comments
0 B
|👍
/👎