dampsinf.m

来自「这是一个 York大学开发的FDTD代码。很好用的工具」· M 代码 · 共 18 行

function v = dampsinf(w,wn,z)
%Generate Damped sine wave spectrum
%dampsin (w,wn,z)
%  w = angular frequency
% wn = Natural frequency
%  z = damping factor

% this gives unit dc component

    v=(wn.^2)' * ones(1,size(w,2))./(ones(size(wn,2),1)*(-(w.^2)) + ((2*z.*wn*j)'*w) + (wn.^2)'*ones(1,size(w,2)));

% this gives unit amplitude sine wave but doesn't like zero z
    v=(abs(2*z.*sqrt(1-z.^2)))'*ones(1,size(w,2)).*v;



end