www.gusucode.com > Matlab源程序 > Matlab源程序/精通Matlab综合辅导与指南-源程序/fsresp.m

    function y=fsresp(num,den,un,wo)
%FSRESP Fourier Series Linear System Response.
% FSRESP(N,D,Un,Wo) returns the complex exponential FS of the
% output of a linear system when the input is given by a FS.
% N and D are the numerator and denominator coefficients
% respectively of the system transfer function.
% Un is the complex exponential Fourier Series of the system input.
% Wo is the fundamental frequency associated with the input.

% D.C. Hanselman, University of Maine, Orono, ME 04469
% 1/12/95
% Copyright (c) 1996 by Prentice-Hall, Inc.

if nargin<4,wo=1;end
N=(length(un)-1)/2;
jNwo=sqrt(-1)*(-N:N)*wo;
y=(polyval(num,jNwo)./polyval(den,jNwo)).*un;