lapulas1.m

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

function lapulas1(b,a,t0)
omega=-5:0.073:6;
sigma=0.05:0.073:1.1;
[sigmagrid,omegagride]=meshgrid(sigma,omega);
sgrid=sigmagrid+j*omegagride;
L1=polyval(b,sgrid)./polyval(a,sgrid);
L1=L1.*exp(-t0.*sgrid);
mesh(sigma,omega,abs(L1));
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,0,10]);