数学实验水仙花数.docx

Matlab基础知识目的:熟悉MATLAB的具体操作与操作键。掌握MATLAB中的常用函数与变量、表达式的定义方法。熟悉MATLAB常用的工作方式M文件的编程工作方式掌握MATLAB语言中的程序结构。任务建立一个M文件,求所有的水仙花数。所谓的水仙花数是指一个三位数,其各位数字的立方和等于该数本身。例如153是一个水仙花数,因为153=13+53+33。用subplot分别在不同的坐标系下画出下列四条曲线,为每幅图形加上标题:概率曲线四叶玫瑰曲线叶形线正弦曲线(1)运用subplot函数可以将多个图画到一个平面上,subplot(m,n,p),其中,m表示是图排成m行,n表示图排成n列,也就是整个figure中有n个图是排成一行的,一共m行。functionsifutux=-5::5;y1=exp(-(x.^2));j=0::2*;%概率曲线p=sin(2*j);%四叶玫瑰线t=0::100;x2=3*t./(1+t.^3)y2=3*t.^2./(1+t.^3)%叶形线a=-2*::2*;y3=sin(a);%正弦曲线subplot(2,2,1);plot(x,y1);subplot(2,2,2)polar(j,p);subplot(2,2,3)plot(x2,y2);subplot(2,2,4)plot(a,y3);问题分析:首先要先设一个三位数s,然后提取它的百位数a,十位数b,个位数c,然后根据条件s=a^3+b^3+c^3,求出s。functionshuifors=100:999;%定义s为三位数a=fix(s/100);%提取百位数b=fix(s/10-10*a);%提取十位数c=s-100*a-10*b;%提取个位数ifs==a.^3+b.^3+c.^3%s=a^3+b^3+c^3disp(s)%显示出sendend然后运行的出水仙花