Math 简单的MATLAB/Octave仿真
对于在这方面有经验的人来说,这应该是一个非常简单的问题,但我对这方面还是新手 我有以下系统(或): 考虑到以下输入:Math 简单的MATLAB/Octave仿真,math,matlab,simulation,octave,Math,Matlab,Simulation,Octave,对于在这方面有经验的人来说,这应该是一个非常简单的问题,但我对这方面还是新手 我有以下系统(或): 考虑到以下输入: u = min(2 - t/7.5, 2*(mod(t, 2) < 1)); 模拟后,xout将有五列。我的问题是:如何知道哪个是y系统的输出? 编辑:好的,简单地说,我想这样绘制解决方案: [tout, xout] = ode23(@func, [0 15], [1.5; 3; -0.5; 0; -1]) a = 1 % what goes here? 1, 2,
u = min(2 - t/7.5, 2*(mod(t, 2) < 1));
模拟后,xout
将有五列。我的问题是:如何知道哪个是y系统的输出?
编辑:好的,简单地说,我想这样绘制解决方案:
[tout, xout] = ode23(@func, [0 15], [1.5; 3; -0.5; 0; -1])
a = 1 % what goes here? 1, 2, 3, 4 or 5?
plot(tout, xout(:,a))
[T,Y,TE,YE,IE]=ode23(odefun,tspan,y0,选项)
下表列出了解算器的输出参数
- 时间点的T列向量李>
- Y解决方案阵列。Y中的每一行对应于T的对应行中返回的时间的解决方案
- TE事件发生的时间李>
- 在事件发生时,您可以选择解决方案李>
- 即消失的事件函数的指数i李>
伊斯滕·菲泽斯 与y相对应的一个,当然是x(1)
如果将代码与方程式进行比较,可以看到x(1)出现在代码中,y出现在方程式中的每个位置。这是我最好的猜测。很清楚;但是对于每一行,我将有五个值(来自五列的值)。哪一个将对应于y系统的输出?原始方程式的项为4y'^2,但您的func的项为4*x(3)^3。我想你希望是4*x(3)^2。是的,那是打字错误。我刚把它修好。
a = 1 % what goes here? 1, 2, 3, 4 or 5?
plot(tout, xout(:,a))