@joseclaudineiferreira

# Ajuste de curvas por mínimos quadrados. # Dados os vetores x e y, supomos que y_i=f(x_i)+erro_i, como erro_i um erro "pequeno". x=seq(-5,5,by=0.5);x n=length(x) y=c(76.8, 66.3, 56.4, 47.2, 38.9, 31.4, 24.4, 17.84, 11.1, 3.7, -0.2, 4.1, 11.21, 1

# # Link para o Colab: https://colab.research.google.com/drive/1owkOikDypwPgOKm3Tg0ae-UgIhWxIF4U?usp=sharing # # Resolver a equação f(x,y)=(f1(x,y),f2(x,y))=(0,0) ou o sistema de equações # f1(x,y)=0 # f2(x,y)=0. # # Primeiro vamos definir as funçõe

# Gráfico em 3 dimensões # Resolver a equação f(x,y)=(f1(x,y),f2(x,y))=(0,0). f1<-function(u){ # Definição de f1 r=u[1] s=u[2] return(exp(r+s)-2) }

# Gráfico em 3 dimensões # Resolver a equação f(x,y)=(f1(x,y),f2(x,y))=(0,0). f1<-function(u){ # Definição de f1 r=u[1] s=u[2] return((cos(r) + s^3) / sqrt(r)) }

# Gráfico em 3 dimensões # Resolver a equação f(x,y)=(f1(x,y),f2(x,y))=(0,0). f1<-function(u){ # Definição de f1 r=u[1] s=u[2] return(exp(r+s)-2) }

# Gráfico em 3 dimensões # Resolver a equação f(x,y)=(f1(x,y),f2(x,y))=(0,0). f1<-function(u){ # Definição de f1 r=u[1] s=u[2] return(exp(r+s)-2) }

# Gráfico em 3 dimensões # Resolver a equação f(x,y)=(f1(x,y),f2(x,y))=(0,0). f1<-function(u){ # Definição de f1 x=u[1] y=u[2] f1=y*x^3+y^4+2*x*y+x+y+1 f1 }

# Gráfico em 3 dimensões # Curva t=seq(0,2.5,by=0.05) u<-function(s){s*cos(s)+2} v<-function(s){(-s^2)*sin(s)/4+1} curve(u,0,2.5) curve(v,0,2.5)

# Gráfico em 3 dimensões # Curva # Superfície z=f(x,y) x <- seq(-1.5,1.5, by=0.05) y <- seq(-2,2,by=0.05) mf<-function(s,r){z=s+r*1i; Mod(z^6+z^4-z^2+8*z+16)^2# f(x,y); auxiliar para gráfico }

c=c(8, 1, 1, 0, 2, 0, 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); A #Definimos a função que calcula a soma dos elementos da diagonal principal, ou o traço, de uma matriz tr<-function(A){ p=0; n=length(A[

# Gráfico em 3 dimensões # Curva # Superfície z=f(x,y) x <- seq(-1.2,0.8, by=0.05) y <- seq(-1.2,1.2,by=0.05) mf<-function(s,r){z=s+r*1i; Mod(z^4+z^3+z^2+z+2)^2# f(x,y); auxiliar para gráfico }

# Gráfico em 3 dimensões # Curva # Superfície z=f(x,y) x <- seq(-1.2,0.8, by=0.05) y <- seq(-1.2,1.2,by=0.05) mf<-function(s,r){z=s+r*1i; Mod(z^4+z^3+z^2+z+2)^2# f(x,y); auxiliar para gráfico }

# Cáclulo de raiz quadrada de s s=5 r=0; r is=1 # Sei que 2^2<5 e 2^3>5 n=10 d=1 for ( i in -is:n){

# Cáclulo de raiz quadrada de s s=5 r=0; r is=1 # Sei que 2^2<5 e 2^3>5 n=10 d=1 for ( i in -is:n){

# Cordic nutiplicar e dividir: esquema simples x=34567.5895 y=0 z=7324567.2511234 psmult<-function(x,z){ y=0 ix=0

# Função a ser maximizada f<-function(s){ x=s[1];y=s[2] p=x*exp(-x^2-y^2);p } # Para fazer grafico a função não pode ser definida como vetorial x <- seq(-2, 2, length = 40) y <- seq(-2,2, length = 40) g <- function(x,y) {x*exp(-x^2-y^2)} # ad

# Máximo de função f<-function(s){ sin(s*log(cos(s))+1)^2 } hh=0.00001 dfn<-function(s){(f(s+hh)-f(s-hh))/(2*hh)} # aproximação para derivada curve(f,-1.3,1)

# Máximo de função f<-function(s){ sin(s*log(cos(s))+1)^2 } hh=0.00001 dfn<-function(s){(f(s+hh)-f(s-hh))/(2*hh)} # aproximação para derivada curve(f,-1.3,1)

#------------ Interpolação f<-function(s){ 1/(1+s^2) } t=c(2,8,12,9,1,13,-2,-8,-5,-3,-4)# x w=f(t) # y n=length(w)

#------------ Interpolação t=c(2,8,12,9,1,13)# x w=c(3,5,7,3,10,7) # y n=length(w) L<-function(i,s){ p=1 for (j in 1:n){