s10fig22.m

来自「Industrial Mathematics」· M 代码 · 共 25 行

% Figure 10.22
%
clear
c=0;
d=0;
dt = 0.1;
dz = 0.1;
for m=1:101         % for t
  for n = 1:9           % for z
    z  = n*dz;
    sz = sqrt(1-z*z);
    t = (m-1)*dt;
    phi = atan(z/sz);
    a(m) = t;
    b(m,n) = 1 - (1/sz)*exp(-z*t)*cos(sz*t + phi);
     end;
  end;
for m=1:101         % for t
  for n = 1:9           % for z
   c = [c; a(m)];
   d= [d;b(m,n)];
    end; end;
%plot(c,d,'k.')
plot(b,'k-')