shiyan1.m

数字信号处理课程的有关实验

%实验一:频谱分析与采样定理
%陈友兴
T=0.0001;               %采样间隔T=0.0002
F=1/T;                 %采样频率为F=1/T
L=0.01;                %记录长度L=0.01   
N=L/T;                 
t=0:T:L;             
a=37;
f1=0:F/N:F;
f2=-F/2:F/N:F/2;
%%%%%%%%%%%%%%%%%%%%%%%%%
x1=cos(100*pi*a*t);
y1=T*abs(fft(x1));
y11=fftshift(y1);
figure(1),
subplot(3,1,1),plot(t,x1);title('正弦信号');
subplot(3,1,2),stem(f1,y1);title('正弦信号频谱');
subplot(3,1,3),stem(f2,y11);title('正弦信号频谱');
%%%%%%%%%%%%%%%%%%%%%
x2=exp(-a*t);
y2=T*abs(fft(x2));
y21=fftshift(y2);
figure(2),
subplot(3,1,1),stem(t,x2);title('指数信号');
subplot(3,1,2),stem(f1,y2);title('指数信号频谱');
subplot(3,1,3),stem(f2,y21);title('指数信号频谱');
%%%%%%%%%%%%%%%%%%%%%
x3=x1.*x2;
y3=T*abs(fft(x3));
y31=fftshift(y3);
figure(3),
subplot(3,1,1),stem(t,x3);title('两信号相乘');
subplot(3,1,2),stem(f1,y3);title('两信号相乘频谱');
subplot(3,1,3),stem(f2,y31);title('两信号相乘频谱');