Tp=
F1=
F2=
Fs=
t=
Sig=exp(li*)
Spectro=
fn=
plot(fn,Spectro)
sig_i=ifft(Spectro)
t_i=
plot(t_i,sig_i)
```
Tp=
F1=
F2=
Fs=
t=
Sig=exp(li*)
Spectro=
fn=
plot(fn,Spectro)
sig_i=ifft(Spectro)
t_i=
plot(t_i,sig_i)
```
你好,代码供参考:
Tp=2;%信号时长
F1=20;% 给定两个频率信号
F2=40;
Fs=200;%采样频率
t=0:1/Fs:Tp;%时间向量
Sig=exp(1i*2*pi*F1*t)+exp(1i*2*pi*F2*t);%信号
Spectro=abs(fft(Sig, 2^nextpow2(numel(t)),2)/numel(Sig));%获取幅值
fn= (0:numel(Spectro)-1)*Fs/numel(Spectro);%获取频率轴
figure(1)%画图用
plot(fn,Spectro)
sig_i=ifft(Spectro);%逆变换
t_i= (1:numel(sig_i))*Tp/numel(sig_i);
figure(2)%画图用
plot(t_i,sig_i)