mathematica吧
关注: 19,842 贴子: 73,929

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

  • 目录:
  • 程序设计及人工智能
  • 0
    大家有兴趣的话可以体验一下,当然了,文档还没得查。正式发布估计不远了。
    xzcyr 17:55
  • 2
    第一次用Mathematica画图,求助大佬们~ 已知f(k)是一个关于k的函数,其中k>1/2。r是其中的一个参数,范围为[0,1] 我想要做的:画出f(x)的图像 我的代码: Assuming[0 <= r <= 1, Plot[-2 (3 - r) r k^3 - (r^3 - 2 r^2 - 7 r + 4) k^2 + (r^3 - 6 r^2 - 4 r) k + r^2 - r, {k, 1/2, 1000}]] 我的结果: 想问下为啥只显示坐标轴,不显示图像呢?
    xzcyr 17:52
  • 5

    广告
    20:15
    不感兴趣
    开通SVIP免广告
  • 4
    func1[f_]:=Module[{}, Table[{f,i},{i,1,3}] ]; func2[f_]:=Module[{i}, Table[{f,i},{i,1,3}] ]; func1[i+j] func2[i+j] 按照文档的说法Table 会局部化变量,但是输入的 i 被视作索引变量了,这是特性吗
    xzcyr 17:50
  • 2
    {DSolve[{t''[s] - (xx'[s] t'[s])/(9000 (1 - xx[s]/9000)) == 0, xx''[s] + xx'[s]^2/( 18000 (1 - xx[s]/9000)) + ((1 - xx[s]/9000) t'[s]^2)/18000 == 0, t[0] == 0, t'[0] == 1, xx[0] == 0, xx'[0] == 0}, {t[s], xx[s]}, {s, 0, 1}], DSolve[{t''[s] - (hh'[s] t'[s])/(9000 (1 - hh[s]/9000)) == 0, hh''[s] + hh'[s]^2/( 18000 (1 - hh[s]/9000)) + ((1 - hh[s]/9000) t'[s]^2)/18000 == 0, t[0] == 0, t'[0] == 1, hh[0] == 0, hh'[0] == 0}, {t[s], hh[s]}, {s, 0, 1}]} (* Out: {{{t[s] -> 18000 ArcTan[s/18000], xx[s] -> -(s^2/36000)}}
    xzcyr 17:44
  • 0
    我想要用mathematica化简一个式子,这个式子占用内存估计要10多G,算过的人告诉我可以要算7,8天。我想知道有没有什么可以远程算的云平台,有没有推荐?
    laivirt 8-1
  • 5
    lz没有系统性的学过积分方程,想用mathematica求解积分方程,查阅帮助文档后发现,mathematica好像没有关于积分方程的文档,百度了一下,一维积分方程可以用Sovle求解,但对于我的问题却失效了,以下是代码: eqn = {Integrate[f[x]/(2\[CenterDot]\[Pi]) 1/(x1 - x), {x, 0, c}] == Subscript[V, \[Infinity]] \[Alpha], f[0] == Infinity, f[c] == 0} Solve[eqn, f[x]] 求问mathematica怎么求解积分方程qvq
  • 3
    贴吧客户端相较于传统的pc网页版弊端颇多(如中括号无法在顶楼正常显示,置顶只能显示一条等),但吧里始终有一小部分用户(访问数约占7%)在使用客户端,故开此帖调查。视调查结果可能会在下个答题日出篇教程。
  • 0
    贴吧客户端相较于传统的pc网页版弊端颇多(如中括号无法在顶楼正常显示,无法自由对帖子文字进行复制,置顶只能显示一条等),因此,如果因为某些原因不得不使用手机访问贴吧,使用网页浏览器访问pc版网页其实是一个不错的选择。这里介绍一下修改手机浏览器设置的方法。我的手机里目前只有小米系统自带浏览器和Opera,就以这俩为例了。其他浏览器的设置修改方法(大概)是大同小异的。先说小米系统自带浏览器。最新版本的小米浏览器设
    xzcyr 6-7
  • 72

    广告
    20:09
    不感兴趣
    开通SVIP免广告
  • 13
    求助有关ListVectorPlot的一些问题,下面的123分别对应着三幅图产生的问题: 1.首先Rainbow最亮的颜色是红色,下面这串代码将大于0.5的部分全部设置为最亮的颜色红色,但其实对于EVPA中的最大值是3.6,这样导致形成的图像大面积的全是红色(存在的问题可能是数据中存在较大的值,而大部分值相对较小,那么图例的颜色范围和实际图像中的颜色会不一致,个人猜测,不一定。)我想要的结果是只在最大值附近才是红色,让红色面积减少。 2.我通过将Vecto
  • 0
    症状: 顶楼中的 [ 和 ] 会被自动转换为 [lbk] 和 [rbk]。此现象在回复帖中不存在。 影响范围: 此bug至少在去年10月16日即已出现,至今(安卓12.58.1.0)未修。IPhone受影响情况不明,欢迎实验补充。 对策: 对于提问者,最简单的对策就是改用电脑访问贴吧网页版,不要使用使用手机客户端发帖。如果确有相关需求(所以说你们到底为什么会有这种需求?电脑没连网?那也可以开手机热点啊,访问网页又花不了什么流量),避免在顶楼张贴代码也是一个
    xzcyr 4-6
  • 1
    刚看见一个相关的问题贴子, 想要回复却不见了. 在Mathematica中"永久的"储存符号定义, 我建议使用Save或者DumpSave 比如 a = 5;b = 10;c = 15; Save[LocalObject["definitions"], "Global`*"]; 再次打开mathematica时, 使用 Get[LocalObject["definitions"]]; 调取存储的a, b, c的数值. 如果符号是物理常量, 还建议使用以下语句直接获取国际单位制下的物理常量数值 比如 First@N@UnitConvert@Quantity[1, "PlanckConstant"] First@N@UnitConvert@Quantity[1, "ReducedPlanckConstant&quot
    物理迷2 7-10
  • 14
    各位前辈好,最近有一个很大的数据集需要用mathematica处理,楼主将之分成了小块希望能够并行化,以充分利用服务器资源。但遇到了许多问题。首先是尝试使用P打头的自动并行方法,失败。然后楼主又试着增加内核数量,手动并行。但在网上查找了一圈(在自带帮助发现了可能有用的math命令,wolfram命令和WolframKernel命令都失败了,谷歌,必应,百度也都没有结果,哭)也没有成功在命令行新建并调用内核。不知道有没有前辈遇到过这个情况,请前辈
    xzcyr 2-4
  • 6
    现在我有一自定义函数是这样的: VectorsOfOnePoint[f_,g_,m_]:=(Module[{s1,s2,s3,r1,r2,x0,y0,z0,s}, s1[s_]:=ReplaceAll[f,{x->x0
    xzcyr 4-7
  • 1
    如图,很容易注意到有大量的类似于thetaL-thetaR的结构,想把它替换成一个值theta,似乎由于指数结构的多样性,没法简单替换QAQ。 代码:expr = Expand[ Denominator[ t1S]/(\[CapitalGamma]^6* Exp[3*I*\[Theta]L + 3*I*\[Theta]R + I*\[CurlyPhi]])]; expr = expr /. Power[E, a_] :> Exp[a]; exprSimplified = expr //. {a_. Exp[ I*p_.*\[Theta]L + I*q_.*\[Theta]R + I*r_.*\[CurlyPhi]] :> a Exp[I*(p - q)*\[Theta] + I*r*\[CurlyPhi]]} expr /. {I *(\[Theta]L - \[Theta]R) -> I*\[Theta]} 结果: {36 + 66 E^(I \[Theta]L - I \[Theta]R + I \[CurlyPhi]) + 66 E^(-
    xzcyr 7-5
  • 32

    广告
    20:03
    不感兴趣
    开通SVIP免广告
  • 9
    刚入职一家规模很小的公司,我需要用mathematica学习一些物理上的东西,做做计算。如果我在公司电脑上装mathematica,用算号器激活,会有风险吗? 目前在使用mathematica eigine+vscode 效果差完整版mathematica太远了
    xzcyr 7-5
  • 7
    请问,从A点到B点的最速降线,假设A点在第二象限,B点在原点,问题是摆线一定是在直线{-Pi t, 2 t}之下的,那如果A点在该直线之上,比如{-1,2},不就无解了吗? ParametricPlot[{3 {-\[Pi] t, 2 t}}~Join~ Table[{-\[Pi] + t - Sin[t], 1 + Cos[t]} r, {r, 9}], {t, 0, \[Pi]}, AspectRatio -> Automatic]
    xzcyr 7-5
  • 4
    Plot[10^-12*x,{x,1.2*10^6,10^8},PlotRange->{{0,1*10^8},{0,5.1*10^-5}},Frame->True,FrameTicks->{{Table[{y,ScientificForm[y]},{y,0,0.00005,0.00001}],Automatic},{Automatic,Automatic}},FrameTicksStyle->Directive[Black,18,FontFamily->"Times"],FrameLabel->{Style["s [GeV^2]",Black,18,FontFamily->"Times"],Style["Subscript[\[Sigma], VBS/VBF]Subscript[d\[ScriptCapitalL], \[Gamma]\[Gamma]]/d\[Tau] [pb]",Black,18,FontFamily->"Times"]},PlotStyle->{ColorData[24][5],Thickness[0.005]},ImageSize->500] 这个图画出来的纵坐标
  • 0
    就是我用的mathematica 11.2,今天突然发现字体变模糊了,正体的字体没问题,斜体变模糊了,如下图 问题是前几天都是正常了,我这几天也没干嘛,字体就变模糊了 我把缩放比例调到150%,字体又正常了,如下图 我用的是win11系统,我在虚拟机上测试了win10系统的情况,虚拟机win10上打开mathematica 11.2,显示正常,斜体不模糊,也就是说只有我的win11才模糊,而且是这几天才出现的状况,前几天都是正常的。 我卸载mathematica 11.2,并删除了C盘的各种mathematic
    247098254 6-29
  • 2
    这段程序就是单纯使用NSolve去求解一个非线性系统的平衡点,但是程序一直运行不出来😢😢,一直显示正在运行。而且CPU基本没运行是怎么回事🥲🥲,刚学这个软件,是不是存在一些操作问题。 NSolve[{ -x0 + 1.16*Tanh[x0] - 1.8*Tanh[x1] + 1.0*Cos[x6]*(x0 - x3) == 0, -x1 + 3.2*Tanh[x0] + 1.1*Tanh[x1] + 1.3*Tanh[x2] == 0, -x2 + 3.0*Tanh[x1] + 1.5*Tanh[x2] == 0, -x3 + 1.16*Tanh[x3] - 1.8*Tanh[x4] + 0.0*Tanh[x5] - 1.0*Cos[x6]*(x0 - x3) == 0, -x4 + 3.2*Tanh[x3] + 1.1*Tanh[x4] + 1.3*Tanh[x5] == 0, -x5 + 3.0*Tanh[x4] + 1.5*Tanh[x5] ==
    Ccclzz 6-29
  • 7
    也算是对这帖(http://tieba.baidu.com/p/2964416898)的2,3,4条的解说吧。有什么意见或建议欢迎提出: http://note.youdao.com/share/?id=abd51087f44c0b6a41ff6022d549dc4
  • 55

    广告
    19:57
    不感兴趣
    开通SVIP免广告
  • 6
    计算的结果可能到10的-20次方这样,但是mathematica好像这么小的数字直接默认为0了,网上找的方法有限还不好用
    xzcyr 7-5
  • 13
    请问下N函数是不是有适用范围的,比如我想对我算出的数值保留两位,结果发现没用,能问下是为什么?谢谢! In[18]:= ({3.34, 3.32, 3.35, 3.32} - 3.24)^2 // Total // Sqrt Out[18]= 0.186815 In[19]:= N[%, 2] Out[19]= 0.186815
    xzcyr 7-6
  • 7
    最近在计算指数很小的EXP的时候,发现了当指数为小数时,计算结果自动归零的问题,有没有大佬给解答一下该怎么办。代码:Exp[N[-784.983,$MachinePrecision]]Exp[N[-800,$MachinePrecision]][图片]
    xzcyr 6-4
  • 7
    一个数用两种类型表示算出的结果为何差距这么大: N[Integrate[ Exp[(10^(-4))*t]*ArcTan[((10^(-4)) - t)^2], {t, 0, 10^(-4)}]/ Integrate[2*t*((10^(-4)) - t), {t, 0, 10^(-4)}], 10] 1.0000000025+ 0.*10^-11i N[Integrate[Exp[0.0001*t]*ArcTan[(0.0001 -t)^2], {t, 0, 0.0001}]/ Integrate[2*t*(0.0001 - t), {t, 0, 0.0001}], 10] -32.7418 请教各位这是怎么回事啊
    xzcyr 7-6
  • 2
    输入:RSolve[{a[n] == \!\( \*UnderoverscriptBox[\(\[Sum]\), \(i = 0\), \(n\)]\((a[i]* a[n - i])\)\), a[0] == a[1] == 1}, a, n] 会报错:RSolve:提供的函数不是给定函数的差分方程 输出:RSolve[{a[n] == \!\( \*UnderoverscriptBox[\(\[Sum]\), \(i = 0\), \(n\)]\(a[ i]\ a[\(-i\) + n]\)\), False == 1/3 == 1}, a, n] 还将a[1]识读成了1/3.将[方程,a,n]的a换成a[n]也没用。求助应该怎么解
  • 0
    提问时请附上引发问题的代码。对标题论断有疑虑的同学请阅读下面这篇对 FindRoot::nlnum 警告的简要解说: https://zhidao.baidu.com/question/810397720290271692.html
    xzcyr 5-2
  • 15

    广告
    19:51
    不感兴趣
    开通SVIP免广告
  • 3
    p1 = {-1, 0}; p2 = {1, 0}; p3 = {0, 1}; curve = BSplineCurve[{p1, p3, p2}]; {RegionQ[curve], RegionDimension[curve]} ContourIntegrate[Sin[z], z \[Element] curve] NContourIntegrate[Sin[z], z \[Element] curve] 除了手动构造NURBS曲线函数以外,有没有直接利用BSplineCurve或BSplineFunction计算积分的方法, 我尝试了DiscretizeRegion,结果不是很满意
    N_update 6-16
  • 5
    大佬们,我生成了一个矩阵,再定义了对矩阵进行一些计算一个的Module,单独计算这个module大概五秒钟可以出结果,但是我使用Parallel一系列的函数程序都会卡死,无法进入计算,不知道是为什么。下面是我的代码
  • 2
    比如我作图像:ContourPlot3D[{z==0,x^2/2+y^2/4==1,z==1+Log[3-x^2/3+y^2/6]},{x,-3,3},{y,-3,3},{z,0,5},Mesh->6],我只能给每个图像规定相同的网格线条数,可不可以第一个为3条,第二个为5条,第三个为2条,该怎么做,请教各位
  • 1
    wolframscript math wolfram 看样子好像是三个命令行。但是具体不知道有啥区别。。。 最近想研究下载latex中直接执行mma的命令,突然发现,wolframscript好像一直卡着。。。软件应该是激活了的,就是不知道啥问题。三台电脑,就一台不行,重装过mma了,也不行。 正常电脑安装完,path配置好, latex直接运行加上shell escape超级快。 这台好像就是卡住了。。。。 单开cmd的话Input [1] 也出不来 谁有遇到过?
  • 3
    大家有用过Wolfram notebook assistant+LLM Kit的吗,感觉怎么样。我刚订阅了一个月的基础版试了一下,感觉好像有点鸡肋,我在vscode里写mma用copilot似乎也能达到差不多的效果。是我不会用吗 还有assistant和用其他模型api的chat相比有什么优势吗,以前没用过chat。
    隨意超 6-8
  • 47

    广告
    19:45
    不感兴趣
    开通SVIP免广告
  • 7
    DA = 1.; r = DA/2; Nx = 1.; Ny = 1.; Ng = 10; SA = Pi/Ng; \[Eta] = 0.5; f = \[Eta] DA; \[CapitalOmega] = Triangle[{{0, 0}, {r Cos[SA], r Sin[SA]}, {r Cos[-SA], r Sin[-SA]}}]; Graphics[\[CapitalOmega]]; \[Beta] = ((1 + D[w[x, y], y]^2)/(1 + D[w[x, y], x]^2))^2; \[Tau] = (Nx/Ny)*\[Beta]; eqn = \[Tau]*Derivative[2, 0][w][x, y] + Derivative[0, 2][w][x, y] == 0; sol = NDSolveValue[{eqn, DirichletCondition[w[x, y] == (x^2 + y^2)/(4 f), y == Tan[SA] x], DirichletCondition[w[x, y] == (x^2 + y^2)/(4 f), y == -Tan[SA] x]}, w[x, y], {x, y} \[Element] \[CapitalOmega]] 解以上非线性偏微分方程会
    xzcyr 6-7
  • 5
    (*如何检查一个表达式是否是广义上的字符串(可能被Style修饰过,也可能由StringJoin或Row拼接而成)*) s1 = StringJoin@{"abc","def"}; s2 = Row@{"abc","def"}; s3 = Row@{"abc",Style["def",Red]}; s4 = Row@{Style["abc",Green],"def"}; s5 = Row@{Style["abc",Blue], Style["def",Pink]}; (*我有(广义上的)字符串s1~s5,我尝试写了RichStringQ,但他工作得不好。请赐教修改*) RichStringQ[expr_]:=StringQ[expr]||MatchQ[expr,_[s_String,___]/;MemberQ[{Style,R
    xzcyr 6-7
  • 4
    反函数的结果多值的,我只要其中一种即可,有没有办法去除这个井号?代码和图片如下 e[a_] =a^2 + (\[CapitalLambda] Subscript[m, 1] Subscript[m, 2] Subscript[m, 3])/(8 a^2 + 2 \[HBar]^2) - (1/( 2048 (a^2 + \[HBar]^2) (4 a^2 + \[HBar]^2)^3))\[CapitalLambda]^2 ((4 \ a^2 + \[HBar]^2)^2 ((4 a^2 + \[HBar]^2) (4 a^2 + 5 \[HBar]^2 - 4 \!\(\*SubsuperscriptBox[\(m\), \(3\), \(2\)]\)) - 4 \!\(\*SubsuperscriptBox[\(m\), \(2\), \(2\)]\) (4 a^2 + \[HBar]^2 + 12 \!\(\*SubsuperscriptBox[\(m\), \(3\), \(2\)]\))) - 4 \!\(\*SubsuperscriptBox[\(m\), \(1\), \(2\)]\) ((4 a^2 + \[HBar
    xzcyr 6-7
  • 9
    我的目标是求解abcd2的值,对于整个问题我的思路是将函数U拟合成D-10乘x的-10次到D10乘x的10次的和,其中D-10到D10都是实部+虚部,然后将φψ函数也表示为洛朗级数,然后比对每个x的n次方前面的系数来求解问题。我的实部部分求解的没有问题,但是虚部部分算出来就明显不对,希望得到帮助
    xzcyr 6-7
  • 1
    如题,请问如何在mathematica中作出这种图?就是画出函数图像后,在坐标轴上标注出点的坐标或者标注一个字母。
  • 11

    广告
    19:39
    不感兴趣
    开通SVIP免广告
  • 1
    \[Integral]1/ 2 Sqrt[-((-4 r^2 - 16 r^3 + 8 r^4 - 8 r^5 + 4 r^6 (1 - \[Omega]^2))/(r^2 (-2 r + r^2)^2))] \[DifferentialD]r 代码太长超过贴吧要求字数了,结果就放了个图
    Repentanze 6-24
  • 2
    我知道一个函数f(x, y)在一个曲线y=x^2的取值, 即f(x, x^2)=g(x)。并且这个函数的各阶偏导数在曲线y=x^2上的取值也知道, 即f^{i, j}(x, x^2)=g_{i, j}(x)。 现在想对f(x, y)的一个表达式(比如说x^2f^2(x, y)+f(x, y)), 求各阶偏导数,然后在曲线y=x^2上取值,该怎么简单操作?因为我不知道f(x, y)的具体表达式,不知道怎么做替换规则。
  • 4
    如何用MMA实现如matlab中的这种局部放大图呢 ,除了用绘图工具,可以用代码实现吗
    xzcyr 4-1
  • 5
    C盘的空间很小,请问怎样将Mathematica的Paclets默认下载地址设置到E盘? In[1]:= $UserBasePacletsDirectory Out[1]= "C:\\Users\\hbkaq\\AppData\\Roaming\\Wolfram\\Paclets"
    冉qx 5-13
  • 5
    比如x^2+y^2+xy+xy^2,处理之后只剩下x^2+y^2+xy,只要有y^2和xy,xy^2这一项就去掉。但如果是x^2+xy^2,处理过后还是x^2+xy^2
  • 5
    RT,是动态图还是平面图啊? ↖( ̄▽ ̄") 反正都不会。求教教~~ 小行星的一般方程和标准方程:&#x2212;0.33781&#x2217;&#xD465;^2+2&#x2217;0.189161&#
    xzcyr 5-3
  • 3

    广告
    19:33
    不感兴趣
    开通SVIP免广告
  • 3
    将复数化为三角表达式与指数表达式,用Mathematica怎么实现?
    xzcyr 5-3
  • 25
    搜了一下,发现吧里好像没有主题帖专门说这个。(尽管我觉得这应该是任何一个有心学习Mathematica的人都能很容易地接触到的常识。)那就专开一帖好了。 为了便于修订,照例放到云笔记上: note.youdao.com/noteshare?id=7334ad886b97de3565d59beea50b5350 ……其实不看云笔记都无所谓,概括起来就那么几句话: 1. For循环写起来很冗长——这其实是其次的; 2. For循环在Mathematica里很慢——这其实也是其次的; 3. For循环的表意逻辑,几乎在任何情况下,都无意义地复
    xzcyr 5-3
  • 9
    OneBounce[k_, ramp_][{t0_, x0_, xp0_, y0_, yp0_}] := Module[{sol, t1, x1, xp1, y1, yp1, gramp, gp}, sol = First[NDSolve[ {x''[t] == 0, x'[t0] == xp0, x[t0] == x0, y''[t] == -9.8 , y'[t0] == yp0, y[t0] == y0}, {x, y}, {t, t0, \[Infinity]}, Method -> {"EventLocator", "Event" :> y[t] - ramp[x[t]]}, MaxStepSize -> 0.01]]; t1 = InterpolatingFunctionDomain[x /. sol][[1, -1]]; {x1, xp1, y1, yp1} = Reflection[k, ramp][{x[t1], x'[t1], y[t1], y'[t1]} /. sol]; Sow[{x[t] /. sol, t0 <= t <= t1}, "X"]; Sow[{ y[t] /. sol, t0 <= t
    xzcyr 5-3

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

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

扫二维码下载贴吧客户端

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

本吧信息 查看详情>>

小吧:小吧主共10

会员: Expr

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

友情贴吧