ex1.m

来自「在数字信号处理中的循环移位程序和门函数的源代码」· M 代码 · 共 25 行

clear all;clc;
wl_delt = -8:0.1:8;
m2 = -8:1:8;
y = zeros(length(m2),length(wl_delt));
 for i = 1:length(wl_delt)
     for j =1:length(m2)
         if (wl_delt(i) - m2(j))~=0
             y(j,i) = (sin(pi*(wl_delt(i) - m2(j)))/(pi*(wl_delt(i) - m2(j))))^2;
         end
     end
 end
 for i = 1:length(m2)
     plot(wl_delt,y(i,:));
     hold on;
 end
 yy = zeros(length(m2),length(m2));
 for j = 1:length(m2)
     for k = 1:length(m2)
         for i = 1:length(wl_delt) 
             yy(j,k) = yy(j,k) + y(j,i)* y(k,i);
         end
     end
 end
 yy;