| 4 | 1/1 | 返回列表 |
| 查看: 1012 | 回復(fù): 3 | ||
陳昌榮新蟲 (初入文壇)
|
[求助]
求助matlab編程
|
|
下面matlab編程中,經(jīng)過(guò)fft2變換后,k和fftk應(yīng)該怎么表示呢?下面中的k和fftk是有問(wèn)題的,但是我不清楚怎么去正確表示,請(qǐng)指教。 n1=100;n2=100; for x=1:2*n1 for y=1:2*n2 test(x,y)=exp(sin(2*(x-1)*pi/n1+3*(y-1)*pi/n2)); testdx(x,y)=2*cos(2*(x-1)*pi/n1).*exp(sin(2*(x-1)*pi/n1+3*(y-1)*pi/n2)); end end testk=fft2(test,2*n1,2*n2); k=[0:n1-1,-n1:-1];%問(wèn)題在此處以及下一行 fftk=repmat(k,200,1); testdxk = fftk.*testk*i; testdx2 = ifft2(testdxk,2*n1,2*n1); figure; mesh(real(testdx2)); title('real part'); figure; mesh(testdx); title('anaylysis part'); |
木蟲 (職業(yè)作家)
獵手

至尊木蟲 (著名寫手)

新蟲 (初入文壇)
|
n1=100;n2=100; for x=1:2*n1 for y=1:2*n2 test(x,y)=exp(sin(2*(x-1)*pi/n1+3*(y-1)*pi/n2)); testdx(x,y)=2*cos(2*(x-1)*pi/n1).*exp(sin(2*(x-1)*pi/n1+3*(y-1)*pi/n2)); end end testk=fft2(test,2*n1,2*n2); k=[0:n1-1,-n1:-1]; fftk=repmat(k,200,1); testdxk = fftk.*testk*i; testdx2 = ifft2(testdxk,2*n1,2*n1); figure; mesh(real(testdx2)); title('real part'); figure; mesh(testdx); title('anaylysis part'); error =abs(testdx-testdx2); figure; mesh(error); 右邊的括號(hào)是我打錯(cuò)了。這個(gè)程序可以運(yùn)行出來(lái),但問(wèn)題在于結(jié)果不對(duì),就是誤差error很大,正確的誤差在10的負(fù)13次方左右,而上面運(yùn)行出來(lái)的結(jié)果在0到10之間,我經(jīng)過(guò)分段測(cè)試后,發(fā)現(xiàn)問(wèn)題在k和fftk上,我不知道怎么去定義這個(gè)200行200列的矩陣fftk,請(qǐng)指教。 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +3 | 匹克i 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +7 | 生物工程調(diào)劑 2026-03-17 | 12/600 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 0854電子信息求調(diào)劑 +3 | α____ 2026-03-22 | 3/150 |
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無(wú)聊 2026-03-21 | 4/200 |
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|