fsresp.m

信号实验常用的simulink模型和mfile,可直接在matlan下运行。

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.
if nargin<4,wo=1;end
N=(length(un)-1)/2;		% highest harmonic
jnWo=sqrt(-1)*(-N:N)*wo;	% frequencies of harmonics
y=(polyval(num,jnWo)./polyval(den,jnWo)).*un;	% output is Yn=H(jnWo)*Un