Floquet-傅里叶方法里面一个积分的验证小程序

%Floquet-傅里叶-方法里面积分的一个验证小程序

function []=main_jifen()
clear;clc
for jj=-100:1:100
    disp(['jj=',num2str(jj),';',newline])
    mm=jifen(jj)
    disp([newline,'--------------'])



end
end

function [out1]=jifen(delta)
syms vac w t k T
%delta=3;
m=2*vac*cos(w*t);
vt=sym(zeros(2));
vt(1,2)=m-m*exp(-1i*k);
vt(2,1)=m-m*exp(1i*k);


M=(1/T)*vt*exp(1i*w*t*delta);

out=int(M,t,0,T);
out1=simplify(subs(out,w,(2*pi/T)));
end

结果表明，当vt等于这个形式时，除了delta为正负1以外，其他积分都是0矩阵，当然主要是看vt的具体形式，不同情况结果不同。