用matlab求解代数式,代数式不能得出具体数值解。

>> syms G T Y P O B Q A S Z X R

WT=(A*G*R*S*Z*log(O/P)*((O/P)^(((log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z))))/log(O/P)) - 1))/(B*Q*(log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z)))) - (A*G*R*S*Z*log(O/P)*((O/P)^(((log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z))))/log(O/P)) - 1))/(B*Q*(log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z)))^2) + (A*G*R*S*Z*log(O/P)^2*((O/P)^(((log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z))))/log(O/P)) - 1))/(B*Q*(log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)^2*(log(O/P) - log((X*Y)/(T*Z)))^3) - (A*G*R*S*T*Z*log(O/P)^2*(O/P)^(((log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z))))/log(O/P))*(1/(T*(log(O/P) - log((X*Y)/(T*Z)))) - (log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)/(T*log(O/P))))/(B*Q*(log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z))))

A=2;G=3;Y=0.1;P=10;O=5;B=50;Q=10;S=11;Z=14;X=45;R=10;T=10;

WT=subs(WT)

WT =

(A*G*R*S*Z*log(O/P)*((O/P)^(((log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z))))/log(O/P)) - 1))/(B*Q*(log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z)))) - (A*G*R*S*Z*log(O/P)*((O/P)^(((log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z))))/log(O/P)) - 1))/(B*Q*(log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z)))^2) + (A*G*R*S*Z*log(O/P)^2*((O/P)^(((log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z))))/log(O/P)) - 1))/(B*Q*(log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)^2*(log(O/P) - log((X*Y)/(T*Z)))^3) - (A*G*R*S*T*Z*log(O/P)^2*(O/P)^(((log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z))))/log(O/P))*(1/(T*(log(O/P) - log((X*Y)/(T*Z)))) - (log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)/(T*log(O/P))))/(B*Q*(log(O/P)/(log(O/P) - log((X*Y)/(T*Z))) - 1)*(log(O/P) - log((X*Y)/(T*Z))))

WT =

(462*log(2)*((1/2)^(((log(2) + log(9/280))*(log(2)/(log(2) + log(9/280)) - 1))/log(2)) - 1))/(25*(log(2) + log(9/280))*(log(2)/(log(2) + log(9/280)) - 1)) + (462*log(2)*((1/2)^(((log(2) + log(9/280))*(log(2)/(log(2) + log(9/280)) - 1))/log(2)) - 1))/(25*(log(2) + log(9/280))^2*(log(2)/(log(2) + log(9/280)) - 1)) - (462*log(2)^2*((1/2)^(((log(2) + log(9/280))*(log(2)/(log(2) + log(9/280)) - 1))/log(2)) - 1))/(25*(log(2) + log(9/280))^3*(log(2)/(log(2) + log(9/280)) - 1)^2) - (924*(1/2)^(((log(2) + log(9/280))*(log(2)/(log(2) + log(9/280)) - 1))/log(2))*log(2)^2*(1/(10*(log(2) + log(9/280))) - (log(2)/(log(2) + log(9/280)) - 1)/(10*log(2))))/(5*(log(2) + log(9/280))*(log(2)/(log(2) + log(9/280)) - 1))

>> vpa(WT, 6)

ans =

-2.67714

>>

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值