ode45

3熱

2答えて

matlabでODE45などの可変時間ステップソルバーを使用している場合、出力の時間幅を定義することになります。つまり、times = [0 50]となります。しかし、RではODEに結果を返す時間点を定義する必要があります。つまり、times = 0:50を指定すると、0,1,2, ... 50に51個の結果が返されます。他の賢明な私は、times = seq(0,50,0.1)のようなシーケン

1熱

1答えて

Ode45 in MATLAB

私はMATLABには新しく、ode45()関数を使用して2つの関数を統合しようとしています。私のコード： [w,x] = ode45(@(w,x) (TSII - TFII) * w , [0 2], 0); TSIIとTFIIは私の2つの関数の値を含む列ベクトル1x151です。 wはまだ価値がありませんが、私はそれをゼロの長さ1x151でいっぱいの列ベクトルにしようとしましたが、それは私の問

1熱

1答えて

私は微分方程式（ローレンツ方程式）のシステムを解決するためにMATLABを使用する方法を学習し、tの関数としての各ソリューションをプロットしようとしていますODE45

使って常微分方程式のシステムを解く X’ = −σx + σy Y’ = ρx − y − xz Z’ = −βz + xy σ = 10、β = 8/3 、およびρ = 28、ならびにx(0) = −8,y(0) = 8およびz(0) = 27である。ここでは、私が使用していたコードです： function xprime = example(t,x) sig = 10; beta

0熱

1答えて

OctaveでODE45シミュレーションでNaN och Infのみを取得するのはなぜですか？

私は、ODE方程式のシステムをシミュレートするためにODE45でOctaveを使用しています。しかし問題は、ODEシミュレーションが間違った値を与えることです。このオクターブのコードを見て：これは与え function dx = dynamik(t, x) b1 = 1000; b2 = 2000; m1 = 10; m2 = 7; M = 2000; g = 9.82; mu =

0熱

1答えて

数学的モデリング - Matlab ode45-forループ

「歴史」を持つ数学モデルのフェーズポートレートを取得するコードを作成する必要があります。私はコードの後に説明します。 close all; clear all; times = 1990:1:2015; hold on b=zeros(1,26); %75-2000 per 5 years b(1:5)=0.0358; b(6:10)=0.0339; b(11:15)=0.0

0熱

1答えて

このMATLABコードをデバッグして、グライダー降下をモデル化するにはどうすればよいですか？

私はこのMatlabプロジェクトを持っていますが、何らかの理由で私はそれを動作させることができなかったので、私はそれについて考えるのをやめることができません。目的：これは万フィートから削除されているグライダーの圧力、温度と密度の変化を計算うMATLABスクリプトです。それが下がるにつれて、計算された新しい値を使用し、ode45とP TとRhoの新しい値だけでなく、すべての点で微分する必要のある

0熱

1答えて

MATLAB ODE45：最後の解を保存するだけです。

ODE45には最後の解を残しておきたいと思います。むしろtspan = [t0 tf]からすべての解を返す。私は返されたベクトルをtfの解にしたいだけです。ドキュメントから Error using horzcat Requested 442368x1828 (6.0GB) array exceeds maximum array size preference. Creation of arrays

0熱

1答えて

MATLAB関数ode45関数 'あまりにも多くの入力エラー'

My機能は次のとおりです。 function [tout, yout] = Lorenz (rho, x0) %define constants sigma = 10; beta = 8/3; %initial conditions y0 = [x0; 1; 0]; f = @(t, y) [sigma*(y(2)-y(1)); (y(1)*(rho-y(3)))-y(2); (y

0熱

1答えて

Matlab ode45ヘルプ未定義変数

function yp = nonlinear (t,y) e=0.2; yp(1)=y(2); yp(2) = (-y(1)-e*y(1)^3); tspan = [0.20]; y0=[0;0] [t,y]=ode45('nonlinear',tpsan,y0) plot (t,y(:,1)) grid xlabel('time') ylabel('u') title (