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)有68人回復(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)有24人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有20人回復(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)驗: +1059
計算強(qiáng)帖: 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 月只藍(lán) at 2012-09-26 10:44:15
感謝您的參與！
您的程序和我的其實是一樣的。
不過，您在算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)驗: +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ā)表

[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 085600材料與化工 +4	安全上岸！ 2026-03-16	4/200	2026-03-17 14:02 by 勇敢太監(jiān)王公公
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9	柴郡貓_ 2026-03-12	9/450	2026-03-17 10:18 by Sammy2
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 341求調(diào)劑 +4	番茄頭--- 2026-03-10	4/200	2026-03-13 23:12 by JourneyLucky
[考研] 0703化學(xué)調(diào)劑 +4	快樂的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 336求調(diào)劑 +6	Iuruoh 2026-03-11	6/300	2026-03-13 22:06 by JourneyLucky
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 【考研調(diào)劑求收留】 +3	Ceciilia 2026-03-11	3/150	2026-03-13 20:18 by JourneyLucky
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao