Message Boards

WOLFRAM COMMUNITY

1427 Views

4 Replies

1 Total Likes

View groups...

Follow this post

Share this post:

GROUPS:

How can I plot a system's response depending on certain condition

Koushik Maharatna

Posted 11 months ago

Hi All, I am working on a system that would generate output with a particular type of square wave as input as long as an input condition (coming from a subsystem) is above a certain threshold. I tried using If statement but it seems not working. Could you please help? Here is my code. Effectively, what I am looking for is, when the output response ('coupling') of the system ('TF1') with the square wave defined as 'inwave' as the input is >= 0.015, the system output ('v1') should be the 'inwave' itself. It Otherwise it should be 0. It is working fine up to the plotting the 'coupling' signal. But the problem is with the If statement. Thanks for your kind help. TF1 = TransferFunctionModel[ Unevaluated[{{(697052. + 25694.6 s)/( 700128. + 25835.5 s + 1. s^2)}}], s, SamplingPeriod ->None, SystemsModelLabels -> None]; PulseWave[duty_, t_, period_, phase_] := UnitBox[Mod[(t/period) + phase, 1]/(2 duty)]; inwave = 11010^-3 PulseWave[0.25, t, 10^-2, 0] - 7010^-3; coupling = OutputResponse[TF1, inwave, {t, 0, 10 10^-2}]; Plot[coupling, {t, 0, 10 10^-2}] v1 = If[coupling >= 0.015, inwave, 0]; Plot[v1, {t, 0, 10 10^-2}, PlotRange -> All]

POSTED BY: Koushik Maharatna

4 Replies

Sort By:

Gianluca Gorni

Gianluca Gorni, University of Udine

Posted 11 months ago

Sorry, I don't know the meaning of your plots, but here a reworking that runs: Clear[TF1, inwave, coupling, v1]; PulseWave[duty_, t_, period_, phase_] := UnitBox[Mod[(t/period) + phase, 1]/(2 duty)]; TF1 = TransferFunctionModel[{{(697052. + 25694.6 s)/(700128. + 25835.5 s + 1. s^2)}}, s, SamplingPeriod -> None, SystemsModelLabels -> None]; inwave[t_] = OutputResponse[TF1, 40*10^-3 PulseWave[0.25, t, 10^-3, 0], {t, 0, 4 10^-3}][[1]]; Plot[inwave[t], {t, 0, 4 10^-3}, PlotLegends -> {"inwave"}] coupling[t_] = OutputResponse[TF1, inwave[t], {t, 0, 4 10^-3}][[1]]; Plot[{coupling[t], inwave[t]}, {t, 0, 4 10^-3}, PlotRange -> All, GridLines -> {{}, {.015}}, PlotLegends -> {"coupling", "inwave"}] v1[t_] = If[coupling[t] >= 0.02, inwave[t], 0]; Plot[v1[t], {t, 0, 4 10^-3}, PlotRange -> {-.01, .045}, GridLines -> {{}, {.015}}, PlotLegends -> {"v1"}]

Sorry, I don't know the meaning of your plots, but here a reworking that runs:

Clear[TF1, inwave, coupling, v1];
PulseWave[duty_, t_, period_, phase_] := 
  UnitBox[Mod[(t/period) + phase, 1]/(2 duty)];
TF1 = TransferFunctionModel[{{(697052. + 25694.6 s)/(700128. + 
        25835.5 s + 1. s^2)}}, s, SamplingPeriod -> None, 
   SystemsModelLabels -> None];
inwave[t_] = 
  OutputResponse[TF1, 
    40*10^-3 PulseWave[0.25, t, 10^-3, 0], {t, 0, 4 10^-3}][[1]];
Plot[inwave[t], {t, 0, 4 10^-3}, PlotLegends -> {"inwave"}]
coupling[t_] = OutputResponse[TF1, inwave[t], {t, 0, 4 10^-3}][[1]];
Plot[{coupling[t], inwave[t]}, {t, 0, 4 10^-3}, PlotRange -> All, 
 GridLines -> {{}, {.015}}, PlotLegends -> {"coupling", "inwave"}]
v1[t_] = If[coupling[t] >= 0.02, inwave[t], 0];
Plot[v1[t], {t, 0, 4 10^-3}, PlotRange -> {-.01, .045}, 
 GridLines -> {{}, {.015}}, PlotLegends -> {"v1"}]

POSTED BY: Gianluca Gorni

Koushik Maharatna

Posted 11 months ago

Thank you so very much! It works perfectly. I am trying to model a biological system of action potential propagation and this is related to that work. Thanks a lot for your kind help.

POSTED BY: Koushik Maharatna

Koushik Maharatna

Posted 11 months ago

Dear Gianluca, Thank you so much. As such it worked in a reverse way, i.e., I am getting the response only when I set coupling <= 0.015, although I can see that coupling for certain times is >=0.015. Not sure why this happens. For coupling >=0.015 it always gives flat line at 0. Any help will be much appreciated. Here is the code again with your suggested modification. PulseWave[duty_, t_, period_, phase_] := UnitBox[Mod[(t/period) + phase, 1]/(2 duty)]; TF1 = TransferFunctionModel[ Unevaluated[{{(697052. + 25694.6 s)/( 700128. + 25835.5 s + 1. s^2)}}], s, SamplingPeriod ->None, SystemsModelLabels -> None]; inwave = OutputResponse[TF1, 40*10^-3 PulseWave[0.25, t, 10^-3, 0], {t, 0, 4 10^-3}]; Plot[inwave, {t, 0, 4 10^-3}]; coupling = OutputResponse[TF1, inwave, {t, 0, 4 10^-3}]; Plot[coupling, {t, 0, 4 10^-3}, PlotRange -> All] v1 = If[coupling[[1]][[1]] >= 0.015, inwave, 0]; Plot[v1, {t, 0, 4 10^-3}, PlotRange -> All]

Dear Gianluca, Thank you so much. As such it worked in a reverse way, i.e., I am getting the response only when I set coupling <= 0.015, although I can see that coupling for certain times is >=0.015. Not sure why this happens. For coupling >=0.015 it always gives flat line at 0. Any help will be much appreciated. Here is the code again with your suggested modification.

PulseWave[duty_, t_, period_, phase_] := 
  UnitBox[Mod[(t/period) + phase, 1]/(2 duty)];
TF1 = TransferFunctionModel[
   Unevaluated[{{(697052. + 25694.6 s)/(
     700128. + 25835.5 s + 1. s^2)}}], s, SamplingPeriod ->None, 
   SystemsModelLabels -> None];
inwave = 
  OutputResponse[TF1, 
   40*10^-3 PulseWave[0.25, t, 10^-3, 0], {t, 0, 4 10^-3}];
Plot[inwave, {t, 0, 4 10^-3}];
coupling = OutputResponse[TF1, inwave, {t, 0, 4 10^-3}];
Plot[coupling, {t, 0, 4 10^-3}, PlotRange -> All]
v1 = If[coupling[[1]][[1]] >= 0.015, inwave, 0];
Plot[v1, {t, 0, 4 10^-3}, PlotRange -> All]

POSTED BY: Koushik Maharatna

Gianluca Gorni

Gianluca Gorni, University of Udine

Posted 11 months ago

Your `coupling` is a list, you cannot compare it to a number. Try this: v1 = If[coupling[[1]] >= 0.015, inwave, 0]

POSTED BY: Gianluca Gorni

Reply to this discussion

Reply Preview

Attachments

Remove Add a file to this post

Follow this discussion

or Discard

Group Abstract

Feedback