function naloga12_1sklop %dvoetažni enoladijski okvir %vnos parametrov naloge m1=10000; %masa 1. etaže (kg) m2=7000; %masa 2. etaže (kg) k=4.5*10^7; %togost etaže (N/m) %omega=46.7; %lastna frekvenca vzbujanja (Hz) omega=70 c=0 a=5 %omega=100; %druga frekvenca vzbujanja ki je daleč stran (Hz) %PARAMETRI čASOVNE INTEGRACIJE tstart = 0; %začetni �as tfinal = 10; %končni �as y0 = [0,0,0,0]'; %za�etne vrednosti vektorja neznank (stolpec vrednosti) %y0 = [ 0.010695 , 0.0033604, -0.009639,0,0,0]'; %za�etne vrednosti vektorja neznank (stolpec vrednosti) Mm=diag([m1,m2]); TK=[-k k 0 -k ]; [vekt,LV]=eig(TK,Mm); %vekt(:,2)' %Omege=sqrt(LV) %MASNA MATRIKA options = odeset('Mass',@masna_matrika); % NUMERI�NO RE�EVANJE PROBLEMA [t,y] = ode23t(@desne_strani,[tstart tfinal],y0,options,m1,m2,c,k,omega,a); %uporabite tudi druge metode dru�ine ODE %GRAFI�NI PRIKAZ REZULTATOV figure(13); clf plot(t,y(:,1:2)); legend('1. et','2. et'); % -------------------------------------------------------------------------- %OPIS DIFERENCIALNIH ENA�B function dydt =desne_strani(t,y,m1,m2,c,k,omega,a) dydt = [ y(3) y(4) -k*y(1)+k*y(2)-m1*a*sin(omega*t) -m2*a*sin(omega*t)-k*y(2) ]; %MASNA MATRIKA SISTEMA ENA�B function M = masna_matrika(t,y,m1,m2,c,k,a,omega) % Mass matrix function M = zeros(4,4); M(1,1) = 1; M(2,2) = 1; M(3,1) = c; M(3,2) = -c; M(3,3) = m1; M(4,2) = c; M(4,3) = m2; M(4,4) = m2; %M = zeros(4,4); % M(1,1) = 1; % M(2,2) = 1; % M(3,3) = m1; % M(4,3) = m2; % M(4,4) = m2;