Зачоть / Untitled2

% разложение матрицы using [L,U]=lu(A) L – нижняя треугольная матрицы с 1 на диагонали;
% U — верхняя треугольная матрица
%разложение матрицы using [Q,R]=qr(A) ,where A=matrix
% QQ^T=E(единичная матрица) и R – верхняя треугольная матрица

N=5; % given size
A=create_matrix(N);

L=lu(A)
Q=qr(A)

d=diag(Q)

% get the the diagonal of of Q or L matrix, using diag(A)let Y=diagonal
% X=1:N X=[1;2;3;4;5;6] for 6x6

X=[1;2;3;4;5]
Y=[-12.2882, -1.7567, 1.6779, 3.5061, -5.6695] % from diag(Q or L) copy the numbers

% line through the points
p=polyfit(X,Y,5);% 1 is for first order polynomial.## TRY DIFFERENT NUMBERS FOR HIGHER ORDER i.e 2,3,(N-1)
lin_fit=polyval(p,X);
k=X(1):X(end);
f=interp1(X,Y,k,'pchip'); %try phchip,spline


figure;
hold on
plot(X,Y,'ro',X,lin_fit)
plot(k,f,'b')

%Aproximation use code from line 20 to 24
x=[1 3 5 7]
y=[-5 -2 1 3]

% operations on matrices A./B A.*B
% extract a column A(:,n) role A(n,:) pick a part B=A(3:5,3:5)