第04章MATLAB程序设计例题源程序.docx

高教社刘卫国《MATLAB程序设计与应用》

高教社刘卫国《MATLAB程序设计与应用》（第二版）例题源程序

第4章：MATLAB程序设计

PAGE10王成优?

PAGE10

王成优?山东大学（威海）MATLAB语言

WangChengyou?ShandongUniversity,Weihai

第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;a

b

然后在MATLAB的命令窗口中输入exch，将会执行该命令文件。程序2：

首先建立函数文件fexch.m：function[a,b]=exch(a,b)c=a;a=b;b=c;

然后在MATLAB的命令窗口调用该函数文件：

clear;x=1:10;

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 计算分段函数：

??cos(x?1)? x2?1, x?10

x x? x

x x? x,

??

x=input(请输入x的值:);

x?10

ifx==10y=cos(x+1)+sqrt(x*x+1);

else

y=x*sqrt(x+sqrt(x));end

y

也可以用单分支if语句来实现：x=input(请输入x的值:);y=cos(x+1)+sqrt(x*x+1);

ifx~=10

y=x*sqrt(x+sqrt(x));

endy

或用以下程序：

x=input(请输入x的值:);ifx==10

y=cos(x+1)+sqrt(x*x+1);

end

ifx~=10

y=x*sqrt(x+sqrt(x));

endy

例4.4 输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则

输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出。

c=input(请输入一个字符,s);ifc=Ac=Z

disp(setstr(abs(c)+abs(a)-abs(A)));elseifc=ac=z

disp(setstr(abs(c)-abs(a)+abs(A)));elseifc=0c=9

disp(abs(c)-abs(0));

else

disp(c);

end

例4.5 某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)：

price200

没有折扣

200≤price500

3%折扣

500≤price1000

5%折扣

1000≤price2500

8%折扣

2500≤price5000

10%折扣

5000≤price

14%折扣

输入所售商品的价格，求其实际销售价格。price=input(请输入商品价格);switchfix(price/100)

case{0,1} %价格小于200rate=0;

case{2,3,4} %价格大于等于200但小于500

rate=3/100;

casenum2cell(5:9) %价格大于等于500但小于1000rate=5/100;

casenum2cell(10:24) %价格大于等于1000但小于2500

rate=8/100;

casenum2cell(25:49) %价格大于等于2500但小于5000rate=10/100;

otherwise %价格大于等于5000

rate=14/100;

end

price=price*(1-rate) %输出商品实际销售价格

例4.6 矩阵乘法运算要求两矩阵的维数相容，否则会出错。先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。

A=[1,2,3;4,5,6];

B=[7,8,9;10,11,12];

try

C=A*B;

catch

C=A.*B;

endC

lasterr %显示出错原因

例4.7 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。

form=100:9