dd.m

来自「很多matlab的源代码」· M 代码 · 共 122 行

switch get (g2,'value')
case 1
subplot(1,1,1);
N=5;
k=0:1:N-1;
wk=2*pi*k/N;
zk=exp(j*wk);
Xk=(zk)./(zk-0.7);
xn=real(idfs(Xk,N));
xtilde=xn'*ones(1,40/N);
xtilde=(xtilde(:))';
subplot(2,2,1);
h=stem(0:39,xtilde);
set(h,'LineWidth',2);
set(h,'MarkerSize',2);
%axis(0,40,-0.1,1.5)
xlabel('n');
ylabel('xtilde(n)');
title('N=5')
e1={'x(n) = (0.7)^n*u(n);'
    'N=5;'
    '例：'
    ' N=5;'
    'k=0:1:N-1;'
    'wk=2*pi*k/N;'
    'zk=exp(j*wk);'
    'Xk=(zk)./(zk-0.7);'
    ' xn=real(idfs(Xk,N));'
    'xtilde=xn‘*ones(1,40/N);'
    'xtilde=(xtilde(:))’'};
set(g3,'string',e1);
case 2
subplot(1,1,1);
N=10;
k=0:1:N-1;
wk=2*pi*k/N;
zk=exp(j*wk);
Xk=(zk)./(zk-0.7);
xn=real(idfs(Xk,N));
xtilde=xn'*ones(1,40/N);
xtilde=(xtilde(:))';
subplot(2,2,1);
h=stem(0:39,xtilde);
set(h,'LineWidth',2);
set(h,'MarkerSize',2);
%axis(0,40,-0.1,1.5)
xlabel('n');
ylabel('xtilde(n)');
title('N=10')
e2={'x(n) = (0.7)^n*u(n);'
    'N=10;'
    '例：'
    ' N=10;'
    'k=0:1:N-1;'
    'wk=2*pi*k/N;'
    'zk=exp(j*wk);'
    'Xk=(zk)./(zk-0.7);'
    'xn=real(idfs(Xk,N));'
    'xtilde=xn‘*ones(1,40/N);'
    'xtilde=(xtilde(:))’'};
set(g3,'string',e2);
case 3
subplot(1,1,1);
N=20;
k=0:1:N-1;
wk=2*pi*k/N;
zk=exp(j*wk);
Xk=(zk)./(zk-0.7);
xn=real(idfs(Xk,N));
xtilde=xn'*ones(1,40/N);
xtilde=(xtilde(:))';
subplot(2,2,1);
h=stem(0:39,xtilde);
set(h,'LineWidth',2);
set(h,'MarkerSize',2);
%axis(0,40,-0.1,1.5)
xlabel('n');
ylabel('xtilde(n)');
title('N=20')
e3={'x(n) = (0.7)^n*u(n);'
    'N=20;'
    '例：'
    ' N=20;'
    'k=0:1:N-1;'
    'wk=2*pi*k/N;'
    'zk=exp(j*wk);'
    'Xk=(zk)./(zk-0.7);'
    'xn=real(idfs(Xk,N));'
    'xtilde=xn‘*ones(1,40/N);'
    'xtilde=(xtilde(:))’'};
set(g3,'string',e3);
case 4
subplot(1,1,1);
N=40;
k=0:1:N-1;
wk=2*pi*k/N;
zk=exp(j*wk);
Xk=(zk)./(zk-0.7);
xn=real(idfs(Xk,N));
xtilde=xn'*ones(1,40/N);
xtilde=(xtilde(:))';
subplot(2,2,1);
h=stem(0:39,xtilde);
set(h,'LineWidth',2);
set(h,'MarkerSize',2);
%axis(0,40,-0.1,1.5)
xlabel('n');
ylabel('xtilde(n)');
title('N=40')
e4={'x(n) = (0.7)^n*u(n);'
    'N=40;'
    '例：'
    ' N=40;'
    'k=0:1:N-1;'
    'wk=2*pi*k/N;'
    'zk=exp(j*wk);'
    '|Xk=(zk)./(zk-0.7);'
    'xn=real(idfs(Xk,N));'
    'xtilde=xn‘*ones(1,40/N);'
    'xtilde=(xtilde(:))’'};
set(g3,'string',e4);
end