exp4_16.m

Matlab复杂分析模拟讲义及例程

% exp4_16.m
clear
close all
clc
%有理传递函数模型
num=[1 1];
den=conv([1 2],conv([1 2],[1 2]));
w=logspace(-1,2);
t=0.5;
%求有理传递函数模型的频率响应
[mag1,pha1]=bode(num,den,w);
%求系统的等效传递函数
[n2,d2]=pade(t,4);
numt=conv(n2,num);
dent=conv(d2,den);
%求系统的频率响应
[mag2,pha2]=bode(numt,dent,w);
%在同一张图上绘制频率响应曲线
subplot(211)
semilogx(w,20*log10(mag1),w,20*log10(mag2),'r--');
title('bode plot')
xlabel('frequency-rad/s');
ylabel('gain db');
grid on
subplot(212)
semilogx(w,pha1,w,pha2,'r--');
xlabel('frequency-rad/s');
ylabel('phase deg');
grid on