«Неограниченное значение функции» при интегрировании функции, которая должна сходиться - PullRequest
0 голосов
/ 17 марта 2020

Я пытаюсь интегрировать функцию (от 0 до бесконечности), которая, как я знаю, должна сходиться, потому что она имеет форму, подобную t, умножающую отрицательную экспоненту, но я получаю следующую ошибку:

Ошибка в интегрировать (c11t, нижний = 0, верхний = Inf, подразделы = 2000): значение конечной функции

Этот код показывает ошибку:

c11<-function(t){
  if(Tr^2>4*D){
   exp(-Tr*abs(t)/2)*sigma^2*(R*cosh(R*abs(t)/2)*(D+gamma22^2)+Tr*(-D+gamma22^2)*sinh(R*abs(t)/2))/(2*D*R*Tr)  
  }else if(Tr^2<4*D){
    exp(-Tr*abs(t)/2)*sigma^2*(R*cos(R*abs(t)/2)*(D+gamma22^2)+Tr*(-D+gamma22^2)*sin(R*abs(t)/2))/(2*D*R*Tr)  
  }else{
    exp(-Tr*abs(t)/2)*sigma^2*(Tr^2*(2-Tr*abs(t))+4*(2+Tr*abs(t))*gamma22^2)/(4*Tr^3)
  }
}

c11abs<-function(t){
  abs(c11(t))
}
c11t<-function(t){
  t*c11abs(t)
}

    gamma22<--10
    Tr<-1.020408
    D<-0.2040816
    sigma<-0.2
    R<-0.4742429


    t11<-integrate(c11t,lower=0,upper=Inf,subdivisions = 2000)$value/integrate(c11abs,lower=0,upper=Inf,subdivisions = 2000)$value

Это выглядело очень странно мне, и я попытался построить функцию c11t, получая кривую, которая начинается с нуля при t = 0, немного растет, а затем быстро уменьшается экспоненциально, вероятно, до нуля, в то время как c11 затухает как экспонента (и всегда больше, чем ноль) предполагая, что интеграция должна быть возможной. Может кто-нибудь, пожалуйста, помогите мне с этой проблемой? Заранее большое спасибо

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...