This preview shows pages 1–3. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.View Full Document
Unformatted text preview: the relevant MATLAB code. clear,clc format long Fs = 5000; d = 0.05; Fo = 350; X = 0.11; k = Fs/d; eta = Fo/(k*X); r1 = 0.5; r2 = 2; X1 = abs((Fo/k)/(1-r1*r1+i*eta)); X2 = abs((Fo/k)/(1-r2*r2+i*eta)); r=0:0.001:2; for j = 1:numel(r) X(j) = (Fo/k)/(1-r(j)*r(j)+i*eta); end plot(real(X),imag(X)) xlabel( 'real part of the complex amplitude [m]' ) ylabel( 'imaginary part of the complex amplitude [m]' ) grid on n = find(r==1) text(real(X(n)),imag(X(n)), '\leftarrow \omega = \omega_n' , 'FontSize' ,18) Assigned on 26.10.2009 Monday Due date 02.11.2009 Monday, beginning of the lecture -0.06-0.04-0.02 0.02 0.04 0.06-0.12-0.1-0.08-0.06-0.04-0.02 real part of the complex amplitude [m] imaginary part of the complex amplitude [m] ← ω = ω n...
View Full Document
This note was uploaded on 03/28/2010 for the course AEROSPACE aero520 taught by Professor Mok during the Spring '10 term at Istanbul Technical University.
- Spring '10