2014 dxdy logo

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

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

Сообщения без ответов | Активные темы | Избранное


» » »




Начать новую тему Ответить на тему На страницу 1, 2  След.
  Печатать страницу | Печатать всю тему Пред. тема | След. тема 
salang 
 не удается построить графики в Mathematica
Сообщение22.10.2014, 13:24 


14/10/12
210
http://zalil.su/?fg=486354
Код:
Clear[c, t, h, eta, r, k0, lambda, del, d, zeta, d0, thetaint, \
thetastart, theta, thetaw, kinc]
s = 1

ns = 10^(-9);
c = 2.99792458 10^8;
h = 720000;
r = 6380000;
eta = 1 + h/r;
lambda = c/( 13.575 10^9);
k0 = 2 Pi/lambda;
del = 7200./18182;
zeta = Pi/(64 k0 del);
d0 = 4096
bandwidth = 320000000.;
res = 1/bandwidth;
sigma = Sqrt[(s/2 c)^2 + 0.513 res^2]
gammabar = 0.012215368000378016`;
gammahat = 0.0381925958945466`;
gamma1 = Sqrt[2/(2/gammabar^2 + 2/gammahat^2)];
gamma2 = Sqrt[2/(2/gammabar^2 - 2/gammahat^2)];

beta = Pi/2;
baseline = 1.1676;
gainsqr[roe_, thetaw_] :=
Exp[-2 ( ((roe Cos[thetaw])/gamma1)^2 + ((roe Sin[thetaw])/gamma2)^2)]

pulse[t_] := If[t == 0, 1, (Sin[Pi t/res]/(Pi t/res))^2]

rough0[t_, sigma_] := 1/(Sqrt[2 Pi] sigma) Exp[-(1/2) (t/sigma)^2]

knrange = 31;
knmid = 3;

istart = -50;
iend = 180;
zetab = 500/h;
npoints = (iend - istart)/0.1;

nsigma = 25;
sigmaint = 0.10;
icre = 1.0;

1

4096

1.49896*10^8


p1 (t) = exp[-t/(4 sigma^2)] ParabolicCylinderD[-0.5, t] exp[-0.0074 t]
Plot[p1 (t), {t, 0, 2.0 10^(-7)}]


Set::write: Tag Times in p1 t is Protected. >>

exp[-0.0074 t] exp[-1.11265*10^-17 t] ParabolicCylinderD[-0.5, t]

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


p2 (t) = Sqrt[
   2 sigma/res] (exp[-0.25 t^2/sigma^2] ParabolicCylinderD[-1.5,
     t/sigma])
Plot[p2 (t), {t, 0, 2.0 10^(-7)}]


Set::write: Tag Times in p2 t is Protected. >>

3.09731*10^8 exp[-1.11265*10^-17 t^2] ParabolicCylinderD[-1.5,
  6.67128*10^-9 t]

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

Plot[ParabolicCylinderD[-1.5, t],
ParabolicCylinderD[-0.5, t], {t, 0, 2.0 10^(-7)}]


Plot::nonopt: Options expected (instead of {t,0,2./10^7}) beyond position 3 in Plot[ParabolicCylinderD[-1.5,t],ParabolicCylinderD[-0.5,t],{t,0,2./10^7}]. An option must be a rule or a list of rules. >>

Plot[ParabolicCylinderD[-1.5, t],
ParabolicCylinderD[-0.5, t], {t, 0, 2./10^7}]
 Профиль  
                  
Vince Diesel 
 Re: не удается построить графики в Mathematica
Сообщение22.10.2014, 17:01 
Заслуженный участник


25/02/11
1797
Если несколько функций, их надо в списке ставить:
Plot[{ParabolicCylinderD[-1.5, t], ParabolicCylinderD[-0.5, t]}, {t, 0, 2./10^7}]
 Профиль  
                  
salang 
 Re: не удается построить графики в Mathematica
Сообщение22.10.2014, 19:58 


14/10/12
210
Vince Diesel в сообщении #921939 писал(а):
Если несколько функций, их надо в списке ставить:$Plot[{ParabolicCylinderD[-1.5, t], ParabolicCylinderD[-0.5, t]}, {t, 0, 2./10^7}]$

я вроде так и сделал
 Профиль  
                  
Aritaborian 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 00:32 
Аватара пользователя


11/06/12
10390
стихия.вздох.мюсли
Нет, вы сделали не так. Сравните
Код:
Plot[f[x], g[x], {x, 0, 1}]
и
Код:
Plot[{f[x], g[x]}, {x, 0, 1}]
 Профиль  
                  
salang 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 08:05 


14/10/12
210
переделал: http://zalil.su/?fg=420543. Но 2 первых графика не строит. Для одной функции фигурные скобки ведь не нужны? А на третьем графике непонятно какой цвет относится к какой функции
Код:
[math]$In[1]:= Clear[c, t, h, eta, r, k0, lambda, del, d, zeta, d0, \
thetaint, thetastart, theta, thetaw, kinc]
s = 1

ns = 10^(-9);
c = 2.99792458 10^8;
h = 720000;
r = 6380000;
eta = 1 + h/r;
lambda = c/( 13.575 10^9);
k0 = 2 Pi/lambda;
del = 7200./18182;
zeta = Pi/(64 k0 del);
d0 = 4096
bandwidth = 320000000.;
res = 1/bandwidth;
sigma = Sqrt[(s/2 c)^2 + 0.513 res^2]
gammabar = 0.012215368000378016`;
gammahat = 0.0381925958945466`;
gamma1 = Sqrt[2/(2/gammabar^2 + 2/gammahat^2)];
gamma2 = Sqrt[2/(2/gammabar^2 - 2/gammahat^2)];

beta = Pi/2;
baseline = 1.1676;
gainsqr[roe_, thetaw_] :=
Exp[-2 ( ((roe Cos[thetaw])/gamma1)^2 + ((roe Sin[thetaw])/gamma2)^2)]

pulse[t_] := If[t == 0, 1, (Sin[Pi t/res]/(Pi t/res))^2]

rough0[t_, sigma_] := 1/(Sqrt[2 Pi] sigma) Exp[-(1/2) (t/sigma)^2]

knrange = 31;
knmid = 3;

istart = -50;
iend = 180;
zetab = 500/h;
npoints = (iend - istart)/0.1;

nsigma = 25;
sigmaint = 0.10;
icre = 1.0;

In[34]:= 
p1 (t) = exp[-t/(4 sigma^2)] ParabolicCylinderD[-0.5, t] exp[-0.0074 t]
Plot[p1 (t), {t, 0, 2.0 10^(-7)}]


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

In[36]:= 
p2 (t) = Sqrt[
   2 sigma/res] (exp[-0.25 t^2/sigma^2] ParabolicCylinderD[-1.5,
     t/sigma])
Plot[p2 (t), {t, 0, 2.0 10^(-7)}]


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

In[41]:= Plot[{ParabolicCylinderD[-1.5, t],
  ParabolicCylinderD[-0.5, t]}, {t, 0, 20 10^(-7)},
ColorFunction -> Automatic]


Out[41]= \!\(\*
GraphicsBox[{{}, {},
{Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwdzH0w0wEcx/EdrlSjzpp0UdappITiyPKNPHXqLqeiQhGjEuUoXYUfMTGO
GWI1T/PY5iGRh9n6FTmlEo4Ll1BqUS7LzEXKtz8+977XPx/G+SseLDUKhcJc
HvaOd3OF/e/tdgNjGS1cvuJgej/PjU2zAgvJXBPaX6SeStCcIS3LpxHtliMR
R9NOgpOrcT1axtlgd5vGggbx0yq0/veI4aC5a+B8Y6gAvZTwur/AnYCWpdFM
dEiUe7QjJxFqxcZstOjBVLePNQe8vH1uoD/m6+i6OmTA5/ZbIegeTUkefwsP
jMySfNFxRvrS5J/Z8JnsO4beFHN0QoObC4laVEf05nDNT7aX7sPHM+ZW6NXF
vIlFnXzYGGBvjPZ+PuVo6lUI7zK5m9Asoi92/7MisG97q4UO/LBVtoYmhPRq
FQW9YaVNMDOuBPL41F8Zy36YvHO0SV4KjOrzE2gLUbz3onU5sG+XvEeH5hU5
CIQVUH6k9xVasN3YRENVCSX9U1J0PBEZ3RkkAj+/mRq0yz5B4VC3GPJ0bIrR
fxNMO+oMq0EceykbfVl4x1OeVgMxP+4moYNqC2fDw2thly7vJpot+VGao/cI
PIJGQtGUKgF18s0jYDbS/dDK0LhIw6g6WHxucxztRO/kJu55DES3qzNa7dti
Dnd82eq51mhmKktApNSDp6LPBP2WbnE6374BFGML+ujhMHPdgskG8E9duw59
atZA4sJ9Am5nA9XRZn5hdAOXRlgwL1WmL5u4r7g6MN0IsbIeOXq3dDBOVtoE
XUcnh9Ae1bO8ZPdmkGdPv0HLu/6kLaiagVhh+Qz9vvui3YOqFhiPYj1GHxKu
mt7iJYH1JfFl6BW+1BO92q1A6+Hkom0d/xTtamuF1G0DKf//5hkzwSFS6LGk
xqBnsiqVvQwZfHHacxUdSWEtGL6WQefaAwHo4Fztb/SIp8BQ43iiz3C0OvTX
k3BE+eIwGupM9QovkHBQ+J2JLv7K1gyJJCGinGGLXmkwOm9FkFAm8rRBv0vK
HOy6R4J2PWmJtpJOvcwRkuDQpNyL5iucWvxrSLjeamKBXtqR/3C3hAQRec4M
HeCr4qs6SBhpzzJF/wO4s+x0
"]]},
{Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData["
1:eJwdzH0wEwAcxvFdrYzU1da8HDKFK3lJznnLz1t0nLqlsoqulJFQXoa6Q/Oe
pqtR3se8lcy7SixrveEabXSVKMKtw8pbNFdUfv3x3Pc+/zwGZy/5MNcQCASb
1WFT/NqqnX8ZO5Jos45jMT+cAgRrs9gUd2hQOb8P7ZUrrE2gHINjMyP2aBFH
0zGewoRykdQGrfsteijoZwy8jNhtif6T2vOulM4Gkkxpgg6Noye4cdJg/Mi4
IVpQrJD623DAU+G1DT1SQtY44HILpBY7tNB9JGFBoX4ODGfEbUEnGep2ZM7e
gXSiQg2tk+gtJ3LzgW4uJaK3RZLG7S8UQXmoxZ/RVauV58iXySVA9SUsof2e
K9zMGHygMWfn0Ez226t2z8rgUYa/Ah34ebtoA6UCQhl2crSmim2wQ1Il1Nhn
D6NrMnd9eTxRBSmuagNoS0Gy37LNPYhXzvehwwvKXHgV1bDylCFB84x3mhCV
94FvbfYKncxmJXQHCUAiShShPax4/EFpLXDDZlrRf1PNOptp9RDXJWtCh1Wk
+E7caAC5tbkAHdTIX4iMbATW4eVKdLrwe1WuVhPwuJMlaEIdT32qtwkOyX3y
0YvhSSxaXDPQvppmo/dTu7lp5i0Q/yuZg14zuZzLHWuBSKulVLRDFpPHvv4A
0rWHE9FvqJYnSpwfghHF+TJ66OIejdKphzDeRo5CH1/QE3pwH4FDGjEMbXHm
IlXPoxV0A6OYaHbRfMT76VbQn6WfRpt2fEwSVT2GmpuC42if+oWcTHob+H6K
9UFPSFZu/Fa2QeHBAG/0gDTEsbiuHUL6u9zRrhWq0/oMITCIpU7o9afUj/Zv
egIyz3lbtL3bStnuF09gtCVw7/+/JYO54NAOIEsOmKLnbt9f7DcQQe/iXSM0
i8D8TesRQXVvrD46OH/TJDX6Kei9ZmqjT3I2dupuFYPOoIyMhmYzLf55Mdjl
1aqjVfS+LFmzxfBBuqiGlmVkf5TkiYFFclJFF87vbw9oEMNml2sq6HOnlIXK
TjHUX+lbh/4HwoHOow==
"]]}},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 1.162},
Method->{},
PlotRange->{{0, 2.*^-6}, {1.1627342014801836`, 1.2162802142574964`}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]\)$[/math]
 Профиль  
                  
Vince Diesel 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 08:22 
Заслуженный участник


25/02/11
1797
У аргументов функций скобки болжны быть квадратными.
 Профиль  
                  
salang 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 09:22 


14/10/12
210
исправил: http://zalil.su/?fg=112234
На первых двух графиках сообщений об ошибках теперь нет, но и самих кривых тоже нет
Как-то возможно узнать какой цвет соответствует какой функции на последнем графике?
Код:
[math]$Clear[c, t, h, eta, r, k0, lambda, del, d, zeta, d0, thetaint, \
thetastart, theta, thetaw, kinc]
s = 1

ns = 10^(-9);
c = 2.99792458 10^8;
h = 720000;
r = 6380000;
eta = 1 + h/r;
lambda = c/( 13.575 10^9);
k0 = 2 Pi/lambda;
del = 7200./18182;
zeta = Pi/(64 k0 del);
d0 = 4096
bandwidth = 320000000.;
res = 1/bandwidth;
sigma = Sqrt[(s/2 c)^2 + 0.513 res^2]
gammabar = 0.012215368000378016`;
gammahat = 0.0381925958945466`;
gamma1 = Sqrt[2/(2/gammabar^2 + 2/gammahat^2)];
gamma2 = Sqrt[2/(2/gammabar^2 - 2/gammahat^2)];

beta = Pi/2;
baseline = 1.1676;
gainsqr[roe_, thetaw_] :=
Exp[-2 ( ((roe Cos[thetaw])/gamma1)^2 + ((roe Sin[thetaw])/gamma2)^2)]

pulse[t_] := If[t == 0, 1, (Sin[Pi t/res]/(Pi t/res))^2]

rough0[t_, sigma_] := 1/(Sqrt[2 Pi] sigma) Exp[-(1/2) (t/sigma)^2]

knrange = 31;
knmid = 3;

istart = -50;
iend = 180;
zetab = 500/h;
npoints = (iend - istart)/0.1;

nsigma = 25;
sigmaint = 0.10;
icre = 1.0;

1

4096

1.49896*10^8


p1[t] = exp[-t/(4 sigma^2)] ParabolicCylinderD[-0.5, t] exp[-0.0074 t]
Plot[p1[t], {t, 0, 2.0 10^(-7)}, PlotRange -> {0, 3 10^(-8)}]


exp[-0.0074 t] exp[-1.11265*10^-17 t] ParabolicCylinderD[-0.5, t]

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


p2[t] = Sqrt[
   2 sigma/res] (exp[-0.25 t^2/sigma^2] ParabolicCylinderD[-1.5,
     t/sigma])
Plot[p2[t], {t, 0, 2.0 10^(-7)}, PlotRange -> {0, 3 10^(-8)}]


3.09731*10^8 exp[-1.11265*10^-17 t^2] ParabolicCylinderD[-1.5,
  6.67128*10^-9 t]

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

Plot[{ParabolicCylinderD[-1.5, t], ParabolicCylinderD[-0.5, t]}, {t,
  0, 2 10^(-6)}, ColorFunction -> Automatic, PlotRange -> {0, 1.4}]


\!\(\*
GraphicsBox[{{}, {},
{Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwdzH0w0wEcx/EdrlSjzpp0UdappITiyPKNPHXqLqeiQhGjEuUoXYUfMTGO
GWI1T/PY5iGRh9n6FTmlEo4Ll1BqUS7LzEXKtz8+977XPx/G+SseLDUKhcJc
HvaOd3OF/e/tdgNjGS1cvuJgej/PjU2zAgvJXBPaX6SeStCcIS3LpxHtliMR
R9NOgpOrcT1axtlgd5vGggbx0yq0/veI4aC5a+B8Y6gAvZTwur/AnYCWpdFM
dEiUe7QjJxFqxcZstOjBVLePNQe8vH1uoD/m6+i6OmTA5/ZbIegeTUkefwsP
jMySfNFxRvrS5J/Z8JnsO4beFHN0QoObC4laVEf05nDNT7aX7sPHM+ZW6NXF
vIlFnXzYGGBvjPZ+PuVo6lUI7zK5m9Asoi92/7MisG97q4UO/LBVtoYmhPRq
FQW9YaVNMDOuBPL41F8Zy36YvHO0SV4KjOrzE2gLUbz3onU5sG+XvEeH5hU5
CIQVUH6k9xVasN3YRENVCSX9U1J0PBEZ3RkkAj+/mRq0yz5B4VC3GPJ0bIrR
fxNMO+oMq0EceykbfVl4x1OeVgMxP+4moYNqC2fDw2thly7vJpot+VGao/cI
PIJGQtGUKgF18s0jYDbS/dDK0LhIw6g6WHxucxztRO/kJu55DES3qzNa7dti
Dnd82eq51mhmKktApNSDp6LPBP2WbnE6374BFGML+ujhMHPdgskG8E9duw59
atZA4sJ9Am5nA9XRZn5hdAOXRlgwL1WmL5u4r7g6MN0IsbIeOXq3dDBOVtoE
XUcnh9Ae1bO8ZPdmkGdPv0HLu/6kLaiagVhh+Qz9vvui3YOqFhiPYj1GHxKu
mt7iJYH1JfFl6BW+1BO92q1A6+Hkom0d/xTtamuF1G0DKf//5hkzwSFS6LGk
xqBnsiqVvQwZfHHacxUdSWEtGL6WQefaAwHo4Fztb/SIp8BQ43iiz3C0OvTX
k3BE+eIwGupM9QovkHBQ+J2JLv7K1gyJJCGinGGLXmkwOm9FkFAm8rRBv0vK
HOy6R4J2PWmJtpJOvcwRkuDQpNyL5iucWvxrSLjeamKBXtqR/3C3hAQRec4M
HeCr4qs6SBhpzzJF/wO4s+x0
"]]},
{Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData["
1:eJwdzH0wEwAcxvFdrYzU1da8HDKFK3lJznnLz1t0nLqlsoqulJFQXoa6Q/Oe
pqtR3se8lcy7SixrveEabXSVKMKtw8pbNFdUfv3x3Pc+/zwGZy/5MNcQCASb
1WFT/NqqnX8ZO5Jos45jMT+cAgRrs9gUd2hQOb8P7ZUrrE2gHINjMyP2aBFH
0zGewoRykdQGrfsteijoZwy8jNhtif6T2vOulM4Gkkxpgg6Noye4cdJg/Mi4
IVpQrJD623DAU+G1DT1SQtY44HILpBY7tNB9JGFBoX4ODGfEbUEnGep2ZM7e
gXSiQg2tk+gtJ3LzgW4uJaK3RZLG7S8UQXmoxZ/RVauV58iXySVA9SUsof2e
K9zMGHygMWfn0Ez226t2z8rgUYa/Ah34ebtoA6UCQhl2crSmim2wQ1Il1Nhn
D6NrMnd9eTxRBSmuagNoS0Gy37LNPYhXzvehwwvKXHgV1bDylCFB84x3mhCV
94FvbfYKncxmJXQHCUAiShShPax4/EFpLXDDZlrRf1PNOptp9RDXJWtCh1Wk
+E7caAC5tbkAHdTIX4iMbATW4eVKdLrwe1WuVhPwuJMlaEIdT32qtwkOyX3y
0YvhSSxaXDPQvppmo/dTu7lp5i0Q/yuZg14zuZzLHWuBSKulVLRDFpPHvv4A
0rWHE9FvqJYnSpwfghHF+TJ66OIejdKphzDeRo5CH1/QE3pwH4FDGjEMbXHm
IlXPoxV0A6OYaHbRfMT76VbQn6WfRpt2fEwSVT2GmpuC42if+oWcTHob+H6K
9UFPSFZu/Fa2QeHBAG/0gDTEsbiuHUL6u9zRrhWq0/oMITCIpU7o9afUj/Zv
egIyz3lbtL3bStnuF09gtCVw7/+/JYO54NAOIEsOmKLnbt9f7DcQQe/iXSM0
i8D8TesRQXVvrD46OH/TJDX6Kei9ZmqjT3I2dupuFYPOoIyMhmYzLf55Mdjl
1aqjVfS+LFmzxfBBuqiGlmVkf5TkiYFFclJFF87vbw9oEMNml2sq6HOnlIXK
TjHUX+lbh/4HwoHOow==
"]]}},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{{0, 2.*^-6}, {0, 1.4}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02], Automatic}]\)$[/math]
 Профиль  
                  
Vince Diesel 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 09:50 
Заслуженный участник


25/02/11
1797
Встроенные функции, в том числе Exp, пишутся с большой буквы. И вообще, я смотрю до первой ошибки, а вам стоит проверить всеь код :-)
 Профиль  
                  
salang 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 11:37 


14/10/12
210
если бы знал как- то проверил бы уже :-(
Исправил, но оба предыдущих вопроса остались:
Код:
[math]$Clear[c, t, h, eta, r, k0, lambda, del, d, zeta, d0, thetaint, \
thetastart, theta, thetaw, kinc]
s = 1

ns = 10^(-9);
c = 2.99792458 10^8;
h = 720000;
r = 6380000;
eta = 1 + h/r;
lambda = c/( 13.575 10^9);
k0 = 2 Pi/lambda;
del = 7200./18182;
zeta = Pi/(64 k0 del);
d0 = 4096
bandwidth = 320000000.;
res = 1/bandwidth;
sigma = Sqrt[(s/2 c)^2 + 0.513 res^2]
gammabar = 0.012215368000378016`;
gammahat = 0.0381925958945466`;
gamma1 = Sqrt[2/(2/gammabar^2 + 2/gammahat^2)];
gamma2 = Sqrt[2/(2/gammabar^2 - 2/gammahat^2)];

beta = Pi/2;
baseline = 1.1676;
gainsqr[roe_, thetaw_] :=
Exp[-2 ( ((roe Cos[thetaw])/gamma1)^2 + ((roe Sin[thetaw])/gamma2)^2)]

pulse[t_] := If[t == 0, 1, (Sin[Pi t/res]/(Pi t/res))^2]

rough0[t_, sigma_] := 1/(Sqrt[2 Pi] sigma) Exp[-(1/2) (t/sigma)^2]

knrange = 31;
knmid = 3;

istart = -50;
iend = 180;
zetab = 500/h;
npoints = (iend - istart)/0.1;

nsigma = 25;
sigmaint = 0.10;
icre = 1.0;

1

4096

1.49896*10^8


p1[t] = Exp[-t/(4 sigma^2)] ParabolicCylinderD[-0.5, t] Exp[-0.0074 t]
Plot[p1[t], {t, 0, 2.0 10^(-7)}, PlotRange -> {0, 1}]


E^(-0.0074 t) ParabolicCylinderD[-0.5, t]

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


p2[t] = Sqrt[
   2 sigma/res] (Exp[-0.25 t^2/sigma^2] ParabolicCylinderD[-1.5,
     t/sigma])
Plot[p2[t], {t, 0, 2.0 10^(-7)}, PlotRange -> {0, 3 10^(-8)}]


3.09731*10^8 E^(-1.11265*10^-17 t^2)
  ParabolicCylinderD[-1.5, 6.67128*10^-9 t]

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

Plot[{ParabolicCylinderD[-1.5, t], ParabolicCylinderD[-0.5, t]}, {t,
  0, 2 10^(-6)}, ColorFunction -> Automatic, PlotRange -> {0, 1.4}]


\!\(\*
GraphicsBox[{{}, {},
{Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwdzH0w0wEcx/EdrlSjzpp0UdappITiyPKNPHXqLqeiQhGjEuUoXYUfMTGO
GWI1T/PY5iGRh9n6FTmlEo4Ll1BqUS7LzEXKtz8+977XPx/G+SseLDUKhcJc
HvaOd3OF/e/tdgNjGS1cvuJgej/PjU2zAgvJXBPaX6SeStCcIS3LpxHtliMR
R9NOgpOrcT1axtlgd5vGggbx0yq0/veI4aC5a+B8Y6gAvZTwur/AnYCWpdFM
dEiUe7QjJxFqxcZstOjBVLePNQe8vH1uoD/m6+i6OmTA5/ZbIegeTUkefwsP
jMySfNFxRvrS5J/Z8JnsO4beFHN0QoObC4laVEf05nDNT7aX7sPHM+ZW6NXF
vIlFnXzYGGBvjPZ+PuVo6lUI7zK5m9Asoi92/7MisG97q4UO/LBVtoYmhPRq
FQW9YaVNMDOuBPL41F8Zy36YvHO0SV4KjOrzE2gLUbz3onU5sG+XvEeH5hU5
CIQVUH6k9xVasN3YRENVCSX9U1J0PBEZ3RkkAj+/mRq0yz5B4VC3GPJ0bIrR
fxNMO+oMq0EceykbfVl4x1OeVgMxP+4moYNqC2fDw2thly7vJpot+VGao/cI
PIJGQtGUKgF18s0jYDbS/dDK0LhIw6g6WHxucxztRO/kJu55DES3qzNa7dti
Dnd82eq51mhmKktApNSDp6LPBP2WbnE6374BFGML+ujhMHPdgskG8E9duw59
atZA4sJ9Am5nA9XRZn5hdAOXRlgwL1WmL5u4r7g6MN0IsbIeOXq3dDBOVtoE
XUcnh9Ae1bO8ZPdmkGdPv0HLu/6kLaiagVhh+Qz9vvui3YOqFhiPYj1GHxKu
mt7iJYH1JfFl6BW+1BO92q1A6+Hkom0d/xTtamuF1G0DKf//5hkzwSFS6LGk
xqBnsiqVvQwZfHHacxUdSWEtGL6WQefaAwHo4Fztb/SIp8BQ43iiz3C0OvTX
k3BE+eIwGupM9QovkHBQ+J2JLv7K1gyJJCGinGGLXmkwOm9FkFAm8rRBv0vK
HOy6R4J2PWmJtpJOvcwRkuDQpNyL5iucWvxrSLjeamKBXtqR/3C3hAQRec4M
HeCr4qs6SBhpzzJF/wO4s+x0
"]]},
{Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData["
1:eJwdzH0wEwAcxvFdrYzU1da8HDKFK3lJznnLz1t0nLqlsoqulJFQXoa6Q/Oe
pqtR3se8lcy7SixrveEabXSVKMKtw8pbNFdUfv3x3Pc+/zwGZy/5MNcQCASb
1WFT/NqqnX8ZO5Jos45jMT+cAgRrs9gUd2hQOb8P7ZUrrE2gHINjMyP2aBFH
0zGewoRykdQGrfsteijoZwy8jNhtif6T2vOulM4Gkkxpgg6Noye4cdJg/Mi4
IVpQrJD623DAU+G1DT1SQtY44HILpBY7tNB9JGFBoX4ODGfEbUEnGep2ZM7e
gXSiQg2tk+gtJ3LzgW4uJaK3RZLG7S8UQXmoxZ/RVauV58iXySVA9SUsof2e
K9zMGHygMWfn0Ez226t2z8rgUYa/Ah34ebtoA6UCQhl2crSmim2wQ1Il1Nhn
D6NrMnd9eTxRBSmuagNoS0Gy37LNPYhXzvehwwvKXHgV1bDylCFB84x3mhCV
94FvbfYKncxmJXQHCUAiShShPax4/EFpLXDDZlrRf1PNOptp9RDXJWtCh1Wk
+E7caAC5tbkAHdTIX4iMbATW4eVKdLrwe1WuVhPwuJMlaEIdT32qtwkOyX3y
0YvhSSxaXDPQvppmo/dTu7lp5i0Q/yuZg14zuZzLHWuBSKulVLRDFpPHvv4A
0rWHE9FvqJYnSpwfghHF+TJ66OIejdKphzDeRo5CH1/QE3pwH4FDGjEMbXHm
IlXPoxV0A6OYaHbRfMT76VbQn6WfRpt2fEwSVT2GmpuC42if+oWcTHob+H6K
9UFPSFZu/Fa2QeHBAG/0gDTEsbiuHUL6u9zRrhWq0/oMITCIpU7o9afUj/Zv
egIyz3lbtL3bStnuF09gtCVw7/+/JYO54NAOIEsOmKLnbt9f7DcQQe/iXSM0
i8D8TesRQXVvrD46OH/TJDX6Kei9ZmqjT3I2dupuFYPOoIyMhmYzLf55Mdjl
1aqjVfS+LFmzxfBBuqiGlmVkf5TkiYFFclJFF87vbw9oEMNml2sq6HOnlIXK
TjHUX+lbh/4HwoHOow==
"]]}},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{{0, 2.*^-6}, {0, 1.4}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02], Automatic}]\)$[/math]
 Профиль  
                  
Vince Diesel 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 16:49 
Заслуженный участник


25/02/11
1797
1. Уберите PlotRange. 2. Можно, рисуя по одному. Или явно указав цвета опцией PlotStyle->{Hue[0],Hue[0.5]}. Первый цвет Hue[0] — красный.
 Профиль  
                  
Aritaborian 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 17:42 
Аватара пользователя


11/06/12
10390
стихия.вздох.мюсли
2. Для некоторых цветов есть свои удобные названия:
Код:
Red, Green, Blue, Cyan, Magenta, Yellow, Black, White, ...
 Профиль  
                  
salang 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 20:17 


14/10/12
210
c цветами я правильно понял, что первая функция красным цветом?
C графиками так и не получилось
Код:
[math]$Clear[c, t, h, eta, r, k0, lambda, del, d, zeta, d0, thetaint, \
thetastart, theta, thetaw, kinc]
s = 1

ns = 10^(-9);
c = 2.99792458 10^8;
h = 720000;
r = 6380000;
eta = 1 + h/r;
lambda = c/( 13.575 10^9);
k0 = 2 Pi/lambda;
del = 7200./18182;
zeta = Pi/(64 k0 del);
d0 = 4096
bandwidth = 320000000.;
res = 1/bandwidth;
sigma = Sqrt[(s/2 c)^2 + 0.513 res^2]
gammabar = 0.012215368000378016`;
gammahat = 0.0381925958945466`;
gamma1 = Sqrt[2/(2/gammabar^2 + 2/gammahat^2)];
gamma2 = Sqrt[2/(2/gammabar^2 - 2/gammahat^2)];

beta = Pi/2;
baseline = 1.1676;
gainsqr[roe_, thetaw_] :=
Exp[-2 ( ((roe Cos[thetaw])/gamma1)^2 + ((roe Sin[thetaw])/gamma2)^2)]

pulse[t_] := If[t == 0, 1, (Sin[Pi t/res]/(Pi t/res))^2]

rough0[t_, sigma_] := 1/(Sqrt[2 Pi] sigma) Exp[-(1/2) (t/sigma)^2]

knrange = 31;
knmid = 3;

istart = -50;
iend = 180;
zetab = 500/h;
npoints = (iend - istart)/0.1;

nsigma = 25;
sigmaint = 0.10;
icre = 1.0;

1

4096

1.49896229`*^8
PlotRange -> {0, 3 10^(-8)}

In[5]:= 
p1[t] = Exp[-t/(4 sigma^2)] ParabolicCylinderD[-0.5, t] Exp[-0.0074 t]
Plot[p1[t], {t, 0, 2.0 10^(-7)}]


Out[5]= E^(-0.0074 t - t/(4 sigma^2)) ParabolicCylinderD[-0.5, t]

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

In[3]:= 
p2[t] = Sqrt[
   2 sigma/res] (Exp[-0.25 t^2/sigma^2] ParabolicCylinderD[-1.5,
     t/sigma])
Plot[p2[t], {t, 0, 2.0 10^(-7)}]


Out[3]= Sqrt[2] E^(-((0.25 t^2)/sigma^2)) Sqrt[sigma/res]
  ParabolicCylinderD[-1.5, t/sigma]

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

In[9]:= Plot[{ParabolicCylinderD[-1.5, t],
  ParabolicCylinderD[-0.5, t]}, {t, 0, 2 10^(-6)},
PlotStyle -> {Red, Black}, PlotRange -> {0, 1.3}]


Out[9]= \!\(\*
GraphicsBox[{{}, {},
{RGBColor[1, 0, 0], LineBox[CompressedData["
1:eJwdzH0w0wEcx/EdrlSjzpp0UdappITiyPKNPHXqLqeiQhGjEuUoXYUfMTGO
GWI1T/PY5iGRh9n6FTmlEo4Ll1BqUS7LzEXKtz8+977XPx/G+SseLDUKhcJc
HvaOd3OF/e/tdgNjGS1cvuJgej/PjU2zAgvJXBPaX6SeStCcIS3LpxHtliMR
R9NOgpOrcT1axtlgd5vGggbx0yq0/veI4aC5a+B8Y6gAvZTwur/AnYCWpdFM
dEiUe7QjJxFqxcZstOjBVLePNQe8vH1uoD/m6+i6OmTA5/ZbIegeTUkefwsP
jMySfNFxRvrS5J/Z8JnsO4beFHN0QoObC4laVEf05nDNT7aX7sPHM+ZW6NXF
vIlFnXzYGGBvjPZ+PuVo6lUI7zK5m9Asoi92/7MisG97q4UO/LBVtoYmhPRq
FQW9YaVNMDOuBPL41F8Zy36YvHO0SV4KjOrzE2gLUbz3onU5sG+XvEeH5hU5
CIQVUH6k9xVasN3YRENVCSX9U1J0PBEZ3RkkAj+/mRq0yz5B4VC3GPJ0bIrR
fxNMO+oMq0EceykbfVl4x1OeVgMxP+4moYNqC2fDw2thly7vJpot+VGao/cI
PIJGQtGUKgF18s0jYDbS/dDK0LhIw6g6WHxucxztRO/kJu55DES3qzNa7dti
Dnd82eq51mhmKktApNSDp6LPBP2WbnE6374BFGML+ujhMHPdgskG8E9duw59
atZA4sJ9Am5nA9XRZn5hdAOXRlgwL1WmL5u4r7g6MN0IsbIeOXq3dDBOVtoE
XUcnh9Ae1bO8ZPdmkGdPv0HLu/6kLaiagVhh+Qz9vvui3YOqFhiPYj1GHxKu
mt7iJYH1JfFl6BW+1BO92q1A6+Hkom0d/xTtamuF1G0DKf//5hkzwSFS6LGk
xqBnsiqVvQwZfHHacxUdSWEtGL6WQefaAwHo4Fztb/SIp8BQ43iiz3C0OvTX
k3BE+eIwGupM9QovkHBQ+J2JLv7K1gyJJCGinGGLXmkwOm9FkFAm8rRBv0vK
HOy6R4J2PWmJtpJOvcwRkuDQpNyL5iucWvxrSLjeamKBXtqR/3C3hAQRec4M
HeCr4qs6SBhpzzJF/wO4s+x0
"]]},
{GrayLevel[0], LineBox[CompressedData["
1:eJwdzH0wEwAcxvFdrYzU1da8HDKFK3lJznnLz1t0nLqlsoqulJFQXoa6Q/Oe
pqtR3se8lcy7SixrveEabXSVKMKtw8pbNFdUfv3x3Pc+/zwGZy/5MNcQCASb
1WFT/NqqnX8ZO5Jos45jMT+cAgRrs9gUd2hQOb8P7ZUrrE2gHINjMyP2aBFH
0zGewoRykdQGrfsteijoZwy8jNhtif6T2vOulM4Gkkxpgg6Noye4cdJg/Mi4
IVpQrJD623DAU+G1DT1SQtY44HILpBY7tNB9JGFBoX4ODGfEbUEnGep2ZM7e
gXSiQg2tk+gtJ3LzgW4uJaK3RZLG7S8UQXmoxZ/RVauV58iXySVA9SUsof2e
K9zMGHygMWfn0Ez226t2z8rgUYa/Ah34ebtoA6UCQhl2crSmim2wQ1Il1Nhn
D6NrMnd9eTxRBSmuagNoS0Gy37LNPYhXzvehwwvKXHgV1bDylCFB84x3mhCV
94FvbfYKncxmJXQHCUAiShShPax4/EFpLXDDZlrRf1PNOptp9RDXJWtCh1Wk
+E7caAC5tbkAHdTIX4iMbATW4eVKdLrwe1WuVhPwuJMlaEIdT32qtwkOyX3y
0YvhSSxaXDPQvppmo/dTu7lp5i0Q/yuZg14zuZzLHWuBSKulVLRDFpPHvv4A
0rWHE9FvqJYnSpwfghHF+TJ66OIejdKphzDeRo5CH1/QE3pwH4FDGjEMbXHm
IlXPoxV0A6OYaHbRfMT76VbQn6WfRpt2fEwSVT2GmpuC42if+oWcTHob+H6K
9UFPSFZu/Fa2QeHBAG/0gDTEsbiuHUL6u9zRrhWq0/oMITCIpU7o9afUj/Zv
egIyz3lbtL3bStnuF09gtCVw7/+/JYO54NAOIEsOmKLnbt9f7DcQQe/iXSM0
i8D8TesRQXVvrD46OH/TJDX6Kei9ZmqjT3I2dupuFYPOoIyMhmYzLf55Mdjl
1aqjVfS+LFmzxfBBuqiGlmVkf5TkiYFFclJFF87vbw9oEMNml2sq6HOnlIXK
TjHUX+lbh/4HwoHOow==
"]]}},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{{0, 2.*^-6}, {0, 1.3}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02], Automatic}]\)$[/math]
 Профиль  
                  
Vince Diesel 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 20:35 
Заслуженный участник


25/02/11
1797
Да. Вместо p1[t]=... должно быть p1[t_]=...
 Профиль  
                  
Aritaborian 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 21:20 
Аватара пользователя


11/06/12
10390
стихия.вздох.мюсли
Может быть, всё-таки p1[t_]:=...?
 Профиль  
                  
salang 
 Re: не удается построить графики в Mathematica
Сообщение24.10.2014, 21:32 


14/10/12
210
Vince Diesel в сообщении #922678 писал(а):
Вместо p1[t]=... должно быть p1[t_]=...

при определении функции или в Plot ?
Как вводится знак $:=$ ? В base typesetting не нашел
 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  Страница 1 из 2
  [ Сообщений: 23 ]  На страницу 1, 2  След.

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

Список форумов » Тематические обсуждения » Computer Science » Околонаучный софт


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

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

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

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