• API
• FAQ
• Tools
• Archive
SHARE
TWEET

# Untitled

a guest Dec 3rd, 2019 81 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. import random
2. import numpy as np
3.
4. def GramSchmidt(*a):
5.     k=len(a[0])
6.     N=len(a)
7.     b = [[0] * k for i in range(N)]
8.     b[0]=a[0]
9.     for i in range(1,N):
10.         sum=a[i]
11.         for j in range(0,i):
12.             scolar_ab=0
13.             scolar_bb=0
14.             proj=[i for i in range(k)]
15.             for n in range(k):
16.                 scolar_ab+=b[j][n]*a[i][n]
17.                 scolar_bb+=b[j][n]*b[j][n]
18.             for n in range(k):
19.                 proj[n]=(scolar_ab/scolar_bb)*b[j][n]
20.             for n in range(k):
21.                 sum[n]-=proj[n]
22.         b[i]=sum
23.     return b
24.
25. l1=np.array([1,-1,1,1])
26. l2=np.array([3,-2,1,2])
27. l3=np.array([-1,-2,1,2])
28. l4=np.array([1,-2,-1,2])
29. print('gs=',GramSchmidt(l1,l2,l3,l4))
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.

Top