Question

matlab知道y求x问题

n=600:10:4000;CDA=2.77;
f=0.013;m=3880;g=9.8;
Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;
h=0.85;%机械效率
r=0.367;i0=5.83;
ig=0.793;
ua=0.377*r.*n/(i0*ig); Ft=Tq*ig*i0*h/r;
Fw=CDA*(ua5).^2/21.15;
D=(Ft-Fw)/(m*g);
上面是程序，我想求一下D=f（D=0.013）时ua的值。请大神伸出援手帮帮小弟，谢谢！！！

Answer 1

应该就在99附近吧

要相对准确的值，可以用插值

interp1(D,ua,0.013)

ans =

   99.0834

但由于你的D(ua)对于函数，所以ua(D)不是一个单质函数

当D>0.03的那一段时，应该是有两个对应的ua值的，但是插值只能返回一个

追问

求F与Ft的交点。用上面哪种方法怎么不行呢？？？求解答。唉，说我字数超限，我只能图片了，不好意思。

追答

这是ua和F（蓝色），Ft（绿色） F-Ft(红色）的关系

交点就是F-Ft等于0时候的点

interp1(F-Ft,ua,0)

ans =99.0834

显然和上一个问题是同一点

interp1(ua,F,99.0834)

ans =1.7801e+003

F和Ft相交于 99.0834 1.7801e+003

Answer 2

既然D=(Ft-Fw)/(m*g); 你可以整理成Fw=.....
进而可以整理成 ua=.... 等号右边的都是已知 就等求解了啊

追问

。。。。我试了这种方法，不行。能帮忙试一下吗？？？

追答

提示什么? 把代码截图贴下 帮你看看

追问

为什么没有结果呢

你说的太麻烦了，我想在程序后面写句命令直接出结果求出来

追答

看看是不是你想要的 如果你想用solve的话 那么需要你把所有的式子都列进去 solve中每一个不是数的变量都被认为是自变量

n=600:10:4000;
CDA=2.77;
f=0.013;m=3880;g=9.8;
Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;
h=0.85;%机械效率
r=0.367;i0=5.83;
ig=0.793;
% ua=0.377*r.*n/(i0*ig);
Ft=Tq*ig*i0*h/r;
% Fw=CDA*(ua).^2/21.15;
% D=(Ft-Fw)/(m*g);
Fw=Ft-f*(m*g);
ua=sqrt(Fw*21.15/CDA)

追问

谢谢！！！