今天我们开始学编程啦~
这是新建新的.m文件，也就是matlab命令文件，
一般命令文件不带输入和输出的参数，只是一些命令的组合，带有输入输出参数的文件是函数文件~
下面我们来编写一个hello world 程序~
输入disp('Hello World!');
然后点击运行
是不是很无聊。。。那就来个好一点的~
输入：
t=pi*(0:0.001:1);
x=sin(t).*sin(10*t);
comet(t,x)
然后运行，，自己研究玩去吧~
m文件就是这样啦，方便多次运行保存和修改，解决复杂的问题更方便~♪(^∇^*)
我们继续学习。。。
接下来看看函数文件的编写
新建一个函数文件
函数的第一行的格式 ：function [输出的参数] =函数名 (输入的参数)
输入和输出的参数可以有多个
保存函数文件的时候，注意文件名要和函数名一样
函数头和函数体之间可以有多行%开始的注释，说明函数的使用方法，当执行 help 函数名 的时候会显示这些注释 ，直到遇到第一个非注释行为止
语法也没什么好讲的，我们直接看一个例子吧 ~(～￣▽￣)～
这是二分法求近似值的程序：
function [x_star,k] = bisect1(fun,a,b,ep)
%二分法求方程的近似值
if nargin<4
ep=1e-5;
end
fa=feval(fun,a);fb=feval(fun,b);
if(fa*fb>0)
disp('error');
x_star=NaN;
k=NaN;
else
k=1;
while abs((b-a)/2)>ep
x=(a+b)/2;
fx=feval(fun,x);
if fa*fx<0
b=x;fb=fx;
else
a=x;fa=fx;
end
k=k+1;
end
x_star=(a+b)/2;
end
(⊙﹏⊙)b看不懂也没关系，看解释~~~
第一行：
function [x_star,k] = bisect1(fun,a,b,ep)
注意输入的参数。x_star, k 都是数字；fun 是函数；a, b, ep都是数字
这个程序目的是输入一个函数fun以及给定的区间[a,b]，求出精确度为ep的零点x_star，k为迭代次数 。
if nargin<4
ep=1e-5;
end
nargin,nargout 是函数文件里的固有变量，意思是输入（输出）参数的个数
这句话就是说，如果只输入前三个参数，那么精确度默认为1e-5
fa=feval(fun,a); fb=feval(fun,b);
feval(fun,a)用于计算函数fun 在点a 处的函数值，这个命令记住就行了。(๑•̀ㅂ•́)و✧
if(fa*fb>0)
disp('error');
x_star=NaN;
k=NaN;
上面这段是无解的情况
else
k=1;
迭代次数初始化
while abs((b-a)/2)>ep
x=(a+b)/2;
fx=feval(fun,x);
区间大于精确度的话就迭代
每循环一次迭代次数加一~~然后没什么可说的了~~~
最后把结果赋值给输出的参数
x_star=(a+b)/2;
写完保存，记得文件名要和函数名相同(⊙v⊙)
然后回到命令窗口
我们来求一个函数f(x)=x^3-x-1 在区间[1,15]上的零点~~
在命令窗口依次输入：
fun=inline('x^3-x-1');
[x_star,k]=bisect1(fun,1,15,0.005)
插一句：关于循环语句 ，和c语言一样的有三种
while
for
do while
注意for(int i=1;i<10;i++)
在matlab写成for i=1:10 ... end
如果是i=i+2,就写成 for i=1:2:10
这个和那个线性生成向量的是一样的啦(～￣▽￣)～
matlab函数也是可以递归调用的，但是效率很差，不推荐使用喔。
哦！我突然想起来c语言里的数组，这在MATLAB里就是用向量，向量用起来还是很方便的，使用前不需要声明，不要动态内存分配，随便用个a(i)就行了
举个例子：求斐波那契数列函数
function [ F ] = Fibonacci( n )
F()=;
F()=;
for i = : n
F(i) = F(i-)+F(i-);
end end
二维同理~
自己多写一些练一练喔 ✧(≖ ◡ ≖✿)