图1-3 抽样信号
【练一练】
用MATLAB信号工具箱中的pulstran函数产生冲激串的信号。
T = 0:1/50E3:10E-3; D = [0:1/1E3:10E-3;.^(0:10)]'; Y = pulstran(T,D,'gauspuls',10E4,; plot(T,Y)
【实验心得】
通过此次试验,首先,让我对MATLAB强大的功能有了进一步的了解。其次,也学会了一些常用信号的表示方法。通过自己动手操作,我知道了pulstran函数的调用方法,可以自行画出冲击串函数。
实验二 信号的Fourier分析
【实验目的】
1) 通过计算周期方波信号的Fourier级数,进一步掌握周期信号Fourier级数的计算方法。
2) 通过求解非周期方波信号的Fourier变换,进一步掌握非周期信号Fourier变换的求解方法。 【实验环境】
装有或以上版本的PC机。 【实验内容】
1. 连续时间周期方波信号及其傅里叶级数计算的程序代码,其结果如图2-1所示。 dt = ; T =2; t =-4:dt:4; x1=rectpuls( ,1); x=0;
for m = -1:1
% 时间变量变化步长 % 定义信号的周期 % 定义信号的时间变化范围
% 产生1个周期的方波信号
w0 = 2*pi/T; % 定义信号的频率
% 扩展1个周期的方波信号
x = x+rectpuls(*T-dt),1); % 产生周期方波信号 end
subplot(221); plot(t,x);
axis([-4 4 0 ]); % 设定坐标变化范围 title('周期方波信号') N=10;
% 定义需要计算的谐波次数为10
% 求得Fourier系数ak
for k=-N : N
ak(N+1+k) = x1*exp(-j*k*w0*t') *dt/T; end k=-N:N; subplot(212);
stem(k,abs(ak),'k.'); title('傅里叶级数');
% 绘制幅度谱
图2-1 连续时间周期方波信号及其Fourier级数
2. 非周期连续时间信号及其Fourier变换的程序代码,其结果如图2-2所示。 width=1; t=-5::5;
y=rectpuls(t,width); subplot(221); plot(t,y); ylim([-1 2]);
% 矩形脉冲信号
% 限定y坐标的范围
title('矩形脉冲信号'); Y=fft(y,1024); Y1=fftshift(Y); subplot(212); plot(abs(Y1));
title('傅里叶变换');
% 快速Fourier变换 % 将频谱分量集中
图2-2 非周期连续时间信号及其Fourier变换
【实验心得】
这次实验是信号的Fourier分析。通过此次实验,我进一步掌握周期信号Fourier级数的计算方法和非周期信号Fourier变换的求解方法。可以通过MATLAB来自己画出要求的图形,对老师的代码也掌握了。