题目:考察具有粘性阻尼的单摆绘制其在区域中的相图。

首先编写一段 MAPLE 子程序 pendphaseplot 如下:

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; s := {};

定义变量

for i from 0 to n1 do a := ui+i*(uf-ui)/n1;

a是u变量变化的步长

for j from 0 to n2 do b := vi+j*(vf-vi)/n2;

b是v变量变化的步长

sys := diff(u(t), t) = v(t), diff(v(t), t) = -k*sin(u(t))-c*v(t);

定义系统

fcns := {u(t), v(t)};

定义方程(带入dsolve中用的,就是个过程变量)

q := dsolve({u(0) = a, v(0) = b, sys}, fcns, type = numeric, method = rkf45);

带入不同的初值求得相轨迹

p := plots[odeplot](q, [u(t), v(t)], ti .. tf, numpoints = n, view = [vu, vv], scaling = CONSTRAINED);

将上一步求得的相轨迹作出相应的图像(通过数值解作图)

s := s union p

union是取并集的操作,其实就是用s保存p

end do

end do;

plots[display]([op(s)])

op的含义就是提取整个s用于作图,从这里感觉maple编程是一个比较笨重的语言,当然了,用op也可以提取固定的一些项,详细看帮助即可

end proc

至此整个程序结束

命令行:

无阻尼情况:

pendphaseplot(0, 2, -4, 4, 10, -3, 3, 6, -5, 5, 100, -10 .. 10, -5 .. 5)

有阻尼的情况:

pendphaseplot(.5, 2, -4, 4, 8, -3, 3, 6, -5, 5, 100, -10 .. 10, -6 .. 6)

求解过程到这里结束,原始的题目和代码出自胡海岩编写的《应用非线性动力学》,csdn不允许上传课本资料。

非线性ode及在工程中的应用第一次大作业相关推荐

  1. 非线性ode及在工程中的应用第四次作业更正版

    1. 这个题先写一个非常简洁的写法(见胡海岩<应用非线性动力学>P49): 不过既然是作业,肯定要自己推一下验证是不是可以完全理解,推导的过程我手写之后扫描存在这里: 当然了结果有一些小问 ...

  2. 南京工程学院C++/VC程序设计大作业

    南京工程学院C++/VC程序设计大作业 源码下载 链接:源码下载链接 提取码:1111 学生成绩输入及统计分析 模拟自动取款机 C++/VC程序设计大作业任务书 一.课外大作业要求: 1.每人从附件1 ...

  3. msflexgrid允许大选择_工程中要用多大的电线电缆?一文教你怎么算

    一.电机功率与配线直径计算 首先要计算100KW负荷的线电流. 对于三相平衡电路而言,三相电路功率的计算公式是:P=1.732IUcosφ. 由三相电路功率公式可推出: 线电流公式:I=P/1.732 ...

  4. 0.5mm的焊锡丝能吃多大电流_工程中要用多大的电线电缆?一文教你怎么算

    一.电机功率与配线直径计算 首先要计算100KW负荷的线电流. 对于三相平衡电路而言,三相电路功率的计算公式是:P=1.732IUcosφ. 由三相电路功率公式可推出: 线电流公式:I=P/1.732 ...

  5. 电动机在计算机控制系统中的角色,机电一体化大作业(二).doc

    机电一体化系统设计结课论文--执行元件系统中电动机的优化 <机电一体化系统设计> 结课论文 论文题目: 电动机的优化节能 学院:机械电气工程学院 专业:工业12级(1)班 姓名: 周-- ...

  6. 机器学习-特征工程中的特征降维

    对于一个机器学习问题,数据和特征决定了机器学习的上限,而模型和算法只是逼近这个上限.由此可见,数据和特征在模型的整个开发过程中是比较重要.特征工程,顾名思义,是对原始数据进行一系列工程处理,将其提炼为 ...

  7. 机器学习-特征工程中的数据预处理

    对于一个机器学习问题,数据和特征决定了机器学习的上限,而模型和算法只是逼近这个上限.由此可见,数据和特征在模型的整个开发过程中是比较重要.特征工程,顾名思义,是对原始数据进行一系列工程处理,将其提炼为 ...

  8. matlab计算海洋浮力频率_水下机器人在海洋石油工程中的应用

    世界海洋总面积约为3.6×108平方千米,占地球总面积的70.8%,在广阔的海洋中,蕴含着非常丰富的生物资源.矿物资源和海洋能源.其中最具经济开采价值的是海底的石油和天然气资源,海洋石油资源储量占世界 ...

  9. 计算机在中学物理的应用分析报告,BGPLUS官方合作 | 加州大学伯克利分校 | 物理:计算机算法及物理在工程中应用...

    原标题:BGPLUS官方合作 | 加州大学伯克利分校 | 物理:计算机算法及物理在工程中应用 课题名称 = 工程学综合课题:计算机算法及物理原理在工程项目优化中的应用= 项目背景 项目集结来自伯克利工 ...

最新文章

  1. Pycharm+Django搭建第一个Python Web程序
  2. java 模块化_Java 9 新特性 - 模块化 - Java 技术驿站-Java 技术驿站
  3. 有些事儿,工程师可能今生仅此一次
  4. Pyhton为什么这么厉害?无意中从一道奥数题发现Python真的爽翻了
  5. oracle的pdb,oracle pdb基本管理
  6. c++简单程序设计-5
  7. 粉刷木板(ybtoj-单调队列)
  8. ApplePay支付使用
  9. 系统集成资质 - 案例分析与常规题答案
  10. 用mysql建立商城数据字典_把mysql数据库生成数据字典,直接可用
  11. 如何寻找互联网红利期产品?
  12. 《深入理解计算机系统(修订版)》读感
  13. 一台空调的容量是多少_大型中央空调工程节能改造方案如何选择
  14. stlink固件版本低,但升级失败的解决办法
  15. Python的一点人生经验
  16. Python之爬取58同城在售楼盘房源信息
  17. 为什么用python写爬虫_老猿为什么写Python爬虫教程
  18. python 通过tif文件生成经度和纬度的二维矩阵
  19. 微信小程序实现音乐播放器(2)
  20. 1. 数据库的四个基本概念

热门文章

  1. java微信普通群发图文消息_java微信发送图文和文字消息
  2. 网络 || icmp协议 || 三层交换机
  3. 淘宝API接口( item_detail - 淘宝商品详情查询)
  4. 以太坊 – 使用truffle console访问智能合约
  5. 酒店会议服务与管理【1】
  6. 三、Win10系统怎么自动清理缓存,无需手动按时清理磁盘,自动释放空间哦
  7. 【Dfinity完整概述】一文带你读懂Dfinity的链上运行
  8. js设置div不可点击
  9. Google 真的「不作恶」吗
  10. 在windows下使用Cygwin模拟unix环境,并安装apt-cyg,repo等工具