| 查看: 2155 | 回复: 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 |
回复此楼» 猜你喜欢
人气不行了
已经有10人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有6人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
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
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 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
