mathematica中解方程组并将结果画出图来
来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/05/08 18:37:45
mathematica中解方程组并将结果画出图来
代码如下
Q = 0
\[Mu] = 0
sol = Solve[
y[a,k,mu]*(a^2 *(1 + r) + r^2 *(-3 + 2 *Q^2/r + r)) ==
a *mu* (1 - r) +
k* Q* (a^2 + Q^2 - r) && (y[a,k,mu]*(r^2 + (a)^2) - mu*(a) -
k*Q*r)^2 == (r^2 - 2*r + a^2 +
Q^2)*(1 + ((a)^2)/2*(1 + mu^2)*y[a,k,mu]^2 -
2*(a)*mu*y[a,k,mu] + (\[Mu]^2)*r^2) &&
r \[GreaterSlantEqual] (1 + (1 - (a^2 + Q^2))^(1/2)),{y[a,k,
mu],r},Reals];
T1 = ParametricPlot[Evaluate[y[1,0,mu] /.sol],{mu,-1,1},
PlotRange -> {0.1,0.5},Frame -> True,AspectRatio -> 1/1]
其中,y=y[a,k,mu].求解出y后画出y随mu变化的图像 a可以取若干个数值,k也可以去若干个数值.但是我的不出图,求大神给看看.
代码如下
Q = 0
\[Mu] = 0
sol = Solve[
y[a,k,mu]*(a^2 *(1 + r) + r^2 *(-3 + 2 *Q^2/r + r)) ==
a *mu* (1 - r) +
k* Q* (a^2 + Q^2 - r) && (y[a,k,mu]*(r^2 + (a)^2) - mu*(a) -
k*Q*r)^2 == (r^2 - 2*r + a^2 +
Q^2)*(1 + ((a)^2)/2*(1 + mu^2)*y[a,k,mu]^2 -
2*(a)*mu*y[a,k,mu] + (\[Mu]^2)*r^2) &&
r \[GreaterSlantEqual] (1 + (1 - (a^2 + Q^2))^(1/2)),{y[a,k,
mu],r},Reals];
T1 = ParametricPlot[Evaluate[y[1,0,mu] /.sol],{mu,-1,1},
PlotRange -> {0.1,0.5},Frame -> True,AspectRatio -> 1/1]
其中,y=y[a,k,mu].求解出y后画出y随mu变化的图像 a可以取若干个数值,k也可以去若干个数值.但是我的不出图,求大神给看看.
可以出的,
再问: 如果我想依次 画出y[1,0,mu] y[0.99,0,mu] y[0.9,0,mu] y[1,1,mu] y[1,2,mu]怎么弄?另外,要是把您画的这个的直线部分的下部分抹掉怎么弄?如果分开弄的话,我知道怎么弄。另外,您画的这个图中直线部分是r=1的时候,亦即方程组中r>=(1 + (1 - (a^2 + Q^2))^(1/2))中的等号部分。求教。
再答: 不好意思,关于直线部分其实也是一个解,ta在列表首First可以获取,Plot[Evaluate[y[1, 0, mu] /. Take[sol,{2,Length@sol}]可以不画这个解;关于你的Solve求解很慢,我现在没想到好办法把它做成函数。
再问: 如果我想依次 画出y[1,0,mu] y[0.99,0,mu] y[0.9,0,mu] y[1,1,mu] y[1,2,mu]怎么弄?另外,要是把您画的这个的直线部分的下部分抹掉怎么弄?如果分开弄的话,我知道怎么弄。另外,您画的这个图中直线部分是r=1的时候,亦即方程组中r>=(1 + (1 - (a^2 + Q^2))^(1/2))中的等号部分。求教。
再答: 不好意思,关于直线部分其实也是一个解,ta在列表首First可以获取,Plot[Evaluate[y[1, 0, mu] /. Take[sol,{2,Length@sol}]可以不画这个解;关于你的Solve求解很慢,我现在没想到好办法把它做成函数。
mathematica中解方程组并将结果画出图来
MATHEMATICA中如何解方程组
mathematica中如何解未知变量方程组
怎么将Mathematica解方程组的解应用到后续的运算中啊?
mathematica解微分方程组
mathematica 解微分方程组
mathematica 解方程组
Mathematica中如何画出映射的函数图
matlab solve 解含参方程组,并画出点来!
Mathematica解四元二次方程组,结果太大了,怎么回事?
mathematica怎么将循环的每个结果放入文件中,打算这样用origin画图.
Mathematica中,如何将FindRoot[]的结果转化成数值形式?