не удается построить графики в Mathematica

[salang]

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]

Если несколько функций, их надо в списке ставить:
Plot[{ParabolicCylinderD[-1.5, t], ParabolicCylinderD[-0.5, t]}, {t, 0, 2./10^7}]

[salang]

Vince Diesel в сообщении #921939 писал(а):

Если несколько функций, их надо в списке ставить:$Plot[{ParabolicCylinderD[-1.5, t], ParabolicCylinderD[-0.5, t]}, {t, 0, 2./10^7}]$

я вроде так и сделал

[Aritaborian]

Нет, вы сделали не так. Сравните

Plot[f[x], g[x], {x, 0, 1}]

и

Plot[{f[x], g[x]}, {x, 0, 1}]

[salang]

переделал: 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]

У аргументов функций скобки болжны быть квадратными.

[salang]

исправил: 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]

Встроенные функции, в том числе Exp, пишутся с большой буквы. И вообще, я смотрю до первой ошибки, а вам стоит проверить всеь код

[salang]

если бы знал как- то проверил бы уже
Исправил, но оба предыдущих вопроса остались:

[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]

1. Уберите PlotRange. 2. Можно, рисуя по одному. Или явно указав цвета опцией PlotStyle->{Hue[0],Hue[0.5]}. Первый цвет Hue[0] — красный.

[Aritaborian]

2. Для некоторых цветов есть свои удобные названия:

Red, Green, Blue, Cyan, Magenta, Yellow, Black, White, ...

[salang]

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]

Да. Вместо p1[t]=... должно быть p1[t_]=...

[Aritaborian]

Может быть, всё-таки p1[t_]:=...?

[salang]

Vince Diesel в сообщении #922678 писал(а):

Вместо p1[t]=... должно быть p1[t_]=...

при определении функции или в Plot ?
Как вводится знак $:=$ ? В base typesetting не нашел