如何通过一个函数句柄作为参数数值吗?

1视图(30天)
考虑下面的代码:
R = 0.5;
l = 1;
g = 9.81;
ω= 0.25;
F = @ (t、y) [y (2), (R *ω^ 2 / l) * cos (y(1) -ω* t) - (g / l) *罪(y (1)));
t0 = 0;
T = 10;
y0 =(π/ 4,π/ 4);
h = 0.05;
[ts,欧美]= VectorRK4 (F, t0, T, y0, h);
[T、Y] =数值(F[0 10][π/ 4π/ 4]);
图;
情节(T Y (: 1),“- - -”、T、Y (:, 2),“-”。);
我想自己RK4模拟算法。错误出现在使用数值命令我获得下列的错误:
错误使用odearguments
@ (T、Y) [Y (2), (R *ω^ 2 / L) * COS (Y(1)ω* T) - (G / L) *罪(Y(1))]必须返回一个列向量。
这是如何固定的?

接受的答案

大卫。马塞罗说道
大卫。马塞罗说道 2022年10月18日
编辑:大卫。马塞罗说道 2022年10月18日
只使用一个分号,数值需要一列输出数组。
R = 0.5;
l = 1;
g = 9.81;
ω= 0.25;
F = @ (t、y) [y (2);(R *ω^ 2 / l) * cos (y(1) -ω* t) - (g / l) *罪(y (1)));
t0 = 0;
T = 10;
y0 =(π/ 4,π/ 4);
h = 0.05;
% (ts,欧美)= VectorRK4 (F, t0, T, y0, h);
[T、Y] =数值(F[0 10][π/ 4π/ 4]);
图;
情节(T Y (: 1),“- - -”、T、Y (:, 2),“-”。)

更多的答案(0)

类别

找到更多的在编程帮助中心文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!