Question

H(s) where, Q-5, and 40-2000, rad./sec. a) Use MATLAB to determine magnitude response and phase response of the filter. b) What type of filter is it? c) What will be the output of this filter if input x,() SCos(1001). Show all calculations step by step as shown in Lecture-20 d) Verify your answer of part (c) by using Simulink model. Attach the snapshot of Simulink model and output. e) What will be the output of this filter if input x(0) SCos(2000t) Show all calculations f) Verify your answer of part (e) by using Simulink model. Attach the snapshot of Simulink model and output. g) what will be the output of this filter if input X3(t) = 5Cos(4000xt). Show all calculations. h) Verify your answer of part (g) by using Simulink model. Attach the snapshot of Simulink model and output. i) Using MATLAB find the poles of H(s). Is the system stable? j) Find and sketch h(t) in MATLAB

All i want is the whole MATLAB code for this nothing else

Answer 1

MATLAB code

clc;

clear all;

syms t l s;

format long

Q=5

w0=2000*pi

a=w0/Q

b=w0*w0

num=[a 0]

den=[1 a b]

H=tf(num,den)

bode(H)

pzmap(H)%pole zero maping

%transfer function in s domain

h=(1257*s)/((s*s)+1257*s+3.9488e07)

%inverse laplace to get impulse responce h(t)

h1=ilaplace(h)

MATLAB results

LIVE EDITOR MEW Run Section Control. %33 Go To BIU M Run and Advance ロゼ」ES Section Run to End Find ▼ Text Refactor ▼ | SECTION RUN NAVIGATE TEXT CODE BodePlotOfDynamicSystemExample.mb clc; clear all: syms t 1 s; format long 4 6.283185307179586e+03 1.256637061435917e+03 7 8 we-2000" pi 3.947841760435743e+07 a-we/ b-we we num1 2 10 11 num-a e1 12 den-[1 a b] 1.256637061435917 13 14 15 16 17 18 19 28 H-tf (num, den) bode (H) pzmapCH)%pole zero maping 0.000000100000000 0.000125663706144 3.947841760435743 transfer function in s domain h- (1257 s)/((s s)+1257 s+3.9488e07) %inverse laplace to get impulse responce h(t) hi ilaplace(h) 1257 s s21257 3.948e07 Continuous-time transfer function.

LIVE EDITOR MEW Run Section Run and Advance Run to Normal Control ▼ % GoTo ▼ BIU M Text Find ▼ Refactor ▼ Run Step Stop NAVIGATE TEXT CODE SECTION RUN BodePlotOfDynamicSystemExample.mbx Bode Diagram 20 30 90 a -45 90 10 10 Frequency (rad/s) 10

LIVE EDITOR INSERT MEW Normal ▼ Run Section Run and Advance Run to GoTo ▼ control ▼ B I U M Texd .. Run Step Stop Find Refactor Section Break NAVIGATE TEXT CODE SECTION RUN BodePlotOfDynamicSystemExample.mlx Pole-Zero Map 8000 6000 4000 8 2000 2000 E -4000 6000 8000 500 100 -700 600 400300 200 0 Real Axis (seconds

FIGURE INSERT MEW Legend Remove L... Colorbar Remove Co...Grid Remove Grid X-Grid X-label Y-Label Y-Grid Text Arrow ANNOTATIONS BodePlotOfDynamicSystemExample.mlbx Pole-Zero Map 8000 6000 X-628.3 Y 6252 4000 8 2000 2000 E -4000 6000 8000 -700 500 400 -300 100 600 200 0 Real Axis (seconds

LIVE EDITOR INSERT MEW 呂Run Section Normal > Control% B IU M Run and Advance Run Step Find ▼ Text Refacfor Break Run to NAVIGATE TEXT CODE SECTION RUN BodePlotOfDynamicSystemExample mk 2000 -4000 -6000 8000 700 600 500 -400 300 200 -100 Real Axis (seconds') 1257 s s2+1257 s +39488000 257156371951 sin 156371951 1257 2 1257e2cos 156371951 156371951

--------------------------------------------------------------------------------------------------------------------

B)answer

from the bode plot magnitude response we can conclude that the given transfer function is a band pass filter

I) answer

from the pole zero map we can derive the 2 poles at

s=-628.3+j6252

s=-628.3-j6252

and one ZERO at

s=0.

answers to C,E,G

,并15)^ 12575 1253 5 (S+ 6281-25aj 1257 Ju we nouw that, od put 125子W 6233 .21 x1o Cos ( Ioo Cos ( 4ooo

Thanks and regards