ex3_10.m

来自「源代码集合」· M 代码 · 共 15 行

A=[0 1; -1 -2]; B=[0;1]; C=[1 2]; D=0;

w=[0.1:0.05:1, 1.5:0.5:10, 15:5:100 150:50:1000]; 
m=[]; I=eye(size(A));
for i=w, m=[m; C*inv(sqrt(-1)*i*I-A)*B+D]; end
length(w)

w1=0.1:.1:1000; m1=[];
for i=w1, m1=[m1; C*inv(sqrt(-1)*i*I-A)*B+D]; end
length(m1)

w2=0.1:1:1000; m2=[];
for i=w2, m2=[m2; C*inv(sqrt(-1)*i*I-A)*B+D]; end
length(m2)