首页 > 其他分享 >用Mathematica和SciPy阐明Jacobi椭圆函数的定义方法

用Mathematica和SciPy阐明Jacobi椭圆函数的定义方法

时间:2023-06-08 19:32:47浏览次数:60  
标签:cn 函数 Mathematica 模数 SciPy 参数 定义方法 Jacobi sn


这,这个,那,那个Jacobi椭圆函数SN和CN类似于三角函数正弦和余弦。它们出现在非线性振动和保形映射等应用中。不幸的是,定义这些函数有多种约定。这篇文章的目的是澄清围绕这些不同公约的混淆。

用Mathematica和SciPy阐明Jacobi椭圆函数的定义方法_Mathematica

上面的图像是函数sn[1]的一个图。

模量、参数和模数角

Jacobi函数有两个输入。我们通常认为Jacobi函数是第一个输入的函数,第二个输入是固定的。这第二个输入是一个“拨号”,你可以转动它来改变他们的行为。

有几种方法可以指定此拨号。我以“拨号”一词开始,而不是“参数”,因为在这个上下文中参数具有技术意义,一种描述刻度盘的方法。除了“参数”,您还可以将Jacobi函数描述为模数模角。这篇文章将是一种罗塞塔石头,展示了描述雅可比椭圆函数的每一种方式是如何相关的。

这,这个,那,那个参数m是[0,1]中的实数。这,这个,那,那个互补参数m' = 1 - m艾布拉莫维茨和斯特贡例如,将Jacobi函数sn和cn编写为sn(u | m)和CN(u | m)他们也用m1=而不是m‘表示互补参数。

这,这个,那,那个模数k的平方根m。如果m代表模数,但这不是传统。相反,m表示参数和k是模数。这,这个,那,那个互补模量k‘是互补参数的平方根。

这,这个,那,那个模角α是由m=SIN 2α。

注意,作为参数m等于零,模数也是零。k以及模块角α。当这三个函数中的任意一个变为零时,Jacobi函数sn和cn收敛到它们对应的正弦和余弦。因此,无论您的刻度盘是参数、模数还是模角,当您将刻度盘转向零时,sn收敛到正弦,cn收敛到余弦。

作为参数m等于1,模数也是1。k,但模块化角度α进入π/2。所以如果你的刻度盘是参数或者模数,它会变成1。但是,如果您认为您的拨号是模块化的角度,它将进入π/2。在这两种情况下,当你将刻度盘向右转时,sn收敛到双曲割线,cn收敛到常量函数1。

季度期

除了参数、模数和模角外,还可以看到Jacobi函数KK“”这些被称为季度期是有充分理由的。函数sn和cn有周期4。K当你沿着真实的轴移动,或者在复杂平面的任何地方水平移动。他们也有第四期IK“”也就是说,当移动距离为4时,函数会重复。K‘垂直[2]。

四分之一周期是模数的函数。季度期K沿着真实的轴线

用Mathematica和SciPy阐明Jacobi椭圆函数的定义方法_ci_02

当我们上面说Jacobi函数有4周期的时候K,这是从变量的角度来看的。u。请注意,当φ=π/2,u = K.

数学中的Jacobi椭圆函数

Mathematica使用u第一个参数的约定和第二个参数的参数约定。

数学函数JacobiSN[u, m]用参数计算函数snu参数m。在A&S的表示法中,sn(u | m).

同样,JacobiCN[u, m]用参数计算函数cnu参数m。在A&S的表示法中,cn(u | m).

到目前为止,我们还没有讨论过Jacobi函数DN,但它是在Mathematica中实现的JacobiDN[u, m].

将振幅φ作为uJacobiAmplitude[um m].

计算季度期间的函数。K从参数mEllipticK[m].

Python中的Jacobi椭圆函数

SciPy库有一个Python函数,它同时计算四个数学函数。功能scipy.special.ellipj有两个论点,um,就像Mathematica一样,并返回sn(u | m),CN(u | m)、DN(u | m),以及振幅φ(um).

功能Km)在Python中实现为scipy.special.ellipk.

相关员额

[1]情节是用JacobiSN[0.5, z]和功能ComplexPlot所述这里.

[2]严格地说,4IK“是”a期间。这是CN的最小垂直周期,但2IK是sn的最小垂直周期。

标签:cn,函数,Mathematica,模数,SciPy,参数,定义方法,Jacobi,sn
From: https://blog.51cto.com/u_16145034/6442459

相关文章

  • 自定义方法,实现按指定属性进行去重
    实例:根据age、sex、name对user对象集合进行去重//构建user对象publicclassUser{privateintid;privateStringname;privateintage;privateStringsex;privatedoubleheight;//构造方法、getter和setter略}//创建对象集合List<User......
  • scipy.signal.butter实现带通滤波器
    fromscipy.signalimportbutter,lfilter#带通滤波器defbutter_bandpass_filter(data,lowcut,highcut,fs,order):low=lowcut*2/fshigh=highcut*2/fsb,a=butter(order,[low,high],btype='bandpass')y=lfilter(b,a,......
  • 自定义方法
          ......
  • LightOJ1007---Mathematically Hard (欧拉函数)
    Mathematicallysomeproblemslookhard.Butwiththehelpofthecomputer,someproblemscanbeeasilysolvable.Inthisproblem,youwillbegiventwointegersaandb.Youhavetofindthesummationofthescoresofthenumbersfromatob(inclusive).T......
  • Makefile知识点-3------makefile定义“空格”变量的定义方法
    下面再介绍两个定义变量时我们需要知道的,请先看一个例子,如果我们要定义一个变量,其值是一个空格,那么我们可以这样来:nullstring:=space:=$(nullstring)#endofthelinenullstring是一个Empty变量,其中什么也没有,而我们的space的值是一个空格。因为在操作符的右边是很难描述......
  • Mathematica学习笔记002-数据导入导出
    如果不能把数据导入导出,Mathematica就只能是个大号计算器了。学会了导入导出,一方面可以把数据、图像结果保存,另一方面也可以将别的程序的中间结果导出成(txt或xls格式),然后交给Mathematica处理,让骑完成高精度计算和绘图。基本操作其实很简单Export["D:\\abc.txt",{{1,2},{3,4......
  • light oj 1007 Mathematically Hard (欧拉函数)
    题目地址:lightoj1007第一发欧拉函数。欧拉函数重要性质:设a为N的质因数,若(N%a==0&&(N/a)%a==0)则有E(N)=E(N/a)*a;若(N%a==0&&(N/a)%a!=0)则有:E(N)=E(N/a)*(a-1)对于这题来说,首先卡MLE。。只能开一个数组。。所以把前缀和也存到......
  • Wolfram Mathematica 大师班:从初学者到专家
    WolframMathematica大师班:从初学者到专家通过实例学习:机器学习、自然科学、统计学、经济、语言学和媒体的巧妙编程课程英文名:WolframMathematicaMasterclassfromBeginnertoExpert此视频教程共8.06GB,中英双语字幕,画质清晰无水印,源码附件全课程地址:https://xueshu.fun/......
  • HDOJ1017 A Mathematical Curiosity
    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1017这个题目其实挺坑的。首先是N,应该挺多人纠结过这个N,N其实是blocks(块),一块有未知个cases。一个块的结束标志是0,0。然后是PE的问题,空格、空行,我也是被坑的好惨。这里应该是每个块之间有一个空行!也就是说,最后一个块是不......
  • 【SciPy】Sparse稀疏矩阵主要存储格式总结(转载)
    原文:【SciPy】Sparse稀疏矩阵主要存储格式总结在数据科学和深度学习等领域常会采用矩阵格式来存储数据,但当矩阵较为庞大且非零元素较少时,运算效率和存储有效率并不高。所以,通常我们采用Sparse稀疏矩阵的方式来存储矩阵,提高存储和运算效率。下面将对SciPy中七种常见的存储方式(COO/......