home *** CD-ROM | disk | FTP | other *** search
/ Borland Programmer's Resource / Borland_Programmers_Resource_CD_1995.iso / utils / mathstud / sinw.m < prev    next >
Encoding:
Text File  |  1995-05-19  |  1.4 KB  |  63 lines
  1. function y=sinw(npoints, f0, df, ph0)
  2. % y=sinw(num_of_points, num_of_cycles, delta_cycles, initial_phase)
  3. % create frequency hopping sine wave composed of segements of linear
  4. % FM sine waves
  5. %
  6. % num_of_points - number of points in each segment
  7. % f0            - inital frequency
  8. % df            - frequency deviation
  9. % ph0           - initial phase
  10. % 
  11. % Example: sinw(256, 2, 3) generates a linear FM signal
  12. % 
  13.  
  14. %       S.Halevy 7/31/92
  15. %       Copyright (c) 1992 by the MathWizards
  16.  
  17. if nargin < 1
  18.   error('insufficient arguments')
  19. end
  20. npoints=npoints(:);   % make sure it is a col vector
  21. [mp,np]=size(npoints);
  22.  
  23. if nargin < 4
  24.    ph0=zeros(mp,np);  % zero initial phase - radians
  25. else
  26.    ph0=ph0(:);
  27.    [m0,n0]=size(ph0);
  28.    if m0 ~= mp | n0 ~= np
  29.       error('initial phase - dimension error')
  30.    end
  31. end
  32.  
  33. if nargin < 3
  34.    df=zeros(mp,np);   % no frequency deviation
  35. else
  36.    df=df(:);
  37.    [m0,n0]=size(df);
  38.    if m0 ~= mp | n0 ~= np
  39.       error('frequency deviation - dimension error')
  40.    end
  41. end
  42.  
  43. if nargin < 2
  44.    f0=ones(mp,np);    % single cycle
  45. else
  46.    f0=f0(:);
  47.    [m0,n0]=size(f0);
  48.    if m0 ~= mp | n0 ~= np
  49.       error('cycles - dimension error')
  50.    end
  51. end
  52.  
  53. y=[];
  54. pi2=atan(1.0)*8;
  55. for i1=1:mp
  56.    n = npoints(i1);
  57.    t = (0:(n-1))/n;                  % time
  58.    ph1 = (f0(i1)+0.5*df(i1)*t).*t;   % phase
  59.    ph1 -= floor(ph1);                % reduce roundoff errors
  60.  
  61.    y=[y sin(pi2*ph1+ph0(i1))];
  62. end
  63.