Finite element method: Difference between revisions

Content deleted Content added
No edit summary
No edit summary
Line 142:
Substituting into (**) and expanding, we obtain:
 
<math>
(***) &sum;''a<sub>k</sub>''&psi;(''e<sub>k</sub>'',''e<sub>j</sub>'')=&sum;''b<sub>k</sub>''&int;<sub>T</sub>''e<sub>k</sub>e<sub>j</sub>'' for all j=1,...,n.
\textrm{(***) } \sum a_k \psi\left(e_k, e_j\right) = \sum b_k \int_T e_k e_j, j=1,\dots,n.
</math>
 
There is now the question of how to invent a suitable ''g<sub>n</sub>'' and there are many approaches, depending on how the set E was chosen. If E is chosen to be some Fourier basis, then ''g<sub>n</sub>'' can be obtained as the projection of ''g'' onto the linear span of F, but other approaches are possible.
Line 158 ⟶ 160:
where ''a'' is the column vector (a1,a2,...,an) and ''b'' is the column vector (b1,b2,...,bn). The matrices P and Q are given by (***):
 
:<math>
:P[j,k]=&psi;(''e<sub>k</sub>'',''e<sub>j</sub>'')
\begin{matrix}
:Q[j,k]=&int;<sub>T</sub>''e<sub>k</sub>e<sub>j</sub>''
\mathrm{P}_{jk} &=& \psi\left(e_k, e_j\right)\\
\mathrm{Q}_{jk} &=& \int_T e_k e_j
\end{matrix}
</math>
 
P is called the ''stiffness matrix'' and Q is called the ''mass matrix''.