Discussion :

[MaODe]

Bonjour,

Je souhaite tracer un graphique de la valeur d'un polynôme de degré i en fonction de x. Chaque terme du polynôme est précédé d'un coefficient qui est issu d'un tirage aléatoire sans remise. Voici le code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
i <- 5 # nombre de valeurs à générer
coef <- data.frame(sample(1:20, i)) # création d'un tableau contenant les coefficients créés
a <- 1:i # degré des polynomes
polynome <- data.frame(coef, a) # Tableau attribuant à chaque coefficient un degré qui lui est associé
names(polynome)=c("i", "a")
 
x <- 2 #Valeur de x quelconque
Vx <- polynome$i*x^polynome$a # Valeur d'un des facteurs ix^a
Polynome <- data.frame(polynome, Vx) 
 
S <- sum(Vx, na.rm=TRUE) # Valeur du polynôme pour une certaine valeur de x

Le problème est que dès que j'essaye de le réaliser pour une séquence de valeurs x et non plus seulement pour une valeur, et que je le mets dans un graphique, j'ai un message d'erreur qui s'affiche :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
x <- seq(from = -5, to = 5, by = 0.01)
Vx <- polynome$i*x^polynome$a
S <- sum(Vx, na.rm=TRUE)
plot(x, S, type = "l")
 
Warning messages:
1: In x^polynome$a :
  longer object length is not a multiple of shorter object length
2: In polynome$i * x^polynome$a :
  longer object length is not a multiple of shorter object length

Quelqu'un a-t-il une idée de comment je peux arriver à mes fins ?

Merci beaucoup

[faubry]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
# i i degre du polynome
# a = vecteur des degres
# x = vecteur des valeurs
# coef = vecteurs des coefficients
a <- 0:i # constante = degre zero
coef <- sample( 1:20, i + 1) # premier coefficient = constante
x <- seq(from = -5, to = 5, by = 0.01)

S <- colSums( kronecker( t( x), a, "^") * co)