mathematica吧
关注: 18,217 贴子: 69,428

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

  • 目录:
  • 程序设计及人工智能
  • 0
    DSolve错误:为避免歧义,因变量的参数应该与自变量完全匹配 什么意思?怎么解决? 急!
    1010__0101 22:08
  • 0
    Exp@# + #^2 &@Range[2, 100, 2] Exp@# + #^2 & // Range[2, 100, 2]
  • 0
    Remove[a]; Do[a[n] = AbsoluteTiming[ Module[{sum = 0, i}, For[i = 1, i <= 10^5, i++, sum += Sin[N@i]]; sum]], {n, 1, 10}]; Do[b[n] = a[n][[1]], {n, 1, 10}]; Table[b[n], {n, 1, 10}]; Mean@%
  • 4
    MinimalPolynomial可以求显式代数数的最小多项式, 如MinimalPolynomial[Sqrt[2] + Sqrt[3] + Sqrt[5], x] 结果是576 - 960 x^2 + 352 x^4 - 40 x^6 + x^8 如果式子带变量,比如Sqrt[a] + Sqrt[b] + Sqrt[c], Sqrt[a + 1] + (b + 1)^(1/3) +a+b +1这样的有没有统一的方法求出最小多项式呢 (前者结果应当是a^4-4 a^3 b+6 a^2 b^2-4 a b^3+b^4-4 a^3 c+4 a^2 b c+4 a b^2 c-4 b^3 c+6 a^2 c^2+4 a b c^2+6 b^2 c^2-4 a c^3-4 b c^3+c^4+(-4 a^3+4 a^2 b+4 a b^2-4 b^3+4 a^2 c-40 a b c+4 b^2 c+4 a c^2+4 b c^2-4 c^3) x^2+(6 a^2+4 a b+6 b^2+4 a c+4 b c+6 c^2) x^4+(-4 a-4 b-4 c) x^6+x^8
  • 5
    b = \[Rho] \[Alpha];p = \[Alpha] - \[Beta] q;w = p - \[Sigma];\[Pi]M[q_] = (w - C) q - b (e q - C) - r b e q;D[FullSimplify[\[Pi]M[q]], q] 得出的结果是这样的 \!\(\*SubscriptBox[\(\[PartialD]\), FractionBox[\(\[Alpha]\ \((\(-1\) + e\ \[Rho] + \[CurlyPhi] + \*FractionBox[\(\(-1\) + e\ \[Rho] + 2\ \[CurlyPhi] - \*SuperscriptBox[\(\[CurlyPhi]\), \(2\)]\), \(\(-2\) + \\[CurlyPhi]\)])\)\), \(2\ \[Beta]\ \((\(-1\) + \[CurlyPhi])\)\)]]\*FractionBox[\(\[Alpha]\ \((2\ C\ \((1 - e\ \[Rho] + 2\ \[Beta]\ \[Rho]\ \((\(-2\) + \[CurlyPhi])\))\)\ \((\(-2\) \+ \[CurlyPhi])\) - \((\(-1\) + e\ \[Rho])\)\
    cxy17506 08:56
  • 10
    如题:想求出sin x=x/6在区间[0, 6]内的所有数值解。请问有什么命令可以一次性把这个区间上的所有所有数值解全部给出来? 求高手解答,谢谢!
  • 1
    如题,在求一个方程的所有零点,可能这个方程比较复杂(包含合流超几何函数),总之用findroot,NSolve等求出来的零点都不全,后来看大佬 @xzcyr 的回复,安装了RootSearch,一开始求解很ok,后来换了个方程找出来的零点也不全了。查了下,发现用RootSearch算完之后,可以用RootSearchSample来找出几乎所有零点,于是打算这两个函数结合起来。 然后新的问题出现了,因为我要求很多个方程的零点,就用到了循环,RootSearch可以循环,循环结束后,不知道怎么
    xzcyr 8-6
  • 1
    heqn=D[u[x,t],t]==D[u[x,t],{x,2}]-D[u[x,t],{x,1}]; ic=u[x,0]==Piecewise[{{x,0<= x<= 0.5},{1-x,0.5<= x<= 1}}]; bc={u[0,t]== 0,u[1,t]== 0}; sol=DSolve[{heqn,bc,ic},u[x,t],{x,t}]
    xzcyr 8-6
  • 11
    有没有大神可以补充的。鼠标使用起来不太方便。
    xzcyr 8-6
  • 13
    如图,想提取图中伸出部分圆柱上边缘上点的坐标,并且以背景的一小格为单位,求问该怎样用代码实现? 我现在用的笨办法是用绘图工具在边缘上点一堆红
    xzcyr 8-6
  • 1
    我这都运行1天了,是写的有问题,还是运行时间不够,再继续运行下去啊! Integrate[ 1/(Sqrt[(x - m)^2 + (y - n)^2 + (z - h)^2])^3, {x, -b, b}, {y, -l, l}, Assumptions -> n > 0 && m > 0 && z > 0 && h > 0 && b > 0 && l > 0]
    xzcyr 8-6
  • 2
    要用欧拉方程算一个最速降线(不是重力场中的),然后因为初始条件是无速度从一点出发,结果数值求解的时候会在出发点碰到1/0。这个可以精确求解吗?还是得加个初速度,但是初速度稍微小一点就会碰到哪个地方步长为零,遇到奇点. 这是代码: Clear[r, \[Theta]]; Remove[r, \[Theta]]; \[Omega] = 3; Needs["VariationalMethods`"]; equ = Numerator[ EulerEquations[(\[Sqrt](r'[\[Theta]]^2 + r[\[Theta]]^2))/(\[Sqrt]( r[\[Theta]]^2 - ( 4/5)^2)), r[\[Theta]], \[Theta]] /. Equal -> Plus]/ r[\[Theta]]
    xzcyr 8-6
  • 4
    求矩阵特征值的时候,用“Eigenvalues[a]”有的矩阵可以输出含有根号的精确解,比如 a = {{0, 1, 2, 0, -2, -1}, {-1, 0, 1, 2, 0, -2}, {-2, -1, 0, 1, 2, 0}, {0, -2, -1, 0, 1, 2}, {2, 0, -2, -1, 0, 1}, {1, 2, 0, -2, -1, 0}} 有的矩阵只能输出含小数点的近似解,比如 b = {{1, 2, 3, 3, 3}, {3, 1, 2, 3, 3}, {3, 3, 1, 2, 3}, {3, 3, 3, 1, 2}, {2, 3, 3, 3, 1}} 请问怎么让他输出含根号的那种解啊?
    xzcyr 8-6
  • 2
    mA = 1; mB = 3; h = 1; l = 1; time = 0.6; \[Theta]0 = ArcTan[h/l]; s = NDSolve[{h/(Sin[\[Theta][t]])^2 \[Theta]'[t] == Sqrt[( 20*h*mB*(1/Sin[\[Theta]0] - 1/Sin[\[Theta][t]]))/( mA + mB*(Cos[\[Theta][t]])^2)], \[Theta][ 0] == \[Theta]0}, \[Theta], {t, 0, time}]; \[Theta] = \[Theta] /. s[[1]]; Plot[{Sqrt[(20*h*mB*(1/Sin[\[Theta]0] - 1/Sin[\[Theta][t]]))/( mA + mB*(Cos[\[Theta][t]])^2)], Sqrt[(20*h*mB*(1/Sin[\[Theta]0] - 1/Sin[\[Theta][t]]))/( mA + mB*(Cos[\[Theta][t]])^2)]*Cos[\[Theta][t]]}, {t, 0, time}] Clear["Global`*"]
    xzcyr 8-6
  • 9
    l=32 dict=Table[2^i,{i,1,l}] value=Range[Length[dict]] X=SparseArray[dict->value] Pos=Table[dict[[RandomInteger[{1,Length[dict]}]]],{i,1,Length[dict]}] SubsetQ[dict,Pos] Y=X[[Pos]]
    xzcyr 8-6
  • 5
    SetDirectory[NotebookDirectory[]]; Import["record.txt"] ReadFile = OpenRead["record.txt"]; ReadLine[ReadFile] ReadLine[ReadFile] Close[ReadFile] ----------------------- 代码如上。请问ReadLine有办法读中文吗?怎么做? 环境:Win11家庭中文版+Mathematica 12.3.1.0中文版
    xzcyr 8-6
  • 1
    使用算回归方程的资源函数的时候,展示函数表达式出了点问题,代码跟官方文档里的例子一样: ResourceFunction[ "RegressionListPlot"][{{0, 10}, {2, 8}, {1, 12}, {6, 1}, {2, 3}, {3, 4}, {4, 6} {1, 1}, {9, 5}, {9, 6}, {12, 2}, {10, 4}, {13, 12}}, "Degree" -> 2, "DisplayEquation" -> True] --------------------------------------------------------------------- 但最后输出的时候函数库貌似出了点问题 这是官方说明文档:https://resources.wolframcloud.com/FunctionRepository/resources/ 不知道是啥情况,
    xzcyr 8-6
  • 4
    有一个m*n大小的数组,想要每行都同时减一个1*n的行向量,请问mma里怎么优雅地实现?
    xzcyr 8-6
  • 5
    假设我已经有个Plot结果存为a了, a = Plot[Sin[x], {x, 0, 6 Pi}] 如何只对a操作,做到把曲线换个颜色再画出来? 谢谢!
    xzcyr 8-6
  • 4
    主要的求解问题如图片所示,代码如下: ClearAll["`*"](*清除变量*) LB=-50; RB=1000; r0=10; b=1; \[Sigma]=0.5; m=1; \[CapitalOmega]h=1; rh=1; tl=30; Ar[r_]:=(m \[CapitalOmega]h)/(1+ProductLog[E^(-1+r/rh)/rh])^2; Br[r_]:=(-4 m^2 rh^2 \[CapitalOmega]h^2+ProductLog[E^(-1+r/rh)/rh] (2+4 m^2+(-1+4 m^2) ProductLog[E^(-1+r/rh)/rh]))/(4 rh^2 (1+ProductLog[E^(-1+r/rh)/rh])^4); kge=D[\[Psi][r,t],{t,2}]==2 I Ar[r] D[\[Psi][r,t],t]+D[\[Psi][r,t],{r,2}]+Br[r] \[Psi][r,t]+NeumannValue[-Derivative[0,1][\[Psi]][r,t],r==LB]+NeumannValue[Derivative[0,1][\[Psi]][r,t],r==RB]; gaussian
  • 7
    考虑到此类问题和“Mathematica已经装好了,但不能用”即“Mathematica一开始就不能用”有一定区别,特开专帖讨论。今后此类问题原则上请不要单开主题。 这里的“不能用”,包括不限于 1. Mathematica突然无法正常打开 2. Mathematica突然变得无法正常工作(1+1也算不了等) 3. Mathematica突然在执行正确的代码时出现莫名其妙的警告。 如果你的“Mathematica一开始就不能用”,请看置顶。 闲话少叙,说一说可能的修复方法: 此类问题似乎和Mathematica的Paclet更新不
  • 3
    程序最后算完结果是 -9.67059*10^-306/-694.319 is too small to represent as a normalized machine number; precision may be lost. 请问这个应该怎么处理就会有结果呀
  • 15
    Solve[-w - \[Delta] + \[Tau] + qs (p + \[Alpha] - \[Tau]) g[qs] == 0, qs] Solve::nsmet: 无法利用 Solve 现有的方法求解该系统. 小白求问,这种式子该怎么解啊
  • 14
    中文的最好,如果英文的也可以 谢谢 不太喜欢看电子的,所以想买一本
  • 43
    1. 帮助文件安装很顺利,没有出现某网友遭遇的尴尬现象。 2. 但是在mathematica中查找某条命令的帮助文件时,会跳到Mathematica官网上的帮助网页。这样,似乎也就无所谓中文版和英文版的区别了。你上的官网是中文版的,网页上帮助文件就是中文版的。 3. 离线无法查找帮助。不知这算不算一个缺点?
  • 5
    代码类似于这样的 p = Permutations[Range[5]] o = Table[ "abcd改坐标的参数为" <> ToString[p[[i]][[5]]] "xyz", {i, 1, 9}] Export["C:\\Users\\8618\\Desktop\\Mathematica document\\In turn point output set\\The ring.txt", o, "Table"] 但是运行出文件里却是这样的 另外,是不是格式的问题,那在mma中什么公式可以显示格式呀
  • 5
    我现在只要打开mma,即便啥也不干,电脑就会偶尔卡一下。我发现只要我重新打开mma,第一次右键桌面都会卡一下。但平时不开mma就不会有这种问题。而且这个卡顿是最近出现的。
  • 3
    q是一个素数且满足:2^q±[2^(q+1)/2]+1是素数,2^q-1无平方因子. 请问用什么方法能尽可能多的找出满足上述条件的q呢?
  • 7
    也算是对这帖(http://tieba.baidu.com/p/2964416898)的2,3,4条的解说吧。有什么意见或建议欢迎提出: http://note.youdao.com/share/?id=abd51087f44c0b6a41ff6022d549dc4
  • 15
    此帖是本吧初代吧主 @mm_酱 发在果壳小组的一篇教程。如今果壳小组已死,mm_酱近期又不太活跃,这里姑且重发一份: https://note.youdao.com/share/?id=058e6037396d925af1f4abe4d54a52 mm_酱你要是想整理下出个v2就吱一声。
  • 5
    https://1drv.ms/u/s!Amd-uedNeUsttDBBDZWPGllTqawk?e=IohgLF
    nameETC 8-12
  • 99
    本文系自此帖(http://tieba.baidu.com/p/2736346260)5楼的回复修订增补而来。吧里的老人们阅读本文时可能会觉得各种眼熟:毕竟重要的经验我都是不断强调的。 我最初会得知Mathematica的存在,是因为《数学建模》老师在课堂上的推荐。我第一眼就被它那极度接近传统数学式书写的使用方法给吸引了——相信这个吧里的许多人也是一样。《数学建模》结课后不久校内就搞了数学建模大赛(目的是为了选拔参加全国大赛的选手),我就参加了,用的当然是Mathem
    Fluid 8-29
  • 2
    ParametricPlot[{200*f, g}, {r, 0.5,10}] 现在我尝试用ParametricPlot[]来画,但面临两个问题:(1)横纵坐标比例尺不能自动调节,比如我的例子中f、g都是r的函数,在上边r的取值范围内,T值都比较小,0.01的量级,G值比较大,10的量级,这样我就不得不把T值放大200倍以观察两者图线关系,有没有更优雅的办法?(2)怎么对f,g分别设置PlotRange?现在只要一用这个命令,就对f,g取了相同的取值范围。 谢谢!
  • 3
    请勿删,评论区有文本 就相当于a=1,b=1,c=a+bt,d=c+at然后求e=c+d
  • 5
    请问如图所示的MeijerG函数乘积的积分是如何用mathematical计算的呢, 我用mathematical算出来的是一个特别复杂的式子 或者说这个式子可以化成MeijerG的形式吗? 新人第一次用这个软件,有点懵逼
    kiginu 5-27
  • 5
    初学mathematica,遇到了点小问题。红框所示的三个矩阵相乘,把两个变量β和γ都替换成了α,输出结果是[Out]651;再运行[In]652给出的结果是替换后的结果。当我把[In]651和[Out]651都删掉后再运行[In]652,发现无法恢复替换之前的结果。我尝试了运行Rdα.Rdβ.Rdγ/.{β->β,γ->γ},发现还是替换后的结果。替换后如果还是想以β和γ输出,应该怎么做?
  • 17
    PikingDistance[list_] := Total[Table[ Max[Table[ SequencePosition[ PData[[m]][[list]], {Data4[[Result[[m]]]][[x]][[n]]}][[1, 1]], {n, 1, Length[Data4[[Result[[m]]]][[x]]]}]] - Min[Table[ SequencePosition[ PData[[m]][[list]], {Data4[[Result[[m]]]][[x]][[n]]}][[1, 1]], {n, 1, Length[Data4[[Result[[m]]]][[x]]]}]], {x, 1, Length[Result[[m]]]}]] 这是我自己定义的一个函数,Data4,PData,Result都是一些数据集(比较大),这个函数实现起来非常低效,对于每个自变量来说要运行0.047s之多
  • 0
    tmr是输出函数,afa是变量,输出曲线为什么不光滑? 代码如下: Clear["`*"] vh = 2*(1.88 - 0.59*Cos[Pi/(1 + d1/(0.65*10^(-9)))]); u = vh*0.89; mu = 4.6; gama = -4*10^(-29)*(1 - Cos[6*afa])*(1 - 1*10^(-8)*I); d1 = 1.95*10^(-9); d2 = 1.95*10^(-9); HBAR = 1.05457266*10^(-34); ME = 9.1093897*10^(-31); ELEC = 1.60217733*10^(-19); Kh = 2.29*10^(10); kc = Sqrt[2*ME*ELEC/HBAR^2]; k := kc*Sqrt[mu] k0 := Sqrt[k^2 - k^2 Sin[x]^2] kh := Sqrt[k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2] khg := Sqrt[ k^2 - (2*Kh*Sin[afa/2]*Sin[afa/2] - k Sin[x] Cos[y])^2 - (2*Kh*Si
    fhn2021 7-26
  • 6
    输入TriangleCenter[{{x1, y1 - 1}, {x1, y1}, {x1 + 1, y1 + 1}}, "Circumcenter"] 提示"{{x1,-1+y1},{x1,y1},{1+x1,1+y1}} 不是一个有效的 Circumsphere 指定." 把"Circumcenter"换成其他的(如"Incenter")都可以进行符号运算。 注意到{{x1, y1 - 1}, {x1, y1}, {x1 + 1, y1 + 1}}是一个非退化的三角形,所以外心的方程是有解的。 求助这是为什么,怎么解决。
  • 4
    比如这种,上面的符号点击了可以翻页,用什么函数添加在面板中啊
    247098254 7-25
  • 7
    需要求这个函数的临界点,求大神指导 h[q, r] = -3/(4*Pi)* Log[2] + (2*Log[1 - r^2] + Log[1 - 2*(J - r^2)])/(4* Pi) - (Log[((r^2)*(1 - Cos[2*q]))/(1 + r^4 - 2*r^2*Cos[2*q])] + 2*Log[(2*J - r^2 - 2*r*Sqrt[2*(J - r^2)]* Cos[q])/(2*((1 + 2*r^2 (J - r^2) - 2*r*Sqrt[2*(J - r^2)]*Cos[q])))])/(4*Pi) /. J -> 0.456; Solve[\!\( \*SubscriptBox[\(\[PartialD]\), \(q\)]\(h[q, r]\)\) == 0 && \!\( \*SubscriptBox[\(\[PartialD]\), \(r\)]\(h[q, r]\)\) == 0, {q, r}]
  • 13
    例如,下式在有理数域是不能分解因式的: Factor[-1 + 8 x^2] 但是在实数域内可以分解成:-(1/2) (Sqrt[2] - 4 x) (
  • 7
    DSolve[ { -y''[t] == (G*m^2)/(x[t]^2 + y[t]^2)*y[t]/Sqrt[x[t]^2 + y[t]^2], -x''[t] == (G*m^2)/(x[t]^2 + y[t]^2)*x[t]/Sqrt[x[t]^2 + y[t]^2] }, {x[t], y[t]}, t]
    jyg_1 7-24
  • 2
    把PolarPlot换成Plot就有图像。
    testzcq 7-24

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

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

扫二维码下载贴吧客户端

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

本吧信息 查看详情>>

小吧:小吧主共10

会员: Expr

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

友情贴吧