mathematica吧
关注: 19,487 贴子: 73,003

自带帮助是最好的教材,提问前请先仔细查查

  • 目录:
  • 程序设计及人工智能
  • 3
    Evaluate在我看来,只能把套在他外面的那个函数的Hold属性去掉 f[x_, y_] := x^2 + y^2 + Evaluate@a; f[x_, y_] := Evaluate[x^2 + y^2 + a]; 在左式中,套在Evaluate外面的函数是Plus而非SetDelayed,因而Evaluate就是废的——Plus没有HoldAttribute。右式中,Evaluate就使:=的HoldAllAttribute消失,使得:=变为了=,输出成功为8。
    江翌日 12-25
  • 3
    ClearAll["Global`*"] Solve[{(4 Subscript[U, dc])/(\[Pi] Sqrt[ 2]) == (1/(I Subscript[C, ft1] \[Omega]) + I Subscript[L, ft1] \[Omega]) Subscript[I, in1] - 1/(I \[Omega] Subscript[C, ft1]) Subscript[I, T1], 0 == -1/(I \[Omega] Subscript[C, ft1]) Subscript[I, in1] + (1/(I \[Omega] Subscript[C, T1]) + I Subscript[L, T1] \[Omega] + 1/( I \[Omega] Subscript[C, ft1]) + I \[Omega] Subscript[M, T1R1] Subscript[I, R1]) Subscript[I, T1], 0 == (1/(I \[Omega] Subscript[C, R1]) + I Subscript[L, R1] \[Omega] + Subscript[R, 1] - I \[Omega] Subscript[M, T1R1] Subscript[I, T1]) Subscript[I, R1], (4 S
  • 4
    Solve[{Subscript[T, 2] == s*(Subscript[T, 1] + d) + a, Subscript[T, 3] == s*Subscript[T, 3] + a, Subscript[T, 4] == Subscript[T, 3] + d, Subscript[T, 6] == s*(Subscript[T, 5] + d) + a}, {s, a}] 请问最上面这个方程,s和a我想只用T表达,在mathematica里怎么表述呢,直接用Solve,s和a解出来是数值解
  • 54
    1. 帮助文件安装很顺利,没有出现某网友遭遇的尴尬现象。 2. 但是在mathematica中查找某条命令的帮助文件时,会跳到Mathematica官网上的帮助网页。这样,似乎也就无所谓中文版和英文版的区别了。你上的官网是中文版的,网页上帮助文件就是中文版的。 3. 离线无法查找帮助。不知这算不算一个缺点?
  • 4
    首先我定义了S函数,我想对w变量进行数值积分 S[w_, l_, m_, mi_, ma_] := { so = NSolve[{r + 2 m*Log[r/(2 m) - 1] == mi}, r]; bd = r /. so; sol = NDSolve[{\[Phi]''[ s] + (w^2 - (R[s] - 2*m)/ R[s]*((l*(l + 1))/R[s]^2 + (2*m)/R[s]^3)) \[Phi][s] == 0, R'[s] == 1 - (2*m)/R[s], \[Phi][mi] == 1, \[Phi]'[mi] == -I*w, R[mi] == bd[[1]]}, {\[Phi], R}, {s, mi, ma}]; Pfun = \[Phi] /. sol[[1, 1]]; Rfun = R /. sol[[1, 2]]; Abs[Exp[-2*I*w*ma]*(I*w*Pfun[ma] + Pfun'[ma])/( I*w*Pfun[ma] - Pfun'[ma])] } NIntegrate[S[x, 0, 1, -20, 20], {x, 0, 1}] 但是我发现NIntegrate
    XzqQAQ 12-19
  • 2
    2.` (0.5` a - 0.5` b - 0.425` E^(-3.` \[Theta]) - 0.7071067811865476` Sqrt[k] \[Theta]) - 1/20 (20 a - 20 b - 17 E^(-3 \[Theta]) - 20 Sqrt[2 b - 2 a b + b^2 + 2 k \[Theta]^2]) // FullSimplify 不知道如何比较,我将两个式子做差了,式子里面的变量的取值范围如下{\[Theta], 0, 1}, {a, 1, 10}, {b, 0, 1}, {k, 0, 1}。 求大佬帮助!!!
  • 3
    r = 0.04; v1 = 0.25; v2 = 0.1; sigma = Sqrt[v1 + v2]; NDSolve[{D[f[t, x], t] + 0.5 (v1 + v2) x^2 D[D[f[t, x], x], x] - r f[t, x] + r x D[f[t, x], x] == 0, f[t, 0.] == Exp[-r (1. - t)], f[t, 1.] - Derivative[0, 1][f][t, 1.] == 0., f[1., x] == 1. - x}, f, {t, 0, 1}, {x, 0, 1}]
    helloapple11 12-17
  • 1
    看似是个简单的事情,但我从帮助文件里面找不到一个函数能干这事。难道数{的个数吗? 比如{1,2,3}是0层,{{1,2},{3,4}}是1层。要是更多的层数,难道真的数花括号的个数吗?
  • 1
    我的wolfram id 已经在两台机器上安装了wolfram engine了,想在第三台机器上安装,前两台机器怎么退出wolffam id,然后在第三台机器注册呢
  • 7
    C1 = -G + (a^2 k (-1 + \[lbk]Gamma[rbk]))/( 2 (4 k \[lbk]Beta[rbk] (-1 + \[lbk]Gamma[rbk]) + \[lbk]Rho[rbk]^2))C2 = -((a^2 k (-1 + \[lbk]Alpha[rbk])^2 (-1 + \[lbk]Gamma[rbk]) (12 k \[lbk]Beta[rbk] (-1 + \\[lbk]Gamma[rbk]) - \[lbk]Rho[rbk]^2))/(2 (4 k \[lbk]Beta[rbk] (-1 + \[lbk]Gamma[rbk]) + \[lbk]Rho[rbk]^2)^2))C3 = C1 - C2 我想定义c1和c2,然后让他俩相加,但是结果却只是简单相加,没有合并运算,该怎么让他俩进行合并运算,如果用Simply的话也是效果不好,有没有大佬知道该怎么办呜呜呜
  • 5
    Plot[lbk]{6 - 6 x}, {x, 0, 10}, PlotRange -> {0, 10}, Frame -> True,PlotStyle -> {{Blue, Dotted}, {Orange, Dotted}, {Blue, Dotted}},PlotLegends -> {"函数1", "函数2", "函数3"}, PlotRangePadding -> 0,FrameLabel -> {{HoldForm[lbk]"左列标签"[rbk],None},{HoldForm[lbk] "下方标签"[rbk],None}},(*MeshStyle\[lbk]Rule[rbk]Directive[lbk]PointSize[lbk]Medium[rbk],\Red[rbk]*)Mesh -> 10[rbk](*Plot[lbk]x,(x,0,1),MeshStyle\[lbk]Rule[rbk]Directive[lbk]PointSize[lbk]Medium[rbk],Red[rbk],\Mesh\[lbk]Rule[rbk]10*)为什么相
  • 3
    代码: ParametricPlot[ Evaluate[{Sin[3 v], Cos[3 v]} + RotationTransform[v][ RotationTransform[u][{8 + Sin[3 u], 3 Cos[3 u]}]]], {u, 0, 2 Pi}, {v, 0, 2 Pi/3}, BoundaryStyle -> None] 结果: 以上应该是定宽曲线旋转平移得到的区域,请问像这种参数区域本身有重叠的情况能不能求出不重叠的区域方程、内外边界方程、内外边界周长、区域面积等
  • 3
    解微分方程,输出的是输入的结果,算不出来,请问是怎么回事?我看网上有的说是下标问题我改了以后还是算不出来,好像是这个定积分的问题,请问该怎么解决? DSolve[f1[x] - 2 \[Lambda]*h*D[f1[x], {x, 2}] + 1/2*\[Lambda]*S*w^2*\!\( \*SubsuperscriptBox[\(\[Integral]\), \(0\), \(l\)]\(f1[ x] \[DifferentialD]x\)\) + \[Lambda]*Subscript[I, y]* D[f1[x], {x, 4}] == 0, f1[x], x] #微分方程求解#
    雪暖冬_ 12-14
  • 3
    明明N[]默认的就是16位精度,凭什么N[,$MachinePrecision] 就慢了500多倍?甚至连N[,5]都比N[]慢了300倍。为什么啊。 L = 50; mat = SparseArray[{Band[{2, 1}] -> RandomPrime[15, L - 1]/RandomPrime[40, L - 1], Band[{1, 2}] -> RandomPrime[15, L - 1]/RandomPrime[40, L - 1]}, {L, L}]; Eigensystem[N[mat]]; // AbsoluteTiming Eigensystem[N[mat, MachinePrecision]]; // AbsoluteTiming Eigensystem[N[mat, $MachinePrecision]]; // AbsoluteTiming Eigensystem[N[mat, 5]]; // AbsoluteTiming
  • 2
    如图,首先约定了v1 v2都是四维向量,这时v1和v2的张量积应该是一个二阶张量,TensorRank也能正确给出2。但是这时候试图将v1⊗v2与给定的4×4矩阵相加,发现v1⊗v2被当做数而非4×4矩阵加在了每个分量上。应该如何让mma将这个当做正确的矩阵加法?
  • 3
    有人遇到过这种问题,程序运行后,在重新打开,所有的代码前面都会添加上"In[1] / Out[1]"之类的字样,而且所有的Cell会合并成一个Cell。如图所示:第一张第二张分别是运行后,和关闭重新打开后的
  • 11
    图一是我想运行的代码,后几张是论文的公式,我的最终目的是想通过这个代码绘制三维图像,就像最后一张图那样
  • 8
    研究了MMA许久了(大约3个月),看了多本书。无奈人老眼花,多个问题无法解决。(对于高手应当不难) 下面程序(意图为每个学生出一道随机题目,并计算结果),有以下疑惑请高手解决: 1,list里面字符串如何对应位置相加?(当前程序,应当是在相乘 ) 2,list里面的公式如何计算它的值? 3,MMA里面的求极限符号,用字符串表达时,位置会混乱。(如何理顺位置?) ClearAll; xuehao = {"01", "02", "03", "04"} xingming = {"
    zglazylazy 12-8
  • 6
    各位大佬mathematica中能否定义一个符号变量是实数啊,类似matlab中的syms a real? 我尝试输入a∈Reals b∈Reals c∈Reals 最后求Re[a*I + b + c*I]得出的还是-Im[a] - Im[c] + Re[b]
    哲e气象 12-7
  • 0
    Solve[(a - \[Tau] - c)^2/4 - (k \[Theta]^2)/2 - (1 - \[Tau] - c)/2*B == 0, c] /. B -> \!\(TraditionalForm\` TagBox[ FractionBox[ RowBox[{"G", "(", "c", ")"}], RowBox[{"g", "(", "c", ")"}]], "MathMLPresentationTag", AutoDelete->True]\) // FullSimplify 上面是具体公式,其中G(C)为概率分布函数,g(c)为概率密度函数。a和k为常数。这种情况下不假设具体的概率分布函数能接出来吗?
  • 1
    比如我有a(t)和b(t),还有关于a,b的函数x(a,b),我想要求出∂x/∂a在t处的取值,该怎么写呢?
    zzzaaasss 12-5
  • 35
    看到有吧友求MMA编程提速技巧,就想到这篇文章了。这篇是Wolfram Blog上的一篇文章http://blog.wolfram.com/2011/12/07/10-tips-for-writing-fast-math,之后在网上看到有人将之翻译过来了,就转载过来给各位吧友看看吧~此翻译文章的原址是http://www.cnblogs.com/mymma/archive/2012/08/27/2659007.html 当人们跟我抱怨说Mathematica不够快时,我一般会叫他们去仔细瞧瞧那些令人恼怒的代码.跑得慢的原因通常不是Mathematica的性能问题,而是使用方式不太合适.我想我应该分享一些我在优化Mathema
    liyajx 12-4
  • 5
    先贴代码。 Subscript[Q, crit] = 4 1/6 (1 + 2 \[Zeta]) CC; ParametricPlot[{{(x^2 + 3 x^4 - (QQ/(4 CC))^2 + 2 x^2 \[Zeta])/( 4 \[Pi] R x^3) /. R -> 1 /. CC -> 1 /. \[Zeta] -> 0 /. QQ -> 0.3, (CC (-x^4 + 3 (QQ/(4 CC))^2 + x^2 (1 + 2 \[Zeta])))/( R x) /. R -> 1 /. CC -> 1 /. \[Zeta] -> 0 /. QQ -> 0.3}, {(x^2 + 3 x^4 - (QQ/(4 CC))^2 + 2 x^2 \[Zeta])/( 4 \[Pi] R x^3) /. R -> 1 /. QQ -> Subscript[Q, crit] /. CC -> 1 /. \[Zeta] -> 0, (CC (-x^4 + 3 (QQ/(4 CC))^2 + x^2 (1 + 2 \[Zeta])))/(R x) /. R -> 1 /. QQ -> Subscript[Q, crit] /. CC -> 1 /. \[Zeta] -&
    liyajx 12-2
  • 3
    求解微分方程 DSolve[2*k[x]*Derivative[2][k][x] - Derivative[1][k][x]^2 == 0, k[x], x] 得到一个通解 {k[x]->x*c[1]+x^2*c[1]^2/(4*c[2])+c[2]}. 可是 k[x]->x^2 的确是方程的一个特解,为什么不包含在通解中呢?
    liyajx 12-2
  • 1
    关于带误差的实验数据,根据模型进行参数拟合,如果模型参数是线性的(model1)就可以得到结果,如果参数是非线性的(model2)就会报错。 求解啊求解 代码如下: data={{1,Around[5.,{0.3829020540280492,0.9823444252245808}]},{2,Around[7.,{0.13859288931953695,0.5062829286140218}]},{3,Around[9.,{0.6690553792750409,0.5611636477117092}]},{4,Around[11.,{0.09811685036353879,0.3115323565071908}]},{5,Around[13.,{0.8971600591178541,0.3276282509045023}]},{6,Around[15.,{0.10957197596386714,0.4900422248350298}]}}; model1=a x+b; model2=a^2 x+b; fit1=Nonlinear
    xzcyr 12-2
  • 2
    Solve[(\!\( \*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(n\)]\(p\ a[i]\)\)) == 0, p, Reals] 但是下面这个把变量提取出来就可以求解? Solve[p (\!\( \*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(n\)]\ \(a[i]\)\)) == 0, p, Reals] 附图如下,我想求解的稍微复杂一些,例如第一种,不把变量提取到外面。那么第一种为什么无法求解呢?
  • 4
    如图,我将一个3维切片图给投影了,但是中间的被投影切片图还在。我想把它去掉,只保留周围的3个投影图,如何操作呀 w = 0.95 \[Theta]0 = Pi/2 m = 1 EOx[r_, \[Theta]_, t_] = (Cos[m*\[Theta] + \[Theta]0]) (1 + I)*(r*Cos[\[Theta]] + t)* Exp[-((r^2 + t^2)/w^2)] EOy[r_, \[Theta]_, t_] = (Sin[m*\[Theta] + \[Theta]0]) (1 + I)*(r*Cos[\[Theta]] + t)* Exp[-((r^2 + t^2)/w^2)] CEOx[r_, \[Theta]_, t_] = Conjugate[EOx[r, \[Theta], t]] CEOy[r_, \[Theta]_, t_] = Conjugate[EOy[r, \[Theta], t]] r = Sqrt[x^2 + y^2] \[Theta] = ArcTan[x, y] INx[x_, y_, t_] = EOx[r, \[
    HXZ2D 12-2
  • 2
    我在学用mathematia推导公式,现在我要把一些在笛卡尔坐标系下的表达式转换到一个自定义的正交曲线坐标系下去表示,请问这个可以如何实现。我在帮助里好像只有球坐标柱坐标之类现成的转换,能够自定义正交曲线坐标系吗?
    xzcyr 12-2
  • 2
    向各位请教下,在绘图的时候,我给出一个x的值,想读取对应一点的y值,应该如何操作? 或者我只想绘制一个十字线,而不是很多线怎么操作? 比如: 我想读取速度为35时候的y值 Manipulate[ ContourPlot[ pitch == (1000*v)/(rpm*25.4/60), {v, 20, 40}, {pitch, 8, 17}, AspectRatio -> 0.3, PlotRange -> {{20, 40}, {8, 16}}], {{rpm, 6000}, 5000, 7000},GridLines -> {Range[20, 40, 2], Range[8,16, 0.5]] 在Mathematica中,可以使用GridLines选项来设置绘图中的网格线。以下是一个简单的例子: Plot[Sin[x], {x,
    xzcyr 12-2
  • 4
    Clear[P, S1, S2, \[Rho]]; c = 3*10^5; neff = 1.467; n = 1.444; Lext = 5/10^5; gB = 2.4/10^18; Aeff = 8/10^17; k = 1.155/10^5; \[Lambda] = 1.55/10^9; Nu = N[c/\[Lambda], 20]; \[CapitalGamma]b = 10^8; v = 2.9824/10^3; Lr = 8; L = 20; \[Alpha] = 0.05; \[Eta] = 1.08/10^4; kb = 1.38/10^23; T = 300; \[Rho]0 = 2.65*10^15; v0 = 2.38; eqns = {(n/c)*D[P[Z, t], t] + D[P[Z, t], Z] == (-(gB/Aeff))*S2[Z, t]*\[Rho][Z, t] - \[Alpha]* P[Z, t], (n/c)*D[S2[Z, t], t] - D[S2[Z, t], Z] == (-\[Alpha])*S2[Z, t] + (gB/Aeff)*P[Z, t]* Conjugate[\[Rho][Z, t]] + \[Eta]*S1[Z, t] - \[Eta]* S2[Z, t], (n/c)*D[S1[Z, t], t] - D
    xzcyr 12-2
  • 3
    例如,在不知道向量维数的情况下, k*r.r-s.r==0&&r!=0向量, 解得 k==s.r/(r.r) 如果不行,有什么软件能做到呢
    xzcyr 12-2
  • 7
    我的思路是将直方图经ImageRotation旋转后,再通过GraphaGrid拼接在一起。 不知道是否存在内置的函数,能直接画出来这样的图?
    xzcyr 12-2
  • 5
    Subscript[\[Pi], dr][p_, w_, g_] = (p - w)*(\[Rho] a - b p + k g) Subscript[\[Pi], dm][p_, w_, g_] = w (\[Rho] a - b p + k g) + p [(1 - \[Rho]) a - b p + k g] - i g^2/2(*利润函数*) Subscript[Subscript[\[Pi], dm], 1][w_, g_] = [w (\[Rho] a - b p + k g) + p [(1 - \[Rho]) a - b p + k g] - i g^2/2] /. {p -> w + x} FullSimplify[D[Subscript[Subscript[\[Pi], dm], 1][w, g], w]]输出结果为什么是这个样子的呀
    xzcyr 12-2
  • 3
    RegionPlot[x^2 + y^3 < 2 && x + y < 1, {x, -2, 2}, {y, -2, 2}] 就是如何完全填充,各个方向不要留空隙,可以保留刻度线和下面的数字,学渣在线感谢解答!!
    xzcyr 12-2
  • 3
    如图: "`*" // Remove; t = 4 Module[{t}, t + 1] t Block[{t = 2}, (t)] t 运行结果是 4 1+t$2275 4 2 4 这么看Block和Module也差不多,对全局变量t的值,都是没有影响的。
    xzcyr 12-2
  • 3
    Z是一个积分,对上面的Π进行积分,F是数量,对F进行求导,解出F,然后代回Z中,得到的结果化简为什么是一个0呢?代码如下: Subscript[n, 0] = F Subscript[n, c][t] = Subscript[n, 0] e^(Subscript[b, c] t) Subscript[P, g][Subscript[Q, g]] = a - b Subscript[Q, g][F] Subscript[Q, g][F] = Subscript[Q, 0] + \[Lambda] *Subscript[n, c][t] Subscript[\[Pi], A][t] = Subscript[P, g][Subscript[Q, g]]*Subscript[Q, g][F] - Subscript[C, g] * Subscript[Q, g][F] Z = Integrate[Subscript[\[Pi], A][t], t] Solve[D[Z] == 0, F] F = (e^(-2 t Subscript[b, c]) (a e^(t Subscrip
    xzcyr 12-2
  • 6
    具体的代码如下: Reduce[(1 - 85/100 E^(-3 \[Theta]) - c)^2 - (2 \[Theta]^2)/11 - 2*(1 - 85/100 E^(-3 \[Theta]) - c)*c >= 0, c]
    xzcyr 12-2
  • 11
    如下图,令4张3D图照此顺序排列,且零平面相接。 用GraphicsGrid排列空白太大了零平面也不相连,用show直接在同一坐标上重叠了。 谢谢解答!
    xzcyr 12-2
  • 6
    解出来的东西不等号左边比右边大?? 代码如下 ToRadicals[ Reduce[-2916*a1^12 + (-432*a1^4 + 486*a1^6 + 432*a1^2*t)^2 >= 0 && -1 <= a1 <= 1 && -1 <= t <= 1, a1]]
    xzcyr 12-2
  • 1
    各位大神帮帮忙啊,我定义了个函数求含参求复共轭,发现对单项式有效,对多项式失效,求帮忙 ClearAll["`*"] conj[x_] := Refine[Conjugate[x], Assumptions -> Element[{v, kx, ky}, Reals]]; conj[I kx v] conj[I ky v] conj[I kx v + I ky v]
  • 1
    请问windows server上多用户运行mathematica,显示一直在运行,就是不出结果,最简单的1+1也是,有什么解决办法吗?
  • 2
    代码如下,sinx的导数为什么不显示为cosx In[lbk]31[rbk]:= D[lbk]sin[lbk]x[rbk], x[rbk] Out[lbk]31[rbk]= Derivative[lbk]1[rbk][lbk]sin[rbk][lbk]x[rbk]
    寒之歌 11-27
  • 4
    在实践中发现对此不甚明了的人意外地多,故开此主题。为方便修订,照例弄个云笔记,请多点一下鼠标: http://note.youdao.com/noteshare?id=f2cdd9e72b66c73f372296800de7
    ZX_10281028 11-30
  • 4
    现在想解一组方程,三个方程三个未知数相互独立,理想情况是得到六组(phi,x,y)的值,但是现在输出的是空白,求帮助!!!!!代码如下:Solve[{c2 == 15.91, c3 == 0.00, d3 == 10.00, L2 == 17.04, L3 == 20.84, theta == 50.57, p1 == 14.98, p2 == 15.38, p3 == 12, p1^2 == x^2 + y^2, p2^2 == (x + L2*cos (phi) - c2)^2 + (y + L2*sin (phi))^2, p3^2 == (x + L3*cos (phi + theta) - c3)^2 + (y + L3*sin (phi + theta) - d3)^2}, {x, y, phi}]
  • 3
    下式: (Fm R V)/(d (R + L s + Cs L R s^2) + Fm (1 + Fv R + Cs R s) V) 想化为如下形式,参考之前的帖子未找到答案,找官网的参考呀没找到答案,请求大佬帮忙,找了mathematica官网仍为找到答案:
  • 1
    题主在做如下积分时发现结果中总是出现复数项: Integrate[Sin[b (-2 h + z)]^2 Sinh[a z]^2 a^2 b^2 x^2, {z, 0, h}] 尝试
    冉qx 5-23
  • 1
    RT,为什么运算时会出现空括号呢,请大家解答一下,谢谢。 l = 550*10^-3; L = 125*10^-3; H = 50*10^-3; h0 = H + 50*10^-6; E1 = 3*10^-6; E2 = 5*10^-6; \[Lambda]1 = 90*10^-3; \[Lambda]2 = 320*10^-3; E3 = 0.000008; \[Lambda]3 = 0.3; \[CurlyPhi]1 = 0; \[CurlyPhi]2 = 0.3*\[Pi]; Z1 = E1*Sin[((2*\[Pi])/\[Lambda]1)*(X - L/2) + \[CurlyPhi]1] - h0/2 - E2*Sin[((2*\[Pi])/\[Lambda]2)*(X - L/2) + \[CurlyPhi]2] Z2 = E1*Sin[((2*\[Pi])/\[Lambda]1)*(X - L/2) + \[CurlyPhi]1] - h0/2 - E3*Sin[((2*\[Pi])/\[Lambda]3)*(X - L/3) + \[CurlyPhi]2]
    傻丶书生 11-20

  • 发贴红色标题
  • 显示红名
  • 签到六倍经验

赠送补签卡1张,获得[经验书购买权]

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频!

本吧信息 查看详情>>

小吧:小吧主共10

会员: Expr

目录: 程序设计及人工智能

友情贴吧