_強(qiáng)鵬

銀蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 395
散金: 100
紅花: 1
帖子: 15
在線: 33.4小時(shí)
蟲號(hào): 2132734
注冊(cè): 2012-11-17
性別: GG
專業(yè): 空間環(huán)境和空間天氣

[求助] 關(guān)于matlab進(jìn)行傅里葉變換和逆變換的程序，求大神指導(dǎo)！

模型為：
k1=1.3;
a=4e7;
b=6e8;
Ep=50e3;
dt=1.5e-9;
T=1.5e-7;
t=0:dt:T;
x_t=k1*Ep*(exp(-a.*t)-exp(-b.*t));

1.做fft后得到頻域；

2.fft得到頻域的值得到乘以系數(shù)（可以是復(fù)數(shù)也可以是實(shí)數(shù)，系數(shù)的模長(zhǎng)小于1）得到另外的一組值；

3將得到的值做ifft得到時(shí)域圖。

求大神寫上程序�。�！

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有88人回復(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ù)
請(qǐng)問(wèn)還有沒(méi)有用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)價(jià)值貼推薦，對(duì)您同樣有幫助:

求教matlab傅里葉變換求解已經(jīng)有5人回復(fù)
【求助MATLAB大神】MATLAB畫多行圖時(shí)，每行的高度不太一樣已經(jīng)有14人回復(fù)
MATLAB微分方程參數(shù)擬合問(wèn)題，求大神已經(jīng)有7人回復(fù)
求助大神：要做參數(shù)回歸，fortran和matlab，該選擇哪個(gè)？已經(jīng)有14人回復(fù)
求大神幫我改一下下面的程序！matlab迭代法的不知道錯(cuò)在哪里！已經(jīng)有7人回復(fù)
MATLAB關(guān)于lsqcurvefit程序幾個(gè)錯(cuò)誤，求大神指導(dǎo)一下已經(jīng)有6人回復(fù)
求助：能不能幫忙寫一個(gè)用matlab求一個(gè)函數(shù)的傅里葉變換的程序已經(jīng)有16人回復(fù)
matlab 傅里葉變換頻譜分析已經(jīng)有12人回復(fù)
如何用matlab進(jìn)行快速傅里葉變換？（請(qǐng)教）已經(jīng)有5人回復(fù)
關(guān)于MATLAB實(shí)現(xiàn)連續(xù)傅里葉變換繪圖的問(wèn)題已經(jīng)有3人回復(fù)

1樓 2013-06-02 20:32:12

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

somomo91

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

專家經(jīng)驗(yàn): +3
計(jì)算強(qiáng)帖: 1
應(yīng)助: 144 (高中生)
金幣: 19353
紅花: 115
帖子: 3854
在線: 1053.2小時(shí)
蟲號(hào): 2399301
注冊(cè): 2013-04-03
專業(yè): 動(dòng)力學(xué)與控制
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+2, 謝謝 2013-06-03 14:44:23

clear all; clc

k1 = 1.3;
a  = 4e7;
b  = 6e8;
Ep = 50e3;
dt = 1.5e-9;
T  = 1.5e-7;
t  = 0:dt:T;
x_t= k1*Ep*(exp(-a.*t)-exp(-b.*t));

X  = fft(x_t);

X2 = X * 0.9;
x_t_2 = ifft(X2);

Fs = 1/dt;             % Sampling frequency
L  = length(t);       % Length of signal
NFFT = 2^nextpow2(L); % Next power of 2 from length of y

f = Fs/2*linspace(0,1,NFFT/2+1);

subplot(3,1,1); plot(t, x_t); xlabel('time')
subplot(3,1,2); plot(f, 2*abs(X(1:NFFT/2+1))) ; xlabel('frequency')
subplot(3,1,3); plot(t, x_t); hold on; plot(t, x_t_2, 'r'); xlabel('time'); legend('original signal', 'reconstructed signal')

贊一下(2人)

回復(fù)此樓

2樓2013-06-03 04:29:43

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

_強(qiáng)鵬

銀蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 395
散金: 100
紅花: 1
帖子: 15
在線: 33.4小時(shí)
蟲號(hào): 2132734
注冊(cè): 2012-11-17
性別: GG
專業(yè): 空間環(huán)境和空間天氣

引用回帖:

2樓: Originally posted by somomo91 at 2013-06-03 04:29:43
clear all; clc

k1 = 1.3;
a  = 4e7;
b  = 6e8;
Ep = 50e3;
dt = 1.5e-9;
T  = 1.5e-7;
t  = 0:dt:T;
x_t= k1*Ep*(exp(-a.*t)-exp(-b.*t));

X  = fft(x_t);

X2 = X * 0.9;
x_t_2 = ifft(X2);

F ...

非常感謝，你這里的取得均勻變換因子0.9，但是因子實(shí)際操作時(shí)是變換的小于1的值，就導(dǎo)致x_t_2 = ifft(X2)為復(fù)數(shù)，那么最后
plot(t, x_t_2, 'r'); 怎么畫，畫abs 還是real，還是直接畫？

贊一下

回復(fù)此樓

3樓2013-06-03 17:59:43

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

somomo91

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

專家經(jīng)驗(yàn): +3
計(jì)算強(qiáng)帖: 1
應(yīng)助: 144 (高中生)
金幣: 19353
紅花: 115
帖子: 3854
在線: 1053.2小時(shí)
蟲號(hào): 2399301
注冊(cè): 2013-04-03
專業(yè): 動(dòng)力學(xué)與控制
管轄: 計(jì)算模擬

引用回帖:

3樓: Originally posted by _強(qiáng)鵬 at 2013-06-03 17:59:43
非常感謝，你這里的取得均勻變換因子0.9，但是因子實(shí)際操作時(shí)是變換的小于1的值，就導(dǎo)致x_t_2 = ifft(X2)為復(fù)數(shù)，那么最后
plot(t, x_t_2, 'r'); 怎么畫，畫abs 還是real，還是直接畫？...

這樣的話，其實(shí)就是個(gè)比較了
你把兩個(gè)都畫出來(lái)，看看那條線更接近理想值，
信號(hào)處理，沒(méi)有絕對(duì)的理論的，

贊一下(1人)

回復(fù)此樓

4樓2013-06-03 18:08:28

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

_強(qiáng)鵬

銀蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 395
散金: 100
紅花: 1
帖子: 15
在線: 33.4小時(shí)
蟲號(hào): 2132734
注冊(cè): 2012-11-17
性別: GG
專業(yè): 空間環(huán)境和空間天氣

引用回帖:

其實(shí)系數(shù)是跟頻率有關(guān)系的，可以定為0.9*f
X2 = X * 0.9.*f;

贊一下

回復(fù)此樓

5樓2013-06-03 18:19:22

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

somomo91

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

專家經(jīng)驗(yàn): +3
計(jì)算強(qiáng)帖: 1
應(yīng)助: 144 (高中生)
金幣: 19353
紅花: 115
帖子: 3854
在線: 1053.2小時(shí)
蟲號(hào): 2399301
注冊(cè): 2013-04-03
專業(yè): 動(dòng)力學(xué)與控制
管轄: 計(jì)算模擬

引用回帖:

5樓: Originally posted by _強(qiáng)鵬 at 2013-06-03 18:19:22
其實(shí)系數(shù)是跟頻率有關(guān)系的，可以定為0.9*f
X2 = X * 0.9.*f;...

不能直接定義為，
X2 = X * 0.9.*f;
否則你可以看看結(jié)果。。。

贊一下

回復(fù)此樓

6樓2013-06-03 18:22:48

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

_強(qiáng)鵬

銀蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 395
散金: 100
紅花: 1
帖子: 15
在線: 33.4小時(shí)
蟲號(hào): 2132734
注冊(cè): 2012-11-17
性別: GG
專業(yè): 空間環(huán)境和空間天氣

引用回帖:

6樓: Originally posted by somomo91 at 2013-06-03 18:22:48
不能直接定義為，
X2 = X * 0.9.*f;
否則你可以看看結(jié)果。。。...

我知道具體的公式太復(fù)雜了，反正系數(shù)是隨頻率變換的小于1的值（可以是實(shí)數(shù)，也可以是復(fù)數(shù)），所以ifft得到的為復(fù)數(shù)，我就想知道復(fù)數(shù)怎么畫圖，畫abs貌似不對(duì)，直接畫，顯示忽略了相位

贊一下

回復(fù)此樓

7樓2013-06-03 19:06:17

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

張石平

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 16.5
帖子: 22
在線: 46.7小時(shí)
蟲號(hào): 2216925
注冊(cè): 2013-01-01
專業(yè): 巖土與基礎(chǔ)工程

引用回帖:

請(qǐng)問(wèn)一下，如先有頻域曲線，想得到時(shí)域曲線，做傅里葉逆變換后，時(shí)間項(xiàng)該如何處理？

贊一下

回復(fù)此樓

8樓2014-08-14 10:57:06

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

活在當(dāng)下bj

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 249.1
散金: 5
帖子: 95
在線: 12.3小時(shí)
蟲號(hào): 3212999
注冊(cè): 2014-05-17
專業(yè): 結(jié)構(gòu)工程

引用回帖:

8樓: Originally posted by 張石平 at 2014-08-14 10:57:06
請(qǐng)問(wèn)一下，如先有頻域曲線，想得到時(shí)域曲線，做傅里葉逆變換后，時(shí)間項(xiàng)該如何處理？...

問(wèn)下，直接把頻域曲線經(jīng)逆傅里葉變換得到時(shí)域曲線，您是怎么處理的？謝謝~

贊一下

回復(fù)此樓

9樓2015-05-31 09:24:26

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

張石平

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 16.5
帖子: 22
在線: 46.7小時(shí)
蟲號(hào): 2216925
注冊(cè): 2013-01-01
專業(yè): 巖土與基礎(chǔ)工程

引用回帖:

9樓: Originally posted by 活在當(dāng)下bj at 2015-05-31 09:24:26
問(wèn)下，直接把頻域曲線經(jīng)逆傅里葉變換得到時(shí)域曲線，您是怎么處理的？謝謝~...

直接用ifft函數(shù)就可以得到，但是我的時(shí)間項(xiàng)一直沒(méi)有得到有意義的處理，惆悵中

贊一下

回復(fù)此樓

10樓2015-06-01 17:49:31

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

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

[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4	share_joy 2026-03-16	4/200	2026-03-17 21:32 by hmn_wj
[考研] 085601材料工程專碩求調(diào)劑 +5	慕寒mio 2026-03-16	5/250	2026-03-17 21:31 by hmn_wj
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 0703化學(xué)調(diào)劑 +6	妮妮ninicgb 2026-03-15	9/450	2026-03-16 16:40 by houyaoxu
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 341求調(diào)劑 +4	番茄頭--- 2026-03-10	4/200	2026-03-13 23:12 by JourneyLucky
[考研] 材料工程調(diào)劑 +4	咪咪空空 2026-03-11	4/200	2026-03-13 19:57 by JourneyLucky
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 考研調(diào)劑 +4	芬達(dá)46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 大連大學(xué)化學(xué)專業(yè)研究生調(diào)劑 +3	琪久. 2026-03-10	8/400	2026-03-11 10:02 by 琪久.