实验2M文件及流程控制语句〔参考解析〕.docVIP

实验M文件流程控制语句 编写一个函数文件，用于求，然后在命令文件中调用。建立函数M文件.m function s=fun(a) s=0; for i=1:10 s=s+a^n; end MATLAB命令： 编写一个函数文件，求小于任意自然数n的斐波那契数列各项。Fibnacci数列定义如下： M文件fib.m function f=fib(n) if n==1 f=1; end if n==2 f(1)=1; f(2)=2; end if(n2) f(1)=1; f(2)=2; for i=3:n f(i)=f(i-1)+f(i-2); end end 编制一个解数论问题的函数文件：取任意整数，若是偶数，则除以2，否则乘3加1，重复此过程，直到整数变为1。M文件collatz.m function c=collatz(n) % collatz % Classic “3n+1” Ploblem from number theory c=n; while n1 if rem(n,2)==0 %返回n/2的余数 n=n/2; else n=3*n+1; end c=[c n]; end 编写一个函数M文件（[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.m function [y1,y2]=fun(x1,x2) if nargin=1nargout=1 y1=x1; elseif nargin=1nargout=2 y1=x1/2; y2=y1; elseif nargin=2nargout=1 y1=x1+x2; elseif nargin=2nargout=2 y1=(x1+x2)/2; y2=y1; end A=rand(3,4)0.7, 编程实现查找矩阵A的每行中第一个0元素所在的列。将结果存放着一个列向量中。 建立脚本M文件 A=rand(3,4)0.7 result=zeros(3,1) for i=1:3 for j=1:4 if A(i,j)==0 result(i)=j break; end end end result （个人所得税纳税问题）根据中华人民共和国个人所得税法规定，公民的工资、薪金所得应该依法缴纳个人所得税。个人所得税计算公式为：在每个人的月收入中超过1600元以上的部分应该纳税，这部分收入称为应纳税所得额。应纳税所得额实行分段累积税率，按如下税率表计算 级数 全月应缴税所得额 税率% 1 2. 3 4 5 6 7 8 9 500 500~2000 2000~5000 5000~20000 20000~40000 40000~60000 60000~80000 80000~100000 100000 5 10 15 20 25 30 35 40 45 设月收入为x元，应缴纳款为y元，求 y和x的函数关系，并编写函数M文件实现此函数。 调用函数，求月收入为6850和12300的应纳税额。 分析：当时不用缴税，y=0 当时，纳税部分是x-1600,税率是5%，所以 当时，500元按5%纳税，应缴纳25元，再多的部分，即x-2100，税率是10%，所以，以此类推可得函数关系式： 0 1 2 3 4 5 6 7 8 9 建立函数M文件tax.m function y=tax(x) %points 收入分段点，t收入分段点应缴纳的税额 rates收入段税率points=[0,1600,2100,3600,6600,21600,41600,61600,81600,1