Run Code
|
API
|
Code Wall
|
Misc
|
Feedback
|
Login
|
Theme
|
Privacy
|
Patreon
Grad-sistemab
# Resolvendo o sistema linear Au=v pelo método dos gradientes. c=c(8, 1, 1, 0, 2, 20, 9, 2, 1, 1, 2, 1, 11,2, 0, 3, 2, 1, -8, 1, 2, 1, 2, 3, 10) A=matrix(c,5,5,,byrow=TRUE); print("matriz A"); A # Entrada da matriz A v=c(1,2,3,4,5); print("Vetor v"); v # Entrada do vetor v B=t(A)%*%A; print("Matriz A*A");B; b=t(A)%*%v; print("vetor A*v"); b # Adaptando o sistema para A*Au=A*v, ou Bu=b, em que B=A*A é positiva definida e b=A*v. n=1000 u=c(1,1,1,1,1) t=0 # Tempo inicial for ( j in 1:n){ w=B%*%u-b h= t(w)%*%w /(t(B%*%w) %*%w) # Escolha de passo ótimo para o u=u-h[1]*w # Método de Euler t=t+h } print("Aproximação para a solução de Au=v"); u # Aproximação para u(t), sendo t a soma dos passos h. print(" tempo t estimado"); t print(" Teste para erro de aproximação"); A%*%u-v # Verificação do erro de aproximação.
run
|
edit
|
history
|
help
0
plot_cars
Kruskal Wallis test
Ex15-12-20
10-08-2020CalcuNum
Sentiment Analysis - using input field
data frame
List of libraries in R
28-09-2020Int-poli-inter-FuncaoComportadabem
31-08-2020-Exemplo AjusteNaoL
Área sob curvas de funções: Integração numérica utilizando a regra do ponto médio.