 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1998 | 回復(fù): 6 | |||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||||
apollosun9283金蟲 (小有名氣)
劍橋公爵
|
[求助]
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ù)
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有278人回復(fù)
-
哈爾濱理工大學物理系招收物理學考研調(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ù)
-
光學工程學碩調(diào)劑信息
已經(jīng)有26人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有34人回復(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ù)
wurongjun
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗: +831 - 應(yīng)助: 791 (博后)
- 貴賓: 0.308
- 金幣: 24609
- 散金: 310
- 紅花: 75
- 帖子: 3004
- 在線: 881.4小時
- 蟲號: 1368482
- 注冊: 2011-08-14
- 性別: GG
- 專業(yè): 計算數(shù)學與科學工程計算
- 管轄: 數(shù)學
【答案】應(yīng)助回帖
★ ★
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 所以你的結(jié)果不對! 應(yīng)該是把每一個小梯形的計算值加起來!步長不一致!只能逐項加!不能用trapz |
7樓2013-11-23 11:28:29
wucheng0107
金蟲 (小有名氣)
開心就好
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1455.3
- 散金: 765
- 紅花: 26
- 帖子: 250
- 在線: 325.2小時
- 蟲號: 1397472
- 注冊: 2011-09-10
- 性別: GG
- 專業(yè): 生物化工與食品化工
2樓2012-09-26 08:35:02
328911052
金蟲 (小有名氣)
- 應(yīng)助: 15 (小學生)
- 金幣: 2701
- 帖子: 190
- 在線: 139.9小時
- 蟲號: 1023892
- 注冊: 2010-05-20
- 性別: GG
- 專業(yè): 化學反應(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
感謝參與，應(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 |
3樓2012-09-26 09:37:04
月只藍
主管區(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è): 宇宙學
- 管轄: 計算模擬區(qū)
4樓2012-09-26 10:44:15
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +6 | yy要上岸呀 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 【0856】化學工程（085602）313 分，本科學科評估A類院�；瘜W工程與工藝，誠求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
|
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進步 2026-03-11 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 一志愿，福州大學材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[基金申請]
有必要更換申報口嗎
 20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學-物理化學 +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] 四川大學085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|
|
|
[考博] 2026年博士申請 +3 | QwQwQW10 2026-03-11 | 3/150 |
|