apollosun9283

金蟲 (小有名氣)

劍橋公爵

應(yīng)助: 0 (幼兒園)
金幣: 831.8
帖子: 181
在線: 104.7小時
蟲號: 699792
注冊: 2009-02-11
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

[求助] MATLAB一個數(shù)值積分的問題

待求參數(shù)an如附件中的公式所示，θ∈[0,2π]，由于公式中的D(θ)是以區(qū)間[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];

  所以只能求得數(shù)值積分。
  我的求解程序如下：
  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

請問我這個程序?qū)�？如果不對，求正確的！
  請大家?guī)蛶兔�，在輸入原始�?shù)據(jù)D的時候，直接在我這里復(fù)制可以了哈！

333333333333333333333333.jpg

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有220人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有29人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有53人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

matlab 積分求助已經(jīng)有8人回復(fù)
有關(guān)MATLAB數(shù)據(jù)的問題已經(jīng)有6人回復(fù)
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復(fù)
求matlab求解一個方程組的問題已經(jīng)有6人回復(fù)
一個簡單的式子，卻算不正確已經(jīng)有10人回復(fù)
各位，在論文里寫明數(shù)值模擬是用matlab進行的，會不會有版權(quán)問題？已經(jīng)有3人回復(fù)
關(guān)于matlab 2010b 的一個小問題已經(jīng)有7人回復(fù)
【求助】matlab求解一個簡單的積分（急！�。。� 已經(jīng)有7人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復(fù)
【求助】問一個關(guān)于matlab中l(wèi)egend的問題已經(jīng)有11人回復(fù)
【求助】請教Matlab高手求定積分近似值問題已經(jīng)有4人回復(fù)

倚樓聽風(fēng)雨，淡看江湖路

1樓 2012-09-25 21:42:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wucheng0107

金蟲 (小有名氣)

開心就好

應(yīng)助: 6 (幼兒園)
金幣: 1455.3
散金: 765
紅花: 26
帖子: 250
在線: 325.2小時
蟲號: 1397472
注冊: 2011-09-10
性別: GG
專業(yè): 生物化工與食品化工

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

來給你添點人氣，我是化工的，學(xué)過Matlab，但是還是有待繼續(xù)努力啊

贊一下

回復(fù)此樓

開心就好！

2樓2012-09-26 08:35:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

328911052

金蟲 (小有名氣)

應(yīng)助: 15 (小學(xué)生)
金幣: 2701
帖子: 190
在線: 139.9小時
蟲號: 1023892
注冊: 2010-05-20
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
dbb627: 金幣+1, 感謝應(yīng)助 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

贊一下

回復(fù)此樓

我不管是否能夠贏得成功，既然選擇了遠方，便只顧風(fēng)雨兼程！

3樓2012-09-26 09:37:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

★
dbb627: 金幣+1, 謝謝應(yīng)助！ 2012-09-26 21:05:49

引用回帖:

3樓: Originally posted by 328911052 at 2012-09-26 09:37:04
clear;clc
D=;

  for n=1:300
   X=linspace(0,2*n*pi,144);
   Y=D.*cos(n*X);
   Z(n)=trapz(X,Y);
  end...

感謝您的參與！
您的程序和我的其實是一樣的。
不過，您在算cos(nθ)的時候，積分上下限取的是[0，2nπ]，是不是取錯了，應(yīng)該是[0，2π]。
此外，最后在調(diào)用梯形積分函數(shù)trapz的時候，您好像漏掉了1/π。
我在您程序上作了上述修改后，我算出的a和您程序算出的Z，吻合得很好。

untitled.jpg

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

4樓2012-09-26 10:44:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huwanpeng2

金蟲 (職業(yè)作家)

應(yīng)助: 421 (碩士)
金幣: 14443.2
紅花: 68
帖子: 4559
在線: 326.2小時
蟲號: 660538
注冊: 2008-11-23
專業(yè): 化學(xué)反應(yīng)工程

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

建議參考黃華江的那本書，化學(xué)工業(yè)出版社，有例子的

贊一下

回復(fù)此樓

skloche.ecust

5樓2012-09-26 11:02:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

328911052

金蟲 (小有名氣)

應(yīng)助: 15 (小學(xué)生)
金幣: 2701
帖子: 190
在線: 139.9小時
蟲號: 1023892
注冊: 2010-05-20
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

4樓: Originally posted by 月只藍 at 2012-09-26 10:44:15
感謝您的參與！
您的程序和我的其實是一樣的。
不過，您在算cos(nθ)的時候，積分上下限取的是，是不是取錯了，應(yīng)該是。
此外，最后在調(diào)用梯形積分函數(shù)trapz的時候，您好像漏掉了1/π。
我在您程序上作了上述修 ...

是的，有些大意，呵呵

贊一下

回復(fù)此樓

我不管是否能夠贏得成功，既然選擇了遠方，便只顧風(fēng)雨兼程！

6樓2012-09-26 12:47:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wurongjun

專家顧問 (職業(yè)作家)

專家經(jīng)驗: +831
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時
蟲號: 1368482
注冊: 2011-08-14
性別: GG
專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
管轄: 數(shù)學(xué)

【答案】應(yīng)助回帖

★ ★
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
所以你的結(jié)果不對!
應(yīng)該是把每一個小梯形的計算值加起來!步長不一致!只能逐項加!不能用trapz

贊一下(1人)

回復(fù)此樓

善惡到頭終有報,人間正道是滄桑.

7樓2013-11-23 11:28:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +10	yangfz 2026-03-17	10/500	2026-03-18 20:14 by walc
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +6	d如愿上岸 2026-03-12	9/450	2026-03-18 19:50 by macy2011
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4	1孫悟空 2026-03-17	4/200	2026-03-18 17:59 by fivewind
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 331求調(diào)劑（0703有機化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 288求調(diào)劑 +4	奇點0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant