24小时热门版块排行榜    

查看: 1314  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

朱十三

铁杆木虫 (正式写手)

[求助] matlab做fft和ifft的问题 已有1人参与

帮我看一下程序,我的程序写得是通过公式算得频域信号(E_THz2),然后给所得频域信号做傅立叶逆变换得到时域信号(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=00.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);
我的问题:
1. 为什么程序里n_THz=sqrt(eps) 和n_THz=real(sqrt(eps))+imag(sqrt(eps))算出来的n_THz结果不一样?eps是复数的时候两种写法应该一样才对
2. 我对傅立叶变换和逆变换不懂,请帮我改一下逆变换的程序

matlab做fft和ifft的问题
2.png


matlab做fft和ifft的问题-1
1.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

话说,以前穷人是不取名的,以出生日期命名,如此算来本来该叫廿三。无奈竟忘记还有廿字,遂减去十,故称十三
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wake_up

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
朱十三: 金币+15, ★★★很有帮助 2015-11-11 21:00:22
1.傅立叶的原理就注定了它没有时间定位的能力(小波可以做到时间定位)。
2. fft在一维的情况下, 你给它一个时间序列(离散数值序列),它给你算出一个频谱。
3. ifft,你给它一个频谱, 它能给你还原一段序列。

因为作FFT时你就没办法给时间信息,只能给序列信息。 所以还原出来的序列,你也只能拿你当给的时间参数算。反变换后的序列应该和你原时间序列一致。
扫地扫地扫心地,心地不扫空扫地,人人都把心地扫,世上无处不净地.
4楼2015-11-11 16:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

wake_up

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
朱十三: 金币+5, 有帮助 2015-11-10 21:23:56
eps在matlab里面是不是保留字啊?你拿来当变量用了?

发自小木虫Android客户端
扫地扫地扫心地,心地不扫空扫地,人人都把心地扫,世上无处不净地.
2楼2015-11-08 23:27:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

朱十三

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by wake_up at 2015-11-08 23:27:54
eps在matlab里面是不是保留字啊?你拿来当变量用了?

应该没有问题,之前也用它,麻烦问一下,我有频域数据,用ifft做傅立叶变换后,频率怎么换算成时间?
话说,以前穷人是不取名的,以出生日期命名,如此算来本来该叫廿三。无奈竟忘记还有廿字,遂减去十,故称十三
3楼2015-11-10 21:26:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

骆驼小枣

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by wake_up at 2015-11-11 16:51:26
1.傅立叶的原理就注定了它没有时间定位的能力(小波可以做到时间定位)。
2. fft在一维的情况下, 你给它一个时间序列(离散数值序列),它给你算出一个频谱。
3. ifft,你给它一个频谱, 它能给你还原一段序列。

...

你好!我将用FFT分析信号,能不能教教我怎磨用软件进行FFT变换 谢谢!可以有偿
5楼2016-03-25 09:20:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 299求调剂 +4 shxchem 2026-03-20 5/250 2026-03-21 00:24 by 星空星月
[考研] 317求调剂 +8 申子申申 2026-03-19 13/650 2026-03-21 00:09 by 刘国森
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 324求调剂 +5 lucky呀呀呀鸭 2026-03-20 5/250 2026-03-20 22:30 by 促天成
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
信息提示
请填处理意见