%
%    sdof_base_accel_transfer.m  ver 1.0  November 7, 2012
%
function[f,H_acc,H_rd]=sdof_base_accel_transfer(iunit,df,num,fn,damp)
%
omegan = 2*pi*fn;
%
f=linspace(0,(num-1)*df,num);
%
omega=2*pi*f;
%
H_acc=zeros(num,1);
H_rd=zeros(num,1);
%
for i=1:num
%
    jterm=(1i)*2*damp*omega(i)*omegan;
%
    den=((omegan^2-omega(i)^2)+jterm);
%
    H_acc(i)=(omegan^2 +jterm )/den; 
    H_rd(i)= -1/den; 
%
end
%
if(iunit==1)
    H_rd=H_rd*386;
end
%
f=fix_size(f);
H_rd=fix_size(H_rd);
H_acc=fix_size(H_acc);