>> a=100:1:999; %定义矩阵
>> b1=rem(a,10); %求取矩阵a中每个元素的个位,并形成和a一样的同型矩阵
>> b2=rem(fix(a/10),10); %求取矩阵a中每个元素的十位,并形成和a一样的同型矩阵
>> b3=fix(a/100); %求取矩阵a中每个元素的百位,并形成和a一样的同型矩阵
>> k=find(a==b1.*b1.*b1+b2.*b2.*b2+b3.*b3.*b3) %个位矩阵每个元素的三次方 加上 十位矩阵每个元素的三次方 加上 百位矩阵每个元素的三次方,获得一个总矩阵与a矩阵比较。将相等的序号为保存到k中。
k =
54 271 272 308
>> c=a(k) %通过k序号位,获得对应数据
c =
153 370 371 407
Matlab:水仙花数
最新推荐文章于 2025-04-05 08:09:51 发布