资源描述：
《chap07 专题1 MATLAB 编程实例.ppt》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、MATLAB编程实例(I)1MATLAB编程实例MATLAB编程实例(I)2目录求素数计算税率猴子选大王有序数组的插入打印杨辉三角打印矩阵的最小值打印鞍点求矩阵的所有K阶子式MATLAB编程实例(I)31.求素数输入一个数，判断它是否是素数，如果不是素数，则至少输出它的一个非1因子MATLAB编程实例(I)4程序%求一个矩阵中的每个元素是否是奇数，相当于MATLAB的函数isprimefunction[resfctr]=is_prime(n)res=ones(size(n));fctr=ones(size(n));forK=1
2、:length(n(:))cur_n=n(K);[res(K)fctr(K)]=is_a_prime(cur_n);end%判断单个数是否是奇数function[resfctr]=is_a_prime(n)res=true;forI=2:floor(sqrt(n))ifmod(n,I)==0fctr=I;%找到一个因子res=false;return;endend%主程序subplot(3,1,1);A=[2:100];stem(A,A.*is_prime(A),'.-')subplot(3,1,2);A=[2:1000];s
3、tem(A,A.*is_prime(A),'.-')subplot(3,1,3);A=[2:5000];stem(A,A.*is_prime(A),'.-')MATLAB编程实例(I)5结果100x100个素数分布MATLAB编程实例(I)6MATLAB编程实例(I)72.计算税率月收入税率1600以下01600-21005%2100-360010%3600-660015%6600-2160020%21600-4160025%41600-6160030%61600-8160035%81600-10160040%101600以上
4、45%MATLAB编程实例(I)8salary=linspace(0,10000,20);plot(salary,get_revenue(salary),'x-');gridon;subplot(2,1,2);salary=linspace(10000,200000,20);plot(salary,get_revenue(salary),'x-');gridon;functionrevenues=get_revenue(salaries)revenues=ones(size(salaries));forI=1:length(s
5、alaries(:))revenues(I)=get_a_revenue(salaries(I));end程序1functionrevenue=get_a_revenue(salary)revenue=0;ifsalary<1600return;end%--1600~2100salary=salary-1600;ifsalary<(2100-1600)revenue=revenue+salary*0.05;return;elserevenue=revenue+(2100-1600)*0.05;end%--2100~3600sa
6、lary=salary-(2100-1600);ifsalary<(3600-2100)revenue=revenue+salary*0.10;return;elserevenue=revenue+(3600-2100)*0.10;end%--3600~6600%--6600~21600%--21600~41600%--41600~61600%--61600~81600%--81600~101600%--81600~101600salary=salary-(81600-61600);ifsalary<(101600-81600
7、)revenue=revenue+salary*0.40;return;elserevenue=revenue+(101600-81600)*0.40;end%--101600~infsalary=salary-(101600-81600);revenue=revenue+salary*0.45;MATLAB编程实例(I)9结果MATLAB编程实例(I)10程序2:查表法functionrevenues=get_a_revenue2(salary)revenues=0;starts=[-inf,1600,2100,3600,6
8、600,21600,41600,61600,81600,101600,inf];pre_fee=[0,0,25,175,625,3625,8625,14625,21625,29625];rates=[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0
