Estoy resolviendo la ecuación diferencial
con condiciones iniciales u (0) = u (1) = 0 , u '' (0) = u '' (1) = 0 . Aquí \ sigma (x) \ geqslant \ sigma_ {0}> 0 es el parámetro. En forma de operador podemos reescribir la ecuación diferencial como Au = f , donde el operador A es definitivo positivo.
(σ2(x)u′′(x))′′=f(x),0⩽x⩽1
u(0)=u(1)=0u′′(0)=u′′(1)=0σ(x)⩾σ0>0Au=fA
Siguiendo el esquema FEM, reduzco mi problema a un problema de optimización
J (u) = (Au, u) - 2 (f, u) \ to \ min_ {u}
Introduzco
J(u)=(Au,u)−2(f,u)→minu
elementos finitos
hk(x) como
vk(x)=⎧⎩⎨1−(x−xkh)2,0,x∈[xk−1,xk+1]otherwise
para cualquier
k=1,…,n−1 , donde
xk=hk ,
h=1n . Los elementos finitos
v0(x) y
vn(x) se introducen de manera similar.
Intento encontrar numéricamente el vector α modo que u(x)=∑nk=0αkvk(x) resuelva el problema de optimización. Tenemos
J(u)=∑i=0n∑j=0nαiαj(Avi,vj)−∑i=0n2αi(vi,f)=αTVα−2αTb→minα,
donde
bi=(f,vi) y
Vi,j=(Avi,vj) . Después de la diferenciación con respecto a
α , recibo
Vα=b,
pero aquí la matriz de rigidez
V es singular. ¿Entonces qué tengo que hacer? ¿Tal vez tengo que elegir otros elementos finitos?