声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2157|回复: 1

[Maple] 单摆画相图程序,请大家帮忙调试一下

[复制链接]
发表于 2008-10-17 11:08 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
> restart;

> pendphaseplot := proc (c, k, ui, uf, n1, vi, vf, n2, ti, tf, n, vu, vv) local s, i, a, j, b, sys, fcns, q, p; end proc;

> s := {};

>
> for i from 0 to n1 do;
> a := ui+i*(uf-ui)/n1;
> for j from 0 to n2 do;
> b := vi+j*(vf-vi)/n2;
> sys := diff(u(t), t) = v(t), diff(v(t), t) = -k*sin(u(t))-c*v(t);
> fcns := {u(t), v(t)};
> q := dsolve({sys, u(0) = a, v(0) = b}, fcns, type = numeric, method = rkf45);
> p := plots[odeplot](q, [u(t), v(t)], ti, tf, numpoints = n, view = [vu, vv], scaling = constrained);
> od:

> od:

> plots[display]([op(s)]);
> end:
> pendphaseplot(0, 2, -4, 4, 10, -3, 3, 6, -5, 5, 100, -10 .. 10, -5 .. 5);
以上是画单摆的相图的程序,总是出错,希望得到高人指教QQ:296620565 请注明:maple


[ 本帖最后由 无水1324 于 2008-10-25 11:59 编辑 ]
回复
分享到:

使用道具 举报

发表于 2008-10-21 14:37 | 显示全部楼层
附件是从http://www.maplesoft.com/applications/找到的,也许有用。

论坛对附件大小和类型有所限制,所以我用Maple另存为word文件,图片和输出结果删除了。用Maple工作表读一下。

pendphaseplot.doc (172.5 KB, 下载次数: 21)

[ 本帖最后由 MathPen001 于 2008-10-21 14:48 编辑 ]

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-24 19:38 , Processed in 0.083433 second(s), 25 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表