Message Boards Message Boards

0
|
5929 Views
|
2 Replies
|
0 Total Likes
View groups...
Share
Share this post:
GROUPS:

How to work out the DC bias of a transfer function to a sine wave input

Posted 4 years ago

For a transfer function of (0.5s +0.1)/s, how can I predict what the response DC bias to a sine wave input may be? I can find this using MATLAB/simulink but is there a way I can predict this theoretically?

Thanks

POSTED BY: Henry Cole
2 Replies
Anonymous User
Anonymous User
Posted 4 years ago

It would be impossible to answer the question without more information - you are assuming a lot.

POSTED BY: Anonymous User
Posted 4 years ago

Perhaps this:

(* transfer function *)
tf = (1/2 s + 1/10)/s // Simplify;

(* the input *)
in = LaplaceTransform[Sin[w t], t, s];

(* output in t domain *)
outt = InverseLaplaceTransform[tf in, s, t];

(* the integral of the output in the t domain *)
(* note division by s to get integral *)
outtIntegral = InverseLaplaceTransform[tf in/s, s, t];

(* the mean value over the first n periods of the input *)
Simplify[outtIntegral/n /. t -> n 2 Pi/w, n \[Element] Integers]

(* \[Pi]/(5 w^2) *)

Edit: Sorry, the integral was over t = 0 to n 2 Pi/w, so that is the normalization period. (the mean value over the first n periods of the input)

Simplify[outtIntegral/(n 2 Pi/w) /. t -> n 2 Pi/w, 
 n \[Element] Integers]

 (* 1/(10 w) *)

Edit: A numerical example with w = 1

dcBias = NIntegrate[outt /. w -> 1, {t, 0, 6 Pi}]/6/Pi

(* 0.100 *)

Plot[{outt /. w -> 1, Callout[dcBias, "DC Bias"]}, {t, 0, 6 Pi}]

enter image description here

POSTED BY: David Keith
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract