应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 请问一下用matlab求封闭曲线的面积，怎么求
【调剂】北京石油化工学院2024年16个专业接受调剂
151/2返回列表12下一页
查看: 8181  |  回复: 14
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

高湘伦

新虫 (初入文坛)

[求助] 请问一下用matlab求封闭曲线的面积，怎么求已有3人参与

x=[-0.0296 -0.0139 0.0017 0.0169 0.0317 0.0460 0.0595 0.0722 0.0839 0.0946 0.1041 0.1125 0.1197 0.1255 0.1300 0.1331 0.1348 0.1351 0.1340 0.1316 0.1277 0.1225 0.1160 0.1082 0.0991 0.0890 0.0777 0.0655 0.0524 0.0385 0.0240 0.0090 -0.0064 -0.0220 -0.0376 -0.0530 -0.0682 -0.0827 -0.0966 -0.1096 -0.1214 -0.1320 -0.1412 -0.1489 -0.1548 -0.1590 -0.1614 -0.1619 -0.1605 -0.1572 -0.1522 -0.1453 -0.1369 -0.1269 -0.1157 -0.1032 -0.0897 -0.0754 -0.0605 -0.0451 -0.0296];
y=[1.0105 1.0208 1.0309 1.0407 1.0500 1.0588 1.0669 1.0743 1.0809 1.0866 1.0914 1.0951 1.0978 1.0995 1.1000 1.0994 1.0978 1.0951 1.0913 1.0866 1.0809 1.0743 1.0669 1.0587 1.0500 1.0406 1.0308 1.0207 1.0104 0.9999 0.9895 0.9792 0.9690 0.9593 0.9499 0.9412 0.9330 0.9256 0.9191 0.9134 0.9086 0.9049 0.9022 0.9005 0.9000 0.9006 0.9022 0.9049 0.9087 0.9134 0.9191 0.9257 0.9331 0.9413 0.9501 0.9594 0.9692 0.9793 0.9896 1.0001 1.0105];
围成的是个椭圆，怎么求面积？是直接用polyarea(x,y)吗？我用c++计算的结果是0.0138542；用polyarea(x,y)计算的结果为0.0132.求问正确结果应该是多少

请问一下用matlab求封闭曲线的面积，怎么求
untitled.jpg
回复此楼

» 猜你喜欢
1楼2019-07-17 12:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lwloveflxgg

禁虫 (知名作家)
★ ★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+4, 感谢指导 2019-07-19 21:21:43
本帖内容被屏蔽
4楼2019-07-19 15:43:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lwloveflxgg

禁虫 (知名作家)
★ ★
高湘伦(月只蓝代发): 金币+2, 感谢指导 2019-07-19 21:21:57
本帖内容被屏蔽
5楼2019-07-19 15:45:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

googon

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
高湘伦(月只蓝代发): 金币+4, 感谢指导 2019-07-19 21:21:36
都是数值结果，看起来差的不是的太多，要是有曲线方程就可精确求解，对于你这个离散点，建议先做插值再进行计算，可能结果更加接近真实值

发自小木虫Android客户端
赞一下(1人) 回复此楼
2楼2019-07-17 15:36:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
高湘伦(月只蓝代发): 金币+4, 感谢指导 2019-07-19 21:21:27
也可以用trapz(x,y),结果是0。0132
赞一下(1人) 回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
3楼2019-07-19 02:45:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美
引用回帖:
8楼: Originally posted by 高湘伦 at 2019-07-25 17:12:26
你能提供一下定积分数值解法的程序代码吗，谢谢
...

如果你一定要自已积分去算，可以用下面代码，按梯形处理。（记得给我花花）

x=[-0.0296 -0.0139 0.0017 0.0169 0.0317 0.0460 0.0595 0.0722 0.0839 0.0946 0.1041 0.1125 0.1197 0.1255 0.1300 0.1331 0.1348 0.1351 0.1340 0.1316 0.1277 0.1225 0.1160 0.1082 0.0991 0.0890 0.0777 0.0655 0.0524 0.0385 0.0240 0.0090 -0.0064 -0.0220 -0.0376 -0.0530 -0.0682 -0.0827 -0.0966 -0.1096 -0.1214 -0.1320 -0.1412 -0.1489 -0.1548 -0.1590 -0.1614 -0.1619 -0.1605 -0.1572 -0.1522 -0.1453 -0.1369 -0.1269 -0.1157 -0.1032 -0.0897 -0.0754 -0.0605 -0.0451 -0.0296]
y=[1.0105 1.0208 1.0309 1.0407 1.0500 1.0588 1.0669 1.0743 1.0809 1.0866 1.0914 1.0951 1.0978 1.0995 1.1000 1.0994 1.0978 1.0951 1.0913 1.0866 1.0809 1.0743 1.0669 1.0587 1.0500 1.0406 1.0308 1.0207 1.0104 0.9999 0.9895 0.9792 0.9690 0.9593 0.9499 0.9412 0.9330 0.9256 0.9191 0.9134 0.9086 0.9049 0.9022 0.9005 0.9000 0.9006 0.9022 0.9049 0.9087 0.9134 0.9191 0.9257 0.9331 0.9413 0.9501 0.9594 0.9692 0.9793 0.9896 1.0001 1.0105]
data=[x' y']
[~,ind1]=max(data(:,1))
[~,ind2]=min(data(:,1))
for ii=1:length(y)-1
    y1(ii)=(y(ii)+y(ii+1))/2
end
x1=diff(x)
data1=[x1',y1']
s1=data1(1:ind1-1,1).*data1(1:ind1-1,2)
S1=sum(s1)
s2=data1(ind1:ind2-1,1).*data1(ind1:ind2-1,2)
S2=sum(s2)
s3=data1(ind2:end,1).*data1(ind2:end,2)
S3=sum(s3)
S=S1+S3+S2
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

高湘伦
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
10楼2019-07-26 04:26:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

匿名

本帖仅楼主可见
6楼2019-07-24 14:41:31
 已阅   申请程序强帖   回复此楼   编辑   查看我的主页

hollya

至尊木虫 (文坛精英)

葛城?{美
引用回帖:
6楼: Originally posted by 在杭州 at 2019-07-24 14:41:31
你怎么这么厉害...

你应该觉得楼上更厉害吖，包括线性变换，积分的基本思想都给你讲啦

发自小木虫IOS客户端
赞一下 回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
7楼2019-07-24 22:58:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

高湘伦

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by lwloveflxgg at 2019-07-19 15:43:32
有两个思路，一个是把它当做一个椭圆（实际上不是椭圆，上半圈向内凹陷，下半圈更像椭圆），用椭圆的参数方程去近似，最后得出长轴a=0.1776，短轴b∈（0.0177，0.03）中值为0.0239 ，利用椭圆面积公式S=pi*a*b
所以 ...

你能提供一下定积分数值解法的程序代码吗，谢谢

发自小木虫Android客户端
赞一下 回复此楼
8楼2019-07-25 17:12:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

高湘伦

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by hollya at 2019-07-19 02:45:55
也可以用trapz(x,y),结果是0。0132

谢谢谢谢

发自小木虫Android客户端
回复此楼
9楼2019-07-25 17:15:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 高湘伦 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] Nature一直在编辑手里，考虑好几天了，是悬了吗 +9 彩虹初见 2024-04-24 9/450 2024-04-25 11:13 by ca0yan9
[考博] 申博求助 +3 dskabdh 2024-04-24 9/450 2024-04-25 08:16 by 湛卢赤霄
[论文投稿] Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？ 10+4 yifeng11 2024-04-20 10/500 2024-04-25 02:14 by yifeng11
[基金申请] 企业博后是否能申请CSC博士后项目？谢谢 66+3 龙包子 2024-04-19 6/300 2024-04-24 21:59 by 3115321
[访问学者] CSC的访问学者申请，没有个评审意见，也不知道怎么改，还有必要申请吗 +4 flyingship 2024-04-20 4/200 2024-04-24 21:09 by 59038mute
[论文投稿] 研二光催化6月底四篇二区什么水平 5+4 wjtab 2024-04-22 10/500 2024-04-24 19:37 by luckyeheng
[有机交流] 紧急求助，有谁用过三甲基硅醇钾TMSOK吗？ 1000+3 genius2008 2024-04-19 10/500 2024-04-24 12:37 by genius2008
[论文投稿] 期刊推荐 20+4 木颜尘ip 2024-04-22 7/350 2024-04-24 10:06 by bobvan
[公派出国] 東京大学 应用化学专业柳田研究室 招收修士/博士生 +5 蝶梦123 2024-04-19 5/250 2024-04-23 22:40 by 90472jvvh
[论文投稿] 无under review，直接DIP，咋么办。 +10 lizhengke06 2024-04-19 16/800 2024-04-23 19:15 by 化学程序员
[教师之家] 发现督导在听课，需要跟TA招呼示意吗？有同学看手机、课件有疏漏算大问题吗？ +11 河西夜郎 2024-04-19 11/550 2024-04-23 15:01 by mylifecg
[电化学] 有偿收购废弃锂硫电池 +3 汤圆喵 2024-04-19 5/250 2024-04-23 10:41 by Happy C
[论文投稿] 一直找不到审稿人 +4 lizhengke06 2024-04-21 5/250 2024-04-22 23:37 by topedit
[教师之家] 河海大学常州校区新进的讲师待遇如何？考核要求如何？ +3 德国小球迷 2024-04-19 18/900 2024-04-22 19:23 by ZHONGWU_U
[高分子] 请问UV灯是365nm的，那么选光引发剂的波长选多少的？要完全一致吗？ +4 engledd2004 2024-04-21 4/200 2024-04-22 16:08 by wangcz23
[考博] 青岛科技大学顾玉超教授课题组招收博士生 +3 yanqingwang 2024-04-19 3/150 2024-04-22 15:21 by 理科生yes
[考研] 问题已经解决。 +13 lekinna 2024-04-19 30/1500 2024-04-22 13:09 by 矛dei到
[论文投稿] 催稿后秒拒 +9 lizhengke06 2024-04-19 11/550 2024-04-21 22:29 by ggjjyyy
[考研] 312求调剂 +3 Leroic000 2024-04-19 3/150 2024-04-20 18:19 by xinaishan
[考研] 考研英语带学 +6 choice666 2024-04-18 6/300 2024-04-19 17:46 by charles-c
信息提示
关闭
请填处理意见
关闭