pronyfreq.m

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

function [F]=pronyfreq(C,Sp,fr,Ts)
% [F]=pronyfreq(C,Sp,fr,Ts)
% Use the method of Prony representation of a time response to determine
% the frequency repsonse or spectrum.
% In:
%  C vector of amplitude coefficients from prony exponential series
%  Sp vector of poles
%  fr=frequency value(s)
%  Ts=time domain sample rate?
% Out:
%  F=complex frequency response
%  J. F. Dawson 27 Oct 1999

%How many poles
Np=rows(Sp)

%How many frequencies
Nf=rows(fr)

% compute jw
jw=2j*pi.*fr;   %column vector

% Compute jw-Sp
W=jw*ones(1,Np);         % Nfx1 * 1xNp => NfxNp
S=ones(Nf,1)*Sp.';       % Nfx1 * 1xNp => NfxNp

%Coefficient matrix
D=ones(Nf,1)*C.';        % Nfx1 * 1xNp => NfxNp

%Do coefficient multiply and summation 
F=D./(W-S); % coefficient/Denominator        % NfxNp  
F=F*ones(Np,1); % Summation                  % NfxNp * Npx1 =>