怎样用MATLAB编程求出原函数,希望给出详细程序: 函数如下:y'=f*(1+m)*(1+(1+N6*N7)^0.5)/N6

MATLAB求解
y'=f*(1+m)*(1+(1+N6*N7)^0.5)/N6;
N6=1+f*m*(p+r*y-f*c);
N7=(1-m*f^2)/(f^2*(1+m)^2);
其中f,m,p,r,c均为常量,把N6、N7代入第一式之后,也就是左边为y',右边是关于y的函数
求高手解答,希望给出具体MATLAB程序

第1个回答  2010-12-23
syms f m p r c y
a=int('f*(1+m)*(1+(1+(1+f*m*(p+r*y-f*c))*((1-m*f^2)/(f^2*(1+m)^2)))^0.5)/(1+f*m*(p+r*y-f*c))','y')

积分结果为

y=(2.0*(1.0 - (1.0*(f^2*m - 1.0)*(f*m*(p - 1.0*c*f + r*y) + 1.0))/(f^2*(m + 1.0)^2))^(1/2))/r + log(2.0*f*m - 2.0*f*m*(1.0 - (1.0*(f^2*m - 1.0)*(f*m*(p - 1.0*c*f + r*y) + 1.0))/(f^2*(m + 1.0)^2))^(1/2))/r - (1.0*log(2.0*f*m + 2.0*f*m*(1.0 - (1.0*(f^2*m - 1.0)*(f*m*(p - 1.0*c*f + r*y) + 1.0))/(f^2*(m + 1.0)^2))^(1/2)))/r + (log(2.0*f - 2.0*f*(1.0 - (1.0*(f^2*m - 1.0)*(f*m*(p - 1.0*c*f + r*y) + 1.0))/(f^2*(m + 1.0)^2))^(1/2)) - 1.0*log(2.0*f + 2.0*f*(1.0 - (1.0*(f^2*m - 1.0)*(f*m*(p - 1.0*c*f + r*y) + 1.0))/(f^2*(m + 1.0)^2))^(1/2)) + 2.0*(1.0 - (1.0*(f^2*m - 1.0)*(f*m*(p - 1.0*c*f + r*y) + 1.0))/(f^2*(m + 1.0)^2))^(1/2))/(m*r) + (log(f*m*(p - 1.0*c*f) + f*m*r*y + 1.0)*(m + 1.0))/(m*r)本回答被网友采纳
第2个回答  2011-01-02
syms m f p r c y;
N6=1+f*m*(p+r*y-f*c);
N7=(1-m*f^2)/(f^2*(1+m)^2);
y1=f*(1+m)*(1+(1+N6*N7)^0.5)/N6;
y2=finverse(y1,y)

运行结果
y2 =

(2*f*m*y-m*f^2-y^2*f*m*p+y^2*f^2*m*c+2*f*y+1-y^2)/r/y^2/f/m