Mirrors_Framework/paparazzi
1
0
Fork 0
mirror of https://github.com/paparazzi/paparazzi.git synced 2026-06-08 02:15:53 +08:00
Code Issues Actions 27 Packages Projects Releases Wiki Activity

*** empty log message ***

Browse Source
This commit is contained in:
Antoine Drouin
2007-02-27 19:42:57 +00:00
parent 5892f23077
commit dba078018c
2 changed files with 219 additions and 265 deletions
Download Patch File Download Diff File Expand all files Collapse all files
sw/logalizer/matlab/plot_prop.m
+13 -6
View File
@@ -15,22 +15,29 @@ function [] = plot_prop()
torque_n = torque_g .* ( 9.81 / 1000 );
K = fminsearch(@lin_model, [1, 1]);
Kt = K(1)
Kq = K(2)
Kt = K(1);
Kq = K(2);
% F = 0.5 * rho * prop_area * ct * prop_rad^2 * omega^2
rho = 1.225;
prop_area = 0.005;
prop_rad = 0.125;
Ct = Kt / (0.5 * rho * prop_area * prop_rad^2)
Cq = Kq / (0.5 * rho * prop_area * prop_rad^2)
plot(omega_square, thrust_n, ...
omega_square, Kt * omega_square, ...
omega_square, torque_n, ...
omega_square, Kq * omega_square )
title('EPP 10*4.5');
title('Propeller caracterisation (EPP 10*4.5)');
legend('thrust', 'fitted thrust', 'torque', 'fitted torque');
xlabel('omega square in rad^2/s^2');
ylabel('forces in N');
tb = annotation('textbox', [0.2 0.65 0.1 0.1]);
set(tb, 'String', ...
[sprintf('Kt = %.2e\n\n', Kt) ...
sprintf('Kq = %.2e', Kq)], ...
[sprintf('Thrust = %.2e * omega^2 ( Ct = %.2e )\n\n', Kt, Ct) ...
sprintf('Torque = %.2e * omega^2 ( Cq = %.2e )', Kq, Cq)], ...
'FitHeightToText', 'on' );
function [sse] = lin_model(params)
sw/logalizer/matlab/quad_sim.mdl
+206 -259
View File
File diff suppressed because it is too large Load Diff