ex1_2.m

数字信号处理 实验 的matlab 源码

%1-2
% Q1.6
clf
c=-(1/12)+(pi/6)*i;
k=2;
n=0:40;
x=k*exp(c*n);
figure(1)
subplot(2,2,1),stem(n,real(x))
xlabel('时间序号 n'),ylabel('振幅'),title('实部 c=-(1/12)+(pi/6)*i')
subplot(2,2,2),stem(n,imag(x))
xlabel('时间序号 n'),ylabel('振幅'),title('虚部 c=-(1/12)+(pi/6)*i')

% Q1.8
c=(1/12)+(pi/6)*i;
k=2;
n=0:40;
x=k*exp(c*n);
subplot(2,2,3),stem(n,real(x))
xlabel('时间序号 n'),ylabel('振幅'),title('实部 c=(1/12)+(pi/6)*i')
subplot(2,2,4),stem(n,imag(x))
xlabel('时间序号 n'),ylabel('振幅'),title('虚部 c=(1/12)+(pi/6)*i')

% Q1.11
figure(2)
n=0:35;
a=1.2;
k=0.2;
x=k*a.^n;
subplot(211),stem(n,x),xlabel('时间序号 n'),ylabel('振幅')
% Q1.16
eng=sum(x.*x);
te=strcat('能量为',num2str(eng));
text(5,100,te)

% Q1.14
a1=0.9;
k1=20;
x1=k1*a1.^n;
subplot(212),stem(n,x1),xlabel('时间序号 n'),ylabel('振幅')
% Q1.16
eng=sum(x.*x);
te=strcat('能量为',num2str(eng));
text(25,10,te)