[2018年最新整理]chap07专题1MATLAB编程实例.pptVIP

下载本文档

1
0
约5.39千字
约 28页
2018-02-20 发布于浙江
举报
版权申诉

[2018年最新整理]chap07专题1MATLAB编程实例.ppt

1、本文档共28页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

[2018年最新整理]chap07专题1MATLAB编程实例

MATLAB 编程实例(I) * MATLAB 编程实例 MATLAB 编程实例(I) * 目录求素数计算税率猴子选大王有序数组的插入打印杨辉三角打印矩阵的最小值打印鞍点求矩阵的所有K阶子式 MATLAB 编程实例(I) * 1. 求素数输入一个数，判断它是否是素数，如果不是素数，则至少输出它的一个非1因子 MATLAB 编程实例(I) * 程序 % 求一个矩阵中的每个元素是否是奇数，相当于MATLAB的函数 isprime function [res fctr]= is_prime(n) res = ones(size(n)); fctr = ones(size(n)); for K=1:length(n(:)) cur_n = n(K); [res(K) fctr(K)] = is_a_prime(cur_n); end % 判断单个数是否是奇数 function [res fctr] = is_a_prime(n) res = true; for I=2:floor(sqrt(n)) if mod(n, I) == 0 fctr = I; % 找到一个因子 res = false; return; end end % 主程序 subplot(3, 1, 1); A=[2:100]; stem(A, A.*is_prime(A), .-) subplot(3, 1, 2); A=[2:1000]; stem(A, A.*is_prime(A), .-) subplot(3, 1, 3); A=[2:5000]; stem(A, A.*is_prime(A), .-) MATLAB 编程实例(I) * 结果 100x100个素数分布 MATLAB 编程实例(I) * MATLAB 编程实例(I) * 2. 计算税率月收入税率 1600以下 0 1600-2100 5% 2100-3600 10% 3600-6600 15% 6600-21600 20% 21600-41600 25% 41600-61600 30% 61600-81600 35% 81600-101600 40% 101600以上 45% MATLAB 编程实例(I) * salary=linspace(0, 10000, 20); plot(salary, get_revenue(salary), x-); grid on; subplot(2, 1, 2); salary=linspace(10000, 200000, 20); plot(salary, get_revenue(salary), x-); grid on; function revenues = get_revenue(salaries) revenues = ones(size(salaries)); for I=1:length(salaries(:)) revenues(I)=get_a_revenue(salaries(I)); end 程序1 function revenue = get_a_revenue(salary) revenue = 0; if salary 1600 return; end % -- 1600 ~ 2100 salary = salary - 1600; if salary (2100-1600) revenue = revenue + salary * 0.05; return; else revenue = revenue + (2100-1600)*0.05; end % -- 2100 ~ 3600 salary = salary - (2100-1600); if salary (3600-2100) revenue = revenue + salary * 0.10; return; else revenue = revenue + (3600-2100)*0.10; end % -- 3600 ~ 6600 % -- 6600 ~ 21600 % -- 21600 ~ 41600 % -- 41600 ~ 61600 % -- 61600 ~ 81600 % -- 81600 ~ 101600 % -- 81600 ~ 101600 salary = salary - (81600-61600); if salary (101600-81600) revenue = revenue + salary * 0.40; return; else