24小时热门版块排行榜

返回列表

njfu3110246

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 2078.9
帖子: 81
在线: 20.4小时
虫号: 1768130

[交流] 该方程无解，求大神帮助

各位大神好，自编了一个程序，求解Rmu，但无解，求助。

clear;clc;
syms beta Rmu;
beta=3.7; %目标beta
KAmu=0.96;KAcv=0.06;KAsigma=KAmu*KAcv;
KPmu=1.00;KPcv=0.05;KPsigma=KPmu*KPcv;
KQ1mu=0.66;KQ1cv=0.19;KQ1sigma=KQ1mu*KQ1cv;
KQ2mu=0.90;KQ2cv=0.04;KQ2sigma=KQ2mu*KQ2cv;
KQ3mu=0.72;KQ3cv=0.12;KQ3sigma=KQ3mu*KQ3cv;
KQ4mu=0.75;KQ4cv=0.07;KQ4sigma=KQ4mu*KQ4cv;
gmu=1.060;gcv=0.070;gsigma=gmu*gcv;
qmu=0.524;qcv=0.288;qsigma=qmu*qcv;
p=1;rG=1.2;rQ=1.4;wc=1.0; KDOL=KQ3mu;
KBmu=1.00;KBcv=0.05;KBsigma=KBmu*KBcv;
fmu=141.7334;fcv=0.1;fsigma=fmu*fcv;fk=fmu-1.645*fsigma;
KB0=KBmu;g0=gmu;q0=qmu;Rmu0=1; %设定初值
Rcv=(KAcv^2+KPcv^2+KQ1cv^2+KQ2cv^2+KQ3cv^2+KQ4cv^2+fcv^2)^0.5;
%%%%%建立方程式并求解%%%%%%
R0=Rmu;Rnorm_sigma=R0*sqrt(log(1+Rcv^2));Rnorm_mu=R0*(1-log(R0)+log(Rmu/(1+Rcv^2)^0.5));
qnorm_alpha=pi/qsigma/6^0.5;qnorm_u=qmu-0.57722/qnorm_alpha;qnorm_F=exp(-exp(-qnorm_alpha*(q0-qnorm_u)));qnorm_f=qnorm_alpha*exp(-qnorm_alpha*(q0-qnorm_u))*exp(-exp(-qnorm_alpha*(q0-qnorm_u)));
qnorm_sigma=normpdf(norminv(qnorm_F,0,1),0,1)/qnorm_f;
qnorm_mu=q0-qnorm_sigma*norminv(qnorm_F,0,1);
Ru=1;
Su1=-(g0+p*qnorm_mu)/(rG+wc*rQ*p);
Su2=-KB0/(rG+wc*rQ*p);
Su3=-KB0*p/(rG+wc*rQ*p);
F=(Rnorm_mu+Su1*(KBmu-KB0)+Su2*gmu+Su3*qnorm_mu)/((Rnorm_sigma^2+Su1^2*KBsigma^2+Su2^2*gsigma^2+Su3^2*qnorm_mu^2)^0.5)-beta
Rmu=solve(F)
Rmu=double(Rmu)

计算结果为方程无解，之后solve、fsolve都用过，同样无解。请各位大神帮忙找找原因。谢谢。
完整程序可见附件。

回复此楼

» 猜你喜欢

球磨粉体时遇到了大的问题，请指教！已经有9人回复
江汉大学解明教授课题组招博士研究生/博士后已经有3人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2019-05-07 21:46:12

已阅关注TA 给TA发消息送TA红花 TA的回帖

njfu3110246

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 2078.9
帖子: 81
在线: 20.4小时
虫号: 1768130

结果显示：

F =

(Rmu*(log((100*219^(1/2)*Rmu)/1533) - log(Rmu) + 1) - 701846099337675961/1170935903116328960)/((22895321347407492655221715902025*Rmu^2)/324518553658426726783156020576256 + 88683260162590427/2305843009213693952)^(1/2) - 37/10

Warning: Explicit solution could not be found.
> In solve at 169
  In rR_cov_norm_exrI_2_5_Failure at 57

Rmu =

[ empty sym ]

Rmu =

   []

Rnorm_mu =

   []

Alpha =

   []

RAlpha =

   []

Error using  /
Matrix dimensions must agree.

Error in rR_cov_norm_exrI_2_5_Failure (line 69)
      SAlpha1=KBsigma*Su1/Alpha;