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個點(diǎn)的形式給出的：
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)濟(jì)功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有292人回復(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)有26人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有31人回復(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進(jìn)行的，會不會有版權(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

來給你添點(diǎn)人氣，我是化工的，學(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ù)此樓

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

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

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

月只藍(lán)

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

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

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

是的，有些大意，呵呵

贊一下

回復(fù)此樓

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

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

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

wurongjun

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

專家經(jīng)驗(yàn): +831
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時
蟲號: 1368482
注冊: 2011-08-14
性別: GG
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
管轄: 數(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)該是把每一個小梯形的計(jì)算值加起來!步長不一致!只能逐項(xiàng)加!不能用trapz

贊一下(1人)

回復(fù)此樓

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

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

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

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

[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 085600材料與化工求調(diào)劑 +5	緒幸與子 2026-03-17	5/250	2026-03-17 16:40 by laoshidan
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎兩個實(shí)用型專利一項(xiàng)省級立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級證書40+數(shù)一英一287 +3	小李0854 2026-03-16	3/150	2026-03-17 13:40 by 熱情沙漠
[考研] 271求調(diào)劑 +12	生如夏花… 2026-03-11	14/700	2026-03-17 10:56 by lovewei0727
[考研] 274求調(diào)劑 +5	時間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[教師之家] 焦慮 +7	水冰月月野兔 2026-03-13	9/450	2026-03-16 10:00 by Quakerbird
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6	步川酷紫123 2026-03-13	6/300	2026-03-13 21:59 by 星空星月
[考研] 材料與化工085600調(diào)劑求老師收留 +9	jiaanl 2026-03-11	9/450	2026-03-13 20:22 by JourneyLucky
[考研] 307求調(diào)劑 +5	超級伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 工科調(diào)劑 +4	Jiang191123！ 2026-03-11	4/200	2026-03-13 15:15 by Miko19
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike