0

我对一组方程有疑问。我有:

x' = f(t, x, u) - it's a set of equations - dimension n
x1' = ..
x2' = ..
x3' = ..

我有 au - 它是一个向量 (u1, u2, u3 ..)

我如何在我的方程组中替换这个 u?例如 :

x1' = sin(t) * u1 + sin(u2)
x2' = u2*x2

u1 = sin(1000t)
u2 = cos(1000t)

我需要

x1' = sin(t) * sin(1000t) + sin(cos(1000t))
           x2' = cos(1000t) * x2

谢谢。

4

1 回答 1

1

好吧,假设您使用的是符号工具箱:

syms t u1 u2 x2;
x1prime = sin(t) * u1 + sin(u2);
x2prime = u2 * x2;

然后你可以使用方法A或B。方法A:

x1prime = subs(x1prime, [u1 u2], [sin(1000*t) cos(1000*t)])
x2prime = subs(x2prime, u2, cos(1000*t))

方法B:

u1 = sin(1000*t);
u2 = cos(1000*t);
x1prime = subs(x1prime)
x2prime = subs(x2prime)

我在 Matlab R2011A 上测试了这两种方法。使用最适合您的。

于 2013-10-06T21:42:27.133 回答