在幅度和相位方面进行频域傅立叶变换后，如何计算传递函数并将其绘制？

Question

在对两个信号进行傅立叶变换之后，我很难按大小和相位来绘制传递函数。

首先，我使用excel xlxs来读取要在时域和频域中绘制的列，然后计算传递函数，就可以了。

但是对于幅度和相位，我很难绘制它们。我很累将它们绘制出来，但这是完全错误的。有人可以帮我吗？

这是代码和excel文件。

solarcell1 = xlsread('solarcell.xlsx','A2:C100005');
t=solarcell1(:,1);
N=length(t);
Channel1V = solarcell1(:,2);
Channel2V = solarcell1(:,3);
sig1=Channel1V;
sig2=Channel2V;
fs=1/((solarcell1(3)-solarcell1(2))*1e-3);
FA1=fs/length(sig1);
FA2=fs/length(sig2);
frange1=-fs/2:FA1:fs/2-FA1;
frange2=-fs/2:FA2:fs/2-FA2;
subplot(3,2,1);
plot(t,sig1);
hold on
plot(t,sig2);
title('Input and Output of Solar cell');
xlabel('Time');ylabel('Amplitude');
subplot(3,2,2);
plot(t,sig2);
title('Output');
xlabel('Time');ylabel('Amplitude');
z1=fftshift(fft(sig1))/length(t);
subplot(3,2,3);
plot(frange1,abs(z1));
title('Input');
xlabel('Freq.');ylabel('Amplitude');
z2=fftshift(fft(sig2))/length(t);
subplot(3,2,4);
plot(frange2,abs(z2));
title('Output');
xlabel('Freq.');ylabel('Amplitude');
TFC=z2./z1;
magnitude=20*log(abs(TFC));
phase=atan2(imag(TFC),real(TFC));
subplot(3,2,5);
bode(frange1(1002:N),magnitude(1002:N));
title('Magnitude');
xlabel('Freq.');ylabel('Magnitude');
subplot(3,2,6);
semilogx(frange1(1002:N),phase(1002:N));
title('Phase');
xlabel('Freq.');ylabel('Phase');

Answer 1

在您提供的代码中，我可以看到三个问题。

首先，您对bode的使用不正确。此函数以dynamic model system作为参数，并且给了两个向量frange1(1002:N)和magnitude(1002:N)。现在，我建议您只使用plot。

subplot(3,2,5);
plot(frange1(1002:N),magnitude(1002:N));
title('Magnitude');
xlabel('Freq.');ylabel('Magnitude');

然后，当x轴的值为负值时，semilogx的使用也会带来风险。使用frange1=-fs/2:FA1:fs/2-FA1;时，您无需索取semilogx。现在，我建议使用plot。

最后，我建议您使用unwrap绘制相位。

subplot(3,2,6);
plot(frange1(1002:N),unwrap(phase(1002:N)));
title('Phase');
xlabel('Freq.');ylabel('Phase');

它绘制：