| 4 | 1/1 | 返回列表 |
| 查看: 1015 | 回復(fù): 3 | ||
陳昌榮新蟲 (初入文壇)
|
[求助]
求助matlab編程
|
|
下面matlab編程中,經(jīng)過fft2變換后,k和fftk應(yīng)該怎么表示呢?下面中的k和fftk是有問題的,但是我不清楚怎么去正確表示,請指教。 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'); |
木蟲 (職業(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); 右邊的括號是我打錯了。這個程序可以運(yùn)行出來,但問題在于結(jié)果不對,就是誤差error很大,正確的誤差在10的負(fù)13次方左右,而上面運(yùn)行出來的結(jié)果在0到10之間,我經(jīng)過分段測試后,發(fā)現(xiàn)問題在k和fftk上,我不知道怎么去定義這個200行200列的矩陣fftk,請指教。 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工學(xué)調(diào)劑 +10 | 用戶573181 2026-03-20 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +7 | holy J 2026-03-21 | 7/350 |
|
|
[考博] 申博26年 +4 | 八6八68 2026-03-19 | 4/200 |
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|