i need to solve I(dw/dt)=MR-MG

3 次查看（过去 30 天）
Casper Stroem
Casper Stroem 2017-10-4
回答： Walter Roberson 2017-10-4
i am a new to matlab and i need to solve I(dw/dt)=MR-MG this is what i have written so far but not sure how to solve numerically for w(t)
clear all
close all
clc
lambda=[14.06;11.71;10.04;8.79;7.81;7.03;6.39;5.86;5.41;5.02;4.69;4.39;4.13;3.90;3.70;3.51;3.35;3.19;3.06;2.93;2.81];
Cp=[.145;.315;.405;.452;.473;.476;.468;.452;.432;.406;.377;.343;.305;.267;.231;.199;.170;.146;.125;.107;.092];
R=30; %meter radius
rho=1.225; %kg/m^3
Vo=10; %m/s
omega=lambda*Vo/R;
I=7.95*10^5 %kgm^2
Mr=0.5*rho*Vo^3*pi*R^2.*Cp./omega; %aero torque Nm
if omega<2.1;
Mg=0;
else
Mg=4.35*10^6*(omega-2.1);
end

请先登录，再回答此问题。

回答（1 个）

Walter Roberson
Walter Roberson 2017-10-4
Change
if omega<2.1;
Mg=0;
else
Mg=4.35*10^6*(omega-2.1);
end
to
Mg = zeros(size(omega));
omega_boundary = 2.1;
mask = omega >= omega_boundary;
Mg(mask) = 4.35*10^6*(omega(mask) - omega_boundary);

类别

在 Help Center 和 File Exchange 中查找有关 Inertias and Loads 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by