Vinnaren i pepparkakshustävlingen!
2009-10-16, 13:57
  #1
Medlem
andreells avatar
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)
Citera
2009-10-16, 14:19
  #2
Medlem
evolutes avatar
Du har ju med 1i i din formel för ztopp och detta är den imaginära enheten.
Citera
2009-10-16, 16:17
  #3
Medlem
andreells avatar
Jo det vet jag. Har dock lyckats lösa detta problem genom att ta absolutbeloppet av de tal jag får fram.. (Reella & imaskinära) Men har däremot stött på ett nytt problem. Hur kan jag plotta två diagram samtidigt? När man kör m-filen vill jag att det ska ploppa up två diagram i samma ruta, ne upptill och ett undertill. tack på förhand!
Citera
2009-10-16, 16:25
  #4
Medlem
Är du ute efter t.ex. subplot så du kan ha två diagram i samma fönster eller hold on så du kan rita en kurva till i samma koordinatsystem?
Citera
2009-10-16, 16:26
  #5
Medlem
evolutes avatar
subplot(2,1,1)
plot(...)
subplot(2,1,2)
plot(...)
Citera

Stöd Flashback

Flashback finansieras genom donationer från våra medlemmar och besökare. Det är med hjälp av dig vi kan fortsätta erbjuda en fri samhällsdebatt. Tack för ditt stöd!

Stöd Flashback