-
3Evaluate在我看来,只能把套在他外面的那个函数的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。
-
3ClearAll["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
-
541. 帮助文件安装很顺利,没有出现某网友遭遇的尴尬现象。 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
-
22.` (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}。 求大佬帮助!!!
-
3r = 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}]
-
1看似是个简单的事情,但我从帮助文件里面找不到一个函数能干这事。难道数{的个数吗? 比如{1,2,3}是0层,{{1,2},{3,4}}是1层。要是更多的层数,难道真的数花括号的个数吗?
-
1我的wolfram id 已经在两台机器上安装了wolfram engine了,想在第三台机器上安装,前两台机器怎么退出wolffam id,然后在第三台机器注册呢
-
7C1 = -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的话也是效果不好,有没有大佬知道该怎么办呜呜呜
-
5Plot[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
-
3
-
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
-
3
-
11
-
8
-
6
-
0Solve[(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处的取值,该怎么写呢?
-
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
-
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] -&
-
3
-
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
-
2
-
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, \[
-
2我在学用mathematia推导公式,现在我要把一些在笛卡尔坐标系下的表达式转换到一个自定义的正交曲线坐标系下去表示,请问这个可以如何实现。我在帮助里好像只有球坐标柱坐标之类现成的转换,能够自定义正交曲线坐标系吗?
-
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,
-
4Clear[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
-
3例如,在不知道向量维数的情况下, k*r.r-s.r==0&&r!=0向量, 解得 k==s.r/(r.r) 如果不行,有什么软件能做到呢
-
7
-
5Subscript[\[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]]输出结果为什么是这个样子的呀
-
3
-
3如图: "`*" // Remove; t = 4 Module[{t}, t + 1] t Block[{t = 2}, (t)] t 运行结果是 4 1+t$2275 4 2 4 这么看Block和Module也差不多,对全局变量t的值,都是没有影响的。
-
3Z是一个积分,对上面的Π进行积分,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
-
6
-
11
-
6
-
1
-
1请问windows server上多用户运行mathematica,显示一直在运行,就是不出结果,最简单的1+1也是,有什么解决办法吗?
-
2
-
4在实践中发现对此不甚明了的人意外地多,故开此主题。为方便修订,照例弄个云笔记,请多点一下鼠标: http://note.youdao.com/noteshare?id=f2cdd9e72b66c73f372296800de7
-
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}] 尝试
-
1RT,为什么运算时会出现空括号呢,请大家解答一下,谢谢。 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]
扫二维码下载贴吧客户端
下载贴吧APP
看高清直播、视频!
看高清直播、视频!