-
-
3我想将向量场 r sin(\[Theta])Overscript[\[CurlyPhi], ^] 从球坐标 {r,\[Theta],\[CurlyPhi]} 转化为直角坐标 {x,y,z},并且保持场点仍以球坐标系表达,而场量以直角坐标系表达: 我可以分两步编程实现,代码如下: In[11]:= TransformedField[ "Cartesian" -> "Spherical", #, {x, y, z} -> {r, \[Theta], \[CurlyPhi]}] & /@ TransformedField[ "Spherical" -> "Cartesian", {0, 0, r Sin[\[Theta]]}, {r, \[Theta], \[CurlyPhi]} -> {x, y, z}] Out[11]= {-r Sin[\[Theta]] Sin[\[CurlyPhi]], r Cos[\[CurlyPhi]] Si
-
214.3版本,使用的是fritx5
-
10今天早上打开 mma 才发现 linux 下的 mma 不能输入中文(亏我还装了快一个月...),折腾了大半天才终于弄好了 下面讲讲解决方案
-
0
-
3贴吧客户端相较于传统的pc网页版弊端颇多(如中括号无法在顶楼正常显示,置顶只能显示一条等),但吧里始终有一小部分用户(访问数约占7%)在使用客户端,故开此帖调查。视调查结果可能会在下个答题日出篇教程。
-
12考虑到此类问题和“Mathematica已经装好了,但不能用”即“Mathematica一开始就不能用”有一定区别,特开专帖讨论。今后此类问题原则上请不要单开主题。 这里的“不能用”,包括不限于 1. Mathematica突然无法正常打开 2. Mathematica突然变得无法正常工作(1+1也算不了等) 3. Mathematica突然在执行正确的代码时出现莫名其妙的警告。 如果你的“Mathematica一开始就不能用”,请看置顶。 闲话少叙,说一说可能的修复方法: 此类问题似乎和Mathematica的Paclet更新不
-
6大家有用过Wolfram notebook assistant+LLM Kit的吗,感觉怎么样。我刚订阅了一个月的基础版试了一下,感觉好像有点鸡肋,我在vscode里写mma用copilot似乎也能达到差不多的效果。是我不会用吗 还有assistant和用其他模型api的chat相比有什么优势吗,以前没用过chat。
-
1被豆包忽悠下了MathpixSnippingTool,但是效果不是很好,容易乱码,符号也容易变成代码 有没有什么好的软件啊?各位大佬
-
6
-
2Mathematica有类似matlab的变量查看器吗? 比如用久了之后,我也不清楚我每个变量值,类型了,有类似matlab的查看器,直接查看现在有哪些变量存在,是
-
10
-
8求助,可能是哪里出问题了、 输入: f1[b_] = GyromagneticRatio["1H"] (1 - δ₁ 10⁻⁶) b; f2[b_] = GyromagneticRatio["1H"] (1 - δ₂ 10⁻⁶) b; f3[b_] = GyromagneticRatio["1H"] (1 - δ₃ 10⁻⁶) b; ave[b_] := (f1[b] + f2[b] + f3[b]) / 3.; Ham0[b_] = (f1[b] - ave[b]) opI[1, "z"] + (f2[b] - ave[b]) opI[2, "z"] + (f3[b] - ave[b]) opI[3, "z"] + 2 π (J12 opI[1].opI[2] + J23 opI[1].opI[3]); {ev1, vecs1} = CorrelatedEigensystems[Ham0[b] /. {J12 -> 25., J23 -> 5, δ₁ -> +0, δ₂ -> 10, δ₃ -> -10},
-
13比如说,我已经在一个.nb文件中写了一个较长的计算过程,这个过程的输入条件是a,b,c这三个变量,输出结果是d,e,一个曲线图。但是我要重用这个过程算很多种工况。也就是说 a=1,b=2,c=3 a=2,b=3,c=4 .... 由于这个过程较长较复杂,我把这个过程再写成一个函数的话太麻烦。 那么,能否另写一个nb文件,在其中设定 a=1,b=2,c=3 a=2,b=3,c=4 .... 我运动这个新的nb文件后,在此nb文件中出现对应的计算结果,即 a=1,b=2,c=3,d,e,一个曲线图 a=2,b=3,c=4,d,e,一个曲线图
-
6感觉还是旧版本的稳定,Mathematica14.3版本已经遇到很多次闪退了!Mathematica是我电脑最常用的软件,基本是开机后第一个打开的软件,我用神经网络算法比较多,经常出现奔溃、闪退,以前用旧版的时候记得从来没出现过。
-
5
-
8我找到Mathprogramming-intro.org这个网站,但并没有在其中找到下载链接,该网站是已经关闭了吗?
-
0嗯……众所周知(好吧,可能只有吧里的老人知道)Leonid Shifrin所著的《Mathematica Programming: An Advanced Introduction》(《Mathematica编程:高级导论》)是学习Mathematica核心语言的重要参考资料之一,尽管个人觉得随着时间的推移——尤其是在老沃自己写的《An Elementary Introduction to the Wolfram Language》(《Wolfram 语言基础入门》)成书之后——其重要性在逐渐下降,它依旧值得一读。本书的原版及(未完成但够看的)翻译版的地址在近期都出现了一些问题,因此
-
23先给各位吧友拜个晚年,各位新年好。
-
3"CUDALink使 Wolfram 语言能够在GPU上使用 CUDA 并行计算架构". 我想要将复数的 List 转换为 CUDAVector, 再使用 CUDADot 计算内积. 但是发现有 CUDATranspose 函数但是没有相应的 CUDA 下的 Conjugate 函数, 因此想要计算内积似乎需要先将 CUDAVector 转换回 List, 再正常地 Conjugate, 最后转换回 CUDAVector. 这个做法有些奇怪, 看上去是绕了一大圈远路, 有谁知道为CUDAVector取复共轭有什么简便的方法吗? Needs["CUDALink`"]; cuvec1 = CUDAVector[Normalize@RandomComplex[{-1 - I, 1 + I}, 10], &
-
6某些工程问题分析时,让y轴向下为正方向,更便于分析问题,比如小球下落问题。那么怎么在plot中打印出这种分析数据呢?我知道可以对数据进行处理,最简单的就是乘以-1,然后再到plot里面画图。
-
4ClearAll["Global`*"] PoissonBracket[f_,g_,coords_List,moms_List]:=Total[MapThread[D[f,#1]*D[g,#2]-D[f,#2]*D[g,#1]&,{coords,moms}]]; q={qx,qy,qz}; p={px,py,pz}; A={Ax[qx,qy,qz],Ay[qx,qy,qz],Az[qx,qy,qz]}; v=(1/m)*(p-e*A); H=1/2*m*(v.v)+e*\[Phi][qx,qy,qz]; rules={D[Az,qy]-D[Ay,qz]->Bx, D[Ax,qz]-D[Az,qx]->By, D[Ay,qx]-D[Ax,qy]->Bz, D[Ay,qz]-D[Az,qy]->-Bx, D[Az,qx]-D[Ax,qz]->-By, D[Ax,qy]-D[Ay,qx]->-Bz}; 我在计算带电粒子在电磁场中运动的泊松括号时,定义了磁矢势与磁场的替换规则,但是在实际计算的时候Mathematica并没有运
-
24近段时间我逐步增加了使用版本11.3和版本12的时间,结果发现,Mathematica近几版以来的某些改变(有的好有的坏)其实并没有写在文档里,故开此帖与大家交流。 先说我发现的几个好的改变: 1. 并行能力比过去强了。比如NDSolve在解偏微分方程时,如果调用的是TensorProductGrid方法,多核机器的CPU会跑满。(版本9时代其实也不是完全没并行,但不会像现在这么高效。)NonlinearModelFit也出现了CPU跑满的情况,不过我在版本9时代对NonlinearModelFit的使用有限,
-
0提问时请附上引发问题的代码。对标题论断有疑虑的同学请阅读下面这篇对 FindRoot::nlnum 警告的简要解说: https://zhidao.baidu.com/question/810397720290271692.html
-
4在实践中发现对此不甚明了的人意外地多,故开此主题。为方便修订,照例弄个云笔记,请多点一下鼠标: http://note.youdao.com/noteshare?id=f2cdd9e72b66c73f372296800de7
-
2也差不多该针对语义式匹配(Semantic Pattern Matching)专开个帖子了。为方便修订,照例开个云笔记,请各位多点一下鼠标: http://note.youdao.com/noteshare?id=77b86cb08aaf29ce8747990c9e71 有什么意见或建议欢迎提出。
-
10data = {{43.68, 1.}, {206.42, 0.76}, {398.11, 0.58}, {1019.47, 0.42}, {1910.44, 0.33}, {2964.39, 0.26}, {4116.84, 0.22}, {5318.17, 0.19}, {6505.9, 0.16}, {7709.48, 0.15}, {8827.55, 0.13}, {9984.18, 0.12}, {11015.08, 0.11}}; model = Sum[ 4/((2*k + 1)*Pi)*1(*c0*)* MittagLefflerE[\[Beta], ((-d)*(2*k + 1)^2*\[Pi]^2*t^\[Beta])]* Sin[(2*k + 1)*Pi*1/2], {k, 0, 10}]; M = NonlinearModelFit[data, model, {{\[Beta], 0.0029}, {d, 0.71}}, t]; General::ovfl: Overflow occurred in computation. General::ovfl: Overflow occurred in computation. General::ovfl: Overflow occurred in computation. General::stop: Furth
-
2
-
0症状: 顶楼中的 [ 和 ] 会被自动转换为 [lbk] 和 [rbk]。此现象在回复帖中不存在。 影响范围: 此bug至少在去年10月16日即已出现,至今(安卓12.58.1.0)未修。IPhone受影响情况不明,欢迎实验补充。 对策: 对于提问者,最简单的对策就是改用电脑访问贴吧网页版,不要使用使用手机客户端发帖。如果确有相关需求(所以说你们到底为什么会有这种需求?电脑没连网?那也可以开手机热点啊,访问网页又花不了什么流量),避免在顶楼张贴代码也是一个
-
3看 史蒂芬大佬演示的 plotgrid 函数 非常好用
-
7随便执行一点有问题的代码即可重现,比如: Plot[lbk]1,1,1[rbk] 此bug于14.2版本引入,截至14.3未修。规避此bug的方法是,双击那三个点打开菜单,再去选“显示堆栈跟踪”。参se帖子《Mathematica 14.2 can not show stack trace window》(编号 310769)。
-
3
-
6请问下各位大神,我用Reduce约化后这个结果,t有一个范围,但是这个b有#号,我不太懂,应该怎么处理呀,有没有办法可以用含t的表达式表示b呢? Reduce[{(256 t + b (-128 + 32 b^2 t - 16 b (-1 + 2 t) (7 + 2 t) + b^3 (-9 + t (14 + t (23 + 4 t)))))/(16 (-8 + b^2 (3 + t))^2) > 0, 0 < b < 1, 0 < t < 1}, {t, b}] (0 < t <= 1/4 && 0 < b < Root[256 t - 128 #1 + (112 - 192 t - 64 t^2) #1^2 + 32 t #1^3 + (-9 + 14 t + 23 t^2 + 4 t^3) #1^4 &, 2]) || (1/4 < t < 1 && 0 < b < 1)
-
14各位大神求解,我想用vb调用mathematica生成的dat或txt来画图,但是vb菜鸟,求各位大神帮忙,小弟感激不尽。
-
6此问题是在阅读SE帖子《Recursive anonymous functions with |-> notation》(编号316940)及《Can someone explain this snippet: (#[#] &)[#[#][#] &]》(1532)中产生的。闲话少叙,谁能给我讲讲怎么理解下面这段代码?: Y[f_] := #[#] &[Function[n, f[#[#]][n]] &] fac[r_] := If[# < 2, 1, #*r[# - 1]] & Y[fac][4] (* 24 *) 啊我并不是想要求教这段代码的语法,这段代码涉及的所有语法我(大概)都懂,但是连在一起我就不明白了……应该说,这段代码(似乎)并没有用到什么高级语
-
12mathematica中我希望定义一个抽象的方向导数求导,这个方向设为v,我不在乎v的坐标。输入d[v, f]得到微分Dt[f]展开式中关于自变量的微分比如Dt[x]能自动替换为d[v, x],比如说d[v, Sin[x y]]=Cos[x y](x d[v, y]+y d[v, x]).然后有一些替换规则把里面得到的d[p, P]替换为其他的不含d[]的表达式。怎么弄?
-
3NSolve[1/2*(Exp[x] + Exp[-x]) == Exp[x^2/8], {x}] NSolve::nsmet 无法利用 NSolve 现有的方法求解该系统. 请尝试使用 Reduce 或 FindInstance. FindRoot[1/2*(Exp[x] + Exp[-x]) == Exp[x^2/8], {x, 8}] 可以得到三个解之一 {x -> 7.23339} 我不太确定吧里是不是有过相关的讨论。
-
3
-
4
-
0鉴于过时教材的“流毒”甚广,特开此标题易懂的主题说此问题。为了便于修订,正文已放在了云笔记上,请大家多点一下鼠标: http://note.youdao.com/noteshare?id=0c2719208239696d61182199327b 有什么意见或建议欢迎提出。
-
2在执行代码 ResourceFunction["FunctionParity"][x*Sin[x], x] 时报错 ResourceFunction::ffail: Failed to find resource function HoldForm[<|Name->FunctionParity,ShortName->FunctionParity,<<6>>,FunctionLocation->CloudObject[https://www.wolframcloud.com/obj/7829e729-2b53-495b-9ad1-2928a85]|>,$Failed]. ,即使更换待求解函数依然报错,请问如何解决?谢谢
-
3为了验证Kirchhoff积分公式,代码如下: In[573]:= \[Lambda] = 10; k = (2 \[Pi])/\[Lambda]; n = {x, y, z}; R = {1000, 0, 0}; \[Psi] = Exp[-I k Sqrt[(x - R[[1]])^2 + y^2 + z^2]]/Norm[R]; \[ScriptCapitalE] = ({1, -1, 0} Exp[(I k (x + y))/Sqrt[2]])/Sqrt[2]; t1 = 1/(4 Pi) (NSurfaceIntegrate[\[ScriptCapitalE] ( Grad[\[Psi],{x, y, z}] . n) - \[Psi] (\!\( \*SubscriptBox[\(\[PartialD]\), \({{x, y, z}}\)]\[ScriptCapitalE]\) . n), {x, y, z} \[Element] Sphere[], Method -> "LocalAdaptive"] // Norm) t2 = 1/(4 Pi) (NSurfaceIntegrate[ Grad[(\[Psi] n . \[ScriptCapitalE]),{x, y, z}] - n
-
58贴吧新弄了个贴条,贴条只能同时存在10个并且只能实际只能显示26字,于是我就在想着是不是该总结些Mathematica新手最常见的问题贴到那上面滚动播出(虽然很多人或许不会看,但有总比没有好嘛……)。试着弄了个草稿(二十戒),大家可以看看,哪些需要修改,哪些需要删节或补充: 二十戒 1 不要为了节约硬盘选择过老的版本。至少要用版本7 2 自带帮助是最好的教材,提问前先打开软件按下F1仔细查查 3 将光标停在不认识的函数前/中/后再
-
2相当古怪的bug,开帖说一下。此bug的重现步骤如下。 1. 执行如下代码: BeginPackage["abc`"]; EndPackage[]; abc`Det; Det[{{1, 2}, {3, 4}}] 执行后会发现Det变红,且最后一行原样返回,这是意料之中的。 2. 执行 Exit 关闭内核。 3. 重启内核,执行 BeginPackage["abc`"]; EndPackage[]; 注意,Det[{{1, 2}, {3, 4}}]和上面两行代码不要放在同一个单元内。执行完上述代码后,你会发现,Det立刻变红了,进一步执行 ?? *`Det 我们会看到 abc`Det 已经冒出来了,尽管我们这一轮根
-
11
-
2求助,输入1+1类似的式子后回车直接输出2了不是换行,shift enter也是输出,有没有什么解决办法
-
6我想问一下这个报错应该如何解决 findBinaryCriticalValueOptimized[start_, end_, precision_] := Module[{left = 1, right = 4, mid, fLeft, fRight, fMid, iter = 0, maxIter = 20, amp},(*初始化函数值*)amp = left; NotebookEvaluate[ "A", InsertResults -> True]; fLeft = ampp; amp = right; NotebookEvaluate[ "A", InsertResults -> True]; fRight = ampp; (*检查区间有效性*) If[fLeft == fRight, Return["错误:函数在区间端点值相同,请调整起始和结束值。f(left) = " <> ToString[fLeft] <> ", f(right) = " <
-
4我想在下面的测试程序中画出来图应该如何做,我使用plot的时候没办法在用图一程序运行的时候跑出来图(意思是我查看保存下来的程序的时候发现保存的程序并没有画图出来),但是使用print的时候会直接在图一中的程序里面输出。 代码: (*定义要测试的变量范围和步长*)startValue = 0.2; (*起始值*) endValue = 1; (*结束值*) stepSize = 0.1; (*步长*) (*循环遍历变量值*) Do[Global`amp = varValue; NotebookEvaluate[ "A"], {varValue, startValue, endValue, stepSize}] Global`amp Plot[Sin
扫二维码下载贴吧客户端
下载贴吧APP
看高清直播、视频!
看高清直播、视频!



