exe0311.m

来自「Matlab程序为清华大学胡广书版信号处理的上机习题答案」· M 代码 · 共 45 行

clear;

xn=[1,2,3,4];
hn=[4,3,2,1];

%yn=xn*hn,卷积
x2nn=[1 0 0 0;
    2 1 0 0;
    3 2 1 0;
    4 3 2 1;
    0 4 3 2;
    0 0 4 3;
    0 0 0 4];
yn=x2nn*hn'


x1=0:1:23;
k=0:1:23;
for kj=0:1:23
    
    xtemp=0;
    for ni=0:1:11
        xtemp=xtemp+(exp(j*pi/6*ni)+exp(-j*pi/6*ni))/2*exp(-j*pi/12*ni*kj);
    end
x1(kj+1)=xtemp
end
x1=abs(x1)
stem(k,x1)

x2=0:1:11;
k2=0:1:11;
for kj=0:1:11
    
    xtemp=0;
    for ni=0:1:11
        xtemp=xtemp+(exp(j*pi/6*ni)+exp(-j*pi/6*ni))/2*exp(-j*pi/6*ni*kj);
    end
x2(kj+1)=xtemp
end
x2=abs(x2)
figure
ylable("我");

stem(k2,x2)