用极坐标画一个五角星Mathematica选读.doc

  1. 1、本文档共9页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
有谁知道怎么用极坐标画一个五角星么,? ParametricPlot[{6/25 Cos[t] + 4/25 Cos[3/2 t], 6/25 Sin[t] - 4/25 Sin[3/2 t]}, {t, 0, 4 Pi}, Axes - False]就知道参数方程有这个 PolarPlot[ 1/Cos[2 \[Pi]/5 - Mod[\[Theta], 4 \[Pi]/5]], {\[Theta], 0, 4 \[Pi]}] 萌新准备做一个迈克尔逊干涉条纹的动态模拟。函数如下I[x,y]=2 + 2 Cos[(4 \[Pi])/(632.8*10^-9)*7*10^-6*Cos[Sqrt[x^2 + y^2]/0.15] DensityPlot[2+2 Cos[(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+2 Cos[(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}],{n,400,760,10}] Manipulate[ DensityPlot[ 2 + 2 Cos[(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 + 2 Cos[(4 \[Pi])/(632.8*10^-9)*7*10^-6* Cos[Sqrt[x^2 + y^2]/0.15]], {x, -0.06, 0.06}, {y, -0.06, 0.06}] DensityPlot[ 2 + 2 Cos[(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 + 2 Cos[(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}]; ListAnimate[movies] 另外由于你这个函数震荡比较严重,所以默认的绘图参数效果不是很好,增加PlotPoints或MaxRecursion选项 可以提高绘图质量,效果如下,但相应的绘图时间也会大大增加,你可以尝试一下,然后自己取舍【图片】 至于换成极坐标我觉得没有太大意义,/questions/3680/how-do-i-draw-a-densityplot-in-polar-coordinates 这里是SE上的一个相关讨论,你可以看下 如果只是针对这个函数的话因为变化只和r有关,所以可能会有一些针对性的优化方法,比如 用ListDensityPlot指定取点位置之类的,但是点的位置指定起来有点费劲,取多了速度太慢,取少了效果又不好,我还没折腾利索所以就先不发了 试着搞了一个出来,大致思路是找到函数变化比较剧烈的地方,然后用ListDensityPlot画出来。【图片】【图片】 可以看到快了5倍左右,而且质量还算能忍了。代码如下: myPlot[n_] := Module[{f, rlist, data}, f[r_] := 2 + 2 Cos[(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, Infinity]~ Join~{0, 0.2} // Sort; data = Table[{r Cos[\[Theta]], r Sin[\[Theta]], f[r]}, {r,

文档评论(0)

希望之星 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档