24小时热门版块排行榜    

查看: 347  |  回复: 0

如丧莫言

新虫 (小有名气)

[求助] 求助用1stopt解一个多元方程组,因为我不会用1stopt

Function
C1*exp(-gama*x0)+C2*exp(gama*x0)=0;
C1+C2-C3-C4-a1=0;
C1-C2-C3+C4=0;
C3*exp(gama*x1)+C4*exp(-gama*x1)+a1-0.002=0;
-0.5*(r-k2)*x1^2+(r-k2)*x1+d-E*A*0.002=0;
E*A*gama*(C3*exp(gama*x1)-C4*exp(-gama*x1))+(r-k2)*x1-0.5*l*(r-k2)=0;
k2*x1+0.5*l*(r-k2)-(k1/gama)*(C1*(1-exp(-gama*x0))+C2*(exp(gama*x0)-1)+C3*(exp(gama*x1)-1)+C4*(1-exp(-gama*x1))+gama*a1*x1)=0;

其中C1 C2 C3 C4 x0 x1 d为未知数,其余常数有下列代码决定:
l=20;
k1=1.2e7;
k2=24000;
E=2.1e11;
A=7745e-6;
gama=sqrt(k1/(E*A));
q1=exp(0.5*l*gama);
q2=exp(-0.5*l*gama);
detaT=30;
detaP=E*A*1.18e-5*detaT;
r=detaP/l;
a1=r/k1;
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 如丧莫言 的主题更新
信息提示
请填处理意见