直接上代码:
%展示vout的Bode图,返回tf函数的返回值,即传递函数模型。
%其中vin为变量,s为拉普拉斯算子;vout为vin和s的函数。
function [Hs_Bode] = showMeBode(vout,vin,s)
[numexpr, denexpr] = numden(vout);
num = fliplr(double(coeffs(expand(numexpr), [s,vin])));
den = fliplr(double(coeffs(expand(denexpr), [s,vin])));
Hs_Bode=tf(num,den)
margin(Hs_Bode)
end
使用实例:
clc;clear;
syms vin s;
vout = vin/(s+1);
figure(1)
Hs = showMeBode(vout,vin,s);
hold on;
输出:
Hs_Bode = 1 ----- s + 1 连续时间传递函数。 模型属性