Построение того же уравнения, но другого графика.L - PullRequest
Новая
       118

Построение того же уравнения, но другого графика.L

1 голос
/ 27 сентября 2019

Я попытался построить разложение Тейлора с помощью экспоненциальной функции f (x) = exp (x) при a = 0 и при n = 1.

Сначала я написал уравнение #### n = 1 иподготовил это.Но он не смог уместить линию на графике Exp(x).И я попытался переместить знак «+» в верхнюю строку, как видно из уравнения # n = 1, и это сработало.

Эти уравнения точно такие же, за исключением расположения знака оператора (+).

В чем разница между уравнениями (#### n = 1, # n = 1) для построения?

Важно ли расположение знаков операторов при построении? 

f1 <- function(x) exp( x )
x <- seq( -1, 1, by = 0.025 )


a <- 0

#### n = 1
f1.10 <- function( x ){
  exp( a ) / factorial( 0 ) * ( x - a )^0 
  + exp( a ) / factorial( 1 ) * ( x - a )^1
}


# n = 1
f1.1 <- function( x ){
  exp( a ) / factorial( 0 ) * ( x - a )^0 +
    exp( a ) / factorial( 1 ) * ( x - a )^1
}



plot( x, exp(x), ylab = "exp(x)", type = "l", lwd = 3 )
curve( f1.1, -0.5, 0.5, add = T, type = "l", lwd = 2, col = "yellow" )
curve( f1.10, -0.5, 0.5, add = T, type = "l", lwd = 2, col = "blue" )

[sample_output] 1

Ответы [ 2 ]

0 голосов
/ 27 сентября 2019

Две функции дают разные результаты, и они не одинаковы.

f1.10(seq(10))
#[1]  1  2  3  4  5  6  7  8  9 10
f1.1(seq(10))
#[1]  2  3  4  5  6  7  8  9 10 11

Когда n равно 1 f1.10, производим 1 и f1.1 производим 2, или мы можем даже произвести случайное число между -0.5 и 0.5 следующим образом, чтобы увидеть, если онипроизвести тот же результат.

set.seed(777)
x <- runif(100, min = -0.5, max = 0.5)
min(f1.10(x))
#[1] -0.4940781
max(f1.10(x))
#[1] 0.4950499
min(f1.1(x))
#[1] 0.5059219
max(f1.1(x))
#[1] 1.49505

Вы можете видеть, что f1.10 ограничено между -0.5 и 0.5, где f1.1 - нет.

0 голосов
/ 27 сентября 2019

В R, + заголовка означает новую строку.Если вы поставите + в конце строки, R распознает следующую строку в скрипте как продолжение текущей строки.
Таким образом, 

f1.10 <- function( x ){
  exp( a ) / factorial( 0 ) * ( x - a )^0 
  + exp( a ) / factorial( 1 ) * ( x - a )^1
}

совпадает с 

f1.11 <- function( x ){
  exp( a ) / factorial( 1 ) * ( x - a )^1
}
...