求教,这个为什么总是报错?13.1版本
mathematica吧
全部回复
仅看楼主
level 6
方程不是原本的方程,我为了debug删得只剩这么点了,居然还报错,删到只剩一边了也报错,到底为什么啊?我看这个报错的官方解释就是自变量的顺序要一致,比如f[x, t]不能写成f[t, x],但我这没有这个错误啊
DSolve[{D[c[t, \[Zeta][t]], t] ==
D[c[t, \[Zeta][t]], \[Zeta][t]]}, c, {t, \[Zeta][t]}]
DSolve[{D[c[t, \[Zeta][t]], t] == 1}, c, {t, \[Zeta][t]}]
DSolve[{1 == D[c[t, \[Zeta][t]], \[Zeta][t]]}, c, {t, \[Zeta][t]}]
2022年12月23日 14点12分 1
level 9
, \[Zeta][t]并不应在D算符内出现,他只应在原函数中被表征,要不然可想而知D[c[t, \[Zeta][t]], \[Zeta][t]]时要求zeta函数的反函数
2022年12月24日 13点12分 4
那怎么写啊?原函数是 \[Zeta][t]和t的“二元”函数
2022年12月25日 08点12分
@加加加藤惠 那原函数就是t的函数,y=x^2+x+1你对x+1求导,链式法则感觉怪怪的,至少我没算过这种奇怪的东西,链式法则
2022年12月28日 08点12分
……看了3遍不知道你在表达啥,别的不说,D[c[t, \[Zeta][t]], \[Zeta][t]]这代码你自己执行过了吗?如果执行过了,你觉得这个输出意味着什么?
2023年01月07日 02点01分
吧务
level 15
……你是不是在搞移动边界问题?是的话把这帖看了,注意里面链的帖也要看:
mathematica.stackexchange.com/q/211080/1871
2023年01月07日 03点01分 5
不明觉厉,看起来好像是要自己构建个函数用来以函数值算出导数值,然后NDSolve就会自己把整个插值函数一步一步爬出来。确实有点厉害,以前我还以为NDSolve不能处理带有这种看起来不能进行符号分析的自定义函数的方程
2023年01月08日 09点01分
@加加加藤惠 ……稍微补两句吧。使用Mathematica解一维动边界问题目前至少有两个可行思路:1. 换坐标系,把动边界变成静边界。2. 直接用移动网格。这两个思路都有效利用了符号计算的优势,具体请看上面链接里的代码。
2023年02月04日 03点02分
@xzcyr ……补一句,后来经过思考,我发现我之前的移动网格实现是错的。(这个思路也不是完全不可行,但是非要写的话,最后只会弄出和换系差不多的东西来。)所以,换系法当前差不多是唯一现实的求解思路。
2024年04月06日 11点04分
1