MATLAB实验指导书（附答案）(4)

f =

1/y*exp(-x*y)

9、对符号方程f = ax2+bx+c 求解。1）对x求解，2）对a求解。 （1）syms a b c x;f=2*x*a+x*b+c;x=solve(f) x =

-c/(2*a+b)

（2）syms a b c x;f=2*x*a+x*b+c;a=solve(f,a) a =

-1/2*(x*b+c)/x

d2ydydy10、求解微分方程2?2?2y?0,y(0)?1,(0)?0。

dxdxdxsyms x y;y=dsolve('D2y+2*Dy+2*y=0','y(0)=1,Dy(0)=0','x') y =

exp(-x)*sin(x)+exp(-x)*cos(x)

15

实验三 MATLAB语言的程序设计

一、实验目的：

1、熟悉MATLAB程序编辑与设计环境

2、掌握各种编程语句语法规则及程序设计方法 3、函数文件的编写和设计

4、了解和熟悉跨空间变量传递和赋值

二、实验基本知识：

1、程序流程控制语句 for循环结构

语法：for i＝初值：增量：终值

语句1 …… 语句n end

说明：1．i＝初值：终值，则增量为1。

2．初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。

while 循环结构 语法：while 逻辑表达式

循环体语句 end

说明：1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。若表达式的

值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。若表达式的值为假，则程序执行end之后的语句。

2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位置加break语句、以便程序能正常执行。（执行循环体的次数不确定；每一次执行循环体后，一定会改变while后面所跟关系式的值。）

16

3、while循环也可以嵌套、其结构如下： while逻辑表达式1

循环体语句1 while逻辑表达式2

循环体语句2 end

循环体语句3 end

if-else-end分支结构 if 表达式 1

语句1

else if 表达式 2(可选)

语句2

else(可选)

语句3

end end

说明：1.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不

满足，则跳出if结构。

2．else if表达式2与else为可选项，这两条语句可依据具体情况取舍。 3．注意：每一个if都对应一个end，即有几个if，记就应有几个end。 switch-case结构 语法：switch表达式

case常量表达式1 语句组1 case常量表达式2 语句组2 …… otherwise

17

语句组n end

说明：1．switch后面的表达式可以是任何类型，如数字、字符串等。

2．当表达式的值与case后面常量表达式的值相等时，就执行这个case后面的语句组如果所有的常量表达式的值都与这个表达式的值不相等时，则执行otherwise后的执行语句。 程序流程控制指令 break、return、pause

2、函数文件的结构和编写方法

函数只能在函数体内对变量进行操作，也就是只能访问函数本身工作空间中的变量。

M函数文件的结构

（1） 函数定义行（function） （2） H1行（函数帮助文本的第一行） （3） 函数帮助文本 （4） 函数体 （5） 注释

函数文件编写后，保存时，其文件名必须与函数名相同。注意：函数名不要与MATLAB自身的函数命令相同。

三、实验内容：

练习A

1、熟悉MATLAB程序编辑与设计环境

18

2、用for循环语句实现求1～100的和。

3、用for循环语句实现编写一个求n阶乘的函数文件。

a4、Fibonacci数组的元素满足Fibonacci 规则： k?2?ak?ak?1(k?1,2,?)且 a 1 ? a 2 ? 1 ，现要求该数组中第一个大于10000的元素。 练习B 思考题

4、已知一维数组A=[2,4,5,8,10],B=[4,6,9,3,4]，用for循环语句实现?AiBn?i?1。

i?1n求和函数可用sum()。.

5、利用for循环求1！+2！+3！+ … +5！的值。.

19