MATLAB 结构化程式与自定义函数

主要内容：

1.script writing 脚本编写

2.structured programming 结构化程序设计

3.structured programming 用户定义函数

一.部分基础操作（脚本编辑器script editor）

新建脚本

for i=1:10
    x=linspace(0,10,101);
    plot(x,sin(x+i));
    print(gcf,'-deps',strcat('plot',num2str(i),'.ps'));
end

1.程序内一部分不运行：

当想让程序内某一行不运行，可以在该行前加符号%，即可将改行变为注解，该行不运行

若想让程序内连续几行不运行，可以将该部分选中，右击选择【注释】，或点击工具栏内【注释】按键

2.程序继续运行:

若想让选中程序继续运行，可选中该部分程序，右击或者在工具栏中单击【取消注释】即可

3.运行部分代码：可用来检测程序的bug，分部分进来检查

在该部分代码前输入%%，即可将代码进行分节，单击选中需要运行部分，点击工具栏内【运行节】按键，即可运行该部分代码

单击【运行或F5】按键，运行整个程序

4.查看程序到某一节点前是否存在错误：

点击第四行，4右侧横线，变为圆点，（中断点）

代表代码在运行到第4行时，就会停止，此时可以查看已执行程序是否符合预期或存在错误

退出：关掉中断点，点击上侧工具栏【继续】按键，即可退出

5.当输入代码未缩进，可全选【ctrl+A】，右击，选择【智能缩进】

二.层次码执行

按照行的顺序，从上到下

更改条件之后，再次输出结果时，需要先清除之前输出的变量

structured programming（结构化程序设计）

if,elseif,else   

if  condition1
    statement1
elseif condition2
      statement2
else
      statement3
end

a=3;
if rem(a,2)==0
    disp('a is even')
else
    disp('a is odd')
end

rem，表示remainder余数

even，表示even number偶数

odd，表示odd number奇数

disp，表示display显示

for

for用来执行一个事情很多次（遍历循环）

for variable=start:increment（增量）:end

     commands

end

for n=1:10
    a(n)=2^n;
end
disp(a)

依次输出2^n,直至2^10

i=1
for n=1:2:10
    a(i)=2^n;
    i=i+1
end
disp(a)

依次输出2的奇数次方，直至n到10

switch,case

switch expression(表达式)

case value 1

          statement1

case value 2

        statement2

...

otherwise

       statement

end

input_num=-1;
switch input_num
    case -1
        disp('negative 1');
    case 0
        disp('zero');
    case 1
        disp('positive 1');
    otherwise 
        disp('other value');
end

otherwise

try,catch

while

（条件循环）

while expression

         statement

end

n=1;
while prod(1:n)<1e100
    n=n+1;
end

prod相乘

prod(1:n)=prod([1 2 3...n])=1*2*3...*n=n!

e100=10^100(10的100次方)

break（与while一起用）

continue

end

pause

return

relational(logical)   operators   关系（逻辑）操作符

<   

<=

>

>=

== equal to

~= not equal to 不等于

&& and

||  or

三.pre-allocating space to variables 为变量预分配空间

为变量预分配空间可以大大减少代码的运算时间

例：未分配

tic
for ii=1:2000
    for jj=1:2000
        A(ii,jj)=ii+jj;
    end
end
toc

历时 1.565193 秒

例：分配

tic
A=zeros(2000,2000);
for ii=1:2000
    for jj=1:2000
        A(ii,jj)=ii+jj;
    end
end
toc

历时 0.021914 秒

tic指令，用来计时

四.一些小Tips

1.【clear all】消除之前的所有变量

2.【close all】关掉所有图形

3.【clc】清空【命令行窗口】

4.使用分号【；】，不会将所有的运行结果都显示在【命令行窗口】，可使页面保持干净

5.【...】换行号      接续上一行

6.【Ctrl+C】在程序执行中，即界面右下角显示busy，可协助跳出程序

五.MATLAB内置函数的内容

1.【>>edit(which('mean.m'))】

funtion---keyword关键字

y---output输出

mean---function name=file name函数名

x---input输入

绿色输出---online help

再往下---MATLAB Code

2.定义函数，进行调用

function x=freebody(x0,v0,t)
x=x0+v0+1/2*9.8*t.*t;

注意点乘

即，定义了freebody函数

将此保存至以freebody命名的对应的文件下

在【命令行窗口】即可以调用此函数

>>freebody(0,0,10)    ans=490

即初始位置与初始时间均为0，经过10s之后，自由落体运动经过的路程

3.定义acc函数

function[a,F]=acc(v2,v1,t2,t1,m)
a=(v2-v1)./(t2-t1);
F=m.*a;

acc=acceleration加速度

4.

nargin---number of function input arguments函数输入参数的个数

nargout---函数输出参数的个数

varargin---variable length input argement list    可变长度输入参数列表

varargout---可变长度输出参数列表

inputname---variable name of function input   函数输入的变量名

mfilename---file name of currently running function   当前运行函数的文件名

5.function handles函数句柄

a way to create anonymous functions

一种创建匿名函数的方式

f=@(x) exp(-2*x);
x=0:0.1:2;
plot(x,f(x));