通信电子线路实验常见错误（A1.2）

1. 容易混淆的字母和数字

字母‘l’和数字‘1’，字母‘o’和数字‘0’，非常容易混淆。编写程序的时候需要仔细核对。

1.1 错误提示信息

1.2 定位出错的代码

2 ‘fftshift’误写成了‘fft’

2.1 实验显现

画出的频谱图如下图所示，和教科书中的图不一致。

2.2 错误原因

2.3 更正

2.3.1 使用fftshift得到的频谱图

2.3.2 不使用fftshift得到的频谱图

3. 完整程序

% 例题 A1.1 对连续时间变化信号进行采样，
% 然后进行fft变换。
clear
clf
f0=2e3;
t0=1/f0;
fs=1e5;
ts=1/fs;
df=1e3;
t=[-t0/2:ts:t0/2];
xn=cos(2*pi*f0*t);
subplot(2,1,1)
stem(t*1000,xn,'fill')
grid on
title('采样时间序列x(n)')
xlabel('T(ms)')
ylabel('Amplitute')
[X,x,dfl]=fftseq(xn,ts,df)
X1=X/fs;
f1=[0:df1:df1*(length(x)-1)]-fs/2;
subplot(2,1,2)
plot(f1/1000,abs(fftshift(X1)))
grid on
title('经补零增长后的FFT结果')
xlabel('f(kHz)')
ylabel('Amplitute')

% 实验中需要使用的fftseq函数

function[M,m,df]=fftseq(m,ts,df)
fs=1/ts;
if nargin==2
    n1=0;
else
    n1=fs/df;
end
n2=length(m);
n=2^(max(nextpow2(n1),nextpow2(n2)));
M=fft(m,n);
m=[m,zeros(1,n-n2)];
df=fs/n;

% 例题 A1.2 对于例题A1.1中的时间序列进行FFT变换，
% 要求频率分辨率为1kHz以下。
clear
clf
f0=2e3;
t0=1/f0;
fs=1e5;
ts=1/fs;
df=1e3;
t=[-t0/2:ts:t0/2];
xn=cos(2*pi*f0*t);

subplot(2,1,1)
stem(t*1000,xn,'fill')
grid on
title('采样时间序列x(n)')
xlabel('T(ms)')
ylabel('Amplitute')


[X,x,df1]=fftseq(xn,ts,df);
X1=X/fs;
f1=[0:df1:df1*(length(x)-1)]-fs/2;
subplot(2,1,2)
plot(f1/1000,abs(fftshift(X1)));                       % plot(f1/1000,abs(fftshift(X1)))
grid on
title('经补零增长后的FFT结果')
xlabel('f(kHz)')
ylabel('Amplitute')