2023年秋季Matlab上机作业.docx

《2023年秋季Matlab上机作业.docx》由会员分享，可在线阅读，更多相关《2023年秋季Matlab上机作业.docx（25页珍藏版）》请在第一文库网上搜索。

1、2023年秋季Mat1ab上机作业1、举例说明m命令文件、m函数文件、m文件三者之间的关系。m文件根据调用方式的不同分为两类：m命令文件和In函数文件。它们的扩展名均为.m,主要区别在于：1) .命令文件没有输入参数，也不返回输出参数，而函数文件可以带输入参数，也可返回输出参数。2) .命令文件对MAT1AB工作空间中的变量进行操作，文件中所有命令的执行结果也完全返回到工作空间中，而函数文件中定义的变量为局部变量，当函数文件执行完毕时，这些变量被清除。3) .命令文件可以直接运行，在MAT1AB命令窗口输入命令文件的名字，就会顺序执行命令文件中的命令，而函数文件不能直接运行，而要以函数调用的方

2、式来调用它。命令文件：c1eara11;c1osea11;c1c;X=1:10;y=2:5;,y=exchange(xzy)函数文件：functiona,b=exchange(a,b)end结果：X=2345y=123456789102、调查现行个人所得税政策，编制一个mat1ab程序计算个人所得税(1、按月扣缴，2、按年扣缴),并测试。按月扣缴程序：c1eara11;c1osea11;c1c;Sa18;ifSa1Tax=0;e1seifSa1=8TaX=(Sd15)*0.03;e1seifSa1=1.7Tax=3*0.03+(Sa1-8)*0.1;e1seifSa1=30Tax=3*0.03

3、+9*0.1+(Sa1-1.7)*0.2;e1seifSa140Tax=3*0.03+9*0.1+13*0.2+(Sa1-30)*0.25;e1seifSa160Tax=3*0.03+9*0.1+13*0.2+10*0.25+(Sa140)*0.3;e1seifSa185Tax=3*0.03+9*0.1+13*0.2+10*0.25+20*0.3+(Sa1-60)*0.35;e1seTax3*0.03+9*0.1+13*0.2+10*0.25+20*0.3+25*0.35+(Sa1-85)*0.45;enddisp(U应缴个人所得税为：Irmm2str(1000*Tax),，元)；结果:应缴个

4、人所得税为：90元按年扣缴程序:c1eara11;c1osea11;c1c;Sa120;ifSa1Tax=0;e1seifSa19.6Tax=(Sa1-6)*0.03;e1seifSa120.4Tax3.6*0.03+(Sa1-9.6)*0.1;e1seifSa136Tax=3.6*0.03+8.4*0.1+(Sa1-20.4)*0.2;e1seifSa148Tax=3.6*0.03+8.4*0.1+15.6*0.2+(Sa136)*0.25;e1seifSa172Tax=3.6*0.03+8.4*0.1+15.6*0.2+12*0.25(Sa1-48)*0.3;e1seifSa1102Tax

5、=3.6*0.03+8.4*0.1+15.6*0.2+12*0.2524*0.3+(Sa172)*0.35;e1seTax=3.6*0.03+8.4*0.1+15.6*0.2+12*0.25+24*0.3+30*0.35(Sa1-102)*0.45;enddisp(1应缴个人所得税为：Inum2str(Tax)z,万元，);结果：应缴个人所得税为：1.148万元3、利用matIab编写一个程序实现三对角矩阵的存储和数据获取，其中对角线上为：1,2,n,第一上对角元素分别为1,1/2,1(nT),第一下对角元素分别为h1/2,，1(n-1),求当n=5,10,30,100时的行列式的值。程序：c

6、1eara11;c1osea11;c1c;n=5f10,30,100;fori=1:4forj=1:n(i)-1A(j,j)=j；A(jrj+1)=1/j;A(j+1,j)=1/j;endA(n(i),n(i)=n(i);det(A)end结果：ans=54.2726ans=1.6367e+06ans=1.1960e+32ans=4.2079e+1574、采用高斯主元消去法解线性方程组，要求自己构造系数矩阵和右端向量，并且未知数个数为3、10、50、100进行试算。程序：c1eara11;c1osea11;c1c;Num=3,10,50,100;fori=1:4A=rand(Num(i)zNum(i);B=rand(Num(i);N,N=size(A);X=zeros(N,1);Y=zeros(N,1);C=zeros(1,N);R=I:N;k=1;whi1ek、一2.25,1).程序：c1eara11;c1osea11;c1c;A=2,1,1；3,1,2；1r2f2;