SYSM 6325
1 0 ) UT % vector fresp contains ( complex-valued) frequency responses % vector w, contains frequencies at which measurements were taken load...
View the step-by-step solution to:

Question

There is one more matlab file to solve this question. Let me know if you can solve this problem as soon as

possible.

Screen Shot 2019-10-14 at 5.33.09 PM.pngex_1_data.m.pngex_1_threemodes.m.png

ex_1_data.m.png

1 0 ) UT
% vector fresp contains ( complex-valued) frequency responses
% vector w, contains frequencies at which measurements were taken
load Kalman_G11QA3. mat
10
-
fresp=02i1(10:10: 650) ;
11
12
frequency vector (in Hz)
13
-
freq=02i1x (10 : 10 :650) ;
14
15
%frequency vector (in rad/sec)
16
-
w = 2*pi*freq;
17
18
%non-parametric model
19
-
G_h = frd( fresp, w) ;
20
-
bode (G_h)
21

ex_1_threemodes.m.png

VOULAWN
load Kalman_G11QA3. mat
fresp=02i1 (20:1:1200) ; %frequency response vector
% frequency vector (in Hz)
freq=02i1x(20:1:1200) ;
8
9
frequency vector (in rad/sec)
10
w = 2*pi*freq;
11
12
%non-parametric model
13
G_h = frd( fresp, w);
% SYS = frd (RESPONSE, FREQS)
14
bode (G_h)
% SYS = frd (RESPONSE, FREQS, TS) creates a discrete-time frd model with
15
sample time TS (a positive value)
16
17
% Frequency Response Data ( frd) models store the frequency response of linear
18
% systems, for example, experimental data collected with a frequency analyzer

Screen Shot 2019-10-14 at 5.33.09 PM.png

Consider the problem of identifying transfer function of a linear time—invariant single—
input single—output system from frequency domain data. Let us assume that the transfer function of the system can be described by transfer function B( )
s
G“) ‘ A(s)
where
B(s) = be + bls + b282 + . .. + bn_1s"_1
A(s) = l+a13+a232+...+ansn We perform experiments that involve applying sine waves to the system at frequencies
£01,002, . . . ,wN and determining the input—output frequency response at those frequencies, A C(jwi). Note that we use G instead of G, as there will be unavoidable errors in the measurements. (i) Show that the problem of finding A(s) and 3(3) can be cast as a least squares
optimization problem. Define a least squares cost function and find the optimal solution. Hint: Note that
AUCUWUW) = BUM), and that the error corresponding to each frequency wi may be defined as
6w: = AUWQGUM) — BUM)
Also note that 6,; is complex-valued. (ii) Download the data file G11QA3.mat and the accompanying Matlab file ex—l—data.m.
Running ex— 1—data.m produces the frequency response function of the first mode of
a cantilever beam. Use the method above to identify a model for the system with
this frequency response data. First, you will need to make an estimate of the order of this system. Explain your rationale. (iii) Evaluate and plot the least squares cost associated with increasing order of the system and use that to determine if your initial estimate was correct. (iv) Can you obtain a better fit to the frequency response data by incorporating a weight— ing function into the least squares cost? Explain how. (v) Download and run ex—l—threemodes.m. This data set contains the frequency response
of the first three modes of the same structure. Repeat steps (ii)—(iv) with this data.
Can you achieve an acceptable outcome by selecting a subset of frequency data points? Explain your observations.

Recently Asked Questions

Why Join Course Hero?

Course Hero has all the homework and study help you need to succeed! We’ve got course-specific notes, study guides, and practice tests along with expert tutors.

  • -

    Study Documents

    Find the best study resources around, tagged to your specific courses. Share your own to gain free Course Hero access.

    Browse Documents
  • -

    Question & Answers

    Get one-on-one homework help from our expert tutors—available online 24/7. Ask your own questions or browse existing Q&A threads. Satisfaction guaranteed!

    Ask a Question
Ask Expert Tutors You can ask You can ask ( soon) You can ask (will expire )
Answers in as fast as 15 minutes