%1 genera matrices especiales, matrices 0, todas las 1 matrices unitarias, matrices de identidad, matrices aleatorias, matrices aleatorias de distribución normal estándar %
zeros();one();eye();rand();randn();
zeros(3 ) ;
ceros(3,2);
A=[1:1:4;1:1:4];
%[1:9];(1:9);1:9;A[1:9]
A( 1 ,2);% Acceso unidimensional a bidimensional, obtenga 1, 1 está organizado en secuencia
B=[1,2,3;1,2,3];
C=A(end-1:end,end- 1: fin);
D=[B([1,3]);A([1,2]),A(7:9)];
ceros(tamaño(A));
x=20+(50-20)*rand(5);%y=a+(ba)x, establece una matriz aleatoria con intervalo de orden k (a, b) y=0.6+sqrt(0.1)*randn(5
);%y=μ+ox, establezca una matriz aleatoria de distribución normal con media μ y varianza de orden 0.1k
% de entero aleatorio y entero de matriz aleatoria
A=fix((90-10+1)*rand(1)+10) ; %rand(1) obtiene un número aleatorio de primer orden, también un número entero aleatorio
% Matriz del cubo de Rubik, cada fila y columna, la suma de las dos diagonales es (n^2+1)n/2
% Llene los números en el rango [101,125] en la matriz de 5x5, la suma es 565, considere 5 El la suma de las matrices de orden es 65, suma 100 a cada número, lo que da 565
M=100+magic(5);
% de matriz de Vandermond, la última columna es 1, la penúltima columna es la matriz especificada y las demás son el producto de la última columna y la penúltima columna
E=vander([1;2;3;5] );
%
formato de matriz de Hilbert rata %producción racional
H=hilb(4);
H=invhilb(4);
formato
%Matriz de Toeplitz , la primera fila y columna las define usted mismo, y cada uno de los demás elementos es igual que el elemento superior izquierdo
toeplitz(1:6)
toeplitz(1:6,1:5)
%Matriz adjunta , la matriz adjunta de polinomios, se compone de vectores de coeficientes polifásicos
%Find x^3-7x+6
p=[1,0,-7,6];
compan(p);
%Matriz de Pascal
Z=pascal(6);
inv(Z)
La matriz %Diagonal extrae el elemento diagonal principal
A=[1:3;4:6];
D1=diag(A); %diag(A,k), la diagonal se puede desplazar hacia arriba o hacia abajo
D=diag ([1, 1,1,1]); % Construya una matriz diagonal, una matriz diagonal de mxm, o puede ir seguida de parámetros (V, k) y el desplazamiento diagonal de otros elementos es 0 diag(diag(A));
% Generar una matriz diagonal con los elementos diagonales principales de A
%Matriz triangular
%Triángulo superior triu(A,k) Triángulo inferior tril(A,k)
%Transposición de matriz
B=[1:3;4:6;7:9];
B1=B.';
B2=B'; %Transposición conjugada, si es un número real, los dos resultados son iguales, equivalente a conjunto (A).'
% Rotación de matriz y giro hacia la izquierda, derecha, arriba y abajo
B3=rot90(B,1); % Girar 90° en sentido antihorario
B4=rot90(B3,-1); % Girar 90° en sentido horario
B5=fliplr(B4); % Voltear izquierda y derecha
B6= flipud(B5);% voltear hacia arriba y hacia abajo, rotar 180 y luego voltear verticalmente, lo que equivale a transponer
%Inversa y pseudoinversa
A=[1:3;4:6];
A1=inv(A); %Encuentra la inversa, es decir, A*A1=A1*A=I
%Usa AX=b para obtener X =A-1b , puedes resolver un sistema de ecuaciones lineales. Es más eficiente
usar
%Determinante, rango, traza
det(A);
rango(A);
traza(A);% es igual a la suma de elementos diagonales o la suma de valores propios, eig, por lo que también puedes juzgar la suma(diag(A) )
%normas vectoriales y matriciales, número de condición de matriz
%norm(V,1) norma(V),norm(V,2) norma(V,inf)vector %norm(
A,1) norma(A) norma(A, inf ) matriz
%cond(A,1) cond(A),cond(A,2) cond(A,inf)%Cuanto más cerca esté el número de condición de 1, mejor será el rendimiento de la matriz
%Valores propios y vectores propios
A=[1:3;4:6;7:9];
res1=eig(A);
[X,D]=eig(A);
[X1,D1]=eig(A, 'nobalance'); %AX=XD,A*inv(D)=I
%Método de valor propio para resolver
la ecuación matriz de coeficientes %p, A=compan(p) y luego eig(A)==roots(p)
% Matriz dispersa
% Modo de almacenamiento completo y modo de almacenamiento disperso, este último almacena (fila, columna) num en columnas, hay dos conceptos %
Uno se refiere a que la matriz tiene más 0 elementos y tiene las características de una matriz dispersa; el otro es Se refiere a almacenar
S = sparse (A) en forma de una matriz dispersa,% se convierte completamente a sparse. Cuando A es sparse, es equivalente a una llamada de asignación S = A. Un vector indefinido sparse (m,
n ) se puede introducir durante la asignación. ;%Generar una matriz dispersa
sparse(u,v,S) en la que todos los elementos de mxn son 0; %Los vectores de igual longitud, las etiquetas de fila, las etiquetas de columna y los elementos distintos de 0 corresponden a matriz dispersa
[u,v,S]=find(A);% Encuentra la posición de los elementos distintos de cero en la matriz A (fila, columna) num full(A);% Devuelve el método de almacenamiento completo correspondiente a la matriz
dispersa matriz A
spconvert(A) %La matriz dispersa se convierte en una matriz de almacenamiento dispersa , donde A está compuesta por la matriz mx3 o mx4 (fila,col) num
spdiags (B,d,row, col);%La matriz original
speye(row,col); %Matriz de unidades de almacenamiento dispersas
% Dos operaciones dispersas aún no están dispersas, las otras están completamente almacenadas