请大虾帮我看下问题在哪？我的QQ：68705890。。不甚感激。。
内容：
[img]http://fm062.img.xiaonei.com/pic001/20070423/1005/L3935888ME.jpg[/img]

[img]http://fm012.img.xiaonei.com/pic001/20070423/1005/L3936074ME.jpg[/img]

我的程序：
clear all

d=8;       %指条宽度
x=32;       %波长
n1=26;            %输入idt指对数
n2=39;            %输出idt指对数
w1=3200;      %输入idt孔径
w2=2550;      %输出idt孔径
f0=122;     %中心频率





for f=0:1:300
    for i=1:26
        x1=(-1)^(i-1)*exp(j*pi*i*(f-f0)/f0);
        y1=sum(x1);
    end
end

for f=0:300    
    for k=1:39
        x2=(-1)^(k-1)*exp(j*pi*k*(f-f0)/f0);
        y2=sum(x2);
    end
    
end

%y3=30*sqrt(-1)*sqrt(1500/2200)/37;

%y4=exp(-sqrt(-1)*2*pi*(32+4522)/x);


z=20*log10(y1*y2);
f=0:1:300;
plot(f,z,'c+')
title('Y型双声路声表面波传感器件Matlab模拟结果')
xlabel('频率(MHz)')
ylabel('插入损耗(dB)')
grid        %网格