在MATLAB中如何将时域方程进行拉普拉斯变换以后，直接绘制BODE图以及NYQUIST曲线

在MATLAB中如何将时域方程进行拉普拉斯变换以后，直接绘制BODE图以及NYQUIST曲线

我们以一个一阶的串联RC低通动态电路的全响应方程为例

syms t s;
Us=15;
T=100;
U0=5;
ft=(Us*(1-exp(-t/T))+U0*exp(-t/T))/Us;
    Fs=laplace(ft,t,s);
    [h,k] = numden(Fs);
    a=sym2poly(h);
    b=sym2poly(k);
    G=tf(a,b);
    subplot(1,2,1) ;bode(G);
    subplot(1,2,2) ;nyquist(G);

我们在其中设置电压源为15V，电容初始电压为5V，时间常数为100.根据电路知识我们得到响应方程为代码中的“ft”。当然这些并不重要，这些可以自由设置。
最重要的是以下部分：
Fs为ft拉氏变换以后的方程
但是Fs不能直接进行BODE图的绘制，因为此时的Fs是最简型，并没有将其有理化，matlab的bode函数无法将其识别，所以我们还需要进行以下处理。

1. 利用numden函数，我们可以得到h与k分别为将Fs进行有理化以后的分子与分母。
2. 利用sym2poly函数，我们将符号多项式的系数依次由高阶到0阶的顺序输出。所以得到了分子h的系数矩阵a，分母k的系数矩阵b。
3. 最后利用tf函数，将分子分母的系数重新组合成一个频域传递函数。

这就是将拉拉氏变换以后的传递函数进行有理化的过程，简而言之这就是将一个分式函数通过分别将分子分母有理化，然后将其系数提取出来，最后又进行组合的过程。
感谢您的观看，希望能够对诸君有所帮助。最后祝阅读者学业有成，工作顺利（^ _ ^）.