0
svar
62
visningar
hejjagbehöverhjälp är nöjd med hjälpen
Fast på en deluppgift av ett projekt. (b). Jag får inte till RK4 metoden och dessutom ingen plot. (ljudvågor)
såhär långt har jag kommit hittills:
b0= 12.3;
c=[5000; tand(b0)]
p1=-20.2090;
p2=17.3368;
p3= 272.9057;
p4=0.7528;
Z = @(z) 4800+ p1 + p2*(z/1000)+ p3*(exp(-p4*z/1000));
Zderv = @(z) p2/1000- (p3*p4*exp(-p4*z/1000)/1000);
q0 =(Z(5000)/cos(12.3).^2);
S= @(x,c) [c(2); -q0.*(Zderv(1)/Z(c(1)).^3)]
xtol= 30*6076;
h =10;
while x < xtol
k1= S(x,c);
k2=S.*[x+0.5.*h, c+0.5*h.*k1]
k3=S.*[x+0.5.*h, c+0.5*h*k2]
k4=S.*[x+h, c+h*k3]
end
skulle verkligen uppskatta lite hjälp!