This is machine translation

Translated by Microsoft
Mouseover text to see original. Click the button below to return to the English version of the page.

Note: This page has been translated by MathWorks. Click here to see
To view all translated materials including this page, select Country from the country navigator on the bottom of this page.

timeresp

Time response for rational function object

Syntax

[y,t] = timeresp(h,u,ts)

Description

[y,t] = timeresp(h,u,ts) computes the output signal, y, that the rational function object, h, produces in response to the given input signal, u.

The input h is the handle of a rational function object returned by rationalfit. ts is a positive scalar value that specifies the sample time of the input signal.

The output y is the output signal. RF Toolbox™ software computes the value of the signal at the time samples in the vector t using the following equation.

Y(n)=sum(C.*X(nDelay/ts))+D*U(nDelay/ts)

where

X(n+1)=F*X(n)+G*U(n)X(1)=0F=exp(A*ts)G=(F1)./A

and A, C, D, and Delay are properties of the rational function object, h.

Examples

collapse all

Define the input signal.

SampleTime = 2e-11;
OverSamplingFactor = 25;
TotalSampleNumber = 2^12;
InputTime = double((1:TotalSampleNumber)')*SampleTime;
InputSignal = sign(randn(1, ...
              ceil(TotalSampleNumber/OverSamplingFactor)));
InputSignal = repmat(InputSignal, [OverSamplingFactor, 1]);
InputSignal = InputSignal(:);

Create a rational function object.

orig_data=read(rfdata.data,'default.s2p');
freq=orig_data.Freq;
data=orig_data.S_Parameters(2,1,:);
fit_data=rationalfit(freq,data);
Warning: Achieved only -32.8 dB accuracy with 48 poles, not -40.0 dB. 

Compute the time response.

[y,t]=timeresp(fit_data,InputSignal,SampleTime);

Introduced in R2007a