Matlab怎么求透过率函数的傅里叶变换

如果我求这个环境的频谱图，已经写出透过率函数了，怎么求它的频谱呢？初学，谢谢您的回复。
clc
clear all
delta_n=5.9e-4;
L1=0.3;%两段PMF长度
L2=0.2;
lambda=[1520:0.1:1580]10^-9;
phi1=pidelta_nL1./lambda;
phi2=pidelta_nL2./lambda;
%夹角90度
theta=pi/4;
theta2=pi/2;%保偏光纤之间的夹角
%透过率公式
T=(sin(theta))^2.(cos(theta2))^2.(cos(phi1+phi2)).^2+...
(cos(theta))^2.(sin(theta2))^2.(cos(phi1-phi2)).^2+...
1/4sin(2*(theta)).sin(2theta2).(cos(2phi1)+cos(2phi2));
subplot(2,2,1)
plot(lambda1e6,T*100)
xlabel('Wavelength / μm'),ylabel('Transmission / 100%'),title('夹角90度')

img

参考链接

用matlab实现傅里叶变换,matlab实现傅里叶变换_匿名悄悄的的博客-CSDN博客 一、傅立叶变化的原理；(1)原理正交级数的展开是其理论基础！将一个在时域收敛的函数展开成一系列不同频率谐波的叠加，从而达到解决周期函数问题的目的。在此基础上进行推广，从而可以对一个非周期函数进行时频变换。从分析的角度看，他是用简单的函数去逼近(或代替)复杂函数，从几何的角度看，它是以一族正交函数为基向量，将函数空间进行正交分解，相应的系数即为坐标。从变幻的角度的看，他建立了周期函数与序列之间的对应... https://blog.csdn.net/weixin_27676369/article/details/116062996?ops_request_misc=&request_id=&biz_id=102&utm_term=Matlab%E6%80%8E%E4%B9%88%E6%B1%82%E9%80%8F%E8%BF%87%E7%8E%87%E5%87%BD%E6%95%B0%E7%9A%84%E5%82%85%E9%87%8C%E5%8F%B6%E5%8F%98%E6%8D%A2&utm_medium=distribute.pc_search_result.none-task-blog-2~blog~sobaiduweb~default-0-116062996.nonecase&spm=1018.2226.3001.4450

希望有用
https://b23.tv/4V8KGCh

请用插入代码片重新插入一下代码，直接贴上的代码乘号不显示