deval
微分方程式の解の構造体を評価
説明
例
定数遅延をもつ DDE の評価
この例では、dde23
を使用して DDE 方程式 y' = ddex1de(t,y,Z)
を解き、解をプロットします。
dde23
を使用して方程式系を解きます。
sol = dde23(@ddex1de, [1 0.2], @ddex1hist, [0 5]);
区間 [0 5]
内の 100 の点で解を評価します。
x = linspace(0,5); y = deval(sol,x);
解をプロットします。
plot(x,y)
ファン デル ポールの方程式の評価
この例では、ode45
を使用して方程式系 y' = vdp1(t,y)
を解き、解の最初の成分をプロットします。
ode45
を使用して方程式系を解きます。
sol = ode45(@vdp1, [0 20], [2 0]);
区間 [0 20]
内の 100 の点で、解の最初の成分を評価します。
x = linspace(0,20,100); y = deval(sol,x,1);
解をプロットします。
plot(x,y)
内挿関数の勾配
ode23
を使用し、区間 で初期条件 y0 = 0
を使用して、単純な ODE である y' = t^2
の解を求めます。
sol = ode23(@(t,y) t^2, [0 3], 0);
7 点で解を評価します。解の構造体 sol
には、deval
がそれらの点で連続解を生成するために使用する内挿関数が含まれます。deval
に 2 番目の出力引数を指定して、指定点における内挿関数の導関数も返します。
x = linspace(0,3,7); [y,yp] = deval(sol,x)
y = 1×7
0 0.0417 0.3333 1.1250 2.6667 5.2083 9.0000
yp = 1×7
0 0.2500 1.0000 2.2500 4.0000 6.2500 9.0000
入力引数
sol
— 解の構造体
構造体
解の構造体。以下の微分方程式ソルバーのいずれかによって返される構造体として指定します。
初期値問題ソルバー (
ode45
、ode23
、ode113
、ode15s
、ode23s
、ode23t
、ode23tb
、ode15i
)遅延微分方程式ソルバー (
dde23
、ddesd
またはddensd
)境界値問題ソルバー (
bvp4c
またはbvp5c
)
例: sol = ode45(@myode,tspan,y0)
データ型: struct
x
— 評価点
ベクトル
評価点。ベクトルとして指定します。x
は、解の値を求める点を指定します。x
の要素は、元の積分区間 [sol.x(1) sol.x(end)]
内に存在しなければなりません。各インデックス i
について、解 y(:,i)
は x(i)
に対応します。
例: 0:0.1:1
例: [2 3 5 8]
データ型: single
| double
idx
— 返される解の成分
ベクトル
返される解の成分。ベクトルとして指定します。解の特定成分にのみ関心がある場合は、この入力を使用します。
例: y = deval(sol,x,[2 3])
は、解の 2 番目の成分と 3 番目の成分のみを返します。
データ型: single
| double
出力引数
y
— 内挿された解
ベクトル | 行列
内挿された解。ベクトルまたは行列として返されます。y
の行数は、返される解の成分の数と等しくなります。
複数の点がある境界値問題の場合、bvp4c
または bvp5c
で取得される解は、インターフェイスで不連続になることがあります。境界点 xc
では、関数 deval
は、xc
の左右の限界の平均値を返します。限界値を取得するには、x
の値を xc
よりわずかに大きい値または小さい値に設定します。
yp
— sol
により生成された連続解の導関数
ベクトル | 行列
sol
により生成された連続解の導関数。ベクトルまたは行列として返されます。yp
は y
と同じサイズであり、y
の個々の解の点で sol
が使用した内挿関数の勾配を示します。
拡張機能
スレッドベースの環境
MATLAB® の backgroundPool
を使用してバックグラウンドでコードを実行するか、Parallel Computing Toolbox™ の ThreadPool
を使用してコードを高速化します。
この関数はスレッドベースの環境を完全にサポートしています。詳細については、スレッドベースの環境での MATLAB 関数の実行を参照してください。
バージョン履歴
R2006a より前に導入
MATLAB コマンド
次の MATLAB コマンドに対応するリンクがクリックされました。
コマンドを MATLAB コマンド ウィンドウに入力して実行してください。Web ブラウザーは MATLAB コマンドをサポートしていません。
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
You can also select a web site from the following list:
How to Get Best Site Performance
Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)