应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 该方程无解，求大神帮助
111/1返回列表
查看: 879  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

njfu3110246

木虫 (小有名气)

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

各位大神好，自编了一个程序，求解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都用过，同样无解。请各位大神帮忙找找原因。谢谢。
完整程序可见附件。
回复此楼

» 猜你喜欢

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

查看全部散金贴
1楼 2019-05-07 21:46:12
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

njfu3110246

木虫 (小有名气)
结果显示：

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;
赞一下 回复此楼
2楼2019-05-07 21:47:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
8771645213楼
2019-05-07 22:13  
njfu3110246(金币+1): 谢谢参与
。 发自小木虫Android客户端
tzynew4楼
2019-05-07 22:14  
njfu3110246(金币+1): 谢谢参与
syhorchid5楼
2019-05-07 22:41  
njfu3110246(金币+1): 谢谢参与
谜团小子6楼
2019-05-07 23:10  
njfu3110246(金币+1): 谢谢参与
祝福 发自小木虫Android客户端
iamluo7楼
2019-05-08 00:03  
njfu3110246(金币+1): 谢谢参与
发自小木虫Android客户端
taro18楼
2019-05-08 00:51  
njfu3110246(金币+1): 谢谢参与
发自小木虫Android客户端
xiaoqiuqiu209楼
2019-05-08 01:07  
njfu3110246(金币+1): 谢谢参与
xhmaohan10楼
2019-05-08 02:22  
njfu3110246(金币+1): 谢谢参与
发自小木虫Android客户端
nono200911楼
2019-05-09 18:16  
njfu3110246(金币+1): 谢谢参与
·
相关版块跳转 我要订阅楼主 njfu3110246 的主题更新
111/1返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 化学专业申博 +3 赵子羊 2026-05-23 4/200 2026-05-24 18:10 by 工大学长
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 1rx34o113h 2026-05-23 3/150 2026-05-24 17:41 by 0i3mu4vkjz
[基金申请] 评审有感 +16 popular289 2026-05-18 27/1350 2026-05-24 17:34 by hhs666
[教师之家] 论文撤稿了 +4 bjvtcliu 2026-05-24 7/350 2026-05-24 17:29 by bjvtcliu
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 hvkbtfonbv 2026-05-23 4/200 2026-05-24 17:21 by 75ui6h7z2t
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:10 by 75ui6h7z2t
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:01 by 75ui6h7z2t
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 5/250 2026-05-24 16:21 by hhx1yx9evi
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 4/200 2026-05-24 16:16 by hhx1yx9evi
[基金申请] 河北省自然科学基金 +6 Peterchao 2026-05-18 9/450 2026-05-24 16:02 by 130067131
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[论文投稿] 投稿求助，期刊 +4 希冀，有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀，有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[有机交流] 反应很差，大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博，可以边工作边读全日制嘛？ 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
信息提示
关闭
请填处理意见
关闭