matlab解第一类不完全椭圆积分方程
t1,t12,t13都已知，可以随便假定一个数，k属于0到1，然后求解k，（里面包含第一类不完全椭圆积分，matlab中有这个函数，好像是ellipticF）怎么用matlab编程啊，里面涉及到复数，对我这个matlab菜鸟来说太难了，跪求各位大神帮帮忙。。。财富不多了 就这点了==
一旦有财富值的话 我都会加上去的 求大神帮帮忙（T.T）
第1个回答 推荐于2016-11-09
假设等式左边为1/2*0.3
f=@(x)-mfun('EllipticF',asin(1),x^2)/imag(mfun('EllipticF',asin(1/x),x^2))-1/2*0.3;
fplot(f,[-10 10])
fsolve(f,0.0001)追问
f=@(x)-mfun('EllipticF',asin(1),x^2)/imag(mfun('EllipticF',asin(1/x),x^2))-1/2*0.3;
fplot(f,[-10 10])
fsolve(f,0.0001)追问
大神~解出来的k值是复数~？
追答fsolve(f,0.0001)
No solution found.
fsolve stopped because the relative size of the current step is less than the
default value of the step size tolerance squared, but the vector of function values
is not near zero as measured by the default value of the function tolerance.
ans =
3.1499e-06 + 1.5627e-05i
解出了复数，回带后为-0.15，不是0
大神~肿么感觉每次运行的结果都不一样啊~初值取不同的值~最后结果也不一样诶~
追答你试试图解吧，Mathematica 画图如下，不知左边等于什么
Plot[-EllipticF[ArcSin[1], x^2]/
Im[EllipticF[ArcSin[1/x], x^2]], {x, -1, 1}]
嗯嗯~好的~谢谢大神了~以后有问题直接请教你了~还希望大神不要嫌我烦....
本回答被提问者采纳
相似回答