| 查看: 1756 | 回复: 21 | |||||
| 本帖产生 1 个 数学EPI ，点击这里进行查看 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
求助——积分求面积遇到问题已有3人参与
|
|||||
|
如图，x2,y2为转子方程，我想求x2,y2曲线围成的面积，用matlab编写程序如下： z1=6; z2=z1-1; R=48.78; r=8.13; Q=R-r; e=7.05; syms theta k=12.6; f=r/e; a=sin(z1*theta); b=f+cos(z1*theta); m=atan(a./b); p=m-theta; %theta c=sin(theta+p); n=asin(f*c); o=n-p; %phi g=(z1/z2)*o; %psi x2t=Q*sin(theta)+e*sin(z2*theta)-k*sin(p); %原始齿形 y2t=Q*cos(theta)-e*cos(z2*theta)+k*cos(p); x3t=x2t.*cos(o-g)-y2t.*sin(o-g)-e*sin(g); %共轭齿形 y3t=x2t.*sin(o-g)+y2t.*cos(o-g)-e*cos(g); d2=diff(x3t); ds2=inline(y3t.*d2); s2=5*quadl(ds2,0,2*pi/5) 但是运行出来有问题，请高手指点，转子的面积该如何求，先谢谢啦  转子方程.png |
回复此楼» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 数学EPI: 1
- 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
【答案】应助回帖
|
明白是怎么回事了. 我把 z2=z1-1 写成 z2=1-z1 了  这样你的问题反而简化了, 因为曲线都是比较简单的. 第一条曲线周期 0-2 Pi 面积: 8768.85621710424, 跟你说的软件的结果很相似. 第二天曲线, 关键问题在于周期不容易表达, 先用牛顿法和合适的初值求出来这个积分区间应该是: 0~5.235987807991032 然后,同样在这个区间上对这个封闭曲线使用格林公式求面积,在我看来用不到任何坐标变换的东西 但是积分的时候所要求的精度不能太高,否则可能报错 6557.98006735277 这个跟你说得软件结果就更加相近了. 看起来软件没什么问题; 至多精度有差异,跟所用算法可能有关系. |
9楼2015-07-20 18:49:32
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24607.5
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.2小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
2楼2015-07-17 11:09:33
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 数学EPI: 2
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
3楼2015-07-17 11:18:14
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 数学EPI: 2
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
4楼2015-07-18 08:32:42