MATLAB 程式设计 Learning Linear Algebra.pptVIP

MATLAB 程式設計 Learning Linear Algebra 方煒 台大生機系 Vector Products, Dot and Cross a=[1,2,3];b=[3,2,1]; C=a*b D=a.*b E=dot(a,b) F=cross(a,b) Ex7_1 Solve a Linear System A=[ 1, 0,1;-1, 1,1;1,-1,1 ]; b=[4;4;2]; % now solve for x x=A\b %we obtain [1;2;3] Ex7_2 Max and Min x=0:.01:5; y=x.*exp(-x.^2); % take a look at the function so we know what it looks like plot(x,y) % find the max and min ymin=min(y) ymax=max(y) % find the max and min along with the array indices imax and imin % where they occur [ymin,imin]=min(y) [ymax,imax]=max(y) Ex7_3 Matrix Inverse A=[1,0,-1;-1,1,1;1,-1,1] % load C with the inverse of A C=inv(A) % verify by matrix multiplication % that A*C is the identity matrix A*C Ex7_4 Transpose and Hermitian Conjugate If your matrices are real, then there is no difference between these two commands %find the transpose of the matrix A A.’ %find the Hermitian conjugate of matrix A A’ [1,2,3] [1,2,3]’ [4;5;6] [4;5;6]’ Ex7_5a Special Matrices % eye: % load I with the 4x4 identity matrix (the programmer who invented this % syntax must have been drunk) I=eye(4,4) % zeros: % load Z with a 5x5 matrix full of zeros Z=zeros(5,5) % ones: % load X with a 3x3 matrix full of ones X=ones(3,3) % rand: % load Y with a 4x6 matrix full of random numbers between 0 and 1 % The random numbers are uniformly distributed on [0,1] Y=rand(4,6) % And to load a single random number just use r=rand % randn: % load Y with a 4x6 matrix full of random numbers with a Gaussian % distribution with zero mean and a variance of 1 Y=randn(4,6) Ex7_6 Determinant and Ex7_7 Norm Ex7_6 Determinant %Find the determinant of a square matrix this way det(A) Ex7_7 Norm of Vector (Magnitude) %Matlab will compute the magnitude of a vector a %(the square root of the sum of the squares of its components) %with the norm command a=[1,2,3] norm(a) Ex7_8 Sum the Elements %For arrays the command sum adds up the elements of the array: % calculate the sum of the squares of the reciprocals of the % integers fr