matlab实时脚本算拉普拉斯反变换和画图

//////////////////////

clc close all clear syms s t num = 80; den = s^4+12*s^3+41*s^2+82*s+80; Fs = num/den res = solve(den) %求解分母多项式等于 0 的根（即特征根） Ft = simplify(ilaplace(Fs, s, t)) % %对 F(s)做拉普拉斯逆变换，将变量从拉普拉斯域的s转换为时域的t，得到时域函数f(t) % ts = 0:0.01:10; Ft = subs(Ft, t, ts); %将符号函数 Ft 中的变量 t 替换为数值向量 ts， % 把符号表达式转换为数值数组 % （MATLAB 的plot函数仅支持数值，无法直接绘制符号函数） plot(ts, Ft, 'LineWidth', 2); %绘制时域函数 f(t) 的曲线： % 横坐标为ts，纵坐标为Ft，线宽设为 2（线条更粗，便于观察）。 hold on; plot(ts, sqrt(2)*exp(-ts), '-.r', ts, -sqrt(2)*exp(-ts), '-.r', 'LineWidth', 1) %绘制两条红色点划线（'-.r'）作为 f(t) 的包络线： % - 上包络线：2​e−t；- 下包络线：−2​e−t；线宽设为 1，与主曲线区分开。 grid on;

///////////// matlab求零极点

clear all %清空工作区 clc format long num=[0 0 0 0 80] %分子 den=[1 12 41 82 80] %分母 G=tf(num, den); %将分子和分母组成传递函数 bode(G, {10,10e6}); %规定bode图横轴弧度/秒的上下限 legend('G') %给当前图形添加图例 grid zpk(G) %将传递函数从「分子 / 分母多项式形式」转换为「零极点增益形式（Zero-Pole-Gain, ZPK）」 %并在命令行窗口输出系统的零点（Zero）、极点（Pole） 和增益（Gain）

////////////////////////////////////////////////

************************************

clear all %清空工作区 format long Vg=10; C=3.3e-6; L=50e-6; R=1; T1=5.2925e-5; K=0.01; % K=0.01*1000; num=[0 Vg*C*R Vg] den=[L*C*R L R] num2=[0 K] den2=[T1 1] G=tf(num, den); Iden=tf(num2, den2); bode(G, Iden, {100,10e6}); legend('G','Iden') grid zpk(G) zpk(Iden) /////////////////// clear all %清空工作区 format long num=[0 4.688e-9 3.256e-4 5.652] %分子 den=[6.284e-11 3.29e-7 0.01106 0.65] %分母 G=tf(num, den); %将分子和分母组成传递函数 bode(G, {100,10e6}); %规定bode图横轴弧度/秒的上下限 legend('G') %给当前图形添加图例 grid zpk(G) %将传递函数从「分子 / 分母多项式形式」转换为「零极点增益形式（Zero-Pole-Gain, ZPK）」 %并在命令行窗口输出系统的零点（Zero）、极点（Pole） 和增益（Gain） //////////////// clear all %清空工作区 clc format long num=[0 7.032e-9 4.884e-4 8.478] %分子 den=[7.5408e-11 3.948e-7 0.013272 0.78] %分母 G=tf(num, den); %将分子和分母组成传递函数 bode(G, {10,10e6}); %规定bode图横轴弧度/秒的上下限 legend('G') %给当前图形添加图例 grid zpk(G) %将传递函数从「分子 / 分母多项式形式」转换为「零极点增益形式（Zero-Pole-Gain, ZPK）」 %并在命令行窗口输出系统的零点（Zero）、极点（Pole） 和增益（Gain） //////////&&&&&&&&&&&&&&&&&&&& clc close all clear syms s t num = 80; den = s^4+12*s^3+41*s^2+82*s+80; Fs = num/den res = solve(den) %求解分母多项式等于 0 的根（即特征根） Ft = simplify(ilaplace(Fs, s, t)) % %对 F(s)做拉普拉斯逆变换，将变量从拉普拉斯域的s转换为时域的t，得到时域函数f(t) % ts = 0:0.01:10; Ft = subs(Ft, t, ts); %将符号函数 Ft 中的变量 t 替换为数值向量 ts， % 把符号表达式转换为数值数组 % （MATLAB 的plot函数仅支持数值，无法直接绘制符号函数） plot(ts, Ft, 'LineWidth', 2); %绘制时域函数 f(t) 的曲线： % 横坐标为ts，纵坐标为Ft，线宽设为 2（线条更粗，便于观察）。 hold on; plot(ts, sqrt(2)*exp(-ts), '-.r', ts, -sqrt(2)*exp(-ts), '-.r', 'LineWidth', 1) %绘制两条红色点划线（'-.r'）作为 f(t) 的包络线： % - 上包络线：2​e−t；- 下包络线：−2​e−t；线宽设为 1，与主曲线区分开。 grid on; &&&&&&&&&&&&*************************