Wolfram mathematica 使用分段绘制Mathematica的问题
我试图在我的一个问题中使用分段绘图,我有两个变量:x和psi。但是,各功能仅在定义的“x”范围内有效,且psi范围相同。我正在尝试绘制这些的3D图——我基本上只有Plot3D[p,{x,0,1},{psi,0.01,1}]-->这些范围是整个绘图范围,我的各个函数的x范围已经在分段函数中定义了 我得到了以下错误:说这个阴谋::排除:。。。。。。必须是等式或\实值函数的列表 谁能帮我一下吗。我尝试遵循以下相同的程序: 但是,我不知道如何处理绘图部分 谢谢 以下是我的代码:Wolfram mathematica 使用分段绘制Mathematica的问题,wolfram-mathematica,plot,piecewise,Wolfram Mathematica,Plot,Piecewise,我试图在我的一个问题中使用分段绘图,我有两个变量:x和psi。但是,各功能仅在定义的“x”范围内有效,且psi范围相同。我正在尝试绘制这些的3D图——我基本上只有Plot3D[p,{x,0,1},{psi,0.01,1}]-->这些范围是整个绘图范围,我的各个函数的x范围已经在分段函数中定义了 我得到了以下错误:说这个阴谋::排除:。。。。。。必须是等式或\实值函数的列表 谁能帮我一下吗。我尝试遵循以下相同的程序: 但是,我不知道如何处理绘图部分 谢谢 以下是我的代码: j = 10; s =
j = 10;
s = 0; r = 0;
K[x_, psi_] :=
Sum[Sin[n*Pi*x]*
Sin[n*Pi*
psi]*(2*Exp[-(n*Pi)^2*
Abs[s + r]] - (Exp[-(n*Pi)^2*Abs[s - r]] -
Exp[-(n*Pi)^2*(s + r)])/(n*Pi)^2 ), {n, 1, j}];
TL[x_, psi_] = Integrate[K[x - y, psi]*y, {y, -10, 10}];
TU[x_, psi_] = Integrate[K[x - y, psi]*(1 - y), {y, -10, 10}];
eq = {TL[x, psi], TU[x, psi]};
cond = {{0 <= x <= 0.5, 0.01 <= psi <= 1}, {0.5 < x <= 1,
0.01 <= psi <= 1}};
p = Piecewise[{eq, cond}];
Plot3D[p, {x, 0, 1}, {psi, 0.01, 1}]
j=10;
s=0;r=0;
K[x,psi]:=
和[Sin[n*Pi*x]*
Sin[n*Pi*
磅/平方英寸]*(2*Exp[-(n*Pi)^2*
Abs[s+r]-(Exp[-(n*Pi)^2*Abs[s-r]]-
Exp[-(n*Pi)^2*(s+r)]/(n*Pi)^2,{n,1,j}];
TL[x_,psi_]=积分[K[x-y,psi]*y,{y,-10,10}];
TU[x_,psi_]=积分[K[x-y,psi]*(1-y),{y,-10,10}];
eq={TL[x,psi],TU[x,psi]};
cond={{0这是一个工作版本:
time = AbsoluteTime[];
j = 10; s = 0; r = 0;
K[x_, psi_] :=
Sum[Sin[n*Pi*x]*Sin[n*Pi*psi]*
(2*Exp[-(n*Pi)^2*Abs[s + r]] -
(Exp[-(n*Pi)^2*Abs[s - r]] -
Exp[-(n*Pi)^2*(s + r)])/(n*Pi)^2), {n, 1, j}];
TL[x_, psi_] := Integrate[K[x - y, psi]*y, {y, -10, 10}];
TU[x_, psi_] := Integrate[K[x - y, psi]*(1 - y), {y, -10, 10}];
Plot3D[Piecewise[
{{TL[x, psi], 0 <= x <= 0.5}, {TU[x, psi], 0.5 < x <= 1}}],
{x, 0, 1}, {psi, 0.01, 1}]
ToString[Round[AbsoluteTime[] - time]] <> " seconds"
time=绝对时间[];
j=10;s=0;r=0;
K[x,psi]:=
和[Sin[n*Pi*x]*Sin[n*Pi*psi]*
(2*Exp[-(n*Pi)^2*Abs[s+r]]-
(Exp[-(n*Pi)^2*Abs[s-r]]-
Exp[-(n*Pi)^2*(s+r)]/(n*Pi)^2,{n,1,j}];
TL[x,psi]:=积分[K[x-y,psi]*y,{y,-10,10}];
TU[x,psi]:=积分[K[x-y,psi]*(1-y),{y,-10,10}];
Plot3D[分段][
{{TL[x,psi],0非常感谢。你太棒了