52 lines
483 B
Matlab
52 lines
483 B
Matlab
% 1
|
|
|
|
v = [2 -1 0 -1 3 -1 -2.5 1]'
|
|
|
|
u = v/sqrt(v'*v)
|
|
|
|
Q = eye(size(v*v')) - 2 * u*u'
|
|
|
|
Q'*Q -eye(length(u))
|
|
Q'-inv(Q)
|
|
inv(Q)-Q
|
|
|
|
|
|
% 2
|
|
a1 = [1 -2 1 0]'
|
|
a2 = [1 2 -1 -2]'
|
|
a3 = [-1 2 -1 1]'
|
|
|
|
A = [a1 a2 a3]
|
|
rref(A)
|
|
rank(A)
|
|
% kolommen zijn l.o
|
|
|
|
b = [-2 -10 5 8]'
|
|
|
|
xhat = inv(A'*A)*A'*b % kleinste kwadraten opl
|
|
|
|
U = rref([A b])
|
|
|
|
% vraag 3
|
|
x = [1;1;1]
|
|
|
|
% x[k] = A^k*x
|
|
A = [0 0 3/4;1 0 -2;0 1 3/2]
|
|
b= [1;3;2]
|
|
c= [0 0 1]
|
|
|
|
T = 20;
|
|
for k = 1:T,
|
|
x = [x A*x(:,k)];
|
|
end
|
|
|
|
plot(0:T, c*x)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|