| 查看: 1283 | 回复: 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=0  0.01*10^12*2*pi)10*10^12*2*pi); % angular frequencyeps=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. 我对傅立叶变换和逆变换不懂，请帮我改一下逆变换的程序  2.png  1.png |
回复此楼» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB做傅立叶逆变换和傅立叶逆变换
已经有1人回复
- Fortran与MATLAB语言
已经有2人回复
- 求matlab做功率谱分析的大神帮助
已经有3人回复
- matlab2014导入txt坐标点后,怎么进行FFT变换并得出变换后的图形?
已经有8人回复
- 关于matlab的fft函数变换产生的相位谱与理论不相符的问题
已经有12人回复
- MATLAB里的FFT。离散快速傅里叶变换
已经有2人回复
- MATLAB里的FFT。离散快速傅里叶变换。
已经有0人回复
- 在MATLAB中重叠相加FFT滤波器
已经有0人回复
- 做快速傅里叶变换的软件FFT
已经有8人回复
- 菜鸟一枚，求matlab相关指导
已经有9人回复
- matlab 频谱分析
已经有12人回复
- 求救, 用Matlab fft作光学仿真
已经有11人回复
- matlab中的fft2问题
已经有0人回复
- MATLAB傅里叶变换求助
已经有7人回复
- 关于matlab进行傅里叶变换和逆变换的程序，求大神指导！
已经有9人回复
- 如何用matlab对图片中曲线做FFT变换
已经有0人回复
- 有没有哪位大牛用FDTD做过二维光子晶体的透射率的？用matlab编程实现的
已经有6人回复
- OFDM中使用 FFT IP核，用Xilinx FPGA
已经有1人回复
- 在matlab中求取均方根值是不是要先做傅里叶变化呢 还是直接在时域里求
已经有0人回复
- fft变换后在图中显示峰值和频率
已经有2人回复
wake_up
木虫 (著名写手)
- 应助: 124 (高中生)
- 金币: 4001.5
- 散金: 1654
- 红花: 24
- 帖子: 1113
- 在线: 424.2小时
- 虫号: 670119
- 注册: 2008-12-08
- 性别: GG
- 专业: 医学图像数据处理与分析
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
朱十三: 金币+5, ★有帮助 2015-11-10 21:23:56
感谢参与，应助指数 +1
朱十三: 金币+5, ★有帮助 2015-11-10 21:23:56
|
eps在matlab里面是不是保留字啊？你拿来当变量用了？ 发自小木虫Android客户端 |
2楼2015-11-08 23:27:54
朱十三
铁杆木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 7680.7
- 散金: 13
- 红花: 19
- 帖子: 909
- 在线: 605.8小时
- 虫号: 1541839
- 注册: 2011-12-18
- 专业: 光学
3楼2015-11-10 21:26:50
wake_up
木虫 (著名写手)
- 应助: 124 (高中生)
- 金币: 4001.5
- 散金: 1654
- 红花: 24
- 帖子: 1113
- 在线: 424.2小时
- 虫号: 670119
- 注册: 2008-12-08
- 性别: GG
- 专业: 医学图像数据处理与分析
4楼2015-11-11 16:51:26
5楼2016-03-25 09:20:16