如何从Matlab中的ode45中获得t的特定值的y和y‘的值？
要从Matlab中的ode45函数中获取特定时间点t的y和y'的值，可以按照以下步骤进行操作：
- 定义微分方程： 首先，需要定义你的微分方程。可以使用一个函数来表示该微分方程，函数的输入参数为时间t和状态变量y，输出为状态变量y'的值。例如，假设你的微分方程是dy/dt = f(t, y)，可以创建一个名为"myODE"的函数文件，其中的代码为：function dydt = myODE(t, y) dydt = f(t, y);
- 调用ode45函数求解微分方程： 在Matlab中，使用ode45函数可以求解常微分方程。调用ode45函数时，需要指定微分方程函数的名称，以及求解时间范围和初始条件。例如，假设求解的时间范围是从t0到tf，初始条件为y0，可以使用以下代码调用ode45函数：[t, y] = ode45(@myODE, [t0 tf], y0);
- 获取特定时间点的y和y'的值：
在得到ode45的解后，可以通过插值方法获取特定时间点t的y和y'的值。Matlab中的interp1函数可以用于插值。具体步骤如下：
- 定义一个数组t_values，包含你想要获取y和y'值的特定时间点。
- 使用interp1函数，将数组t_values作为输入，从t数组中获取对应时间点t的y和y'的值。例如，可以使用以下代码获取t_values时间点的y和y'的值：y_values = interp1(t, y, t_values); y_prime_values = interp1(t, y_prime, t_values);
- 注意：在上述代码中，y_prime代表微分方程中的y'。
这样，你就可以从Matlab中的ode45函数中获得特定时间点t的y和y'的值了。
请注意，以上答案是以一般性的角度给出的，针对具体问题的答案可能会有所不同。另外，由于问题中要求不涉及特定云计算品牌商的信息，因此没有提供与腾讯云相关的产品链接。
相关·内容
我正在尝试使用ode45来解决一个ODE系统：
[X,Y]= ode45(@sys,[0, T],y0);
哪里,
function dy = sys(t,y)
dy(1) = f_1(y)
dy(2) = f_2(y)
dy(3) = f_3(y)
end
问题是，函数ode45要求y0是初始值[y_1(0), y_2(0), y_3(0)]，而在我的系统中，我只有值[y_2(0), y_3(0), y_3(T)]可用。
从数学上讲，这组初始/终止条件应该足以约束系统，但是有什么方法可以通过ode45或MATLAB中的任何其他函数来使用它吗
浏览 3提问于2013-06-14得票数 7
回答已采纳
我正在用Matlab来解一个微分方程。我想强制ode45采取恒定的步长，所以在求解方程时，它总是在T轴上递增0.01。我该怎么做呢？
ode45一直在采取优化的、随机的步长，而我似乎想不出如何让它采用一致的、小的0.01步长。代码如下：
options= odeset('Reltol',0.001,'Stats','on');
%figure(1);
%clf;
init=[xo yo zo]';
tspan=[to tf];
%tspan = t0:0.01:tf;
[T,Y]=ode45(name,tspan,init,optio
浏览 0提问于2010-03-31得票数 6
回答已采纳
1回答
我试图弄清楚如何将外部变量(随时间-荷尔蒙水平变化)输入到ode中。我在考虑使用事件。我尝试做优化(在Matlab中使用内置事件优化微分方程)
为了学习如何使用它们，我尝试在中运行代码，然而，我做错了一些事情。我猜问题出在我的函数myfun上，但我不知道为什么。
我试着运行以下命令：
function ydot=myfun(t,y,a)
S1 = 2 *a* 2 * y(2) - 2 * 3 * y(1)*y(1);
S2 = 3 * y(1)*y(1) + a - 5 * y(2);
ydot =
浏览 2提问于2017-02-27得票数 0
1回答
我在解一首歌
y''(t) +a y(t) +b=0
用Matlab的ode45解算器。它迭代直到找到一个解决方案，这符合错误准则。为了进行比较，我希望在每个步骤中阅读y的前一个值。
做这件事最好的方法是什么？
我的主要文件是
[t,y] = ode45(@odefnc,tspan,[0 0]',options);
而odefnc是
function dx=odefnc(t,x)
...
if history(end)<13
dx=[x(2),-a*x(1)-b]'
else
dx=[x(2),-c*x(1)-d]'
end
浏览 0提问于2013-08-06得票数 1
对我来说，不清楚如何在matlab中使用ode45。我有一个向量u，有人知道这个错误是什么吗？
function [v, t, d] x =deep(l)
x = max(u(:,1));
v = u;
t = 1;
d = 1;
options = odeset('Events', @deep)
但是当我调用ode45时
ode45(@odefun, [0 85], [0 0], options)
我得到了这个错误：
eventValue = feval(eventFcn,t0,y0,eventArgs{:});
Error in ode45 (line 148)
o
浏览 3提问于2020-11-10得票数 0
我需要解一个二阶微分方程，我特别需要，在最后一个时间点，y的一阶导数值。我的代码如下： [T Y]=ode45(@(t y)vdp4(t,y,0.3),[0 1],[0.3/4,((3*0.3)^0.5)/2]); 我知道输出将包含ode45计算函数时使用的值。为了得到特定时间的y值，建议在MATLAB documentation中给出两个以上的时间点。我也这么做了。 tspan=[0:0.01:1]
[T Y]=ode45(@(t y)vdp4(t,y,0.3),tspan,[0.3/4,((3*0.3)^0.5)/2]); T向量仍然没有从0到1的所有值(最后一个值是0.39)。这特别是
浏览 116提问于2020-01-12得票数 0
我有一个带有句柄@func的函数，它的初始条件为y0，我想在t = [0, tspan]上进行测试。确切地说，我需要做什么来增加ode45使用的时间步骤的数量，而不改变tspan
我看到了ode45的MATLAB文档，并看到我需要以某种方式将options输入更改为ode45。但是，我并不真正理解如何做到这一点，因为没有提供示例代码。
浏览 11提问于2015-11-04得票数 2
回答已采纳
我想从python调用MATLAB的ode45函数，就像下面的python代码中调用MATLAB的isprime()函数一样
import matlab.engine
eng = matlab.engine.start_matlab()
tf = eng.isprime(37)
print(tf)
但我并不知道如何在python中实现代码。我试过了
import matlab.engine
def dydt(t,y):
dydt= 2*y
return dydt
eng = matlab.engine.start
浏览 5提问于2019-10-17得票数 1
回答已采纳
3回答
我想知道如何将MATLAB函数ode45导出为python。根据文档，应该如下所示：
MATLAB: [t,y]=ode45(@vdp1,[0 20],[2 0]);
Python: import numpy as np
def vdp1(t,y):
dydt= np.array([y[1], (1-y[0]**2)*y[1]-y[0]])
return dydt
import scipy integrate
l=scipy.integrate.ode(vdp1([
浏览 1提问于2018-01-25得票数 6
回答已采纳
2回答
我有从流行病传播中导出的微分方程。我想得到数值解。这是方程式，
T是一个自变量，取值范围为[0,100]。初始值为
y1 = 0.99; y2 = 0.01; y3 = 0;
一开始，我打算在matlab中用ode45函数来处理这些问题，但是，我不知道如何表示级数和组合。所以我在这里请求帮助。
**
问题是如何将方程的右侧表示为odefun，它是ode45函数中的一个参数。
**
浏览 3提问于2011-05-25得票数 1
回答已采纳
1回答
我想用Simulink作用域绘制一些数据，这是我简化的Simulink模型：
在MATLAB函数块中有我的代码：
1. function pCO2_a = Oxy(Qg)
2.
3. coder.extrinsic('testWQ');
4.
5. pCO2_a = zeros(41,1);
6. pCO2_a = testWQ();
7. end
我尝试调用的是另一个ode函数：
1. function y = testWQ
2.
3. Qg = 3;
4. tspan = [0 5];
5. y0 = 0;
6. [t, Y] = ode45(@f,tspan,y0);
浏览 7提问于2020-02-20得票数 0
回答已采纳
如果我有一首颂歌，并以两种方式写成，就像下面这样：
function re=rabdab()
x=linspace(0,2000,2000)';
tic;
[T,Y] = ode45(@fun,[x],[0 1 1]);
[T,Y] = ode45(@fun,[x],[0 1 1]);
[T,Y] = ode45(@fun,[x],[0 1 1]);
toc;
tic;
[A,B] = ode45(@fun2,[x],[0 1 1]);
[A,B] = ode45(@fun2,[x],[0 1 1]);
[A,B] = ode45(@fun2,[x],[0 1 1]);
toc;
浏览 1提问于2013-05-27得票数 3
回答已采纳
1回答
当我尝试使用这个Matlab代码时，它会变成一个无限循环。我试图在ode45内部执行集成
clear
clc
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'+',T,Y(:,2),'*',T,Y(:,3),'.')
function dy = rigid(t,y)
dy = zeros(3,1); % a col
浏览 0提问于2014-02-01得票数 0
回答已采纳
1回答
我需要使用ode45在t=1、1.5和3处找到解决方案的近似值，然后在0.5、4上绘制解决方案 %% 7) exp(y)+(t*(exp(y)-sin(y))*y'=0, y(2)=1.5
% 7a) Find approximate values of the solution at t=1, 1.5, and 3 then plot
% the solution on [0.5,4].
[t,y]=ode45(@(t,y) -exp(y)./(t.*(exp(y))-sin(y)),0.5:.2:4,1.5) 如上所示，初始条件从t =2开始，而不是0。如何在从t=2开
浏览 61提问于2019-03-14得票数 0
回答已采纳
我使用MATLAB R2016b -学生版开发了一个m文件，该文件以一个符号微分方程f(t，y)为输入，并根据初始条件输出一个斜率场和解曲线。代码是
prompt={'dy/dt =','tspan','y0 ='};
title='Slope Field Calculator';
answer=inputdlg(prompt,title);
tspan = str2num(answer{2}); %#ok<*ST2NM>
y0 = str2double(answer{3});
syms t y
f = symfun
浏览 3提问于2017-06-23得票数 1
回答已采纳
我正在尝试用matlab的ode45函数来解决耦合的常微分方程：
下面是我的名为'Rot‘的函数，用来描述这些使用matlab ode45的ODE。
function omega= Rot(t,y)
omega(2,1)=(0.03*sin(3*t)*((cos(Y(1)))^2)+0.002*t^3*sin(y(1)))...
/-((cos(Y(1)))^2)+((sin(Y(1)))^2);
omega(1,1)=((0.002*t^2-omega(2,1)*sin(y(1)))...
/-cos(y(3))*sin(y(2)))*cos(y(2))+0.01*t^2+
浏览 3提问于2018-11-15得票数 1
1回答
我首先定义了dy/dt=y和dy/dt=t的函数：
function dy=d(y):
dy=y
end
function ddy=dd(t):
ddy=t
end
然后我分别使用了ode45：
[t,y]=ode45('d',[1 10],1)
[t,y]=ode45('dd',[1 10],1)
它返回以下错误:使用d太多输入参数的错误。
我的问题是：
我哪里出错了？
Matlab是如何知道y还是t是自变量的？当我定义第一个函数时，它可以被合理地解释为dt/dy=y，而不是dy/dt=y。
浏览 3提问于2016-01-03得票数 1
回答已采纳
1回答
在了解如何在八度/matlab中实现事件时，在微分方程的解算方面，我遇到了一些困难。
例如，考虑求解微分方程y‘=-y的简单代码：
function dy = odefun(t,y)
dy = -y;
endfunction
options = odeset('RelTol',1e-4,'AbsTol',[1e-4]);
[T,Y] = ode45(@odefun,[0 12],[1],options);
现在我想介绍一个随机事件。例如，在某个固定的时间步长上，我想随机改变y的值，然后根据微分方程继续演化。我该怎么做？
浏览 2提问于2014-02-01得票数 2
回答已采纳
1回答
我理解数值求解的逻辑，即matlab从未知值开始循环，直到方程收敛到一个值为止。然而，我不明白的是如何正确地进入我的方程式。我认为使用ode45函数是最好的方法。我有下面的等式U^(n+1) = U^n - (t'*3250/10)-(t'/2)(.004(v^n)^2/10。我怀疑我的方程需要不同的形式，但是我不确定让matlab求解方程的正确方法。
我试着在matlab中按原样输入方程，但是它抱怨v和n是未知变量，我不确定如何处理它们。求解这个方程的最终目标是找到v的值。
C=@(t,v) u^n-(3250*t'/10)-(t'/2)*((.004*(v^
浏览 1提问于2014-09-16得票数 0
2回答
我是Matlab的新手，很抱歉问了一个平庸的问题，可能已经回答了问题，但我没能解决这个问题。
我有这首歌要解
g''(x)+f(x)g'(x)
g(0)=0，g'(0)=0.5
f(x)是已知的向量，所以我尝试
xspan = linspace(0,10,100);
y0 = [0,0.5];
%f=f(xspan), known function, for example f=1./xspan;
[x,y] = ode45('odefun',xspan,y0);
对于odefun
function dy=odefun(x,y)
global
浏览 3提问于2014-05-14得票数 0
回答已采纳
我是Matlab的新手。我希望你能帮助我。我必须解决一个使用ODE45函数的ODE系统。这是描述我的等式的函数。
function dNdt = rateEquations(t, y)
%populations of corresponding state
Ng = y(1);
Ns = y(2);
Nt = y(3);
%All constants used are dropped for the sake of easy reading.
请注意参数F。
%rate equations
dNs = s0 * Ng * F - Ns/ t_S1;
dNt
浏览 2提问于2011-10-07得票数 6
回答已采纳
作为Matlab语言的初学者，我正在使用包含要求解的常微分方程的test.m来探索ode45函数。 在编辑器中，我通过两个选项分别调用这个函数。在一种情况下，我定义了选项1 t=[0 50];
y0=[0 2];
[t,y]=ode45(@(t,y)test(t,y),t,y0); 在选项2中，我要求提供一个结构输出 t=[0 50];
y0=[0 2];
sol=ode45(@(t,y)test(t,y),t,y0); 然而，与选项1相比，结构选项2的结果时间步长更少，因此我的图更“粗糙”。 我找不到一种方法来增加步骤数，以完善option 2....Any想法的解决方案？
浏览 17提问于2021-03-16得票数 0
回答已采纳
1回答
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 ('u vs. t')
hold on;
对不起，我绝对是matlab的新手，当我试图执行代码时，它写着“变量t的未定义函数”。我试图用ode45来解一个微分方
浏览 1提问于2017-07-11得票数 0
回答已采纳
1回答
我正在Matlab中使用ode45进行实验。我已经学习了如何将参数传递给ode函数，但我仍然有一个问题。假设我想要计算一辆车的轨迹(速度剖面)，我有一个函数，例如getAcceleration，它给出了汽车的加速度，也给出了正确的齿轮：[acceleration, gear] = getAcceleration(speed,modelStructure)，其中modelStructure代表了汽车的模型。
ode的功能是：
function [dy] = car(t,y,modelStructure)
dy = zeros(2,1);
dy(1) = y(2
浏览 2提问于2012-11-26得票数 3
回答已采纳
y'''' + (a1 + a2*v(t)^2)*y'' + (a3 + a4*v(t)^2)*y = 0
y(0) = 2, Dy(0) = 0, D2y(0) = 0, D3y(0) = 0
我试着用Matlab的dsolve来求解上面的四阶常微分方程，但是我得到的结果非常大。正如我现在所理解的，dsolve在这里不能工作，我将不得不使用ode45来解决这个方程。
我试着在上阅读Matlab的文档，但我不能完全理解它。此外，在我的情况下，系数也是时间相关的。
如何将此公式转换为ode45可理解的格式
浏览 3提问于2013-03-07得票数 2
回答已采纳
我正在尝试解决我的微分方程式课的这个作业，我遇到了一些问题。
给定这个ODEs系统：
dy(1) = y(2) − y(3);
dy(2) = −3*y(2) − 5*sin(ω*y(1));
dy(3) = y(1)*y(2);
和omega=3在一起。初始值：y(1)=0; y(2)=4; y(3)=1;
在t=0和t=20之间绘制f图
f(t) = |2*cos(y1(t)) + y2(t)|
我首先尝试用下面的代码用ODE45模拟方程系统：
关于编辑器：
function dy=modelo10(t,y)
global omega
dy = zeros(3,1);
浏览 14提问于2019-05-07得票数 1
2回答
如果有人能帮忙解决以下问题，我将不胜感激。我有以下颂歌：
dr/dt = 4*exp(0.8*t) - 0.5*r ,r(0)=2, t[0,1] (1)
我用两种不同的方法解决了问题。采用Runge方法(4阶)和Matlab中的ode45。我将这两种结果与解析解进行了比较，分析结果如下：
r(t) = 4/1.3 (exp(0.8*t) - exp(-0.5*t)) + 2*exp(-0.5*t)
当我绘制每种方法相对于精确解的绝对误差时，我得到以下结果：
对于RK-方法，我的代码是：
h=1/50;
浏览 6提问于2014-02-18得票数 9
回答已采纳
我有一个二阶微分方程 function Xdot =second_order(t,X)
y=X(1);
y_d=X(2);
Xdot=zeros(2,1);
Xdot(1)=y_d;
Xdot(2)= -5*y_d-7;
y_dd=Xdot(2);
end 使用MATLAB的ode45命令，我能够模拟它，我能够以这种方式捕获向量中的两个状态变量y和y_d，作为时间的函数。 我想在矢量中捕获y_dd作为时间的函数，但我无法保存它。如何将其视为输出变量？我尝试了对我的函数文件function [Xdot,y_dd] =second_order(t,X)进行以下更改，但是我不知道如何使用o
浏览 82提问于2020-01-02得票数 1
回答已采纳
因此，我尝试使用Matlab中的ode45求解器来模拟简单的抛射体运动(无空气阻力等)。这是我到目前为止的代码： function [x,y] = trajectory_without_AR_45(v0,theta, dt)
%Path of mortar without air resistance using ode45
g = 9.81;
t_start = 0;
t_end = 100;
%Initial Conditions
y01 = 0; %initial x
y02 = v0 * cos(theta); %finding initial velocity in x dir
浏览 37提问于2018-12-19得票数 0
在Matlab中可以使用匿名函数作为事件函数吗？我想要做的基本上是
opt = odeset('Events', @(t,y) (deal(y(end)-t^2,1,0)));
[T,Y] = ode45(@odefun,[tstart tend],y0,opt);
但是，这将返回一个错误，报告要处理的输出数量必须完全匹配。有没有其他方法可以让匿名函数返回多个参数？
浏览 4提问于2013-02-19得票数 5
回答已采纳
所以我需要解x''(t) = -x(t)^p，初始条件是x(0)= 0，v(0) = x'(0) = v_o = 1，参数p的值是1。
这就是我所拥有的：
function [t, velocity, x] = ode_oscilation(p)
y=[0;0;0];
% transform system to the canonical form
function y = oscilation_equation(x,p)
y=zeros(2,1);
y(1)=y(2);
y(2)=-(x)^p;
浏览 5提问于2012-11-05得票数 0
回答已采纳
1回答
我试图用Matlab中的函数ode45来求解三个微分方程组。我并不真正理解我所犯的错误，我需要一些帮助来理解我做错了什么。
微分方程如下：
F1 = -k1y1+k2(y2-y1)
F2 = -k2(y2-y1)+k3(y3-y2)
F3 = -k3(y3-y2)
我在Matlab中的代码是：
function dz = kopplad(t, z)
global m1 m2 m3 k1 k2 k3
dz = [z(2)
-k1*z(1)/m1 + k2*(z(2)-z(1))/m1
z(4)
-k2*(z(2)-z(1))+k3(z(3)-z(2))/m2
z(6)
-k3(z(3)-z(2
浏览 1提问于2013-05-08得票数 1
我正在尝试在Matlab中模拟Kuramoto的隐含。我试着用ode45来解决这个问题。我还看到其他人使用Runge-kutta方法。我知道ode45使用龙格-库塔方法，然而，我从每个方法中获得的值有可疑的不同。
kuramoto= @(x,K,N,Omega)Omega+(K/N)*sum(sin(x*ones(1,N)-(ones(N,1)*x')))'
%Kuramoto is a model of N coupled ocilators (such as multiple radiowaves)
%The solution to the model is the pha
浏览 1提问于2018-08-29得票数 1
2回答
matlab中的ODE45函数接受参数：
(function,[tinitial tfinal],yinitial)
但在这里，我相信，时间跨度是预先确定的。我如何为它赋值一个向量？我的意思是如何解决域1:0.1:5的ODE
谢谢。
浏览 1提问于2013-04-06得票数 0
回答已采纳
1回答
我想用数值的方法来计算带有奇点的微分。
例如，y'=y^2+1，y(0)=0。解析解不难找到，y=tan(x)。
当我应用数值方法时，会出现问题，请参阅Matlab中的代码。
tspan = [0 6];
y0 = [0; 0];
ode = @(t, y) y.^2+1;
[t, y] = ode45(ode, tspan, y0);
plot(t, y(:,1))
axis([0 6 -20 20])
它给出了下面的图，即在第一个奇异点之后，它忽略了解的部分。
我的问题:如何找到具有奇点的微分方程的全数值解？
提前感谢！
浏览 3提问于2021-12-16得票数 0
回答已采纳
因此，在MATLAB中，当使用ODE45处理多个ODEs时，您可以像我的教授所说的那样，将术语“打包”到一个术语中，如下所示。
function [dYdt] = HW3_P2_ftn(t,Y)
%Unpack the state variables
D = Y(1);
CA = Y(2);
所以现在Y是由两个不同变量的向量组成的。
在我的代码末尾，在ODE45之后，我只想绘制Y(2)而不是Y(1)。
有人能给我解释一下我该怎么做吗？
浏览 2提问于2020-02-05得票数 0
1回答
我要在Matlab中求解这个二阶ODE：
(a + f(t))·(dx/dt)·(d²x/dt²) + g(t) + ((h(t) + i(t)·(d²x/dt² > b·(c-x)))·(dx/dt) + j(t))·(dx/dt)² + k(t)·(t > d) = 0
哪里
a，b，c，d是已知的常量
f(t)、g(t)、h(t)、i(t)、j(t)、k(t)是依赖于t的已知函数。
x是这个职位
dx/dt是速度
d²x/dt²是加速度
注意到这两个条件
i(t)在(d²x/dt² > b·(c-x))方程中的引入
浏览 2提问于2017-03-04得票数 1
回答已采纳
标准的Matlab ode45调用是：
[T,Y] = solver(odefun,tspan,y0,options);
但是我也想要每个状态向量Y的odefun的值。似乎应该有一个选项来返回像这样的内容
[T,Y, DYDT] = solver(odefun,tspan,y0,options);
因为这将是最有效的方法，因为例程已经在内部计算了每个点的导数。显而易见的工作是这样做
DYDT = odefun(T,Y);
但是，我需要编写一个修改过的odefun版本，因为传递给solver的应该是一个表示单个状态的列向量。这里很明显，我想要计算solver返回的所有状态下的DYDT。有什么建议
浏览 0提问于2014-10-10得票数 0
定义像下面这样的微分方程在MATLAB中没有问题，我可以用ODE45函数来解决它。
# example.m
x = pi / 2;
x_span=[0 pi/2];
ic=[0 1];
[X OUT] = ode45(@fun,x_span,ic)
#fun.m
function out=fun(x,s)
y1 = s(1)
y2 = s(2)
out=[y2;x*y1];
但是，如果我添加基于x的新函数，则ODE45无法工作。
# example.m
x = pi / 2;
A1 = sin(x);
A2 = sin(x)+1;
x_span=[0 pi/2];
ic=[0 1];
[X
浏览 2提问于2015-07-10得票数 1
回答已采纳
1回答
我试图积分一组常微分方程，并在给定的点上计算解及其一阶导数。Matlab提供了求解方程的"ode45“和计算解及其一阶导数的"deval”，但是在我的情况下，我无法使用它们。
Matlab给出了以下示例：
sol =ode45(@sol 1，0 20,2 0)；
X=直线空间(0，20,100)；
Y= deval(sol，1)；
地块(x，y)；
我猜"vdp1“是由Matlab定义的，这个例子运行得很好。但是，当我尝试使用以下代码时：
sol =ode45(‘vdp1 1’，0.20，2.0)；
X=直线空间(0，20,100)；
Y= deval(sol
浏览 4提问于2014-10-12得票数 1
回答已采纳
1回答
我正在尝试在Matlab中模拟方程y“+ 4y = 0，初始条件y(0) =1和y'(0) =0。下面的代码片段显示了上述常微分方程的位置作为时间函数的工作近似值：
clear
syms y(t)
%Differential Equation Conversion
P = odeToVectorField(diff(y, 2) == -4*y);
M = matlabFunction(P,'vars',{'t','Y'});
%Position Approximation
pos_solution = ode45(M,[0
浏览 3提问于2016-05-24得票数 0
我正试着跟踪这个。我试图复制第20页图7中的解决方案的图表。
首先，我想重新创建左边的图片。问题中的系统就是我作为dX所拥有的。下面是我在一个m文件中的内容：
function dX = CompetitionModel(t,X)
bs = 8*10^(-3);
bl = 4*10^(-3);
bh = 6.4*10^(-3);
N = bs + bl + bh;
K = 10^8;
m1 = 2*10^(-5);
m2 = 9*10^(-9);
p = 5*10^(-13);
I = 10^(-3);
T = 0;
浏览 1提问于2014-07-18得票数 1
回答已采纳
1回答
我正在尝试在Matlab中解决一个ODE，我有以下问题:我的代码如下：
xinit=[0.19;25;0;7];
t=0:768:76800; %% 101 cells
[t1,y]=ode45(@Model_Bio,t,xinit);
在函数@Model_Bio中，我有一个参数，我需要读取每个时间步的相应值(101个值)！所以，我的@Model_Bio在某种程度上是这样的：
load 'mydata'
a=mydata;
.....
问题是，当我执行ode45(@Model_Bio，t，xinit)时，它使用自动时间步长调用我的函数(例如50次！)我的问题是我不能为我
浏览 1提问于2015-07-03得票数 0
1回答
我有一个具有离散值的向量，需要传递到我的ODE系统中，我想使用ode45命令。当我使用这个向量时，它需要在求解器内进行插值。有什么办法可以做到吗？
我有一个线性的耦合系统。
dxdt = f(t)*x + g(t)*y
dydt = g(t)*x + h(t)*y
我有三个向量f( t )，g(t)和h(t)作为t的函数，我需要把它们传递给求解者。
我可以用C或C++编写解算器的代码。有人建议，在Matlab中这样做会更快。有人能建议一个办法吗？
浏览 10提问于2014-10-29得票数 0
回答已采纳
1回答
我正在使用ode45。当我声明我的函数:f= @(x，y) (cos(x-y) - x*y -y)/(2+2x)时，我得到一个错误：“意外的MATLAB表达式”
为什么会发生这种情况？
谢谢。
浏览 0提问于2014-09-03得票数 0
2回答
我在C中实现了一个同等版本的Matlab脚本。为了运行ode45，我选择了GNU科学库。但是ode45为每个版本产生不同的输出。我已经工作了一段时间了，我找不到问题。我使用gsl_odeiv2_driver_apply_fixed_step进行与Matlab中相同的步骤。
Matlab脚本
function ExpoGrowthEqn
% code to solve the exponential growth equation
% dN/dt = r*N, N(0)=N0
% parameter values
r=0.2;
N0=10;
% numerical parameters
s
浏览 5提问于2013-04-04得票数 0
我要在Matlab中求解这个二阶ODE：
(a + f(t))·(dx/dt)·(d²x/dt²) + g(t) + ((h(t) + i(t)·(d²x/dt² > b·(c-x)))·(dx/dt) + j(t))·(dx/dt)² + k(t)·(t > d) = 0
哪里
a，b，c，d是已知的常量
f(t)、g(t)、h(t)、i(t)、j(t)、k(t)是依赖于t的已知函数。
x是这个职位
dx/dt是速度
d²x/dt²是加速度
注意到这两个条件
i(t)在(d²x/dt² > b·(c-x))方程中的引入
浏览 0提问于2017-09-24得票数 1
回答已采纳
1回答
如何计算曲线到达特定x坐标所需的时间(在Matlab中)。假设我们有: dx/dt = x^2 + y^2和dy/dt = 5.x.y，曲线从点(a，b)开始。在ode45的帮助下，我能够得到曲线的图形。我需要计算曲线到达x= c，(c>a)所需的时间。有人告诉我，这可以通过插值来完成，但我不知道如何编写代码。
浏览 1提问于2013-12-07得票数 0
我试图在ode45中使用MAtlab，并希望修复MAtlab使用的点数(时间步骤数)。在'refine'中使用ode45选项似乎没有帮助。例如，如果我将'refine'设置为10，Matlab将返回一个101数组。改变'RelTol'和'AbsTol'也没有帮助。我知道将tspan写成[0,t1,t2,t3,...,tn]是可能的，这解决了这个问题，但是我想通过'refine'选项来修正点数。
浏览 2提问于2014-01-03得票数 2
2回答
我编写了一个函数和一个脚本来计算和绘制一个1自由度系统的响应，使用matlab中的4阶runge方法。下面是我想要绘制的系统对不同输入值的响应。我有零初始条件，输入如下：
F(t) = 5*sin(w*t)
其中w=2 5 10
我不知道如何在函数中使用循环来形成w的所有值的状态空间，并求解我的方程，并将其绘制为w的不同值，从而为每次迭代生成不同的输入。注意，一种方法是为所有三种情况编写三个不同的函数并计算响应，我不想这样做。我只想要一个函数和我的脚本，它使用ode45函数来求解我的微分方程，但是画出w的所有值的响应。你建议我怎么做。以下是我的脚本和功能：
function xDot = My
浏览 4提问于2020-10-06得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
