数学实验第五讲-函数的建立解析.ppt

Mathematics Laboratory 阮小娥博士 数学实验 阮小娥教授 Fall, 2014 应数31,信计31,统计31 Email:wruanxe@xjtu.edu.cn Tel1.函数的建立 建立函数命令1： function 1 function y=f1(x) %声明建立一个名为f1的函数 y=x^2-3; %建立函数，x可以为向量 function [输出变量列表]=函数名（输入变量列表） 建立M-函数由function语句引导，基本格式： 适用范围： 适用于多次调用且比较复杂的程序 函数文件不能独立运行！ 保存时文件名必须与函数名相同！ 996 +699=1695 1695 +5961=7656 7656 +6567=14223 14223 +32241=46464 例1. 回文数 对于一个自然数，若将各位数字倒序排出，加到原来的数字上，反复这样多次后，若能得到一个从左到右读与从右到左读完全一样的数，则称该自然数能产生回文数或者对称数。 例如，996 可以产生一个回文数 46464，如下 通过编程计算，找出102-199之间之间的数及其产生的回文数，显示产生过程，并计算这样的数占总数个数的比例。 2 解: 将一个数各位数字倒序排出，并且原数字与其倒序数均要参与运算，其结果又可能被多次调用。因此，将求一个数的倒序数建立为某一函数，然后调用即可。 1695字长为4， 1695=1*10^(3)+6*10^2+9*10+5， 5961=5*10^(3)+9*10^2+6*10+1 5=mod(1695,10)，即1695除以10后的余数； 9=mod(1695-5)/10,10),即169除以10后的余数； 依次类推，可用循环语句求各位的数； 从而可求出1695的倒序数为： Huiwenshu1.m 3 分析如下： % bnumber.m %This is a function to reach a backward number function [yy]=bnumber(y) i=1; z(i)=y; m(i)=mod(z(i),10); n(i)=floor(z(i)-m(i))/10; y1=m(i); y2=m(i); minus=y-y1; yy=y2; while minus~=0i=10 i=i+1; z(i)= n(i-1); m(i)=mod(z(i),10); n(i)=floor((z(i)-m(i))/10); y1=m(i)*10^(i-1)+y1; minus=y-y1; end y2=m(i); for j=1:i-1 y2=m(i-j)*10^(j)+y2; end yy=y2; 4 %This is to pick up and display symmetric numbers N1=0; N2=0; N3=0;Nsyn=[]; ssyn=[]; ymin=102; ymax=199;kmax=10; for y=ymin:ymax k=1; y1(1)=y; y2(1)=bnumber(y1(1)); minus=y1(1)-y2(1); if minus==0 ssyn=[ssyn,y]; N1=N1+1; else while minus~=0k=kmax k=k+1; y1(k)=y1(k-1)+y2(k-1); y2(k)=bnumber(y1(k)); minus=y1(k)-y2(k); end % End of pickup 5 %Begin of printing symmetric-oriented numbers if kkmax N2= N2+1; k=1; y1(1)=y; y2(1)=bnumber(y1(1)); minus=y1(1)-y2(1); while minus~=0 k=k+1; y1(k)=y1(k-1)+y2(k-1);