| 查看: 2036 | 回复: 6 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
apollosun9283金虫 (小有名气)
剑桥公爵
|
[求助]
MATLAB一个数值积分的问题
|
||||
|
待求参数an如附件中的公式所示，θ∈[0,2π]，由于公式中的D(θ)是以区间[0,2π]上均匀的离散的144个点的形式给出的： D=[2.524,2.522,2.518,2.512,2.504,2.495,2.484,2.472,2.459,2.444,2.427,2.41,2.391,2.372,... 2.351,2.329,2.305,2.281,2.256,2.229,2.201,2.172,2.142,2.111,2.079,... 2.046,2.011,1.975,1.938,1.9,1.86,1.82,1.778,1.735,1.69,1.645,1.598,... 1.55,1.501,1.451,1.399,1.347,1.294,1.239,1.184,1.129,1.072,1.015,... 0.958,0.901,0.844,0.786,0.729,0.673,0.618,0.563,0.51,0.459,0.409,... 0.361,0.315,0.272,0.231,0.193,0.158,0.127,0.098,0.074,0.052,0.034,... 0.02,0.01,0.003,0,0,0.004,0.012,0.023,0.037,0.054,0.074,0.097,0.123,... 0.152,0.183,0.216,0.251,0.289,0.328,0.369,0.412,0.456,0.501,0.548,... 0.596,0.645,0.694,0.745,0.797,0.849,0.901,0.954,1.008,1.062,1.116,... 1.17,1.225,1.279,1.334,1.388,1.442,1.496,1.55,1.603,1.655,1.707,1.758,... 1.808,1.858,1.906,1.953,1.999,2.043,2.086,2.127,2.167,2.205,2.241,... 2.275,2.307,2.337,2.365,2.391,2.414,2.435,2.454,2.471,2.485,2.497,... 2.507,2.515,2.52,2.524,2.525]; 所以只能求得数值积分。 我的求解程序如下： b=linspace(0,2*pi,144); for n=1:300 for i=1:144 c(i)=D(i)*cos(n*(i-1)*2*pi/143); end a(n)=1/pi*trapz(b,c) end 请问我这个程序对吗？如果不对，求正确的！ 请大家帮帮忙，在输入原始数据D的时候，直接在我这里复制可以了哈！  333333333333333333333333.jpg |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有18人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有229人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分求助
已经有8人回复
- 有关MATLAB数据的问题
已经有6人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 一个简单的式子，却算不正确
已经有10人回复
- 各位，在论文里写明数值模拟是用matlab进行的，会不会有版权问题？
已经有3人回复
- 关于matlab 2010b 的一个小问题
已经有7人回复
- 【求助】matlab求解一个简单的积分（急！！！）
已经有7人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】问一个关于matlab中legend的问题
已经有11人回复
- 【求助】请教Matlab高手求定积分近似值问题
已经有4人回复
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
【答案】应助回帖
★ ★
dbb627: 金币+2, 感谢参与讨论 2013-11-23 13:45:54
dbb627: 金币+2, 感谢参与讨论 2013-11-23 13:45:54
|
你的D值不是均匀的! 比如取第一行: D=[2.524,2.522,2.518,2.512,2.504,2.495,2.484,2.472,2.459,2.444,2.427,2.41,2.391,2.372]: >> diff(D) ans = Columns 1 through 7 -0.0020 -0.0040 -0.0060 -0.0080 -0.0090 -0.0110 -0.0120 Columns 8 through 13 -0.0130 -0.0150 -0.0170 -0.0170 -0.0190 -0.0190 所以你的结果不对! 应该是把每一个小梯形的计算值加起来!步长不一致!只能逐项加!不能用trapz |
7楼2013-11-23 11:28:29
wucheng0107
金虫 (小有名气)
开心就好
- 应助: 6 (幼儿园)
- 金币: 1455.3
- 散金: 765
- 红花: 26
- 帖子: 250
- 在线: 325.4小时
- 虫号: 1397472
- 注册: 2011-09-10
- 性别: GG
- 专业: 生物化工与食品化工
2楼2012-09-26 08:35:02
328911052
金虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 2701
- 帖子: 190
- 在线: 139.9小时
- 虫号: 1023892
- 注册: 2010-05-20
- 性别: GG
- 专业: 化学反应工程
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
dbb627: 金币+1, 感谢应助 2012-09-26 21:06:30
apollosun9283: 金币+30, ★★★很有帮助 2012-10-15 10:11:14
感谢参与，应助指数 +1
dbb627: 金币+1, 感谢应助 2012-09-26 21:06:30
apollosun9283: 金币+30, ★★★很有帮助 2012-10-15 10:11:14
|
clear;clc D=[2.524,2.522,2.518,2.512,2.504,2.495,2.484,2.472,2.459,2.444,2.427,2.41,2.391,2.372,... 2.351,2.329,2.305,2.281,2.256,2.229,2.201,2.172,2.142,2.111,2.079,... 2.046,2.011,1.975,1.938,1.9,1.86,1.82,1.778,1.735,1.69,1.645,1.598,... 1.55,1.501,1.451,1.399,1.347,1.294,1.239,1.184,1.129,1.072,1.015,... 0.958,0.901,0.844,0.786,0.729,0.673,0.618,0.563,0.51,0.459,0.409,... 0.361,0.315,0.272,0.231,0.193,0.158,0.127,0.098,0.074,0.052,0.034,... 0.02,0.01,0.003,0,0,0.004,0.012,0.023,0.037,0.054,0.074,0.097,0.123,... 0.152,0.183,0.216,0.251,0.289,0.328,0.369,0.412,0.456,0.501,0.548,... 0.596,0.645,0.694,0.745,0.797,0.849,0.901,0.954,1.008,1.062,1.116,... 1.17,1.225,1.279,1.334,1.388,1.442,1.496,1.55,1.603,1.655,1.707,1.758,... 1.808,1.858,1.906,1.953,1.999,2.043,2.086,2.127,2.167,2.205,2.241,... 2.275,2.307,2.337,2.365,2.391,2.414,2.435,2.454,2.471,2.485,2.497,... 2.507,2.515,2.52,2.524,2.525]; for n=1:300 X=linspace(0,2*n*pi,144); Y=D.*cos(n*X); Z(n)=trapz(X,Y); end |
3楼2012-09-26 09:37:04
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2012-09-26 10:44:15
