为什么两个地块不相同

作者: 电动少女
发布时间: 2023-12-22 12:13:28 (4月前)
转自：

2 条回复

0#
回复此人
满目山河 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> <P> 功能 <code> gradient </code> ，使用一个输出参数，计算x方向的偏导数。因此，您在第一位代码中计算的是d / dx r1，而不是d / dt r1。使用两个输出参数，它计算x和y导数。 </p> <P> 在第二位代码中，计算的导数被转置，导致交换x和t轴。 </p> <P> 以下代码修复了这两个错误： </p> <pre class="lang-matlab prettyprint-override"> <code> X=-1:.05:1; T=-1:.05:1; m=1+1*1i; [x,t]=meshgrid(X,T); r1=m.*sin(x+4.*t); [dr1dx,dr1dt]=gradient(r1); dr1dx=dr1dx/mean(diff(X)); dr1dt=dr1dt/mean(diff(T)); dr1dt_true=4.*m.*cos(x+4.*t); figure subplot(1,2,1) surf(x,t,abs(dr1dt)); xlabel('x') ylabel('t') title('discrete approximation') subplot(1,2,2) surf(x,t,abs(dr1dt_true)); xlabel('x') ylabel('t') title('analytical') </code> </pre> <P> <a href="https://i.stack.imgur.com/gHAhA.png" rel="nofollow noreferrer"> <img src =“https://i.stack.imgur.com/gHAhA.png”alt =“以上代码输出”/> </A> </p> </DIV>

编辑

登录后才能参与评论