傅里叶变换的线性性质

简介

\qquad 傅里叶变换的线性性质是其最基本的性质之一。假设信号$g(t)$和$h(t)$是两个不同的信号，那么它们在时域作线性变换组合之后，等效在频域作线性变换组合，用公式可以表示为

$$\Im [A\ast g(t)+B\ast h(t)]=A\ast \Im [g(t)]+B\ast \Im [h(t)]\text{\hspace{1em}(1)}$$

其中，A和B都是常量，$\Im$是傅里叶变换算符，

代码

%%%%%%%%%%%%验证傅里叶变换的线性性质%%%%%%%%%%%%%%
clc;
clear;
close all;

%% 产生余弦信号g
gA0 = 1;%      幅值是1
gFre = 1e3;%    频率1kHz
gSampleFre = 1e4;% 采样频率10kHz
M = 100;%     采样点数1000
gt = 0:1/gSampleFre:(M-1)/gSampleFre;%时间序列
g = gA0*cos(2*pi*gFre*gt);%信号g
g_FFT = fftshift(fft(g))/M;%信号g的频谱

%% 产生余弦信号h
hB0 = 1;%      幅值是1
hFre = 2e3;%    频率2kHz
hSampleFre = 1e4;% 采样频率10kHz
ht = 0:1/hSampleFre:(M-1)/hSampleFre;%时间序列
h = hB0*cos(2*pi*hFre*ht);%信号h
h_FFT = fftshift(fft(h))/M;%信号h的频谱

%% 添加线性因子
A = 5;
B = 8;
Ag = A*g;
Bh = B*h;
AgBh = Ag + Bh;
AgBh_FFT = fftshift(fft(AgBh))/M;%线性变换后信号的频谱
AgFFTBhFFT = A*g_FFT+B*h_FFT;%频谱的线性变换

%% 频域坐标
Ft = -gSampleFre/2:gSampleFre/M:gSampleFre/2-gSampleFre/M;

%% 验证线性定理
figure('Name','原始信号');
plot(gt,g,'r',ht,h,'b');
xlabel('时间(S)');ylabel('振幅');legend('g','h');

figure('Name','原始信号频谱');
plot(Ft,abs(g_FFT),'r',Ft,abs(h_FFT),'b');
xlabel('频率(Hz)');ylabel('振幅');legend('G','H');

figure('Name','线性变换的信号频谱');
plot(Ft,abs(AgBh_FFT),'r-o',Ft,abs(AgFFTBhFFT),'g-d');
xlabel('频率(Hz)');ylabel('振幅');legend('线性变换的频谱','频谱的线性变换');

结果

1.原始信号

2.原始信号频谱

3.性质验证