- 积分
- 4
- 注册时间
- 2006-9-12
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2009-2-25 18:28:02
|
显示全部楼层
来自 北京
本帖最后由 tianben 于 2009-2-26 16:22 编辑
功夫不负有心人,还是找到了,而且问题也搞定了。
matlab中有这个命令:EllipticPi。现在我把练习椭圆积分的小程序分享给大家,相关命令的使用可以参考附件说明。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 关于椭圆积分的学习, 参考 'mfunlist' and 'mhelp function'
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %
k=0.4;
func1=(1-k^2*sin(phi)^2)*sqrt(1-k^2*sin(phi)^2); % 相当于n=k^2
result1=int(1/func1,phi,0,pi/2) % =25/21*EllipticE(2/5),当n=k^2时,第三类椭圆积分(EllipticPi)的形式可以转化为第二类(EllipticE)。
% result11=vpa(result1,5)
result11=vpa(result1)
result2=mfun('EllipticPi',k^2,k)
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %
k=0.4; % 椭圆积分的模
n=0.7; % 椭圆积分的特征数
func1=(1-n*sin(phi)^2)*sqrt(1-k^2*sin(phi)^2);
result1=int(1/func1,phi,0,pi/2) % =EllipticPi(7/10,2/5),result1与result2结果相同
% result11=vpa(result1,5)
result11=vpa(result1)
result2=mfun('EllipticPi',n,k)
% % % % % % % % % % % % % % % % % % % % % % % % %第一类椭圆积分EllipticK
k_ellip=0.4;
ellipint=mfun('EllipticK',k_ellip)
ellipint=ellipke(k_ellip^2)
[K2,E2] =ellipke(k_ellip^2)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
评分
-
2
查看全部评分
-
|