];y=[];k=1;Рfor q=2:nР for p=1:q-1Р if gcd(q,p)==1Р x(k)=p/q;Р y(k)=1/q;Р k=k+1;Р endР endРendРplot(x,y,'.b');Рaxis([0,1,0,1])Рn=100是图像如下:Р异常情况说明:这道题只是在书中已有的程序上附加n=input('Enter n,negative quits:')这一语句,可以改变不同的n值来得到图像,只是n的数值不能太大,否则计算机运行的速度会很慢。在这里,还可以很直观的感受黎曼函数,分析它的众多解析性质。Р实验结果报告及实验总结:在这几个程序中,其实都还有别的作图法,只是符号作图简单快捷,便于操作。Р第一题中,把两个三维图形放在一个坐标系里画出,不仅美观,也告诉我们作图时一定要注意参数的取值范围;Р第二题墨西哥帽子花费时间较长,尤其是第二个。这里用到了subplot函数,它把一个窗口划分成几个子区域,便于在一个窗口下分析比较。在这里,同样体现了不同的参数范围会影响所绘出的图形,其次我也明白了合作的意义,这道题就很能说明问题;Р第三题稍微简单,在这里可以更好地理解解析几何中的空间曲面,matlab的强大之处仅限于此;Р六题是练习书本上已有的知识,但在编写程序的过程中还不能完全脱离书本,在以后的学习中仍需要更加努力。Р思考与深入:Р这几道题,虽然是做完了但有些是和同学交流之后才有的解决方案,课后要加强复习这方面的工作;另外,matlab把几何中难以绘制的图形给直观的展现了出来,让我对于这几个函数图形有了较为深刻的理解。以后,随着学习的深入,我可以在几何方面不能理解的函数或图形编写成新的函数用于学习,这样,对我几何和matlab这个软件的掌握都大有裨益。对于自己不能解决的问题,首先要查阅相关知识,让不能解决的话,可以向老师或同学寻求帮助。Р教师评语:
全文预览
