求助

求助
这是一段四个级联法珀腔的matlab程序，要求转换到频谱，再滤掉高频信号，剩下的再转换回时域

syms r1 r2 r3  r4 r5 lambda l1 l2 l3 l4;

t1 = sqrt(1-r1^2);

t2 = sqrt(1-r2^2);

t3 = sqrt(1-r3^2);

t4 = sqrt(1-r4^2);

t5 = sqrt(1-r5^2);

theta1 = 2*pi*l1/lambda;

theta2 = 2*pi*l2/lambda;

theta3 = 2*pi*l3/lambda;

theta4 = 2*pi*l4/lambda;

Cav1 = [exp(-1i*theta1), -r1*exp(1i*theta1);-r1*exp(-1i*theta1),exp(+1i*theta1)]/t1;

Cav2 = [exp(-1i*theta2), -r2*exp(1i*theta2);-r2*exp(-1i*theta2),exp(+1i*theta2)]/t2;

Cav3 = [exp(-1i*theta3), -r3*exp(1i*theta3);-r3*exp(-1i*theta3),exp(+1i*theta3)]/t3;

Cav4 = [exp(-1i*theta4), -r4*exp(1i*theta4);-r4*exp(-1i*theta4),exp(+1i*theta4)]/t4;

X = Cav1*Cav2*Cav3*Cav4*t1*t2*t3*t4;

t = t1*t2*t3*t4*t5/(X(1,1)-r5*X(1,2));

xlam = (1500:0.1:1600)*1E-9;

tnum = subs(t,{r1,r2,r3,r4,r5,l1,l2,l3,l4,lambda},{0.05,0.05,0.333,0.306,0.306,5E-4,5E-4,26E-6,3E-4,xlam});

tnumdbl = double(tnum);

plot(xlam*1E9, abs(tnumdbl))

谢谢