实验五M文件和MATLAB程序设计.docxVIP

实验五 M文件和MATLAB程序设计 一、 实验目的 matlab作为一种高级计算机语言，不仅可以命令行方式完成操作，也 具有数据结构、控制流、输入输出等能力，本次实验通过熟悉和掌握m文 件的建立与使用方法，以及函数与控制程序流程语句的使用，使学牛具备 一定的编程和程序调试能力。 掌握M文件的使用方法。 掌握讦语句和switch语句的使用 掌握循环语句的使用 通过练习理解MATLAB编程方法。 二、 实验原理 m文件 用matlab语言编写的程序，称为m文件。M文件根据调用方式的不 同分为两类，命令文件(Script file)和函数文件(Function file)o区别？ 程序控制结构 顺序结构 选择结构 讦语句 a)单分支if语句b)双分支讦语句 c)多分支讦 语句 switch 语句 try语句 循环结构 for语句 while 语句 break 语句、continue 语句、return 使用,区别？ 函数文件 function输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 三、 实验要求 首先上机练习PPT屮各种流程控制语句的有关实例。 然后上机练习下面的实验习题。 四、 实验习题 1.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除 之，若为奇数，则与3相乘再加上重复此过程，最终得到的结果为1。 如： 2?1 3?10?5?16?8?4?2?1 6?3?10?5?16?8?4?2?1 运行下面的程序，按程序提示输入21,235,7等数來验证这一结论。 %classic ”3n+1 problem from number theory? while 1 n 二 in put(#39;E nter n’n egative quits:#39;); if nlt;=O break end a 二 n; while ngt;l 讦 rem(n?2)==0 n=n/2; else n=3*n+l; end a=[a,n]; end a end Enter n’negative quits:3 3 10 5 16 8 4 2 1 i2?100002? 编程求满足?i?lm的最小m值。 a=0; i=l; while (alt;10000) a=a+pow2(i); end m=i-l; m 13 3.编写一个函数，计算下面函数的值，给岀x的值, i=i+l; 调用该函数后, 返回y的值。 function [y]二myfunl(x) ?sinx,?y(x)??x, ??x?6,? x?00?x?3x?3y=sin(x); x?00?x?3x?3 y=sin(x); else讦 xgt;3 y=-x+6; end -1 gt;gt; y=myfunl(O)y = function y 二 myfun l(x) if xlt;=0 else 讦 xgt;0amp;xlt;二 3 y=x; 运行结果:gt;gt; y=myfunl(-pi/2) y = 0 gt;gt; y=myfunl(2) y = 2 gt;gt; y=myfun 1 ⑷ y 二 2 用如下迭代公式求a, a的值分别为：3, 17, 113o迭代的终止 条件为 100次。分别对迭代结果和xn?l?xn?10?5,迭代初值x0?1.0,迭代次数 不超过 xn?la?准确值进行比较，并统计迭代次数。 迭代公式：xn?l?22xn?l function [x, n]=sqrt_a(a) x=1.0; for k=l:100 m=x; x=x/2+a/(2*x); if abs(x-m)lt;=10A(-5) break end end x; n=k; s=(x-sqrt(a)); if slt;=10A(-5) disp(#39;正确#39;); else disp(#39;错误#39;); end下面调用举例： 运行结果： gt；gt; [x ,n]=sqrt_a(3)正确 x = 1.7321 n 二 5 gt;gt; [x, n]=sqrt_a(l刀 正确 x 4.1231 n 二 6 gt;gt; [x,n]=sqrt_a(113)正确 10.6301 8 定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并 在命令文件中调用该函数文件。 程序设计： function [e In s c]=num(x) e=exp(x) ln=log(x) s=si n(x) c=cos(x) end 运行结果： gt；gt; num(5i) 0.2837-0.9589i In = 1.6094 + 1.5708i 0 +74.2032i 74.2099 ans = 0.2837-0.9589i 设 f(x)=l (x?2)2??0.11,编写一个