exa6_2.m

来自「MATLAB 的一些编程,基本的程序,为出学MATLAB的初学者」· M 代码 · 共 12 行

function y=exa6_2(a,y0)
n=length(a)-1;
p=roots(a);
v=rot90(vander(p));
c=inv(v)*y0;
t=0:0.1:10;y=zeros(size(t));
for ii=1:n
    y=c(ii)*exp(p(ii)*t)+y;
end
plot(t,y);grid
s=['y0=[',num2str(y0'),']']; 
legend(s);