Мне нужно найти числовой корень, используя метод деления пополам, и распечатать значения переменных, задействованных на каждой итерации, пока они не достигнут определенного значения.

 bisection <- function(x1, x2){
  l <- vector(mode="integer")
    l[1] <- x1
  r <- vector(mode="integer")
    r[1] <- x2
  m <- vector(mode="integer")
  gl <- vector(mode="integer")
  gr <- vector(mode="integer")
  gm <- vector(mode="integer")
  
  root <- 5e-8  
  i <- 1
  repeat{
    m[i] <- (l[i]+r[i])/2
    gl[i] <- gx(l[i])
    gr[i] <- gx(r[i])
    gm[i] <- gx(m[i])
        
    if (isTRUE(abs(gm[i]) > root) && isTRUE(gl[i]*gm[i] < 0)){
      l[i+1] <- l[i]
      r[i+1] <- m[i]
    }
    if (isTRUE(gm[i] > root) && isTRUE(gr[i]*gm[i] < 0)){
      l[i+1] <- m[i]
      r[i+1] <- r[i]
    }
    else if (isTRUE(abs(gm[i]) <= root)){
      j <- c(0:(length(gm)-1))
      df <- data.frame(j, l,r,m,gl,gr,gm)
      names(df) <- c("i", "xl","xr","xm", "gxl","gxr", "gxm")
      print(df)
      break
    }
  }
}

Когда я пытаюсь запустить эту функцию с bisection(1,1.5), на выходе получается только одна строка итерации, даже если решение для нее вручную приведет как минимум к 12 итерациям. Тоже виснет (?).

Я не знаю, где я ошибаюсь. Пожалуйста помоги.

Отредактировано, чтобы сказать, что функция gx такова: gx <- function(x){x^3-x-1}

r function root bisection
0
Aza Amihan 1 Окт 2020 в 14:03

1 ответ

Лучший ответ

В вашем коде есть две проблемы:

  1. Во втором заявлении if, которое вы забыли abs, оно должно быть if (isTRUE(abs(gm[i]) > root)
  2. Вы забыли увеличить значение счетчика, т.е. i остается равным 1.

Кроме того, я внес небольшую поправку в ваш код. Вместо использования if-if-else-if сначала проверьте, найден ли корень, и, если нет, обновите l и r с помощью одного if-else. Облегчает чтение и понимание.

bisection <- function(x1, x2){
  l <- vector(mode="integer")
  l[1] <- x1
  r <- vector(mode="integer")
  r[1] <- x2
  m <- vector(mode="integer")
  gl <- vector(mode="integer")
  gr <- vector(mode="integer")
  gm <- vector(mode="integer")
  
  root <- 5e-8  
  i <- 1
  repeat{
    m[i] <- (l[i]+r[i])/2
    gl[i] <- gx(l[i])
    gr[i] <- gx(r[i])
    gm[i] <- gx(m[i])
    
    if (isTRUE(abs(gm[i]) <= root)){
      j <- c(0:(length(gm)-1))
      df <- data.frame(j, l,r,m,gl,gr,gm)
      names(df) <- c("i", "xl","xr","xm", "gxl","gxr", "gxm")
      print(df)
      break
    }
    
    if (isTRUE(abs(gm[i]) > root) && isTRUE(gl[i]*gm[i] < 0)){
      l[i+1] <- l[i]
      r[i+1] <- m[i]
    } else {
      l[i+1] <- m[i]
      r[i+1] <- r[i]
    }
    
    i <- i +1
  }
}

gx <- function(x) x^2 - 2

bisection(1, 1.5)
#>     i       xl       xr       xm           gxl          gxr           gxm
#> 1   0 1.000000 1.500000 1.250000 -1.000000e+00 2.500000e-01 -4.375000e-01
#> 2   1 1.250000 1.500000 1.375000 -4.375000e-01 2.500000e-01 -1.093750e-01
#> 3   2 1.375000 1.500000 1.437500 -1.093750e-01 2.500000e-01  6.640625e-02
#> 4   3 1.375000 1.437500 1.406250 -1.093750e-01 6.640625e-02 -2.246094e-02
#> 5   4 1.406250 1.437500 1.421875 -2.246094e-02 6.640625e-02  2.172852e-02
#> 6   5 1.406250 1.421875 1.414062 -2.246094e-02 2.172852e-02 -4.272461e-04
#> 7   6 1.414062 1.421875 1.417969 -4.272461e-04 2.172852e-02  1.063538e-02
#> 8   7 1.414062 1.417969 1.416016 -4.272461e-04 1.063538e-02  5.100250e-03
#> 9   8 1.414062 1.416016 1.415039 -4.272461e-04 5.100250e-03  2.335548e-03
#> 10  9 1.414062 1.415039 1.414551 -4.272461e-04 2.335548e-03  9.539127e-04
#> 11 10 1.414062 1.414551 1.414307 -4.272461e-04 9.539127e-04  2.632737e-04
#> 12 11 1.414062 1.414307 1.414185 -4.272461e-04 2.632737e-04 -8.200109e-05
#> 13 12 1.414185 1.414307 1.414246 -8.200109e-05 2.632737e-04  9.063259e-05
#> 14 13 1.414185 1.414246 1.414215 -8.200109e-05 9.063259e-05  4.314817e-06
#> 15 14 1.414185 1.414215 1.414200 -8.200109e-05 4.314817e-06 -3.884337e-05
#> 16 15 1.414200 1.414215 1.414207 -3.884337e-05 4.314817e-06 -1.726433e-05
#> 17 16 1.414207 1.414215 1.414211 -1.726433e-05 4.314817e-06 -6.474773e-06
#> 18 17 1.414211 1.414215 1.414213 -6.474773e-06 4.314817e-06 -1.079981e-06
#> 19 18 1.414213 1.414215 1.414214 -1.079981e-06 4.314817e-06  1.617417e-06
#> 20 19 1.414213 1.414214 1.414214 -1.079981e-06 1.617417e-06  2.687177e-07
#> 21 20 1.414213 1.414214 1.414213 -1.079981e-06 2.687177e-07 -4.056319e-07
#> 22 21 1.414213 1.414214 1.414214 -4.056319e-07 2.687177e-07 -6.845708e-08
#> 23 22 1.414214 1.414214 1.414214 -6.845708e-08 2.687177e-07  1.001303e-07
#> 24 23 1.414214 1.414214 1.414214 -6.845708e-08 1.001303e-07  1.583661e-08
0
stefan 1 Окт 2020 в 11:42

Похожие вопросы

31 В чем разница между методом и функцией?
18 Бинарный поиск (бисекция) в Python
25 Есть ли способ для некорневых процессов привязаться к «привилегированным» портам в Linux?
28 Определите, работает ли он на корневом устройстве
23 Какая польза от JavaScript-метода bind?
11 Как преобразовать множитель в целое \ числовое без потери информации?
4 Полномиальное нахождение корня пополам
2 Метод бисекции (численный анализ)
21 Как проверить, запущен ли в сценарии bash от имени пользователя root
2 Метод Ньютона в R
1 Метод двумерного деления пополам - поиск корней
1 Как в R накапливать данные из цикла в числовой вектор над циклом
2 Нахождение корня функции методом деления пополам в R
2 Поиск корня VBA через деление пополам
1 Метод деления пополам нахождения корня в R
1 Нахождение квадратного корня из числа методом деления пополам
💵 Получи $100 на хостинг на 60 дней.
Регистрируйся!

Новые вопросы

2 как создать табличку из списка с именованными векторами
3 Как я могу вычесть значения каждого столбца моего фрейма данных из всех других столбцов?
1 Как перебирать многоугольники, которые перекрываются, и случайным образом выбирать точки из каждого в r
2 Можно ли добавить новую строку после нескольких строк таким образом, чтобы вычислялось среднее значение каждого столбца в R?
1 Как я могу использовать один и тот же цвет заливки для одной и той же категории на разных графиках столбцов ggplot?
1 используйте функцию lapply для изменения нескольких матриц в R
1 Как создать максимальную переменную для группы, исключая текущее наблюдение
1 Замените определенные текстовые значения в строке фрейма данных на основе истинного условия в строке перед
1 Использование памяти Polars по сравнению с {data.table}
1 добавить «n» по оси Y в скобках на гистограмме пропорций при просмотре одной категориальной переменной
Все вопросы по теме r >

r

R — это бесплатный язык программирования с открытым исходным кодом и программная среда для статистических вычислений, биоинформатики, визуализации и общих вычислений. Пожалуйста, используйте минимально воспроизводимые примеры, которые другие могут запустить с помощью копирования и вставки. Показать желаемый результат полностью. Используйте dput() для данных и укажите все небазовые пакеты с помощью library(). Не вставляйте изображения для данных или кода, вместо этого используйте блоки кода с отступом. Для вопросов по статистике используйте https://stats.stackexchange.com.

Подробнее про r...