如何提取黑白照片中封闭曲线的数据（即，x和y坐标）？ 在stackexchange也发帖问了，但没有得到我想要的结果。 我想用FindGeometricTransform获得坐标变换。 下面是我的粗糙的代码： 下面的网址前加上：https: 我如果加上，就不能发帖，不知为何。 imgImport = Import["//i.stack.imgur.com/c0KCa.png"];img = ImageTake[imgImport, {55, 286}, {125, 469}]cor = ImageValuePositions[img, {0., 0., 0.}];{o, x, y} = {{6.56169201921564`, 22.11360816415005`}, {28.934949660676992`, 5.2325141690732515`}, {6.679185088841892`, 6.5546

代码如下: m = 512; freq = {1, 2^(2/12), 2^(4/12), 2^(5/12), 2^(7/12), 2^(9/12), 2^(11/12), 2}; Playmusic[lbk]song_[rbk] := Do[lbk]Play[lbk]Sin[lbk]2*Pi*m*t*freq[lbk][lbk]song[lbk][lbk]i[rbk][rbk][rbk][rbk][rbk], {t, 0, 0.8}, PlayRange -> {0, 1}[rbk], {i, 1, Length[lbk]song[rbk]}[rbk]; music = {1, 2, 3, 4, 5, 6, 7, 8, 7, 6, 5, 4, 3, 2, 1}; Playmusic[lbk]music[rbk] 所示代码来源于数学实验课本，但版本较旧，存在一些问题，实际运行并没有发出音乐，向各位大佬求助。

大佬们 我现在用UnsameQ这样 h[lbk]x_，y_[rbk] ：＝q[lbk]x，y[rbk] /； y＝！＝x 但是这样规定的话后面替换h[lbk]a，b[rbk]这样的就会输出成q[lbk]_，_[rbk]的形式 怎么能不改变h并规定y和x不相等啊

怎么用Cases，Pattern和FullForm把列表里面的这种分数提取出来，例：{1/x，2x/3，1/2，3，√2x/3y}里面分母带有未知数的分数怎么提取出来。看了半天不知道怎么弄

在算多项式的时候想到，如果我想将一个多项式使用（x+y）来合并同类项，就是合并之后的式子是a(x+y)^3+b(x+y)^2+c(x+y)这种形式，其中a，b，c也不一定是常数。

一直被吞，所以我发截图描述了

Subscript[\[CurlyEpsilon], \[ScriptR]] = 1 - (\[Omega] - Subscript[\[Omega], \[ScriptC]])/\[Omega]* \!\(\*SubsuperscriptBox[\(\[Omega]\), \(\[ScriptP]\), \(2\)]\)/((\ \[Omega] - Subscript[\[Omega], \[ScriptC]])^2 + \[Nu]^2) - I*\[Nu]/\[Omega]* \!\(\*SubsuperscriptBox[\(\[Omega]\), \(\[ScriptP]\), \(2\)]\)/((\ \[Omega] - Subscript[\[Omega], \[ScriptC]])^2 + \[Nu]^2); \ (*冷磁化等离子体相对介电常数*) Subscript[\[Eta], 2] = Sqrt[ Subscript[\[Mu], 0]/(Subscript[\[CurlyEpsilon], 0]* Subscript[\[CurlyEpsilon], \[ScriptR]])];(*媒质2中的波阻抗*) Subscript[\[ScriptK], 1] = \[Omega]*

Solve[lbk]((2^\[lbk]Delta[rbk]*\[lbk]Alpha[rbk]*k)/(1-\[lbk]Alpha[rbk]))^\[lbk]Alpha[rbk]-(1-\[lbk]Beta[rbk])*((k*\[lbk]Beta[rbk])/(1-\[lbk]Beta[rbk]))^\[lbk]Beta[rbk]==0,k[rbk] 请问为啥简单方程出来的解这么奇怪

rt setprecision要一个个设定太累了

求大佬指点

物理意义上是没问题的，但每次算一半就报错出现无限，咋办呢

，同时我试了 NetTrain 指定GPU训练，没问题

示例如下： (ParallelDo[ Sow@i, {i,1,10^5}]//Reap)[[2,1]]//First 此时Reap似乎收集不到数据

Solve[-w - \[Delta] + \[Tau] + qs (p + \[Alpha] - \[Tau]) g[qs] == 0, qs] Solve::nsmet: 无法利用 Solve 现有的方法求解该系统. 小白求问，这种式子该怎么解啊

求各位大佬们帮帮忙，下面是我写的代码： Kang = (210*10^9)/(1 - 0.09)*1/12*0.2*0.0032^3; rou = 7800; zukang = 1668*200; lamana = 0.1; Td = 0.004; P0 = 124234.645; aef = 0.8; l = 0.2; A = 0.2*0.0032; wn = Sqrt[Kang/(rou*A)]*((n*Pi)/l)^2; Wn = Sqrt[2/(rou*A*l)]*Sin[(n*Pi*0.1)/l]; kesin = zukang/(rou*A)/(2*wn); Dn1 = (1 - (-1)^n)*l/(n*Pi)*Sqrt[2/(rou*A*l)]*P0*aef; Dn2 = (-kesin + Sqrt[kesin^2 - 1])*wn; Dn3 = (-kesin - Sqrt[kesin^2 - 1])*wn; Dn4 = Sqrt[1 - kesin^2]*wn; En1 = -1*(Dn1*(2*kesin*Dn3 + wn))/(wn^4*lamana*Td*Sqrt[kesin^2 - 1]); Fn1 = (Dn1*(2*kesin*Dn2 + wn))/(wn^4*lamana

请问为什么这里清除了变量，还是会代入之前的值呢？（之前M赋值为6）

这是直接设定参数代码： \[Theta] = 90 Degree \[Phi] = 0 Degree a = 0.1 b = 0.1 ParametricPlot[{x = Cos[t] (1/ 2 (Cos[\[Phi]]^2/b^2 + (Cos[\[Phi]]^2 Csc[\[Theta]]^2)/a^2 + Sin[\[Phi]]^2/a^2 + (Csc[\[Theta]]^2 Sin[\[Phi]]^2)/b^2 - Sqrt[4 (-((Cos[\[Phi]] Csc[\[Theta]] Sin[\[Phi]])/ a^2) + (Cos[\[Phi]] Csc[\[Theta]] Sin[\[Phi]])/ b^2)^2 + (-(Cos[\[Phi]]^2/ b^2) + (Cos[\[Phi]]^2 Csc[\[Theta]]^2)/a^2 - Sin[\[Phi]]^2/a^2 + (Csc[\[Theta]]^2 Sin[\[Phi]]^2)/ b^2)^2]))^(-1/2), y = Sin[t] (1/ 2 (Cos[\[Phi]]^2/b^2 + (Cos[\[Phi]]^2 Csc[\[Theta]]^2)/a^2 + Sin[\[Phi]]^2/a^2 + (Csc[\[Theta]]^2 Sin[\[

其他的地方输入一个v-pn 也会这样现实无效

请教各位大神，网上等高线的范例： ContourPlot[ Evaluate[ BesselJ[1, r/2]*Cos[3 t ] /. {r -> Sqrt[x^2 + y^2], t -> ArcTan[x, y]}], {x, -5, 5}, {y, -5, 5}, PlotPoints -> 50, ColorFunction -> ColorData["DarkRainbow"], Contours -> 25]。 上图是案例，下面是我绘制的，左边多出的白线请问该如何让其消失。

请问大佬 怎么用mathematica画出跳棋棋盘？

如图所示，不知道错在哪里了 求助各位大神，小女在此跪谢了！！

……总之算是把自带的教程稀里糊涂啃了一遍，然后开始试验，可是第一个例子就碰壁了： Reap@Do[Sow[i++], {i, 10}] Reap@ParallelDo[Sow[i++], {i, 10}] 期望中的结果当然是两者一样，但是第2个只能获得一个空列表。这里的i好像也不存在要分配变量的问题。（实际上分配了之后也没变化……）有人知道怎么回事吗？

目的是基于mathematics仿真杨氏双缝干涉实验的干涉图样，以下是鄙人所写的代码： (*定义参量，单位均为m*) \[Lambda] = 632.8*10^-9; (*入射光波长*) d = 1*10^-3; (*双缝间距*) L = 1; (*双缝到屏的垂直距离*) (*定义干涉函数T*) Subscript[r, 1] = Sqrt[L^2 + (x - d/2)^2]; (*S1到观察点P的距离*) Subscript[r, 2] = Sqrt[L^2 + (x + d/2)^2]; (*S2到观察点P的距离*) \[Delta] = Subscript[r, 2] - Subscript[r, 1] = ( Subscript[r, 2]^2 - Subscript[r, 1]^2)/( Subscript[r, 2] + Subscript[r, 1]) = (2 xd)/( Sqrt[L^2 + (x - d/2)^2] + Sqrt[L^2 + (x +

这个其实也是说过多次的事情，索性开个主题，方便查找。吧里的老人应该都知道我只在公开的网站（贴吧/知道/知乎/Stackexchange…）回答Mathematica相关的问题。我也没有加入任何与Mathematica相关的qq群/微信群（确切地说我没微信）。这么做的原因主要有2个： 1. 封闭的讨论对知识的传承非常不利。 2. 我想尽可能地对抗“社交蒸发冷却效应”。 顺便，这套做法倒也不是一朝形成的。尽管我很早就意识到了封闭式讨论的缺点，但在早期，面对对方主动奉上

请教一个问题，比如0的α次方为什么输出不显示0而是显示0的α次方呢？

用“龙格库塔 mathematica”在百度翻了10页，只看到2个基于Mathematica的经典龙格库塔法实现。一个是一篇论文，用的是For循环。（顺便摘要里还说经典龙格库塔是常微分方程近似解法里“最常用的一种”，这句话即使放在论文发表的06年也是大有问题。）还有一个是我在知道上的回答，引用的是Stackexchange上使用MapThread和Rule的代码，这个答案如今看来写得也不好。 自编RK4也算是个常见的作业题，如果目前中文网络上的相关实现真的只有上面两个，那是不

代码如下： Subscript[h, 0] = 5; theta0 = Pi/4; Subscript[c, 0] = 1; Subscript[l, 0] = Subscript[h, 0]/Sin[theta0]; Subscript[l, 1] = Sqrt[((Subscript[l, 0]/2)^2) + (Subscript[c, 0]^2) ]; a = Subscript[h, 0] /Tan[theta0]; h = Subscript[h, 0] - x; theta01 = ArcTan[h/a]; Subscript[l, 01] = h/Sin[theta01]; thetamid = ArcCos[Subscript[l, 01]/(2*Subscript[l, 1] )]

请问如果定义矩阵或者向量中的元素都是复数，而且不是具体的数字，比如a+bi这种，我想对矩阵或者向量共轭转置，怎么才能让结果中不出现Conjugate[a]-iConjugate[b]这种形式，出现的是a-bi这种啊？

这个也算是月经问题了，但是吧里好像一直没有标题足够明了的帖子出现，故开此帖。百度目前彻底限死了外链，只好贴截图了：

我想作出一个函数的图像，但是运行后mma提示说”1/7^960.765. 太小，不适合被表示为归一化机器数；可能无法保持原来的精度“。请问这种情况该怎么解决？有办法提高精度让mma不再报错么？ 我写的代码如下： Plot[(27 2^((-2 + a)/(-1 + a)) 5^(a/(-1 + a)) 7^(1/(1 - a)) a^2)/((-1 + a) (4 + 9 a)), {a, 1.001, 3}] 运行结果如下： 请吧里高手赐教，谢谢！！

Solve[lbk]{-(3/2) Log[lbk]a[rbk]^3 Subscript[lbk]a, 1[rbk] Subscript[lbk]k, 1[rbk] + 27/2 Log[lbk]a[rbk]^3 Subscript[lbk]a, 2[rbk] Subscript[lbk]k, 1[rbk] + 9 Log[lbk]a[rbk] Subscript[lbk]a, 1[rbk] Subscript[lbk]a, 2[rbk] Subscript[lbk]k, 1[rbk] - 6 Log[lbk]a[rbk]\!\(\*SubsuperscriptBox[lbk]\(a\), \(2\), \(2\)[rbk]\) Subscript[lbk]k, 1[rbk] - 18 Log[lbk]a[rbk] Subscript[lbk]b, 1[rbk] Subscript[lbk]b, 2[rbk] Subscript[lbk]k, 1[rbk] + 12 Log[lbk]a[rbk]\!\(\*SubsuperscriptBox[lbk]\(b\), \(2\), \(2\)[rbk]\) Subscript[lbk]k, 1[rbk] == 0, 3 Log[lbk]a[rbk]^3 Subscript[lbk]a, 1[rbk] Subscript[lbk]k, 1

就如 @jiaoeyushushu 所说，现在是大会前夕，所以犹豫了很久要不要在这个时候开帖说这事，但想来想去还是要说几句，毕竟之前的几次讨论展现出的“风向”实在是让我有些，失望。 ……今天也快过了，这里只捡重点说吧。总之，大家还记不记得前段时间（好吧，快一年了……）的 如龙6 事件？（忘了的或者不知道的请自行搜索 “如龙6 台湾”。） 你们觉得为什么最后世嘉肯改？ 你们觉得国人买了多少世嘉的游戏？ Wolfram公司说到底也是家商用公司。

如图，现在想把图a和c合并在一起，但是两者纵坐标不同，按照一般的Show没法2个都正常显示， 下面分别是a和c图的代码， （a） J1 = 2; E0 = -3*B/2; E1 = (-B - J2)/2; E2 = (B - J2)/2; E3 = -B/2 + 1/4*(J2 - Sqrt[J2^2 + 8*J1^2]); E4 = B/2 + 1/4*(J2 - Sqrt[J2^2 + 8*J1^2]); E5 = -B/2 + 1/4*(J2 + Sqrt[J2^2 + 8*J1^2]); E6 = B/2 + 1/4*(J2 + Sqrt[J2^2 + 8*J1^2]); E7 = 3*B/2; a3 = a4 = (-J2 - Sqrt[J2^2 + 8*J1^2])/(2*J1); a5 = a6 = Sqrt[2 + a4^2]; N3 = N4 = Sqrt[2 + a3^2]; N5 = N6 = Sqrt[2 + a6^2]; T = 0.05; u = Exp[-E0/T] + ((a3^2)*Exp[-E3/T]/(N3^2)) + (a5^2)*Ex

请问一下这个节开启符号能不能默认开启，我在选项设置找了很长时间没发现这个选项，望好心人告知。

我想用Mathematica计算出下图的积分近似值，但总是无法做到，问题出在哪儿？ 代码NIntegrate[lbk]Abs[lbk]Sin[lbk]x^4[rbk]/(Sqrt[lbk]x[rbk]+x^2),{x,0,\[lbk]Infinity[rbk]}[rbk][rbk]

Table[(#^n) &, {n, 1, 3}] (*我想遍历n得到 {#^1,#^2,#^3} 的纯函数列表，但是这句话里的n却没被赋值。输出为：{#1^n&,#1^n&,#1^n&}*) (*我要怎么做才能从1到3遍历指数n得到 {#^1,#^2,#^3}？*) 然后我想到了以前在吧里咨询时学到的Evaluate Table[Evaluate[#^n] &, {n, 1, 3}] (*这样就能得到 {#1 &, #1^2 &, #1^3 &} 了*) 发出来记录下。

如题。在计算中针对相同的微分方程采用两种列写方法，得到的结果不同，如图所示： 为什么会有误差，应该怎么解决，请各位dalao解惑

我是想计算一个由A的特征向量组成的矩阵T，由T组成的对角矩阵，我使用for循环但是他一直报错，麻烦各位大佬帮我看看 fem[lbk]\[lbk]CurlyEpsilon[rbk]1_, \[lbk]CurlyEpsilon[rbk]2_, \[lbk]CurlyEpsilon[rbk]3_, \\[lbk]CurlyEpsilon[rbk]4_, \[lbk]CurlyEpsilon[rbk]5_, \[lbk]CurlyEpsilon[rbk]6_, \\[lbk]CurlyEpsilon[rbk]7_, t_[rbk] :=Module[lbk]{A, \[lbk]CapitalLambda[rbk], T, T1, TT, TT1, SICO, C\[lbk]CapitalLambda[rbk], S\[lbk]CapitalLambda[rbk]}, A = {{0, \[lbk]CurlyEpsilon[rbk]1, \[lbk]CurlyEpsilon[rbk]2, \[lbk]CurlyEpsilon[rbk]3, 0, 0}, {\[lbk]CurlyEpsilon

为什么在Mathematica中Series定义的多元函数展开方式与一般数学中定义的泰勒展开式不同？ Mathematica中，Series是对f的分量分别进行展开，如下图： 而数学中定义的多元函数泰勒展开式为 这样造成，同一个函数如Sin[x^2+y^2]，在{0,0}处展开时，这两种方式产生的结果不一样。 代码： x0 = {0, 0}; p = {x, y}; f[{x_, y_}] := Sin[x^2 + y^2]^2 df[{x_, y_}] = D[f[{x, y}], {{x, y}}]; ddf[{x_, y_}] = D[f[{x, y}], {{x, y}, 2}]; Mathematica中Series方式： In[32]:= Series[f[{x, y}], {x, 0, 2}, {y, 0, 2}] // Normal Out[32

\[Lambda]p0 = 0.788;(*Centerwavelength of pump laser*) \[Lambda]0 = 2*\[Lambda]p0;(*Centerwavelength of laser*) c = 2.997*10^14;(*Speed of light*) vp = c/\[Lambda]p0;(* Pump freq.*) wp = 2*\[Pi]*vp;(* Angular freq.of pump*) z0 = 20*L;(* Depth of focus for pump*) A = Sqrt[z0*\[Lambda]p0/\[Pi]];(* Waist size for pump*) (*dlp=5*10^-2;*)(* bandwidth for pump*) dlp = 2.2*10^-2;(* bandwidth for pump*) dvp = (dlp*c)/\[Lambda]p0^2; dwp = 2*\[Pi]*dvp; v = c/\[Lambda]0;(* Center freq.of laser*) w = 2*\[Pi]*v;(* Angular freq.of laser*) L = 20*10^3;(* Crystal length*) (*ordinary refractive index is n2 fr

前面的式子计算有点长，然后系统一直显示正在运行中，然后就一直卡住运行，重启了也没有用，请问有没有什么解决方法？

最近在做一个东西，需要用到多元泰勒展开~然后现在我查了Series的用法，发现它只能进行一元泰勒展开。

\[Lambda]p0 = 0.788;(*Centerwavelength of pump laser*) \[Lambda]0 = 2*\[Lambda]p0;(*Centerwavelength of laser*) c = 2.997*10^14;(*Speed of light*) vp = c/\[Lambda]p0;(* Pump freq.*) wp = 2*\[Pi]*vp;(* Angular freq.of pump*) z0 = 20*L;(* Depth of focus for pump*) A = Sqrt[z0*\[Lambda]p0/\[Pi]];(* Waist size for pump*) (*dlp=5*10^-2;*)(* bandwidth for pump*) dlp = 2.2*10^-2;(* bandwidth for pump*) dvp = (dlp*c)/\[Lambda]p0^2; dwp = 2*\[Pi]*dvp; v = c/\[Lambda]0;(* Center freq.of laser*) w = 2*\[Pi]*v;(* Angular freq.of laser*) L = 20*10^3;(* Crystal length*) (*ordinary refractive index is n2 fr

Graphics[{PatternFilling[{"Checkerboard", Red}], Rectangle[]}] U2 = RegionPlot[\[Pi]rYR < \[Pi]rNR && \[Pi]rNR > \[Pi]rYA && \ \[Pi]rNR > \[Pi]rNA && \[Pi]rNR > 0 && P2NR > 0 && P2YR > 0 && P1NR > 0 && P1YR > 0 && P1YR > W1YR && P1NR > W1NR && P2YR > W2YR && P2NR > W2NR && X1YR > 0 && X2YR > 0 && X1NR > 0 && X2NR > 0, {q, 0, 2}, {\[Eta], 1, 3}, PlotStyle -> Directive[Opacity[0.5], Black], FrameLabel -> {q, \[Eta]

比如有 a={{1,2},{11,22},{111,222}} a[[0,25]] (*会报：Part::partd: 部分指定 {{1,2},{11,22},{111,222}}[[0,25]] 比对象深度更长.*) a[[3,3]] (*会报：Part::partw: {111,222} 的部分 3 不存在.*) 有没有办法对于未知形状的列表，在运行 a[[3,3]] 之前就能知道 3,3 这个位置有没有元素？

Reduce[{\[Alpha]*(( A^2 b - 2 A b (-1 + cr) \[Eta] + \[Eta] ((1 + b (-2 + f))^2 + b (-1 + cr)^2 \[Eta] - 4 b \[Beta]1^2 \[Theta]))/( 4 b \[Eta])) - (( A^2*b + \[Eta] + \[Eta]*b^2*(-2 + f) + \[Eta]^2* b*(-1 + cr)^2 + \[Beta]1^2*\[Theta]^2*b + 2 A*b*(\[Eta] - \[Eta]*cr + \[Beta]1*\[Theta]) + 2*\[Eta]* b*(-2 + f + \[Beta]1*\[Theta] - \[Beta]1*\[Theta]*cr))/( 16 \[Eta]*b)) >= 0; (1 - \[Alpha])*(( A^2 b - 2 A b (-1 + cr) \[Eta] + \[Eta] ((1 + b (-2 + f))^2 + b (-1 + cr)^2 \[Eta] - 4 b \[Beta]1^2 \[Theta]))/( 4 b \[Eta])) - (( A^2*b + \[Eta] + \[Eta]*b^2*(-2 + f) + \[Eta]^2*b*(-1 + cr)^2 + 9 \[Be