第04章-matlab程序设计-例题源程序

《第04章-matlab程序设计-例题源程序》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、高教社刘卫国《MATLAB程序设计与应用》（第二版）例题源程序第4章：MATLAB程序设计第4章MATLAB程序设计例4.1建立一个命令文件将变量a,b的值互换，然后运行该命令文件。程序1：首先建立命令文件并以文件名exch.m存盘：clear;a=1:10;b=[11,12,13,14;15,16,17,18];c=a;a=b;b=c;ab然后在MATLAB的命令窗口中输入exch，将会执行该命令文件。程序2：首先建立函数文件fexch.m：function[a,b]=exch(a,b)c=a;a=b;b=c;然后在MATLAB的命令窗口调用该函数文件：clear;x=1:10

2、;y=[11,12,13,14;15,16,17,18];[x,y]=fexch(x,y)例4.2求一元二次方程ax2+bx+c=0的根。a=input('a=?');b=input('b=?');c=input('c=?');d=b*b-4*a*c;x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);例4.3计算分段函数：x=input('请输入x的值:');10王成优©山东大学（威海）MATLAB语言WangChengyou©ShandongUni

3、versity,Weihai高教社刘卫国《MATLAB程序设计与应用》（第二版）例题源程序第4章：MATLAB程序设计ifx==10y=cos(x+1)+sqrt(x*x+1);elsey=x*sqrt(x+sqrt(x));endy也可以用单分支if语句来实现：x=input('请输入x的值:');y=cos(x+1)+sqrt(x*x+1);ifx~=10y=x*sqrt(x+sqrt(x));endy或用以下程序：x=input('请输入x的值:');ifx==10y=cos(x+1)+sqrt(x*x+1);endifx~=10y=x*sqrt(x+sqrt(x));en

4、dy例4.4输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出。c=input('请输入一个字符','s');ifc>='A'&c<='Z'disp(setstr(abs(c)+abs('a')-abs('A')));elseifc>='a'&c<='z'disp(setstr(abs(c)-abs('a')+abs('A')));elseifc>='0'&c<='9'disp(abs(c)-abs('0'));else10王成优©山东大学（威海）MATLAB语言WangChengyou

5、©ShandongUniversity,Weihai高教社刘卫国《MATLAB程序设计与应用》（第二版）例题源程序第4章：MATLAB程序设计disp(c);end例4.5某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)：price<200没有折扣200≤price<5003%折扣500≤price<10005%折扣1000≤price<25008%折扣2500≤price<500010%折扣5000≤price14%折扣输入所售商品的价格，求其实际销售价格。price=input('请输入商品价格');switchfix(price/100)case

6、{0,1}%价格小于200rate=0;case{2,3,4}%价格大于等于200但小于500rate=3/100;casenum2cell(5:9)%价格大于等于500但小于1000rate=5/100;casenum2cell(10:24)%价格大于等于1000但小于2500rate=8/100;casenum2cell(25:49)%价格大于等于2500但小于5000rate=10/100;otherwise%价格大于等于5000rate=14/100;endprice=price*(1-rate)%输出商品实际销售价格例4.6矩阵乘法运算要求两矩阵的维数相容，否则会出错。

7、先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。A=[1,2,3;4,5,6];B=[7,8,9;10,11,12];tryC=A*B;10王成优©山东大学（威海）MATLAB语言WangChengyou©ShandongUniversity,Weihai高教社刘卫国《MATLAB程序设计与应用》（第二版）例题源程序第4章：MATLAB程序设计catchC=A.*B;endClasterr%显示出错原因例4.7一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全