What am I doing wrong?
1 次查看(过去 30 天)
显示 更早的评论
Hello!
I am trying to solve an equation but I get this error:
Error using sym>notimplemented (line 2682)
Function 'ge' is not implemented for MuPAD symbolic objects.
Error in sym/ge (line 823)
notimplemented('ge');
Error in Programa_drag_torque_paper (line 41)
a=double(s{i}(1)>=r_1)&&double(s{i}(1)<=r_2);
I think that the problem is related with how are determined the variables (sym, double,...) but I don't know how to fix it.
Here is the code:
clear all
clc
r_1=178*10^-3;
r_2=218*10^-3;
h=0.2*10^-3;
N=10;
Q_o=5.83*10^-5;
rho=820;
mu=0.0082
sigma=0.030
omega=[700:100:2200]*(2*pi/60);
omega_rpm=[700:100:2200];
theta=0*(2*pi/60);
r_m=(r_1+r_2)/2;
Re_h=(rho*omega*r_1*h)/mu;
m=zeros(16,1);
G_r=m;
f=m;
syms r_o_1
syms r
s=cell(16,1);
r_o_final_1=m;
T_1=m;
t_final_1=m;
for i=1:16
if Re_h(i)>=500
f(i)=0.885*Re_h(i)^-0.367;
else
f(i)=0.09;
end
G_r(i)=1/12*(1+(0.00069.*Re_h(i).^0.88));
s{i}=solve(((rho*(omega(i)^2/2))*(f(i)+(1/4))*(r_o_1^2))-((mu*Q_o*r_o_1)/(2*pi*r_m*h^3*G_r(i)))-((rho*(omega(i)^2/2))*(f(i)+(1/4))*r_1^2)-((2*sigma*cos(theta))/h)+((mu*Q_o*r_1)/(2*pi*r_m*h^3*G_r(i))),r_o_1);
a=double(s{i}(1)>=r_1)&&double(s{i}(1)<=r_2);
b=double(s{i}(2)>=r_1)&&double(s{i}(2)<=r_2);
if a
r_o_final_1(i)=s{i}(1);
elseif b
r_o_final_1(i)=s{i}(2);
else
r_o_final_1(i)=r_2;
end
T_1(i)=int(((mu*omega(i)*r^3)/h)*(1+(0.0012*Re_h(i)^0.94)),r_1,r_o_final_1(i));
T_final_1(i)=2*N*pi*T_1(i);
end
It would be great if you could find a solution.
Thank you so much for your answers!!!
0 个评论
采纳的回答
Walter Roberson
2016-10-13
a=double(s{i}(1)) >=r_1&&double(s{i}(1)) <=r_2;
b=double(s{i}(2))>=r_1&&double(s{i}(2)) <=r_2;
Notice I moved the )
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!