clear all, close all [v_a,interp_power] = load_aircraft_eff() % % %Interpolierte Effizienzkurve plotten % plot(measured_airspeed, measured_power, 'linewidth', 3) % hold on % plot(v_a, interp_power, 'linewidth', 2) % hold off % legend('Rohdaten','Interpolierte Daten','location','southeast' ) % xlabel('airspeed/ km/h') % ylabel('power/ Watt') % ylim([0 270]) % xlim([-40 120]) v_w=40; %windspeed in km/h WD=45; %windangle in ° %Track points typischer Flug P_x = [0 30 30.5 31 30.5 32 32 33 33 34 34 35 35 36 36 39 20]; P_y = [0 40 35 35.6 43 40 35 34 42 40 35 37 45 44 34 35 5]; % %Track points im Kreis % P_x = [0 1 2 2 1 0 -1 0]; % P_y = [0 1 1 0 -1 -1 0 0]; %Track points einmal hin und her P_x = [0 1 0 ]; P_y = [0 0 0 ]; trip_length=0; for track=1:(numel(P_x)-1) %optional zum debuggen % figure(4) % h = plot([P_x(track) P_x(track+1)] , [P_y(track) P_y(track+1)], 'linewidth', 3); % hold on TC=(atan2d((P_x(track+1)-P_x(track)),(P_y(track+1)-P_y(track)))); track_length(track)= sqrt((P_y(track+1)-P_y(track))^2+(P_x(track+1)-P_x(track))^2); [v_best(track),res_eff_best(track)] = eff_calc(interp_power,v_a, TC, v_w, WD ); trip_length=trip_length+track_length(track); end figure(2) colorbar_eff=jet(200); trip_eff=0; for track=1:(numel(P_x)-1) trip_eff = trip_eff + track_length(track)/trip_length*res_eff_best(track); h = plot([P_x(track) P_x(track+1)] , [P_y(track) P_y(track+1)], 'linewidth', 3); hold on set(h, 'color', squeeze((colorbar_eff(ceil(200-res_eff_best(track)),:)))); end trip_eff legendCell = cellstr(num2str(ceil(v_best)', 'v\\_a\\_best=%-d km/h')); legend([legendCell],'location', 'northwest');