查看: 101|回复: 0|关注: 0

[未答复] fsolve 解带参数的方程

[复制链接]

新手

7 麦片

财富积分


050


2

主题

4

帖子

0

最佳答案
发表于 5 天前 | 显示全部楼层 |阅读模式
本帖最后由 dtspslsh 于 2020-3-26 07:18 编辑

相关代码如下:
syms a p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p11 p12 p13 p14 p15 p16
mu1=1.1094
mu2=0.1931
mu3=-0.5
sigma1=1
sigma2=1
sigma3=1
rv=exp(mu3+(sigma3^2)/2)
A=1-0.5*erfc(-(log(a)-mu1)/(sigma1*2^0.5))
B=1-0.5*erfc(-(log(a)-mu1)/(sigma1*2^0.5))
e1=rv/(rv+vpaintegral(A*B^(p3 + p12 + p7+p6),a,[0 inf]))-p1==0
e2=rv/(rv+vpaintegral(A*B^(p4 + p13 + p7),a,[0 inf]))-p2==0
e3=rv/(rv+vpaintegral(A*B^(p1 +  p6),a,[0 inf]))-p3==0
e4=rv/(rv+vpaintegral(A*B^(p16 + p2 + p9),a,[0 inf]))-p4==0
e5=rv/(rv+vpaintegral(A*B^(p7 + p14 + p9),a,[0 inf]))-p5==0
e6=rv/(rv+vpaintegral(A*B^(p3 + p13 + p8+p1),a,[0 inf]))-p6==0
e7=rv/(rv+vpaintegral(A*B^(p1 + p2 + p5),a,[0 inf]))-p7==0
e8=rv/(rv+vpaintegral(A*B^(p9 + p6 + p16),a,[0 inf]))-p8==0
e9=rv/(rv+vpaintegral(A*B^(p5 + p8 + p4),a,[0 inf]))-p9==0
e10=rv/(rv+vpaintegral(A*B^(p11 + p12 + p15),a,[0 inf]))-p10==0
e11=rv/(rv+vpaintegral(A*B^(p10 + p14),a,[0 inf]))-p11==0
e12=rv/(rv+vpaintegral(A*B^(p1 + p15 + p10),a,[0 inf]))-p12==0
e13=rv/(rv+vpaintegral(A*B^(p6 + p2 + p16),a,[0 inf]))-p13==0
e14=rv/(rv+vpaintegral(A*B^(p11 + p15 + p5),a,[0 inf]))-p14==0
e15=rv/(rv+vpaintegral(A*B^(p10 + p12 + p14),a,[0 inf]))-p15==0
e16=rv/(rv+vpaintegral(A*B^(p13 + p8 + p4),a,[0 inf]))-p16==0

在这里我只想得到关于p1-p16的值,a只是用来定积分的参数
请问下接下来的应该怎么算呢? 用vpasolve会卡死,似乎只能用fsolve...

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表