matlab解线性方程组log方程组

查看: 1671|回复: 5|关注: 0
怎么解含参数k的非线性方程组解析解
-x(1)^2+exp(-x(2))=0
k1*x(1)exp(x(2))-k2*x(1)exp(x(2))+k2-2*x(1)+2*x(1)^3*exp(x(2))=0
要使结果用k1,k2表示
下面是用Mathematica软件的结果,有两组
{{x1 -& 1/4 (k2 - Sqrt[8 k1 + (-8 + k2) k2]),
&&x2 -& Log[(4 k1 + k2 (-4 + k2 + Sqrt[8 k1 + (-8 + k2) k2]))/(2 (k1 - k2)^2)]},
{x1 -& 1/4 (k2 + Sqrt[8 k1 + (-8 + k2) k2]),
&&x2 -& Log[(4 k1 + k2 (-4 + k2 - Sqrt[8 k1 + (-8 + k2) k2]))/(2 (k1 - k2)^2)]}}
xmdalong 发表于
下面是用Mathematica软件的结果,有两组
{{x1 -& 1/4 (k2 - Sqrt[8 k1 + (-8 + k2) k2]),
&&x2 -& Log[(4&&...
能用matlab解吗?
xmdalong 发表于
下面是用Mathematica软件的结果,有两组
{{x1 -& 1/4 (k2 - Sqrt[8 k1 + (-8 + k2) k2]),
&&x2 -& Log[(4&&...
太谢谢你了,能不能帮我再解一个方程组
-x(1)^2+exp(-x(2))=0
k1*x(1)exp(x(2))-k2*x(1)exp(x(2))+k2-2*x(1)+2*x(1)^3*exp(x(2))=0
x(1)-x(3)^2+x4=0
-k3*x(3)-k4*x(1)+k4*x(3)^2-k4*x(4)+x(1)^2-exp(-x(2))+2x(1)*x(3)-2x(3)^3+2x(3)*x(4)=0
结果也是要用k表示
关注者: 6
syms k1 k2 x1 x2
eq1=-x1^2+exp(-x2);
eq2=k1*x1*exp(x2)-k2*x1*exp(x2)+k2-2*x1+2*x1^3*exp(x2);
[x1 x2]=solve(eq1,eq2,'x1','x2')
关注者: 6
tititata 发表于
太谢谢你了,能不能帮我再解一个方程组
-x(1)^2+exp(-x(2))=0
k1*x(1)exp(x(2))-k2*x(1)exp(x(2))+k2-2*x ...
syms k1 k2 k3 k4 x1 x2 x3 x4
eq1=-x1^2+exp(-x2);
eq2=k1*x1*exp(x2)-k2*x1*exp(x2)+k2-2*x1+2*x1^3*exp(x2);
eq3=x1-x3^2+x4;
eq4=-k3*x3-k4*x1+k4*x3^2-k4*x1+x1^2-exp(-x2)+2*x1*x3-2*x3^3+2*x3*x4;
[x1 x2 x3 x4]=solve(eq1,eq2,eq3,eq4,'x1','x2','x3','x4')
站长推荐 /2
为其他会员撰写推荐,表示您认可他们对于论坛的出色贡献;
向其他坛友请求推荐,令您的声望页面更有说服力。
赶快学习新功能怎么用吧,推荐写起来!~
Powered by求高手matlab非线性方程组解法!!1)E=A*normcdf(d1,0,1)-B*exp(-r*t)*normcdf(d2,0,1)2) sigE=A*sigA*normcdf(d1,0,1)/Ed1=(log(A/B)+(r+0.5*sigA^2)*t)/(sigA*sqrt(t))d2=d1-sigA*sqrt(t)其中E、sigE、B、t、r已知数值需要求A、sigA的_百度作业帮
求高手matlab非线性方程组解法!!1)E=A*normcdf(d1,0,1)-B*exp(-r*t)*normcdf(d2,0,1)2) sigE=A*sigA*normcdf(d1,0,1)/Ed1=(log(A/B)+(r+0.5*sigA^2)*t)/(sigA*sqrt(t))d2=d1-sigA*sqrt(t)其中E、sigE、B、t、r已知数值需要求A、sigA的
求高手matlab非线性方程组解法!!1)E=A*normcdf(d1,0,1)-B*exp(-r*t)*normcdf(d2,0,1)2) sigE=A*sigA*normcdf(d1,0,1)/Ed1=(log(A/B)+(r+0.5*sigA^2)*t)/(sigA*sqrt(t))d2=d1-sigA*sqrt(t)其中E、sigE、B、t、r已知数值需要求A、sigA的数值解!请教高手!!解方程全部过程谢谢。。我不是学数学的所以一窍不通。。非常感谢!原始方程:
sigA(1)=2;设定迭代初值A(1)=2;设定迭代初值E=1;sigE=1;B=1;t=1;r=1;for&i=1:100设定迭代的步数d1=(log(A(i))-log(B)+(r+0.5*sigA(i)^2)*t)/(sigA(i)*sqrt(t));d2=d1-sigA(i)*sqrt(t)&;A(i+1)=(E+B*exp(-r*t)*normcdf(d2,0,1))&/normcdf(d1,0,1);sigA(i+1)=sigE*E/(A(i)*normcdf(d1,0,1));endx=1:101;subplot(1,2,1);plot(x(1,50:101),A(1,50:101));subplot(1,2,2);plot(x(1,50:101),sigA(1,50:101));画出图形检验是否是收敛的结果。如果图形震荡那么就不是收敛的结果要么改变迭代初值,要么说明这样的迭代解法有问题。在给定的初值和条件下的到的是收敛的结果。为:sigA=0.416A=1.101
原始方程组什么样的?求高手matlab非线性方程组解法!!1)E=A*normcdf(d1,0,1)-B*exp(-r*t)*normcdf(d2,0,1)2) sigE=A*sigA*normcdf(d1,0,1)/Ed1=(log(A/B)+(r+0.5*sigA^2)*t)/(sigA*sqrt(t))d2=d1-sigA*sqrt(t)其中E、sigE、B、t、r已知数值需要求A、sigA的_百度作业帮
求高手matlab非线性方程组解法!!1)E=A*normcdf(d1,0,1)-B*exp(-r*t)*normcdf(d2,0,1)2) sigE=A*sigA*normcdf(d1,0,1)/Ed1=(log(A/B)+(r+0.5*sigA^2)*t)/(sigA*sqrt(t))d2=d1-sigA*sqrt(t)其中E、sigE、B、t、r已知数值需要求A、sigA的
求高手matlab非线性方程组解法!!1)E=A*normcdf(d1,0,1)-B*exp(-r*t)*normcdf(d2,0,1)2) sigE=A*sigA*normcdf(d1,0,1)/Ed1=(log(A/B)+(r+0.5*sigA^2)*t)/(sigA*sqrt(t))d2=d1-sigA*sqrt(t)其中E、sigE、B、t、r已知数值需要求A、sigA的数值解!请教高手!!解方程全部过程谢谢。。我不是学数学的所以一窍不通。。非常感谢!原始方程:
sigA(1)=2;设定迭代初值A(1)=2;设定迭代初值E=1;sigE=1;B=1;t=1;r=1;for&i=1:100设定迭代的步数d1=(log(A(i))-log(B)+(r+0.5*sigA(i)^2)*t)/(sigA(i)*sqrt(t));d2=d1-sigA(i)*sqrt(t)&;A(i+1)=(E+B*exp(-r*t)*normcdf(d2,0,1))&/normcdf(d1,0,1);sigA(i+1)=sigE*E/(A(i)*normcdf(d1,0,1));endx=1:101;subplot(1,2,1);plot(x(1,50:101),A(1,50:101));subplot(1,2,2);plot(x(1,50:101),sigA(1,50:101));画出图形检验是否是收敛的结果。如果图形震荡那么就不是收敛的结果要么改变迭代初值,要么说明这样的迭代解法有问题。在给定的初值和条件下的到的是收敛的结果。为:sigA=0.416A=1.101
原始方程组什么样的?matlab解对数方程组_百度知道
matlab解对数方程组
[x,y,z]=solve('0.29744*x/(3+0.1*x+y+z)=(log10((3+x+y+z)/(3+0.1*x+y+z)))^2','0.79*x/(3+0.1*y+z)=(log10((3+x+y+z)/(3+0.1*x+y+z)*(3+y+z)/(3+0.1*y+z)))^2','1.25*x/(3+0.1*z)=(log10((3+x+y+z)/(3+0.1*x+y+z)*(3+y+z)/(3+0.*1y+z)*(3+z)/(3+0.1*z)))^2') returns出现Error: Missing MATLAB operator.!!哪里有问题啊!!
提问者采纳
所有的等号换成减号
其他类似问题
为您推荐:
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁

我要回帖

更多关于 matlab解偏微分方程组 的文章

 

随机推荐