- 积分
- 33
- 注册时间
- 2005-6-2
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2009-11-22 09:54:28
|
显示全部楼层
来自 美国
本帖最后由 smarten 于 2009-11-21 20:04 编辑
- Manipulate[
- solab = FindRoot[{x[0.015] == -0.0033764, x'[0.015] == -0.430055} /.
- sol2, {{a, a0}, {b, b0}}];
- Plot[x[t] /. sol2 /. solab // Evaluate, {t, 0, 10}, PlotRange -> All,
- PlotLabel -> {{a, b} /. solab}], {{a0, 1}, 0, 5}, {{b0, 1}, 0, 5},
- Initialization :> {sol1 =
- DSolve[x''[t] + 25/2 x'[t] + 256 x[t] == a Sin[b t], x, t];
- solc = Solve[{x[Pi/b + 15/1000] == 0, x'[Pi/b + 15/1000] == 0} /.
- sol1[[1]], {C[1], C[2]}]; sol2 = sol1[[1]] /. solc[[1]]}]
复制代码 小图。
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
评分
-
1
查看全部评分
-
|