MATLAB编程.doc

MATLAB程序设计一、M文件 包含MATLAB语言代码的文件称为M文件，其后缀为M，M文件分为函数程序和脚本程序两种。函数程序用于扩展MATLAB语言的应用、脚本程序用于重复执行一系列指令。编辑M文件可以使用各种文本编辑器，也可以使用MATLAB自带的M文件编辑器。1、 脚本M文件的编写 FILE菜单NEW选项M-FILE选项M文件编辑器（程序编辑器） 例1：x=-pi:0.1:pi; y=sin(x); plot(x,y)按F5运行例2：x=-2:0.1:2; （直方图） y=exp(-x.*x); bar(x,y)保存为MF2。M， mf2 可运行。例3：x=4,3,7,2,1,6,5; （饼图） ex=0,0,0,0,0,1,0; pie(x,ex);2、 函数M文件的编写 (注意保存时函数名与文件名相同)例1：function c=myfile(a,b);c=sqrt(a.2+b.2); x=3; y=4; myfile(x,y)ans = 5例2：function f=fact(n);% fact(n) retuerns the factorial of n,usually denoteed by N!f=prod(1:n); m=5; fact(m)ans = 120 m=10; fact(m)ans = 3628800 help fact fact(n) retuerns the factorial of n,usually denoteed by N!例3：function s=ts(a,b,c); p=0.5*(a+b+c); s=sqrt(p*(p-a)*(p-b)*(p-c);保存时以函数名ts .m , area=ts(3,4,5)调用二、流程控制1、 FOR循环格式：for x=array command end例1： for n=1:10 x(n)=sin(n*pi/10); end xx = Columns 1 through 7 0.3090 0.5878 0.8090 0.9511 1.0000 0.9511 0.8090 Columns 8 through 10 0.5878 0.3090 0.0000例2： data=3,9,45,6;7 16 -1 5data = 3 9 45 6 7 16 -1 5 for n=data x=n(1)-n(2) endx = -4x = -7x = 46x = 1例3： for n=1:5 for m=5:-1:1 a(n,m)=n2+m2; end disp(n) end 1 2 3 4 5 aa = 2 5 10 17 26 5 8 13 20 29 10 13 18 25 34 17 20 25 32 41 26 29 34 41 502、 WHILE 循环 格式：while expression command end例1：a0=0;a1=1;disp(a0);disp(a1);n=2;while n petals help petals 用于创建花瓣图案的脚本程序petals.m查看工作空间变量情况。曲线拟合 t=0 0.3 0.8 1.1 1.6 2.3; y=0.5 0.82 1.14 1.25 1.35 1.40; plot(t,y,o),grid on根据图形，可用一个多项式函数y=a0+a1x+a2x2 对其进行线性拟合。1 t1 t1 2 y1 1 t2 t2 2 a0 y2 1 t3 t3 2 a1 = y3 1 t4 t4 2 a2 y4 1 t5 t5 2 y5 1 t6 t6 2 y6解这个过限定性系统的线性方程组可得未知系数a0、a1、a2。 x=ones(size(t) t t.2x = 1.0000 0 0 1.0000 0.3000 0.0900 1.0000 0.8000 0.6400 1.0000 1.1000 1.2100 1.0000 1.6000 2.5600 1.0000 2.3000 5.2900 a=xya = 0.5318 0.9191 -0.2387因此这个二次多项式函数为