lapulas2.m

来自「the code of the book come form the book 」· M 代码 · 共 17 行

function lapulas2(b,a,t0)
omega=-5:0.073:6;
sigma=0.05:0.073:1.1;

%omega=0:0.16:1.5;
%sigma=0.05:0.09:10;
[sigmagrid,omegagride]=meshgrid(sigma,omega);
sgrid=sigmagrid+j*omegagride;
L2=polyval(b,sgrid)./polyval(a,sgrid);
L2=L2.*exp(-t0*sgrid);
mesh(sigma,omega,angle(L2));
set(gca,'color',[1 1 1]);
set(gca,'XColor',[0 0 0]);                                                                                          
set(gca,'YColor',[0 0 0]);                                                                                          
set(gca,'ZColor',[0 0 0]);
axis([0,1.1,-5,6,-4,4]);