| 查看: 862 | 回复: 10 | |||
[交流]
该方程无解，求大神帮助
|
|||
|
各位大神好，自编了一个程序，求解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都用过，同样无解。请各位大神帮忙找找原因。谢谢。 完整程序可见附件。 |
回复此楼» 猜你喜欢
材料考研调剂
已经有3人回复
-
材料调剂
已经有12人回复
-
英一数一408，总分284，二战真诚求调剂
已经有14人回复
-
085410 一志愿211 22408分数359求调剂
已经有4人回复
-
271求调剂
已经有19人回复
-
385分 生物学（071000）求调剂
已经有3人回复
-
一志愿安徽大学计算机科学与技术学硕，331分求调剂
已经有3人回复
-
318求调剂，计算材料方向
已经有8人回复
-
291求调剂
已经有25人回复
-
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂
已经有6人回复
» 抢金币啦！回帖就可以得到:
-
天津大学智能人机交互团队招收机械非全硕士带薪培养
+1/497
-
广东工业大学-木质纤维素高值化利用团队招博士研究生
+1/87
-
天津科技大学-（新晋院士团队）先进纤维与纸基功能材料团队B+专业
+1/86
-
焊接专业应届毕业生招聘
+1/81
-
河南工业大学化学化工学院 2026 年硕士调剂正在进行！
+1/43
-
大连工业大学纺织与材料工程学院超临界流体技术课题组 研究生调剂
+2/42
-
汕头大学土木与智慧建设工程系调剂信息
+1/37
-
中山大学地理科学与规划学院退化土壤环境修复课题组招聘博士后及科研助理
+1/35
-
招聘 博士教师有编制 ，材料加工工程，金属材料方向 月薪1.8W起， 地点陕西渭南
+1/32
-
西南科技大学杨文彬教授课题组接收材料、化工、化学类专业研究生调剂 。
+1/26
-
材料/化学/环境 5-10人专硕, A类线，武汉纺织大学材料学院，能源高分子方向
+1/17
-
福建师范大学环境与资源学院2026级学硕专硕调剂（网络远程面试）
+1/10
-
西安工程大学能源化工新材料研发项目组接收2026年计算机类调剂研究生
+1/9
-
天津科技大学经管学院姚伟教授团队：知识生态研究中心接受管理科学与工程2026调剂生
+1/9
-
浙江农林大学-桂金山课题组（植物细胞壁-生物质与生物能源）-2026年硕士研究生调剂
+1/9
-
南京林业大学化工院有机合成方向硕士招调剂
+1/8
-
西北农林科技大学水产养殖专业学术型硕士调剂名额1个
+1/5
-
0831生医工调剂本9一志愿北航生医工，338分求调剂
+1/4
-
西工大海洋研究院招收微尺度传热方向博士后1名
+1/4
-
2026年江西科技师范大学Tide课题组招收化学、生物、药学、材料相关学硕或专硕调剂生
+1/1
|
结果显示： 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
简单回复
2019-05-07 22:13
njfu3110246(金币+1): 谢谢参与
。 发自小木虫Android客户端
tzynew4楼
2019-05-07 22:14
njfu3110246(金币+1): 谢谢参与
2019-05-07 22:41
njfu3110246(金币+1): 谢谢参与
谜团小子6楼
2019-05-07 23:10
njfu3110246(金币+1): 谢谢参与
祝福 发自小木虫Android客户端
iamluo7楼
2019-05-08 00:03
taro18楼
2019-05-08 00:51
2019-05-08 01:07
njfu3110246(金币+1): 谢谢参与
xhmaohan10楼
2019-05-08 02:22
nono200911楼
2019-05-09 18:16
njfu3110246(金币+1): 谢谢参与
·
