- 积分
- 33
- 注册时间
- 2002-7-6
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2006-10-20 16:02:58
|
显示全部楼层
来自 湖南长沙
不仿把下面的回答作为精华贴(7)
注意这里sin对应有很简单的反函数arccos,对于反函数不简单的情况呢?
u[r_,θ_]:=3/8-1/2 r^2 Cos[2 θ]+1/8 r^4Cos[4 θ];
v[x_,y_]:=Module[{r,θ},r=Sqrt[x^2+y^2];θ=If[ArcTan[x,y]>0,ArcTan[x,y],2π+ArcTan[x,y]];u[r,θ]//Evaluate];
a=Table[(i-10)/10,{i,0,20}]//N;
c={-.8,-.5,-0.1,0.1,.5,.8};
x1=Table[#,{i,Length[a]}]&/@c//Flatten;
y1=Nest[Join[#,a]&,a,Length[c]-1];
b=y/.MapThread[FindRoot[v[#2,y]==#1,{y,-1,1}]&,{y1,x1}];
z1=ContourPlot[v[x,y]//Evaluate,{x,-1.1,1.1},{y,-1.1,1.1},PlotPoints->1000,Contours->a,ContourShading->False,Epilog->MapThread[Text[ToString[#1],{#3,#2},Background->GrayLevel[1.0]]&,{y1,b,x1}]];
得到:
[ 本帖最后由 scqxaie 于 2006-10-20 16:19 编辑 ] |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|