Mam napisać program w MatLabie który opisze drgania takiego oscylatora gdy działa na niego siła wymuszająca i siła oporu.
Czy siła oporu to siła tarcia?
Bo np. w literaturze znalazłem takie dwa wzory: oraz i mam problem ze skryptem, zresztą nie tylko ja bo koledzy też, stworzyliśmy coś takiego:
t(1)=0; % początek pomiaru czasu dt=0.002; %krok czasowy N = 15000; %liczba kroków czasowych x(1)=1; v(1)=0; %warunki początkowe v(2)=v(1); %dla poprawy dokładności obliczeń %Parametry oscylatora k = %stała sprężystości sprężyny (N/m) m= %masa doczepiona do sprężyny gamma= %współczynnik tłumienia (tarcia) (jaką to ma jednostkę?) Fm = %amplituda siły wymuszającej harmonicznej W= %częstość siły wymuszającej F(1)= %wartość siły zewnetrznej w chwili początkowej w0=sqrt(k/m) p=10^(-7); for i=1:N %v(i+1)=v(i)-(k/m)*x(i)*dt-gamma*v(i)*dt-F(i)/m*dt; v(i+1)=v(i)-(k/m)*x(i)*dt-gamma*v(i)/abs(v(i)+p)*dt-F(i)/m*dt; %z siłą tarcia x(i+1)=x(i)+v(i+1)*dt; t(i+1)=i*dt; F(i+1)=Fm*cos(W*i*dt); end
I jest mały spór, czy w pętli na końcu nie powinien stać sinus?
I po drugie, co wiąże się z moim pytaniem z posta, jak jestem na początku pętli to która linijkę do obliczeń wybrać obliczając prędkość w danej chwili?