5t.m

来自「包含50个数字信号处理方面的matlab程序」· M 代码 · 共 27 行

x11=[1 1 1 1];
n=0:5;
x12=cos(pi*n/4);
y11=circonvt(x11,x12,8)
y12=conv(x11,x12)
y13=[y11(1:1:8),zeros(1,1)]
e1=y13-y12
x21=[1 -1 1 -1];
x22=[1 0 -1 0];
y21=circonvt(x21,x22,5)
y22=conv(x21,x22)
y23=[y21(1:1:5),zeros(1,2)]
e2=y23-y22
n=0:15;
x31=cos(2*pi*n/32);
x32=sin(2*pi*n/32);
y31=circonvt(x31,x32,32)
y32=conv(x31,x32)
y33=[y32(1:1:31),zeros(1,1)]
e3=y31-y33
n=0:9;
x41=(0.8).^n;
x42=(-0.8).^n;
y41=circonvt(x41,x42,15)
y42=conv(x41,x42)
y43=[y41(1:1:15),zeros(1,4)]
e4=y43-y42