-
-
4
-
5
-
1现在是手动二分法,输入两个自变量x,y,寻找输出的两个量为特定的值时与之对应的xy,现在想能不能在mma里或者python中二分法求解x,y。已知函数单调,和xy大致范围,不会出现不同xy对应同一个值
-
2c = 0.5 \[Alpha] = 0.9 d = 4 k = 4 \[Beta] = ( d (4 - \[Alpha]) (\[Alpha] - 2) + 4 [3 w + c (\[Alpha] - k)] - \[Alpha] (c k + w - \[Alpha]))/( c k (4 - 3 \[Alpha])) Plot[(2 w + \[Alpha] (1 + c) + 2 k c (\[Beta] - 1))/(\[Alpha] (4 - \[Alpha])), {w, 3.7, 4.8}]
-
23此帖是本吧初代吧主 @mm_酱 发在果壳小组的一篇教程。如今果壳小组已死,mm_酱近期又不太活跃,这里姑且重发一份: https://note.youdao.com/share/?id=058e6037396d925af1f4abe4d54a52 mm_酱你要是想整理下出个v2就吱一声。
-
3
-
1吧里关心微分方程的人不多,但这个例子确实很独特,姑且提一句。我们知道(好吧,相当多的人不知道,想着自己写经典龙格库塔的人时不时就会冒出来),NDSolve默认的常微分方程求解器(细节不明——甚至可以说是NDSolve中最不透明的部分——似乎一般是LSODA方法,有时似乎也会调用StiffnessSwitching方法)性能是非常好的,在NDSolve求解出错的时候,常微分方程求解器永远应该放在最后检查。而当常微分方程求解器出现性能问题时,“使用MaxStepSize限
-
4下图是运算结果,我希望更改频率和线径来得出一个量来进行比较,假如如下图输入所示fs为3*1的矩阵,线径为1*3的矩阵,最后不应该α是一个3*3的矩阵,为啥是现在这种情况,请各位大佬帮忙看看 代码如下 fs = {70 Quantity[lbk]1, "Kilohertz"[rbk], 140 Quantity[lbk]1, "Kilohertz"[rbk], 150 Quantity[lbk]1, "Kilohertz"[rbk]};(*开关频率*) WireDs = {{0.1 Quantity[lbk]1, "Millimeters"[rbk], 0.07 Quantity[lbk]1, "Millimeters"[rbk], 0.05 Quantity[lbk]1, "Millimeters"[
-
0重积分我想用那个boole函数表示结果很多简单的都积不出来,第二形曲面积分完全不知道该怎么用mma来算
-
2
-
1
-
0
-
0症状: 顶楼中的 [ 和 ] 会被自动转换为 [lbk] 和 [rbk]。此现象在回复帖中不存在。 影响范围: 此bug至少在去年10月16日即已出现,至今(安卓12.58.1.0)未修。IPhone受影响情况不明,欢迎实验补充。 对策: 对于提问者,最简单的对策就是改用电脑访问贴吧网页版,不要使用使用手机客户端发帖。如果确有相关需求(所以说你们到底为什么会有这种需求?电脑没连网?那也可以开手机热点啊,访问网页又花不了什么流量),避免在顶楼张贴代码也是一个
-
1
-
1
-
11
-
0
-
1,求结微分方程出现这种结果,怎么处理,这是代码,{{\[lbk]Alpha[rbk][lbk]t[rbk] -> InverseFunction[lbk]-(( g^2 \[lbk]CapitalDelta[rbk] ArcTan[lbk](2 Sqrt[lbk]k[rbk] Sqrt[lbk]#1[rbk])/ Sqrt[lbk]-2 \[lbk]Delta[rbk] + g^2 \[lbk]CapitalDelta[rbk][rbk][rbk])/( 2 Sqrt[lbk]k[rbk] Sqrt[lbk]-2 \[lbk]Delta[rbk] + g^2 \[lbk]CapitalDelta[rbk][rbk])) + Sqrt[lbk]#1[rbk] &[rbk][lbk]-t \[lbk]CapitalOmega[rbk][rbk]}}
-
22链接不发了。楼下简述更新内容
-
5分别画两幅图,一幅图是实线另一幅是虚线,为什么虚线图例线段长度比实线长,之前我都是直接调整图例框让虚线看起来和实线线段长度一样
-
4
-
0我想将nb文件另存为pdf,但是这个pdf的字体怎么选呢?以前都自动是黑体,但是有一次突然变成别的体了一直没变回来,我想要黑体,应该怎么设置呢
-
13
-
14代码和运行结果如下: 在吧内尝试着找了一下没有找到相关的帖子,我自己觉得应该是与替换规则的使用或者代码的运算顺序有关系,但是不知道应该去找哪些资料来看,所以如果存在现有的资料的话就拜托大家帮忙指个路了 mKc = 0.493677; mpi0 = 0.1349768; f = 0.093; \[Epsilon] = 0.001; omega[m_, q_] := Sqrt[m^2 + q^2]; fcut[p_, cut_] := Exp[-(p/cut)^2]; g[s_, cut_] := NIntegrate[((4*Pi)/(2*Pi)^3)*((omega[mKc, q] + omega[mpi0, q])/(2* omega[mKc, q]*omega[mpi0, q]))* ((q^2* fcut[q, cut])/(s - (omega[mKc, q] + omega[mpi0
-
2明明是收敛的结果我的MMA非要说它不收敛
-
4\[CapitalOmega] = ImplicitRegion[True, {{z, 0, 200}}]; \[CapitalGamma]2 = NeumannValue[0., x == 0 || x == 200] s = NDSolveValue[{D[n[t, z], t] - D[n[t, z], {z, 2}] - \[Tau]/ nb \[Alpha] (1 - R) Pdensity/(Sqrt[2 \[Pi] ] \[Tau]L hv ) gauss[t, 0, \[Tau]L/\[Tau]] Exp[-\[Alpha] Ldi z] == -Ldi e^2 nb Ldi/(k \ T \[Epsilon] \[Epsilon]b) (0 + eF0O[z]) D[n[t, z] , z] - Ldi e^2 nb Ldi/(k T \[Epsilon] \[Epsilon]b) (p0[z] - n0[z] 0 + 1 - E^(-2 u)) n[t, z] + \[CapitalGamma]2, n[-10 \[Tau]L/\[Tau], z] == 0 }, n, {t, -10 \[Tau]L/\[Tau], 5}, {z} \[Element] \[CapitalOmega], PrecisionGoal -> 35, Method ->
-
3
-
1
-
3这个是代码 Plot[0.0172[1 + 3.93 (10^-3)[x - 20]]/10^6, {x, 20, 320}] 我尝试加入PlotRange给因变量范围在1.72*10^-8到4*10^8也不行,查了个帖子说用plot[re]指令,也不行
-
10
-
1mma 同样作为 lisp-like 语言有 Throw/Catch 有 Dialog 甚至有 Goto 却没有 call/cc 的官方实现也是够奇怪的
-
21昨天@cckinger 同学跟我提到了遗传算法,我觉得很有意思,便尝试用MMTC实现了一下.第一次玩,题目简单,高手莫笑.
-
4请问刚接触mathematica应该如何学习应用,看书还是看网课,本人大一想接触数学建模相关的比赛,求大佬指点🙏
-
1
-
2已经把公式中其他内容带入数值,运行不出结果,而且单独计算里面的积分也运行不出来 DSolve[P[\[Omega]] == \!\( \*SubsuperscriptBox[\(\[Integral]\), \(0\), \(Pi/2 - 10^\((\(-3\))\)\)]\(\((\(-\((\((66416.66137691677`\ \ \((\(-3.6551383444493544`*^14\) + \*SuperscriptBox[\(E\), \(38.683223086147535`\ \*SqrtBox[\(0.7514187629` + 0.005167396856733549`\ Cos[\[Theta]]\)]\)])\)\ \ \((\(-2\)\ \(( \*FractionBox[\(9\), \(1600\)] + \*FractionBox[\(5.618`*^21\), SuperscriptBox[\(\[Omega]\), \(2\)]])\)\ \*SuperscriptBox[\(Cos[\[Theta]]\), \(2\)] + 3.0056750516`\ Cos[ 2\ \[T
-
5题一: 已知 f(x)+f(1/x)+3f(-x)=x , 求f(x)的解析式 题二: 已知 f(f(x))=x^2-3x+4 , 求f(2)的值和f(x)的解析式
-
8
-
10
-
2Subscript[\[Pi], m]=(2 a c (t+\[Beta])-c^2 (t+\[Beta])^2+2 t T (-4+2 \[Alpha]^2+(t+\[Beta])^2 \[Gamma]^2)+((-1+\[Theta]) Overscript[a, _]-\[Theta] Q) (2 a+(-1+\[Theta]) Overscript[a, _]-\[Theta] Q))/(-8+4 \[Alpha]^2+2 (t+\[Beta])^2 \[Gamma]^2); Subscript[\[Pi], r]=-(((-2+\[Alpha]^2) (2 a-c (t+\[Beta]))+a (t+\[Beta])^2 \[Gamma]^2+(-2+\[Alpha]^2+(t+\[Beta])^2 \[Gamma]^2) ((-1+\[Theta]) Overscript[a, _]-\[Theta] Q))^2/(2 (-2+\[Alpha]^2) (-4+2 \[Alpha]^2+(t+\[Beta])^2 \[Gamma]^2)^2));A=Q*\[Theta]+(1-\[Theta])*Overscript[a, _];optw=((-2+\[Alpha]^2) (A+c (t-\[Beta]))+A t (t+\[Beta]) \[Gamma]^2)/(-4+
-
8
-
2假设我有一个关系图,以整数之间的因数关系举例: ns = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10}; as = {1 \[DirectedEdge] 2, 1 \[DirectedEdge] 3, 1 \[DirectedEdge] 4, 1 \[DirectedEdge] 5, 1 \[DirectedEdge] 6, 1 \[DirectedEdge] 7, 1 \[DirectedEdge] 8, 1 \[DirectedEdge] 9, 1 \[DirectedEdge] 10, 2 \[DirectedEdge] 4, 2 \[DirectedEdge] 6, 2 \[DirectedEdge] 8, 2 \[DirectedEdge] 10, 3 \[DirectedEdge] 6, 3 \[DirectedEdge] 9, 4 \[DirectedEdge] 8, 5 \[DirectedEdge] 10}; GraphPlot[as, PlotTheme -> "ClassicDiagram", DataRange -> {{0, 1}, {0, 1}}, GraphLayout -> "CircularE
-
15我暂时无法将这个错误缩减为一段简单的代码, 完全的代码在二楼. 问题完全出在被命名为problem的函数上, 所以只看problem这一块就好了. problem函数计算了两个积分w0和w1, 第一个积分w1在数学上的结果是1. 尝试运行测试函数 test[1] 第一个积分结果为1, 很好, test[2] 第一个积分结果就不正确了. 但! 如果退出内核, 先运行 test[2] 第一个积分结果为1, 反而后运行的 test[1] 不正确了, 也就是只能正确计算首次接受的参数. 错误可能是Δ[index_] := Δ[index] = … 这样的“有
-
7也算是对这帖(http://tieba.baidu.com/p/2964416898)的2,3,4条的解说吧。有什么意见或建议欢迎提出: http://note.youdao.com/share/?id=abd51087f44c0b6a41ff6022d549dc4
-
4下面是代码: n = 4; Fig = Table[1, {i, 1, n}]; For[l = 1, l <= n, l++, p1 = 1; p2 = p1 (1 + 10^(-l - 13)); Clear[{s1, s2}]; s1 = y /. NDSolve[{y'[x] == y[x] Cos[p1 x + y[x]], y[0] == 1}, y, {x, 0, 31}, AccuracyGoal -> 20, PrecisionGoal -> 20, WorkingPrecision -> 64][[1]]; s2 = y /. NDSolve[{y'[x] == y[x] Cos[p2 x + y[x]], y[0] == 1}, y, {x, 0, 31}, AccuracyGoal -> 20, PrecisionGoal -> 20, WorkingPrecision -> 64][[1]]; Fig[[l]] = Plot[{s1[x] - s2[x]}, {x, 0, 30}, PlotRange -> Full]; ] Fig 试了很多方法, 但是, 只要p1 和p2 的区别小于10^(-17)左
扫二维码下载贴吧客户端
下载贴吧APP
看高清直播、视频!
看高清直播、视频!