Ekspert i MatLab

Hvad er dit problem?

Det hele omhandler at jeg skal plotte en en planet masse M som funktion af en parameter m_x . problemet er nu at værdierne for dette M bliver så små at matlab giver mig de samme tal ud for hver eneste værdi af m_x. Jeg bruger endda format double, men dette kan ikke hjælpe mig. Nu ved jeg at det er svært at forholde dig til dette da du næsten ikke ved hvad der foregår. Måske du kan kigge på det jeg skal plotte (formlerne er sat helt rigtigt ind, det har jeg checket for)?

Okay jeg kan da prøve.

Koden står neden under. Du kan eventuelt prøve at sætte kopiere dne ind i dit matlab, også ser du at grafen for M bliver en konstant linje, som det ikke skal være.

format long

mx = 1 : 10^(27) : 2*10^(31); // dette er mine værdier for mx

Rr = ((2*10^9)./((mx).^(-(1/3))));
Mr = ((2*10.^(31)./(mx)));

C0 = (3*10^(12).*Rr)./(6.67*10^(-8)*5*Mr);

C1 = (4*pi*5*Rr.^3)./(3*Mr);

M = 2*10.^(31).*(1+C1.*mx.^(-2).*(((1-C0.*mx.^(-4/3)).^(-3))-1));

R = 2*10^9.*(1-C0.*mx.^(-4/3)).^(-1);

(mx.^(-4/3).*C0)   // problemet ligger her. Denne del kommer fra udtrykket for M. Dette led giver mig ens værdier.

plot(mx,M)
 

Jeg håber virkelig at du kan hjælpe mig

#4

Mener du ikke

mx = 1*10^(27) : 2*10^(31); // dette er mine værdier for mx

med et gangetegn i stedet for et : efter det første 1  ?

#5

Nej, det jeg vil have frem er værdier fra 1 til 2*10^(31) i steps af 10^(27). Disse er vilkårligt valgt.

Hvad så wuy123 kan du komme frem til noget fornuftigt?

#6

Men mener du ikke værdier fra 1·10³¹ til 2·10³¹ med den steplængde?

Du regner med alt for små tal og derfor opstår der en masse afrundingsfejl. Matlab kan f.eks. ikke kende forskel på tallene 1-1e-17 og 1. Du bliver nødt til at undgå den slags operationer.

#6

jo det kan jeg sagtens bruge, men det vil ikke ændre på de små tal.

#9

Det kan jeg ikke på nogen måde. Jeg aner virkelig ikke hvad der skal gøres her

Nu har jeg sku fundet ud af det. Når jeg sætter værdierne for C0 og C1 ind i M så ryger al afhængigheden af mx. Så min formel for M afhænger slet ikke af mx

#7

Er du sikker på, at dine formler er korrekte? Her er en tabel over de indgående størrelser beregnet i Excel med dine udtryk:

#7

Hej Andersen, jeg sætter meget pris på at du vil hjælpe. Jeg har fundet ud af hvad problemet er. Som skrevet oven over hvis jeg smider mine udtryk for C0 og C1 ind i M så vil M ikke mere afhænge af mx. Dette er pga at jeg har plottet min M som funktion af 1/M_ref hvor M_ref = M₀/mx hvor M₀ bare er det første udtryk i fomlen for M. Altså er M₀ = 2*10^(31). Men i stedet for at jeg bruger den eksakte værdi for M₀ så bruger jeg bare at M₀ = mx * M_ref hvor jeg nu pludselig få ren lineær afhængighed af mx i M. Også plotter jeg bare M som funktion af M_ref. Jeg håber min forklaring var tydelig nok.

#13

I den slags parametriske studier er det ofte en fordel at danne dimensionsløse størrelser, hvis forløb man ønsker kortlagt. Det giver som regel størrelser inden for et par dekader, og kan også give bedre stabilitet, rent numerisk.