當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

朱十三

應(yīng)助: 8 (幼兒園)
金幣: 7643.7
散金: 13
紅花: 19
帖子: 909
在線: 606小時(shí)
蟲(chóng)號(hào): 1541839
注冊(cè): 2011-12-18
專業(yè): 光學(xué)

[求助] matlab做fft和ifft的問(wèn)題已有1人參與

幫我看一下程序,我的程序?qū)懙檬峭ㄟ^(guò)公式算得頻域信號(hào)（E_THz2），然后給所得頻域信號(hào)做傅立葉逆變換得到時(shí)域信號(hào)（E_t）
我程序如下：

clear all;
clc;

n_g=3.24;                               % optical index
tao_F=100*10^(-15);             % full width at half maximum
tao=tao_F/(2*(log(2))^(1/2));  % full width at 1/e maximum
c=3*10^8;                            % light velocity

omega=0

0.01*10^12*2*pi)

10*10^12*2*pi);       % angular frequency
eps=6+((3.92*((5.32*2*pi*10^12)^2)))./((5.32*2*pi*10^12)^2-(omega).^2+(2*i*2*pi*10^12).*omega); % dielecreic function at THz range

freq=omega./(2*pi*10^12);                                  % frequency

kai=eps-1;                                                             % second susceptibility

%%%%%%%
%n_THz=sqrt(eps);
n_THz=real(sqrt(eps))+imag(sqrt(eps));                % index of ZnTe at THz range
%%%%%%
z=500*10^(-6);                                                    % crystal thickness

E_THz=sqrt(2)*pi.*kai./(n_THz.^2-n_g^2)*tao.*exp((-1/4).*omega.^2.*tao^2).*...
(0.5.*(1+n_g./n_THz).* exp(i.*omega.*n_THz.*z./c)+0.5.*(1-n_g./n_THz).*exp(-i.*omega.*n_THz.*z./c)-exp(i*omega.*n_g*z./c));
% THz electric field

E_THz2=abs(E_THz).^2; % THz complex electric field modular

E=[freq', E_THz2'];
save('E_THz(omega).dat','E','-ascii');

figure(1)
plot(freq,E_THz2);
xlabel('frequency (THz)');
ylabel('E_T_H_z (arb.u.)');

E_t=ifft(E_THz2); %inverse fourier transform
figure(2)
plot(E_t)；
我的問(wèn)題：
1. 為什么程序里n_THz=sqrt(eps) 和n_THz=real(sqrt(eps))+imag(sqrt(eps))算出來(lái)的n_THz結(jié)果不一樣？eps是復(fù)數(shù)的時(shí)候兩種寫法應(yīng)該一樣才對(duì)
2. 我對(duì)傅立葉變換和逆變換不懂，請(qǐng)幫我改一下逆變換的程序

2.png

1.png

回復(fù)此樓

» 猜你喜歡

329求調(diào)劑已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有11人回復(fù)
一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一）已經(jīng)有3人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有14人回復(fù)
320求調(diào)劑0856 已經(jīng)有3人回復(fù)
307求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑已經(jīng)有15人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有14人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有18人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

MATLAB做傅立葉逆變換和傅立葉逆變換已經(jīng)有1人回復(fù)
Fortran與MATLAB語(yǔ)言已經(jīng)有2人回復(fù)
求matlab做功率譜分析的大神幫助已經(jīng)有3人回復(fù)
matlab2014導(dǎo)入txt坐標(biāo)點(diǎn)后,怎么進(jìn)行FFT變換并得出變換后的圖形? 已經(jīng)有8人回復(fù)
關(guān)于matlab的fft函數(shù)變換產(chǎn)生的相位譜與理論不相符的問(wèn)題已經(jīng)有12人回復(fù)
MATLAB里的FFT。離散快速傅里葉變換已經(jīng)有2人回復(fù)
MATLAB里的FFT。離散快速傅里葉變換。已經(jīng)有0人回復(fù)
在MATLAB中重疊相加FFT濾波器已經(jīng)有0人回復(fù)
做快速傅里葉變換的軟件FFT 已經(jīng)有8人回復(fù)
菜鳥(niǎo)一枚，求matlab相關(guān)指導(dǎo) 已經(jīng)有9人回復(fù)
matlab 頻譜分析已經(jīng)有12人回復(fù)
求救, 用Matlab fft作光學(xué)仿真已經(jīng)有11人回復(fù)
matlab中的fft2問(wèn)題已經(jīng)有0人回復(fù)
MATLAB傅里葉變換求助已經(jīng)有7人回復(fù)
關(guān)于matlab進(jìn)行傅里葉變換和逆變換的程序，求大神指導(dǎo)！已經(jīng)有9人回復(fù)
如何用matlab對(duì)圖片中曲線做FFT變換已經(jīng)有0人回復(fù)
有沒(méi)有哪位大牛用FDTD做過(guò)二維光子晶體的透射率的？用matlab編程實(shí)現(xiàn)的已經(jīng)有6人回復(fù)
OFDM中使用 FFT IP核，用Xilinx FPGA 已經(jīng)有1人回復(fù)
在matlab中求取均方根值是不是要先做傅里葉變化呢還是直接在時(shí)域里求已經(jīng)有0人回復(fù)
fft變換后在圖中顯示峰值和頻率已經(jīng)有2人回復(fù)

話說(shuō)，以前窮人是不取名的，以出生日期命名，如此算來(lái)本來(lái)該叫廿三。無(wú)奈竟忘記還有廿字，遂減去十，故稱十三

1樓 2015-11-08 21:22:14

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

駱駝小棗

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 46
帖子: 24
在線: 11.9小時(shí)
蟲(chóng)號(hào): 4074190
注冊(cè): 2015-09-15
專業(yè): 膠體與界面化學(xué)

引用回帖:

4樓: Originally posted by wake_up at 2015-11-11 16:51:26
1.傅立葉的原理就注定了它沒(méi)有時(shí)間定位的能力(小波可以做到時(shí)間定位)。
2. fft在一維的情況下，你給它一個(gè)時(shí)間序列(離散數(shù)值序列)，它給你算出一個(gè)頻譜。
3. ifft，你給它一個(gè)頻譜，它能給你還原一段序列。

...

你好！我將用FFT分析信號(hào)，能不能教教我怎磨用軟件進(jìn)行FFT變換謝謝！可以有償

贊一下

回復(fù)此樓

5樓2016-03-25 09:20:16

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

查看全部 5 個(gè)回答

wake_up

木蟲(chóng) (著名寫手)

應(yīng)助: 124 (高中生)
金幣: 4011.5
散金: 1654
紅花: 24
帖子: 1113
在線: 424.2小時(shí)
蟲(chóng)號(hào): 670119
注冊(cè): 2008-12-08
性別: GG
專業(yè): 醫(yī)學(xué)圖像數(shù)據(jù)處理與分析

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
朱十三: 金幣+5, ★有幫助 2015-11-10 21:23:56

eps在matlab里面是不是保留字��？你拿來(lái)當(dāng)變量用了？

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

掃地掃地掃心地，心地不掃空掃地，人人都把心地掃，世上無(wú)處不凈地．

2樓2015-11-08 23:27:54

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

朱十三

鐵桿木蟲(chóng) (正式寫手)

應(yīng)助: 8 (幼兒園)
金幣: 7643.7
散金: 13
紅花: 19
帖子: 909
在線: 606小時(shí)
蟲(chóng)號(hào): 1541839
注冊(cè): 2011-12-18
專業(yè): 光學(xué)

引用回帖:

2樓: Originally posted by wake_up at 2015-11-08 23:27:54
eps在matlab里面是不是保留字�。磕隳脕�(lái)當(dāng)變量用了？

應(yīng)該沒(méi)有問(wèn)題，之前也用它，麻煩問(wèn)一下，我有頻域數(shù)據(jù)，用ifft做傅立葉變換后，頻率怎么換算成時(shí)間？

贊一下

回復(fù)此樓

話說(shuō)，以前窮人是不取名的，以出生日期命名，如此算來(lái)本來(lái)該叫廿三。無(wú)奈竟忘記還有廿字，遂減去十，故稱十三

3樓2015-11-10 21:26:50

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

wake_up

木蟲(chóng) (著名寫手)

應(yīng)助: 124 (高中生)
金幣: 4011.5
散金: 1654
紅花: 24
帖子: 1113
在線: 424.2小時(shí)
蟲(chóng)號(hào): 670119
注冊(cè): 2008-12-08
性別: GG
專業(yè): 醫(yī)學(xué)圖像數(shù)據(jù)處理與分析

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
朱十三: 金幣+15, ★★★很有幫助 2015-11-11 21:00:22

1.傅立葉的原理就注定了它沒(méi)有時(shí)間定位的能力(小波可以做到時(shí)間定位)。
2. fft在一維的情況下，你給它一個(gè)時(shí)間序列(離散數(shù)值序列)，它給你算出一個(gè)頻譜。
3. ifft，你給它一個(gè)頻譜，它能給你還原一段序列。

因?yàn)樽鱂FT時(shí)你就沒(méi)辦法給時(shí)間信息，只能給序列信息。所以還原出來(lái)的序列，你也只能拿你當(dāng)給的時(shí)間參數(shù)算。反變換后的序列應(yīng)該和你原時(shí)間序列一致。

贊一下

回復(fù)此樓

掃地掃地掃心地，心地不掃空掃地，人人都把心地掃，世上無(wú)處不凈地．

4樓2015-11-11 16:51:26

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

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

[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	5/250	2026-03-19 18:03 by zcl123
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 324分 085600材料化工求調(diào)劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 化學(xué)求調(diào)劑 +3	臨澤境llllll 2026-03-17	4/200	2026-03-19 13:59 by houyaoxu
[考研] 281求調(diào)劑（0805） +9	煙汐憶海 2026-03-16	19/950	2026-03-19 11:42 by laoshidan
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無(wú)際的草原
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽(yáng)CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 085600材料與化工 +5	安全上岸！ 2026-03-16	5/250	2026-03-18 15:33 by cmz0325
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢(mèng)人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo