Vibrating Beam Inverse Problem梁振动反问题.pptVIP

Vibrating Beam Inverse Problem Team K.E.Y Scott Clark ● Asya Monds ● Hanh Pham SAMSI Undergraduate Workshop 2007 Outline First Model (spring) Potential Problems How to improve Second Model (beam) Results The first model: Spring Model We have observations : (t1; y1); … ; (tm; ym). The goal is to estimate the unknown parameters C and K. Checking the assumptions Homoscedasticity Assumption: Normality assumption Independence Assumption Beam Model Cost Function Minimize it. How? 7 parameters, YI, CI, ρ, etc Extremum may be dense in parameter space Find “reasonable” values Set beam to same as patch, search near given data, try to minimize a new cost function A new cost function Needs to take into account spatial variations as well as frequency variations from the model and the data So we use a weighted least squares cost function minimized a simplex method (fminsearch). This doesn’t work. Phase change too much to overcome. What now then? Limit the search. Fewer parameters, smaller variations. And then, it works! (kind of) The data The parameters found gamma 0.17916724273807 air damping YI_beam ** 0.00000669769791 beam -- Youngs modulus CI_beam 0.01013780051727 beam -- internal damping Kp 0.00036441688104 Kp for beam rho_patch ** 0.08291057427864 linear density of patch YI_patch ** 0.20912090251615 patch -- Youngs modulus CI_patch 0.00009783931434 patch -- internal damping Thank You Any Questions? * * Our cost function: Now we need to minimize the cost function. After running the script, we get: C=0.7284; K=1537.8 *