Har problem med att få en vettig graf utav detta. Jag får komplexa tal som "ztopp" vilket jag inte vet om det stämmer? Uppgiften är iallafall att plotta en graf för "ztopp", jag antar att det är funktion utav hastigheten i detta fall. Någon som kan hjälpa? Tack på förhand..
% kvartsbilsmodell :
% Values: Vägbana :
lambda=16; %m
zeta=0.04; %m
A = dlmread('bil246.dat');
v=A(:,1); V=v(1:2:end); %Bilens hastighet (km/h)
m=325; %kg
mt=45; %kg
kt=500000; %N/m
k=16000; %N/m
ct=0; %Ns/m
c=200; %Ns/m
f=V./lambda ; %Hz
w=2*pi.*f; %vinkelhastighet
%rörelseamplituden:
ztopp=(((k+1i.*w.*c).*(kt+1i.*w.*ct)).*zeta)./((k+1i.*w.*c)-((w.^2).*m).*...
[((k+kt)+1i.*w.*(c+ct)-(w.^2).*mt)]-(k+1i.*w.*c).^2)
plot(V,ztopp)