Esempio n. 1
0
 def test_projection_matrix(self):
     M = n.diag(n.arange(1,9))
     Z = n.zeros((8,2))
     Z[0,0] = 1; Z[1,1]=1
     PP = uf.projection_matrix(Z)
     #print M.shape, PP.shape
     M_tilde = n.dot(PP.H,n.dot(M,PP))
     self.assertEqual(M_tilde[0,0],0)
     self.assertEqual(M_tilde[1,1],0)
Esempio n. 2
0
 def test_projection_matrix(self):
     M = n.diag(n.arange(1, 9))
     Z = n.zeros((8, 2))
     Z[0, 0] = 1
     Z[1, 1] = 1
     PP = uf.projection_matrix(Z)
     #print M.shape, PP.shape
     M_tilde = n.dot(PP.H, n.dot(M, PP))
     self.assertEqual(M_tilde[0, 0], 0)
     self.assertEqual(M_tilde[1, 1], 0)
Esempio n. 3
0
 def test_pseudo_inverse(self):
     M = n.diag(n.arange(1,9))
     Z = n.zeros((8,2))
     Z[0,0] = 1; Z[1,1]=1
     PP = uf.projection_matrix(Z)
     M_tilde = n.dot(PP.H,n.dot(M,PP))
     Minv = uf.pseudo_inverse(M,num_remov=2)
     M_id = n.dot(Minv,M_tilde)
     proj_id = n.dot(PP.H,n.dot(n.identity(8),PP))
     self.assertTrue(n.all(M_id==proj_id))
Esempio n. 4
0
 def test_pseudo_inverse(self):
     M = n.diag(n.arange(1, 9))
     Z = n.zeros((8, 2))
     Z[0, 0] = 1
     Z[1, 1] = 1
     PP = uf.projection_matrix(Z)
     M_tilde = n.dot(PP.H, n.dot(M, PP))
     Minv = uf.pseudo_inverse(M, num_remov=2)
     M_id = n.dot(Minv, M_tilde)
     proj_id = n.dot(PP.H, n.dot(n.identity(8), PP))
     self.assertTrue(n.all(M_id == proj_id))