0

私はD方程式を解こうとしていますが、わかりませんがy[0]、知っていy[x1]=y1ます。

関連するxrangeでのみDSolveを解決したいx=[x1, infinitny]

どのように機能しますか?

動作しない例を添付

dsolv2 = DSolve[{y'[x] == c*0.5*t12[x, l2]^2 - alpha*y[x], y[twhenrcomesin] == zwhenrcomesin, x >= twhenrcomesin}, y[x], x]
dsolv2 = Flatten[dsolv2]
zsecondphase[x_] = y[x] /. dsolv2[[1]]

DSolveが不等式条件を許可しないことは承知していますが、私が探しているものを説明するためにそれを入れました(t12[x,l2]l2がわかっているので、xにのみ依存する値を与えます)。

編集

t12[j24_, lambda242_] := (cinv1 - cinv2)/(cop2 - cop1 + (h2*lambda242)*E^(p*j24));
cinv1 = 30; cinv2 = 4; cinv3 = 3; h2 = 1.4; h3 = 1.2; alpha = 0.04; z = 50; p = 0.06; cop1 = 0; cop2 = 1; cop3 = 1.3; teta2 = 0.19; teta3 =0.1; co2 = -0.6; z0 = 10;l2 = 0.1;
4

2 に答える 2

1

方程式は一次で線形であるため、非常に一般的な解を得ることができます。

generic = DSolve[{y'[x] == f[x] - alpha*y[x], y[x0] == y0}, y[x], x]

次に、特定の用語に置き換えることができます:

c = 1;
x0 = 1;
y0 = 1;
solution[x_] = generic[[1, 1, 2]] /. {f[x_] -> c*0.5*t12[x, l2]^2}   


Plot[solution[x], {x, x0, 100}]

例

于 2012-06-25T10:20:31.460 に答える
0

この例の何が問題になっていますか?

t12[x_] := Exp[-x .01] Sin[x];
dsolv2 = Chop@DSolve[{y'[x] == c*0.5*t12[x]^2 - alpha*y[x], y[1] == 1}, y[x], x];
Plot[y[x] /. dsolv2[[1]] /. {alpha -> 1, c -> 1}, {x, 1, 100}, PlotRange -> Full]

ここに画像の説明を入力してください

編集

あなたの解説について:

ドメインを制限するために区分的関数を使用してみてください。

t12[x_] := Piecewise[{{ Exp[-x .01] Sin[x], x >= 1}, {Indeterminate, True}}] ;
dsolv2 = Chop@DSolve[{y'[x] == c*0.5*t12[x]^2 - alpha*y[x], y[1] == 1}, y[x], x];
Plot[y[x] /. dsolv2[[1]] /. {alpha -> 1, c -> 1}, {x, 1, 100}, PlotRange -> Full]
于 2012-06-24T19:12:10.293 に答える