内容发布更新时间 : 2024/9/30 12:33:08星期一 下面是文章的全部内容请认真阅读。

《数学实验》报告一

学号 2013530250243 号 姓名 刘金花 成绩 实验内容： 实验1.3 M文件与程序设计

一、实验目的 通过本实验了解MATLAB中两种形式的M文件，掌握利用MATLAB的三种程序结构编写简单程序，处理实际问题的基本方法。

二、实验内容与要求

1.建立函数

的M文件，并计算f(10)和f(1000).

结果 f(10)=0.0102 f(1000)=1.0715e+286 Matlab命令 x=input('请输入x的值:'); f(x)=2^x/(x^5+1); disp(f(x)) 3.给定如下的数列：1,1,2,3,5,8,13,21,34,55,89，…，每一项都是前两项的和，其递推公式为：F1?1，F2?1，Fn?2?Fn?1?Fn，该数列称为斐波那契（Fibonacci）数列. 请你编写一个程序，求这个数列的每一项，并验证其通项公式为：

nn?????11?51?5??????. ??Fn????5??2??2?????Matlab命令 首先建立一个函数文件fib.c 程序如下： function [f] = fib(n) 结果 Uc.m Input x please:7 f= if n==1 f=1; elseif n==2 f=[1 1]; else b=fib(n-1); f=[b,b(end-1)+b(end)]; end end 然后建立一个脚本文件uc.m 程序如下： n=input('Input x please:'); a=fib(n) disp(a) 然后在matlab命令窗口调用该函数文件，得到斐波那契数列的每一项 在建立一个脚本文件u1c.m，验证斐波那契数列的通项公式 程序如下： n=input('Input n please:'); a=(1/sqrt(5))*[((1+sqrt(5))/2)^n-((1-sqrt(5))/2)^n] disp(a) 最后在matlab命令窗口调用该函数文件，观察所得结果与前面所得结果的最后一项是否一致，如果一致，则证明该通项公式就是斐波那契数列的通项公式。 1 1 2 3 5 8 13 1 1 2 3 5 8 13 U1c.m Input n please:7 F= 13.0000 13.0000 因为所得结果13与前面所得结果的最后一项一致，所以该该通项公式就是斐波那契数列的通项公式。 4.《中华人民共和国个人所得税法》规定：公民的个人工资、薪金应依法缴纳个人所得税. 所得税计算办法为：在每个人的月收入中超过3500元的部分应该纳税，这部分收入称为应纳税所得额. 应纳税计算公式为

工资、薪金所得个人 =应纳税所得额×适用税率—速算扣除数

所得税应纳税额 适用税率按表1.7计算：

等 级 1 2 3 4 5 6 7 表1.7 个人所得税税率表（适用于工资、薪金） 全月应纳税所得额 税率（%） 速算扣除数 3 0 不超过1500元的部分 10 105 超过1500元，不到4500元的部分 20 555 超过4500元，不到9000元的部分 25 1005 超过9000元，不到35000元的部分 30 2755 超过35000元，不到55000元的部分 35 5505 超过55000元，不到80000元的部分 45 13505 超过80000元的部分 若某人的月收入为x，试建立税款y与收入x之间的函数M文件.

Matlab命令 程序如下： 首先建立一个脚本文件fc.m x=input('请输入x的值:'); if x<=1500 y=x*0.03-0; elseif 1500