实验2M文件和流程控制语句[参考答案].doc

《实验2M文件和流程控制语句[参考答案].doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验2M文件和流程控制语句参考答案1、编写一个函数文件fun.m，用于求，然后在命令文件中给定a的值，调用函数fun。建立函数M文件fun.mfunctions=fun(a)s=0;fori=1:10s=s+a^n;endMATLAB命令：>>a=3;>>s=fun(a)2、编写一个函数文件，求小于任意自然数n的斐波那契数列各项。Fibnacci数列定义如下：建立函数M文件fib.mfunctionf=fib(n)ifn==1f=1;endifn==2f(1)=1;f(2)=2;endif(n>2)f(

2、1)=1;f(2)=2;fori=3:nf(i)=f(i-1)+f(i-2);endend3、编制一个解数论问题的函数文件：取任意整数，若是偶数，则除以2，否则乘3加1，重复此过程，直到整数变为1。建立函数M文件collatz.mfunctionc=collatz(n)%collatz%Classic“3n+1”Ploblemfromnumbertheoryc=n;whilen>1ifrem(n,2)==0%返回n/2的余数n=n/2;elsen=3*n+1;endc=[cn];end1、编写一个函数M

3、文件（[y1,y2]=fun(x1,x2)），使之可以处理1个或两个输入参数，一个或两个输出参数，满足如下条件：当只有一个输入参数x1时：如果只有一个输出参数y1，则y1=x1;如果有两个输出参数y1,y2，则y1=y2=x1/2;当有两个输入参数x1,x2时：如果只有一个输出参数y1，则y1=x1+x2;如果有两个输出参数y1,y2，则y1=y2=(x1+x2)/2;建立函数M文件fun.mfunction[y1,y2]=fun(x1,x2)ifnargin=1&nargout=1y1=x1;else

4、ifnargin=1&nargout=2y1=x1/2;y2=y1;elseifnargin=2&nargout=1y1=x1+x2;elseifnargin=2&nargout=2y1=(x1+x2)/2;y2=y1;end2、A=rand(3,4)<0.7,编程实现查找矩阵A的每行中第一个0元素所在的列。将结果存放着一个列向量中。建立脚本M文件A=rand(3,4)<0.7result=zeros(3,1)fori=1:3forj=1:4ifA(i,j)==0result(i)=jbreak;end

5、endendresult1、（个人所得税纳税问题）根据中华人民共和国个人所得税法规定，公民的工资、薪金所得应该依法缴纳个人所得税。个人所得税计算公式为：在每个人的月收入中超过1600元以上的部分应该纳税，这部分收入称为应纳税所得额。应纳税所得额实行分段累积税率，按如下税率表计算级数全月应缴税所得额税率%12.3456789<500500~20002000~50005000~2000020000~4000040000~6000060000~8000080000~100000>100000510152025

6、30354045设月收入为x元，应缴纳款为y元，求1）y和x的函数关系，并编写函数M文件实现此函数。2）调用函数，求月收入为6850和12300的应纳税额。分析：当时不用缴税，y=0当时，纳税部分是x-1600,税率是5%，所以当时，500元按5%纳税，应缴纳25元，再多的部分，即x-2100，税率是10%，所以，以此类推可得函数关系式：0123456789建立函数M文件tax.mfunctiony=tax(x)%points收入分段点，t收入分段点应缴纳的税额rates收入段税率points=[0,1

7、600,2100,3600,6600,21600,41600,61600,81600,101600,inf];t=[0,0,25,175,625,3625,8625,14625,21625,29625];rates=[0,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45];n=length(points);p=1;%收入x在points中的位置fori=1:n-1ifx>points(i)&x<=points(i+1)p=i;break;endendy=t(p)+(x-p

8、oints(p))*rates(p);y=tax(6850)