用极坐标画一个五角星 mathematica

用极坐标画一个五角星 mathematica

ID:34482201

大小:4.12 MB

页数:9页

时间:2019-03-06

用极坐标画一个五角星 mathematica_第1页
用极坐标画一个五角星 mathematica_第2页
用极坐标画一个五角星 mathematica_第3页
用极坐标画一个五角星 mathematica_第4页
用极坐标画一个五角星 mathematica_第5页
资源描述:

《用极坐标画一个五角星 mathematica》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库

1、实用文案有谁知道怎么用极坐标画一个五角星么,?ParametricPlot[{6/25Cos[t]+4/25Cos[3/2t],6/25Sin[t]-4/25Sin[3/2t]},{t,0,4Pi},Axes->False]就知道参数方程有这个PolarPlot[1/Cos[2[Pi]/5-Mod[[Theta],4[Pi]/5]],{[Theta],0,4[Pi]}]萌新准备做一个迈克尔逊干涉条纹的动态模拟。函数如下I[x,y]=2+2Cos[(4[Pi])/(632.8*10^-9

2、)*7*10^-6*Cos[Sqrt[x^2+y^2]/0.15]DensityPlot[2+2Cos[(4Pi/(632.8*10^(-9)))*7*10^(-6)*Cos[(Sqrt[x^2+y^2])]/0.15],{x,-0.2,0.2},{y,-0.2,0.2}]Manipulate[DensityPlot[2+2Cos[(4Pi/(n*10^(-9)))*7*10^(-6)*Cos[(Sqrt[x^2+y^2])]/0.15],{x,-0.06,0.06},{y,-0.06,0.06}

3、],{n,400,760,10}]标准实用文案Manipulate[DensityPlot[2+2Cos[(4[Pi])/(n*10^-9)*7*10^-6*Cos[Sqrt[x^2+y^2]/0.15]],{x,-0.06,0.06},{y,-0.06,0.06}],{n,400,760,10}]标准实用文案DensityPlot[2+2Cos[(4[Pi])/(632.8*10^-9)*7*10^-6*Cos[Sqrt[x^2+y^2]/0.15]],{x,-0.06,0.06},{y,-

4、0.06,0.06}]DensityPlot[2+2Cos[(4[Pi])/(632.8*10^-9)*7*10^-6*标准实用文案Cos[Sqrt[x^2+y^2]/0.15]],{x,-0.2,0.2},{y,-0.2,0.2}]播放中图像错误是因为在播放的两帧之间时间太短,电脑来不及进行精确的计算,所以图像模糊,播放停止后电脑进行较为精确的计算,图像就恢复了解决方法是提前绘制好图像然后用ListAnimate播放,代码如下movies=Table[DensityPlot[2+2Cos[(4

5、[Pi])/(n*10^-9)*7*10^-6*Cos[Sqrt[x^2+y^2]/0.15]],{x,-0.06,0.06},{y,-0.06,0.06}],{n,400,760,10}];ListAnimate[movies]另外由于你这个函数震荡比较严重,所以默认的绘图参数效果不是很好,增加PlotPoints或MaxRecursion选项可以提高绘图质量,效果如下,但相应的绘图时间也会大大增加,你可以尝试一下,然后自己取舍【图片】标准实用文案至于换成极坐标我觉得没有太大意义,http:/

6、/mathematica.stackexchange.com/questions/3680/how-do-i-draw-a-densityplot-in-polar-coordinates这里是SE上的一个相关讨论,你可以看下如果只是针对这个函数的话因为变化只和r有关,所以可能会有一些针对性的优化方法,比如用ListDensityPlot指定取点位置之类的,但是点的位置指定起来有点费劲,取多了速度太慢,取少了效果又不好,我还没折腾利索所以就先不发了试着搞了一个出来,大致思路是找到函数变化比较剧烈的

7、地方,然后用ListDensityPlot画出来。【图片】【图片】标准实用文案可以看到快了5倍左右,而且质量还算能忍了。代码如下:myPlot[n_]:=Module[{f,rlist,data},f[r_]:=2+2Cos[(4[Pi])/(n*10^-9)*7*10^-6*Cos[r/0.15]];rlist=Cases[Normal@Plot[f'[r]f''[r],{r,0,0.2},MeshFunctions->(#2&),Mesh->{{0}}],Point[{x_,y_}]:>x,

8、Infinity]~Join~{0,0.2}//Sort;data=Table[{rCos[[Theta]],rSin[[Theta]],f[r]},{r,标准实用文案rlist},{[Theta],0,2[Pi],[Pi]/30.}]~Flatten~1;ListDensityPlot[data]];myPlot[632.8]//AbsoluteTiming大致思路就是这样,如果有进一步的需求比如提高绘图质量或者更改绘图范围之类的也可以在这基础上继续修改。顺便体格问题,

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。