2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему На страницу Пред.  1, 2, 3, 4, 5 ... 7  След.
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение29.12.2013, 11:20 
Аватара пользователя


15/01/06
200
1. Куда пропало определение t?
2. Так и не избавились от P(t).

Код:
c = 3*10^8;
H = 7.2*10^5;
\[Tau] = 47*10^-6;
Fd = 7200;
Tp = 50*10^-6;
f = 3.5*10^8;
v = 7500;
Q = 18000;
L = 100;
s2 = 0.003;
s1 = 2;
\[CapitalTheta]0 = 0.018;
m = 256;
\[Lambda] = 0.022;
n = 3;
lm = 25;
d = 1.2;
\[Lambda] = 0.022;
n = 3;
t = List[ -10^-8, -10^-9, 0, 10^-16, 10^-13, 10^-9, 10^-8, 20*10^-8];
p = Sum[NIntegrate[
   Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2)/
       H^2 - (0.00005*n^2*\[Pi]^2*d^2*
        y^2)/(\[CapitalTheta]0^2*\[Lambda]^2*H^2) - (100*lm^2*x^2*
        s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*s1^2)/(H^2*
        L^2*\[Pi]^2*
        s2^2) - (5.55*(y^2 +
          x^2))/(H^2*\[CapitalTheta]0^2) - \[Pi]*(f^2*(t - k*Tp -
            x^2/(c*H) - y^2/(c*H))^2 +
        2*Q*(v*k*Tp)/(H*
            d) + \[Tau]^2*((v*k*Tp)/(H*d))^2)], {x, -Infinity,
    Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40,
   AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255,
   255}]
ListPlot[Transpose[{t, p}]]

У меня этот код работает без ошибок. Правда с таким масштабом по t и p лучше все же, наверное, строить логарифмические графики. И еще одно замечание - я-то привел метод Монте-Карло как тот, который проинтегрирует без ошибок, но с ним надо быть аккуратным, он может для одной и той же функции выдавать различные результаты.

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение29.12.2013, 16:47 


14/10/12
210
Спасибо, получилось, но гарфик корявый.
Как сделать плавную, т.е. непрерывную кривую на графике?
Вообще возможно задание построение графиков вида y=x(t)?

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение29.12.2013, 17:59 
Аватара пользователя


15/01/06
200
Ну какая дискретность, такой в общем-то и график. Сделайте разбиение по t чаще, получите более плавную кривую (наверное). Графики функции конечно же строить можно, например, если переделать в функцию
Код:
pp[tt_] :=
Sum[NIntegrate[
   Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2)/
       H^2 - (0.00005*n^2*\[Pi]^2*d^2*
        y^2)/(\[CapitalTheta]0^2*\[Lambda]^2*H^2) - (100*lm^2*x^2*
        s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*s1^2)/(H^2*
        L^2*\[Pi]^2*
        s2^2) - (5.55*(y^2 +
          x^2))/(H^2*\[CapitalTheta]0^2) - \[Pi]*(f^2*(tt - k*Tp -
            x^2/(c*H) - y^2/(c*H))^2 +
        2*Q*(v*k*Tp)/(H*
            d) + \[Tau]^2*((v*k*Tp)/(H*d))^2)], {x, -Infinity,
    Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40,
   AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255,
   255}]
Plot[pp[t], {t, -10^(-8), 20*10^(-8)}]

Но для такой функции график строиться будет так долго, что не факт, что и дождетесь результата.

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение29.12.2013, 18:19 


14/10/12
210
а как сделать логарифмический масштаб по вертикали?

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение29.12.2013, 18:23 
Аватара пользователя


11/06/12
10390
стихия.вздох.мюсли
LogPlot.

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение29.12.2013, 18:54 


14/10/12
210
Leierkastenmann в сообщении #807641 писал(а):
график строиться будет так долго

как-то можно задать шаг аргумента для графика?

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение29.12.2013, 19:10 
Аватара пользователя


11/06/12
10390
стихия.вздох.мюсли
Опцией PlotPoints.

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение29.12.2013, 21:43 


14/10/12
210
это правильно, что в начале объявлено pp[tt_], а построение идет Plot[pp[t] ?

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение29.12.2013, 21:57 
Аватара пользователя


11/06/12
10390
стихия.вздох.мюсли
А что вас смущает?

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение30.12.2013, 08:02 


14/10/12
210
отличие tt_ от t

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение30.12.2013, 18:54 
Аватара пользователя


11/06/12
10390
стихия.вздох.мюсли
salang, вы, я смотрю, не очень сильны в Mathematica.
Объясняю. Вот мы определяем функцию: f[x_] := x^2.
А вот мы её вычисляем: f[3 a]. Как вы думаете, что получим? Правильно, 9 а^2 ;-) Если в определении стоит, скажем, икс, это не значит, что и в вычислении должен фигурировать икс. Можете, впрочем, в своём коде заменить tt на t или наоборот, всё по-прежнему будет работать (или не работать ;-)
А знак подчёркивания _ в определении функции как раз указывает, что на месте икса может быть что угодно.
Загляните в эти разделы справки: Defining Functions, Functions as Procedures, How to Work with Variables and Functions.

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение04.01.2014, 15:03 


14/10/12
210
PlotPoints задает число выборок, а есть возможность задать абсолютный шаг?
В ListLinePlot не получается вставить PlotPoints. Это ошибка?
Cколько примерно для моих данных занимает расчет непрерывной кривой на средней машине?
Код:
c = 3*10^8;
H = 7.2*10^5;
\[Tau] = 47*10^-6;
Fd = 7200;
Tp = 50*10^-6;
f = 3.5*10^8;
v = 7500;
Q = 18000;
L = 100;
s2 = 0.003;
s1 = 2;
\[CapitalTheta]0 = 0.018;
m = 256;
\[Lambda] = 0.022;
n = 3;
lm = 25;
d = 1.2;
\[Lambda] = 0.022;
n = 3;
t = List[ -10^-8, -10^-9, 0, 10^-11, 10^-10, 10^-9, 10^-8, 20*10^-7];

pp[tt_] = Sum[
  NIntegrate [
   Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2 )/H^2 - (
     0.00005*n^2*\[Pi]^2*d^2*y^2 )/(\[CapitalTheta]0^2*\[Lambda]^2 *H^2) - (
     100*lm^2*x^2*s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*s1^2)/(
     H^2*L^2*\[Pi]^2*s2^2) - (5.55*(y^2 + x^2))/(
     H^2*\[CapitalTheta]0^2) - \[Pi] * (f^2 *(t - k*Tp - x^2/(c*H) - y^2/(
           c*H))^2 +
        2*Q*(v*k*Tp)/(H*d) + \[Tau]^2*((v*k*Tp)/(H*d))^2 )], {x, -Infinity,
    Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40,
   AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255, 255}]


In[1]:= Plot[pp[t], {t, -10^(-8), 20*10^(-8)}, PlotPoints -> 20]
ListLogPlotPlot[Transpose[{t, p}]]

Out[1]= \!\(\*
GraphicsBox[{},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{{-1.*^-8, 2.*^-7}, {0., 0.}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]\)

During evaluation of In[1]:= Transpose::nmtx: The first two levels of the one-dimensional list {t,p} cannot be transposed. >>

Out[2]= ListLogPlotPlot[Transpose[{t, p}]]

Sum[NIntegrate[
  Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2)/
      H^2 - (0.00005*n^2*\[Pi]^2*d^2*y^2)/(\[CapitalTheta]0^2*\[Lambda]^2*
       H^2) - (100*lm^2*x^2*s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*
       s1^2)/(H^2*L^2*\[Pi]^2*
       s2^2) - (5.55*(y^2 +
         x^2))/(H^2*\[CapitalTheta]0^2) - \[Pi]*(f^2*(t - k*Tp - x^2/(c*H) -
           y^2/(c*H))^2 +
       2*Q*(v*k*Tp)/(H*d) + \[Tau]^2*((v*k*Tp)/(H*d))^2)], {x, -Infinity,
   Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40,
  AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255, 255}]

{1.5983*10^-12, 335063., 836068., 936968., 90404., 0., 0., 702430.}

ListLinePlot[%, AxesLabel -> {"t", "P(t)"}, PlotStyle -> PointSize[0.01]]

\!\(\*
GraphicsBox[{{}, {{}, {},
{RGBColor[0.24720000000000017`, 0.24, 0.6], PointSize[0.01],
      LineBox[{{1., 1.598299392712959*^-12}, {2., 335062.66072612564`}, {3.,
       836067.9196629806}, {4., 936968.342828458}, {5., 90403.98299123366}, {
       6., 0.}, {7., 0.}, {8., 702429.7247402478}}]}}, {}},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{
FormBox["\"t\"", TraditionalForm],
FormBox["\"P(t)\"", TraditionalForm]},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{{0, 8.}, {0, 936968.342828458}},
PlotRangeClipping->True,
PlotRangePadding->{{0.16, 0.16}, {18739.36685656916, 18739.36685656916}}]\)

t = List[ -10^-8, -10^-9, 0, 10^-9, 10^-8, 10^-7, 10^-6, 20*10^-5];

Файл: http://zalil.ru/34881989/1310a13d.52c89978/6.nb

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение05.01.2014, 02:43 
Аватара пользователя


11/06/12
10390
стихия.вздох.мюсли
salang в сообщении #809443 писал(а):
есть возможность задать абсолютный шаг?
Используйте дополнительно опцию MaxRecursion -> 0.
salang в сообщении #809443 писал(а):
В ListLinePlot не получается вставить PlotPoints. Это ошибка?
Разумеется, ошибка. ListLinePlot строит график по уже имеющемуся набору точек.
salang в сообщении #809443 писал(а):
Сколько примерно для моих данных занимает расчет непрерывной кривой на средней машине?
А что такое средняя машина?

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение05.01.2014, 19:07 


14/10/12
210
MaxRecursion вставляет без стрелки-указателя :-(
Под средней машиной я имел в виду обычный Prescott 2,8+4Гб RAM
Код:
c = 3*10^8;
H = 7.2*10^5;
\[Tau] = 47*10^-6;
Fd = 7200;
Tp = 50*10^-6;
f = 3.5*10^8;
v = 7500;
Q = 18000;
L = 100;
s2 = 0.003;
s1 = 2;
\[CapitalTheta]0 = 0.018;
m = 256;
\[Lambda] = 0.022;
n = 3;
lm = 25;
d = 1.2;
\[Lambda] = 0.022;
n = 3;
t = List[ -10^-8, -10^-9, 0, 10^-11, 10^-10, 10^-9, 10^-8, 20*10^-7];

In[1]:= pp[tt_] = Sum[
  NIntegrate [
   Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2 )/H^2 - (
     0.00005*n^2*\[Pi]^2*d^2*y^2 )/(\[CapitalTheta]0^2*\[Lambda]^2 *H^2) - (
     100*lm^2*x^2*s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*s1^2)/(
     H^2*L^2*\[Pi]^2*s2^2) - (5.55*(y^2 + x^2))/(
     H^2*\[CapitalTheta]0^2) - \[Pi] * (f^2 *(t - k*Tp - x^2/(c*H) - y^2/(
           c*H))^2 +
        2*Q*(v*k*Tp)/(H*d) + \[Tau]^2*((v*k*Tp)/(H*d))^2 )], {x, -Infinity,
    Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40,
   AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255, 255}]


During evaluation of In[1]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((510 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(65025 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[1]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((508 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(64516 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[1]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((506 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(64009 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[1]:= General::stop: Further output of NIntegrate::inumr will be suppressed during this calculation. >>

Out[1]= 511 NIntegrate[
  Exp[(x^2 (\[Pi] Fd^2 \[Tau]^2 (m - Abs[k])^2))/H^2 - (
    0.00005 n^2 \[Pi]^2 d^2 y^2)/(\[CapitalTheta]0^2 \[Lambda]^2 H^2) - (
    100 lm^2 x^2 s1^2)/(H^2 L^2 \[Pi]^2 s2^2) - (100 lm^2 y^2 s1^2)/(
    H^2 L^2 \[Pi]^2 s2^2) - (5.55 (y^2 + x^2))/(
    H^2 \[CapitalTheta]0^2) - \[Pi] (f^2 (t - k Tp - x^2/(c H) - y^2/(
          c H))^2 + (2 Q (v k Tp))/(
       H d) + \[Tau]^2 ((v k Tp)/(
         H d))^2)], {x, -\[Infinity], \[Infinity]}, {y, -\[Infinity], \
\[Infinity]}, MaxRecursion -> 40, AccuracyGoal -> 60,
  Method -> "AdaptiveMonteCarlo"]

In[2]:= Plot[pp[t], {t, -10^(-8), 20*10^(-8)}, PlotPoints -> 20]
ListLogPlotPlot[Transpose[{t, p}]]

During evaluation of In[2]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-(5.55 (<<1>>))/(H^2 \[CapitalTheta]0^2)-<<1>>/<<1>>-\[Pi] ((2 k Q Tp v)/(d H)+f^2 (<<1>>)^2+(k^2 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))+(Fd^2 \[Pi] x^2 \[Tau]^2 (m-Abs[k])^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[2]:= NIntegrate::inumr: The integrand E^(<<9>>+(3.14159 Fd^2 x^2 \[Tau]^2 (m-1. <<1>>)^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[2]:= NIntegrate::inumr: The integrand E^(<<9>>+(3.14159 Fd^2 x^2 \[Tau]^2 (m-1. <<1>>)^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[2]:= General::stop: Further output of NIntegrate::inumr will be suppressed during this calculation. >>

Out[2]= \!\(\*
GraphicsBox[{},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{{-1.*^-8, 2.*^-7}, {0., 0.}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]\)

During evaluation of In[2]:= Transpose::nmtx: The first two levels of the one-dimensional list {t,p} cannot be transposed. >>

ListLogPlotPlot[Transpose[{t, p}], MaxRecursion]

In[4]:= p[tt_] = Sum[
  NIntegrate[
   Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2)/
       H^2 - (0.00005*n^2*\[Pi]^2*d^2*y^2)/(\[CapitalTheta]0^2*\[Lambda]^2*
        H^2) - (100*lm^2*x^2*s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*
        s1^2)/(H^2*L^2*\[Pi]^2*
        s2^2) - (5.55*(y^2 +
          x^2))/(H^2*\[CapitalTheta]0^2) - \[Pi]*(f^2*(t - k*Tp - x^2/(c*H) -
            y^2/(c*H))^2 +
        2*Q*(v*k*Tp)/(H*d) + \[Tau]^2*((v*k*Tp)/(H*d))^2)], {x, -Infinity,
    Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40,
   AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255, 255}]

During evaluation of In[4]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((510 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(65025 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[4]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((508 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(64516 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[4]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((506 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(64009 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[4]:= General::stop: Further output of NIntegrate::inumr will be suppressed during this calculation. >>

Out[4]= 511 NIntegrate[
  Exp[(x^2 (\[Pi] Fd^2 \[Tau]^2 (m - Abs[k])^2))/H^2 - (
    0.00005 n^2 \[Pi]^2 d^2 y^2)/(\[CapitalTheta]0^2 \[Lambda]^2 H^2) - (
    100 lm^2 x^2 s1^2)/(H^2 L^2 \[Pi]^2 s2^2) - (100 lm^2 y^2 s1^2)/(
    H^2 L^2 \[Pi]^2 s2^2) - (5.55 (y^2 + x^2))/(
    H^2 \[CapitalTheta]0^2) - \[Pi] (f^2 (t - k Tp - x^2/(c H) - y^2/(
          c H))^2 + (2 Q (v k Tp))/(
       H d) + \[Tau]^2 ((v k Tp)/(
         H d))^2)], {x, -\[Infinity], \[Infinity]}, {y, -\[Infinity], \
\[Infinity]}, MaxRecursion -> 40, AccuracyGoal -> 60,
  Method -> "AdaptiveMonteCarlo"]

ListLinePlot[%, AxesLabel -> {"t", "P(t)"}, PlotStyle -> PointSize[0.01]]

NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-(5.55 (<<1>>))/(H^2 \[CapitalTheta]0^2)-<<1>>/<<1>>-\[Pi] ((2 k Q Tp v)/(d H)+f^2 (t+<<3>>)^2+(k^2 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))+(Fd^2 \[Pi] x^2 \[Tau]^2 (m-Abs[k])^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-(5.55 (<<1>>))/(H^2 \[CapitalTheta]0^2)-<<1>>/<<1>>-\[Pi] ((2 k Q Tp v)/(d H)+f^2 (t+<<3>>)^2+(k^2 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))+(Fd^2 \[Pi] x^2 \[Tau]^2 (m-Abs[k])^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-(5.55 (<<1>>))/(H^2 \[CapitalTheta]0^2)-<<1>>/<<1>>-\[Pi] ((2 k Q Tp v)/(d H)+f^2 (t+<<3>>)^2+(k^2 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))+(Fd^2 \[Pi] x^2 \[Tau]^2 (m-Abs[k])^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

General::stop: Further output of NIntegrate::inumr will be suppressed during this calculation. >>

ListLinePlot::lpn: 511 NIntegrate[Exp[(x^2 (\[Pi] Power[<<2>>] Power[<<2>>] Power[<<2>>]))/Power[<<2>>]-(0.00005 Power[<<2>>] Power[<<2>>] Power[<<2>>] Power[<<2>>])/Times[<<3>>]-(100 Power[<<2>>] Power[<<2>>] Power[<<2>>])/Times[<<4>>]-(100 Power[<<2>>] Power[<<2>>] Power[<<2>>])/Times[<<4>>]-(5.55 Plus[<<2>>])/Times[<<2>>]-\[Pi] (Times[<<2>>]+Times[<<3>>]+Times[<<2>>])],{x,-\[Infinity],\[Infinity]},<<3>>,Method->AdaptiveMonteCarlo] is not a list of numbers or pairs of numbers. >>

 Профиль  
                  
 
 Re: численное решение интеграла от экспоненты в Mathematica
Сообщение05.01.2014, 19:52 
Аватара пользователя


11/06/12
10390
стихия.вздох.мюсли
salang в сообщении #809852 писал(а):
MaxRecursion вставляет без стрелки-указателя
Простите, не понял вас.
И прячьте код под спойлеры, что ли, а то как-то раздражает листать эти простыни вверх-вниз.

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 93 ]  На страницу Пред.  1, 2, 3, 4, 5 ... 7  След.

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group