首页 > 代码库 > MATLAB入门笔记

MATLAB入门笔记

% MATLAB:MATLAB 为 Matrix Laboratory ,用来处理矩阵可编程可实现算法逻辑的计算工具
%
% eg:1 绘制正弦和余弦曲线
x=[0:0.1:4*pi]; %建立角度向量
plot(x,sin(x),x,cos(x)) %画图
axis([0,2*pi,-1,1]) %设定坐标轴范围
title(′正弦和余弦曲线′); %图形名称
xlabel(′θ(弧度)′) %x轴名
ylabel(′y′) %y轴名
legend(′sinθ′,′cosθ′); %图例
grid on; %网格线


% MATLAB指令系统
% Who或whos 显示在当前工作区中的所有变量名,前者显示变量名,后者还显示变量的大小、字节数和类型
% disp(x) 显示x的内容,它可以是矩阵或字符串
% which test 显示 M 文件 test.m 所在的目录
% type test 在命令窗口下显示 test.m 的内容
% cd,chdir,pwd 显示目前的工作目录
% Load 文件名 调出mat文件中的数据。也可以调出文本文件,但是文本文件只能是由数字组成的矩阵形式
% diary
% 建立一个文本文件,记录在MATLAB中输入的所有命令和它们的输出,但是不能包括图形。如果想把你的输入存入一个特定的文件中,可使用 diary
% filename 建立文件。使用diary off↙命令可以停止记录
%
% what 返回目前目录下M,MAT,MEX文件的列表
% echo 控制是否显示M文件执行的每一条命令
% clc 擦除MATLAB工作区中所有显示的内容
% clf 擦除MATLAB图形 工作区中的图形
% hold 控制当前图形窗口对象是否被刷新(在图表中保持现有图形)
% dir,ls 列出指定目录下的文件和子目录清单
% path 显示目前的搜索路径,可以用File菜单中的 Set path 观察和修改路径
% quit 退出工作区可以用,也可选择File菜单中Exit命令
%
% 键盘操作部分
% ↑/Ctrl+p 重调前一行 (用于调出前面的命令进行修改,重新计算)
% ↓/Ctrl+n 重调下一行
% →/Ctrl+b 向前移一个字符
% ←/Ctrl+f 向后移一个字符
% Ctrl+→ 右移一个字
% Ctrl+← 左移一个字
% Home 移动到行首
% End 移动到行尾
% Esc 清除一行
% Del 删除光标后的字符
% Backspace 删除光标前的一个字符
% Ctrl+k 删除到行尾
% help 查看当前帮助系统
% look for 关键字查找,若是要全文搜索,需要加上 -all
% 模糊查询: 输入关键字,然后按tab键,系统会列出这几个关键字开头的命令
%
% 三种基本数据类型:
% 1.数值型数据:输入的数据为数值数据,包括实数和负数
% 2.字符串数据: 用英文格式单引号加以界定的数字,字符,各种符号,表达式,方程式和汉字等
% 3.符号型数据: 用sym或syms把字符、表达式、方程、矩阵等定义成数学符号,运算结果为数学表达式
%
% 变量符号的命名规则基本类似与JAVA,大小写是区分的,长度不能大于63个。
%
% 表达式使用示例:
% x = sin(pi/4) (若省略“=”,系统自动建立一个名为ans的变量)
% x = sin(pi/4); 若以分号结束,则只进行计算不会马上显示结果,但是不写分号,则会马上显示结果
% 一行中若是写几个语句,他们之间要用逗号或分号隔开
% 表达式的变量需要先定义后使用!
%
% 符号变量的声明: syms sym
% syms a1 a2 a3
% a1,a2,a3是需要定义为符号变量的标识符,不能是数字、函数表达式或方程式。
% 变量a1,a2,a3,不能用引号界定,而是用“空格”隔开
%
% x = sym(1/2),y = sym(‘1/2‘), z = sym(‘sin(x)+exp(x)‘)
% 显然,sym定义的字符参量可以是数字,字符串,也可以是字符串变量名,字符表达式或字符方程
%
% eg 2: 求极限
% lim [x^(1/m)-a^(1/a)]/x-a , x->a

syms x a;
y = sym(‘[x^(1/m)-a^(1/a)]/x-a‘)
limit(y,x,a)

% MATLAB的固定变量:
% ans 用于结果的缺省变量名
% i j 虚数单位,定义为根号-1
% pi 圆周率
% Inf 无穷大,如1/0
% NaN 不确定量 如0/0
% realmin 最小正实数 2.225073858507201×10-308
% realmax 最大正实数 1.797693134862316×10+308
%
% 对变量操作的命令
% who 显示当前工作空间中所有变量的一个简单列表
% whos 列出变量的大小、数据格式等详细信息
% clear 清除指定变量
% save 文件名 变量名 将变量列中所列出的变量保存到磁盘文件中,变量列表中,各个不同的变量之间不能用逗号,只能用空格来分隔
% load 文件名 变量名 将以前用save命令保存的变量从磁盘文件中调入MATLAB工作空间
%
%
% MATLAB主要操作之一:矩阵
% 常用操作矩阵的变量的方法:
% 1.在命令窗口中输入
% 2.语句生成
% 3.函数生成
% 4.小矩阵生成大矩阵
%
% 使用原则:
% 1.矩阵的元素排列在方括号内
% 2.按行输入,每行内的元素使用空格或逗号分隔
% 3.行与行之间用分号,或回车键隔开
% 4.矩阵生成不但可以使用纯数字(含复数),也可以使用变量和表达式
%
% 方括号表达式示例: A = [ 1,2,3;4,5,6;7,8,9]
% 语句生成示例: B = a:h:b 或 B = [a:h:b]
% 其中a为初始值,h为步长,b为终止值。当步长为1时可省略h参数,另外h也可以取负数
% C = linspace(a,b,n) 等价于 C = a:(b-a)/(n-1):b
% 意义是:在线性空间上,行矢量的值从a到b,数据个数为n,缺省则n默认为100.
%
% 常用创建矩阵的函数:
% ones(m,n) 产生m x n的全1矩阵
% zeros(m,n) 产生m x n的全0矩阵
% rand(m,n) 产生m x n的随机矩阵
% eye(m,n) 产生m x n的对角线全1矩阵
% length(X) 返回矩阵最长的维的维度
% size(X) 返回矩阵每一维的长度
%
% 4.小矩阵生成大矩阵示例:
% a = [ 1,2;3 4]
% c = [a a+2 3*a a-2]
% c =
% 1 2 3 4
% 3 4 5 6
% 3 6 -1 0
% 9 12 1 2
%
%
% 矩阵的引用:
% 一.引用矩阵元素
% 1.通过下表引用矩阵元素
% 2.矩阵元素的序号来引用矩阵元素。按列排序,m x n 矩阵A,矩阵元素A(i,j)的序号为(j-1)*m + i (m为列的长度)
% 二.矩阵的拆分
% 1.A(:,j)
% A(:,j)表示取A矩阵的第j列全部元素;
% A(i:i+m,:)表示取A矩阵第i~i+m行的全部元素;
% A(:,k:k+m)表示取A矩阵第k~k+m列的全部元素;
% A(i:i+m,k:k+m)表示取A矩阵第i~i+m行内,并在第k~k+m列中的所有元素
%
%
% MATLAB对数值矩阵的两种不同的运算方法:
% 矩阵算法
% 把矩阵看作一个整体,各种运算完全按照线性代数代表的矩阵运算法则进行,运算的书写形式和运算符号都与矩阵理论完全相同。
% 数组算法
% 把矩阵看作由其元素构成的一组数据(数组),各种运算是在参与运算矩阵的对应元素之间进行的数与数的运算,这种运算方便对大批数据的处理和一次求出多个函数值。
% 数组算法的运算符主要有 .* ./ .\ .^
%
% MATLAB运算符:
% + 加 A+B A B必须大小相同,或一个是标量
% - 减 A-B A B必须大小相同,或一个是标量
% * 矩阵乘 A*B A 的列数等于B的行数
% .* 数组乘 A.*B A B必须大小相同
% \ 矩阵左除 A\B =A-1* B 等效于A*X=B求X A-1=inv(A)
% .\ 数组左除 A.\B Bij/Aij
% / 矩阵右除 A/B =A*B-1 等效于X*B=A求X
% ./ 数组右除 A./B Aij/Bij
% ^ 矩阵乘方 A^m A自乘m次
% .^ 数组乘方 A.^m A中每个元素的m次方
%
% 总结:1.直接使用运算符是矩阵运算,前面带小数点为数组计算。
% 2.矩阵运算的结果是由原本的矩阵得到一个新的矩阵,而数组运算是对原本矩阵的矩阵中的元素做对应的运算
%
%
% 关系运算符:
% ==
% ~=
% >
% <
% >=
% <=
%
% 逻辑运算符:
% &
% |
% !
% 说明:逻辑运算认定任何非零的元素都为真“1”,任何零元素都为假“0”.
%
% 转置:
% A‘ —— 共轭转置(将矩阵中的每个元素取共轭复数,再转置)
% A.‘ —— 普通转置(直接对矩阵做转置)
%
% 逆矩阵与行列式计算:
% 求 逆: inv(A)
% 求行列式: det(A)
% 注意:这两种运算都要求矩阵A为方阵!
%
% 矩阵下标说明:MATLAB的矩阵中不同于其他语言,矩阵的下标从1开始,而不是从0开始!
%
% 重点1:四则运算与幂运算
% 1.只有维数相同的矩阵才能进行加减运算
% 2.只有当两个矩阵中前一个矩阵的列数和后一个矩阵的行数相同时,才可以进行乘法运算。
% 3.a\b运算等效于求a*X=b的解,而a/b等效于求x*b=a的解。(这一条可以用来求方程的解)
% 4.只有仿真才可以求幂
% 5.点运算是两个维数相同矩阵对应元素之间的运算。
%
% 特殊矩阵之对角矩阵:V = [a1,a2,a3···,an]
% A = diag(V)
% A = a1 0 0·
% 0 a2 0 0
% 0 0 a3·
% · · ·an
%
% 数值显示格式控制:Format命令
% format(short) : 短格式(5位定点书)99.1253
% format long : 长格式(15位定点书)99.123456789000000
% format short e : 短格式e方式 9.91e+001
% format long e : 长格式e方式 9.912345678900000e+001
% format bank : 2位十进制(小数部分)99.12
% format hex : 十六进制格式
%
% 重点一:求线性方程的解:
% x1 + 4*x2 - 7*x3 + 6*x4 = 0
% 2*x2 + x3 + x4 = -8
% x2 + x3 + 3*x4 = -2
% x1 + x3 - x4 = 1
%
% 实际代码:
A = [1 4 -7 6;0 2 1 1;0 1 1 3;1 0 1 -1]
B = [0;-1;-2;1]
x = A\B

% MATLAB绘图篇:
% help graph2d可以得到所用画二维图形的命令
% help graph3d可以得到所用画三维图形的命令
%
% 一:绘制二维图形
% 格式:
% plot(X,‘S‘)
% plot(X,Y,‘S‘)
% plot(X1,Y1,‘S1‘,X2,Y2,‘S2‘···,Xn,Yn,‘Sn‘)
% 说明:
% 参数‘S’控制数据点的 标记 曲线类型 和 曲线色彩 , 三者置于一对单引号内。
%
% 常用的绘图选项:
% - 实线 * 用星号标出数据点
% -- 虚线 . 用点号标出数据点
% : 点线 。 用圆圈号标出数据点
% -. 点划线 x 用叉号标出数据点
% b 蓝色 + 用加号标出数据点
% g 绿色 s 用正方形标出数据点
% r 红色 D 用菱形出数据点
% c 青色 V 用下三角标出数据点
% m 洋红 ^ 用上三角标出数据点
% y 黄色 < 用左三角标出数据点
% k 黑色 > 用右三角标出数据点
% w 白色 H 用六角形标出数据点
% P 用五角形标出数据点
%

% 格式一:plot(X,‘S‘)
% 说明:
% 当X是实向量时,以该向量元素的下标为横坐标,元素值为纵坐标画出一条连续曲线。
% 当X是实矩阵时,则按列绘制每列元素值相对其下标的曲线,曲线条数等于输入参数矩阵的列数。
% 示例:
x=0:0.1:2*pi;%length(x)=63
y1=sin(x);
y2=cos(x);
y=[y1‘,y2‘];
z=[y1;y2];
figure(1)
plot(y1 ,‘--r+‘)
figure(2)
plot(y)
figure(3)
plot(z)

% 格式二:plot(X,Y,‘S‘)
% 说明:
% 当X和Y为向量时,长度必须相同,分别用于存储x坐标和y坐标数据。
% 当X是向量,Y是有一维与X同维的矩阵时,则绘制出多根不同颜色的曲线。曲线条数等于Y矩阵的另一维数,X被作为这些曲线共同的横坐标。
% 示例:
x=0:0.1:2*pi;
y=[sin(x);cos(x)];
figure(1)
plot(x,y(1,:))
figure(2)
plot(x,y)

% 格式三:
% plot(X1,Y1,′S1′,X2,Y2,′S2′ ,......,X3,Y3,′S3′)
% 说明:
% 当输入参数都为向量时,X1和Y1,X2和Y2,…,Xn和Yn分别组成一组向量对,
% 每一组向量对的长度可以不同。每一向量对可以绘制出一条曲线,这样可以
% 在同一坐标内绘制出多条曲线。
% 示例:
x=0:0.1:2*pi;
y=[sin(x);cos(x)];
plot(x,y(1,:),′-r′,x(10:60),y(2,10:60),′--ko′)

% 格式四:绘制 y = f(X)图形
% fplot(fname,lims,‘S‘)
% 说明:
% 其中fname为函数名或单引号界定的函数表达式,
% lims为x,y的取值范围,′S′定义与plot函数相同。
% 示例:
fplot(‘cos(tan(pi*x))‘,[ 0,1],‘-ro‘)

% 格式五:绘制隐函数图形 f(x,y)= 0
% 方式一:ezplot(f,[xmin,xmax,ymin,ymax])
% 说明:
% 在区间xmin<x<xmax和ymin<y<ymax绘制f(x,y) = 0的图形。,默认区间-2π<x<2π和-2π<y<2π
% 方式二: x=X(t) , y = Y(t)
% ezplot(X,Y,[tmin,tmax])
% 说明:
% 在区间tmin < t < tmax 绘制 x = X(t)和y = Y(t)的图形,默认区间0< t < 2PI
% 方式三: y = f(x)
% ezplot(f,[a,b])
% 说明:
% 在区间a<x<b绘制y= f(x)的图形,默认区间-2π<x<2π
% 示例:
ezplot(‘x^2+y^2-9‘,[-3 3 -3 3]),axis equal
ezplot(‘8*cos(t)‘,‘4*sqrt(2)*sin(t)‘,[0,2*pi])
ezplot(‘cos(tan(pi*x))‘,[ 0,1])

% 特殊坐标二位图形
% 1.极坐标曲线
% 格式: polar(theta,rho,‘S‘)
% theta:角度向量,rho:幅值向量,‘S’控制参数
% 示例:
theta=0:0.1:8*pi;
r=cos(4*theta)+1/4;
polar(theta,r)

% 对角坐标曲线 (用法与plot相同)
% 函数名 功能
% semilogx x轴对数坐标,y轴线性坐标
% semilogy x轴对数坐标,y轴线性坐标
% loglog x y轴均为对数坐标
% 示例:
x=0:0.01:5;
y=10.^x;
plot(x,y),grid on

x=0:0.01:5;
y=10 .^x;
semilogy(x,y)
grid on

% 二维图形处理
% 标注方式:
% title(′图形名称′)
% xlabel(′x轴说明′)
% ylabel(′y轴说明′)
% text(x,y,′图形说明′)
% legend(′图例1′,′图例2′,...)
%
% 函数中的说明文字,除使用标准的ASCII字符外,还可使用LaTeX格式的控制字符,这样就可以在图形上添加希腊字母、数学符号及公式等内容。例如,text(0.3,0.5, ′sin({\omega}t+{\beta}) ′)将得到标注效果sin(ωt+β)。
% 上述函数除legend外,均可以用于三维函数。
%
% 示例:
x=[0:0.1:2*pi];
plot(x,sin(x),x,cos(x))
title(′正弦和余弦曲线′);
xlabel(′θ(弧度)′)
ylabel(′y′)
text(3,0.5,′sinθ cosθ 曲线′);
legend(′sinθ′,′cosθ′);
%
% 坐标控制:axis函数
% 主要格式
% axis([xmin xmax ymin ymax zmin zmax])
% axis equal:纵、横坐标轴采用等长刻度。
% axis auto:使用缺省设置。
% axis off:取消坐标轴。
% axis on:显示坐标轴。
%
% grid on/off:控制是否画网格线。
% hold on/off 控制是否刷新当前轴及图形
%
% 示例:
t=0:0.1:2*pi;
plot(sin(t),cos(t))
hold on
plot(t/2/pi,sin(t),‘r‘)
plot(-t/2/pi,sin(t),‘r‘)
grid on;
axis equal
%
% 二.三维图形的基本函数
% 1.三维曲线图
% plot3函数
% 格式:Plot3(x1,y1,z1, ′S1′, x2,y2,z2, ′S2′, …xn,yn,zn, ′Sn′)
% 示例:画画螺旋线
t=-pi:0.1:8*pi;
x=sin(t);
y=cos(t);
plot3(x,y,t,′-r′)
xlabel(′sin(t)′);
ylabel(′cos(t)′);
zlabel(′t′)
%
% 2.三维网格图
% meshgrid函数:产生平面区域内的网格坐标矩阵
% 格式: [X,Y]=meshgrid(A,B)
% 说明:
% 语句执行后,矩阵X的每一行都是向量A,行数等于向量B的元素的个数,矩阵Y的每一列都是向量B,列数等于向量A的元素的个数。
%
% mesh函数
% 格式:mesh(x,y,z)
% 说明:一般情况下,x,y,z是维数相同的矩阵。x,y是网格坐标矩阵,z是网格点上的高度矩阵。
% 示例:三维网格图
[x,y]=meshgrid(0:0.08:2*pi);
z=sin(x).*cos(y);
figure(1)
mesh(x,y,z)
xlabel(‘x‘),ylabel(‘y‘)
zlabel(‘sin(x)cos(x)‘)
grid on,box on
figure(2)
mesh(z),box on
%
% 3.三维曲面图
% 格式:surf(x,y,z)
% 示例:
[x,y]=meshgrid(0:0.08:2*pi);
z=sin(x).*cos(y);
surf(x,y,z)
zlabel(′sin(x)cos(x)′)
grid on
box on
% 示例:显示夫琅禾费衍射图样
%
w=3e-5;
lamda=6.328e-7;
z=1e5;
A=pi*w^2;
k=2*pi/lamda;
a=linspace(-4000,4000,401);
[x,y]=meshgrid(a);
r=sqrt(x.^2+y.^2);
m=k*w.*r/z;
I=(A/lamda/z)^2 * ( 2*besselj(1,m)./m ).^2;
J=max(max(I));
I=I./J;
figure(1)
surf(x,y,I)
xlabel(‘x‘),ylabel(‘y‘),zlabel(‘I‘)
shading(‘interp’)%色彩效果
colormap(jet),colorbar
axis([-4000 4000 -4000 4000 0 1])

figure(2)
plot(x(1,:),I(201,:))
xlabel(‘x‘),ylabel(‘I‘)
axis([ -4000 4000 0 1])
%
% 4.图形窗口和分割
% 在实际应用中,有时需要在不同图形窗口或一个图形窗口中绘制若干个独立的图形,
% 这就需要选取不同的图形窗口或对图形窗口分割。
% figure函数
% 格式: figure(n)
% 说明:
% 该函数打开不同的图形窗口。
% n为图形窗口排序号。
% 默认时打开的是1号图形窗,即当前窗。
%
% 示例:
t=linspace(0,2*pi,50);
figure(1)
plot(t,sin(t))
figure(2)
plot(t,cos(t))

% 在实际应用中,有时需要在一个图形窗口中绘制若干个独立的图形,这就需要对图形窗口分割。
% subplot函数
% 格式:subplot(m,n,p)
% 说明:
% 该函数将当前图形窗口分成m×n个绘图区,即每行n个,共m行。
% 区号按行优先编号,且选定第p个区为当前活动区。
% 在每一个绘图区允许以不同的坐标系单独绘制图形。
%
% 示例:
t=0:0.1:2*pi;
subplot(2,2,1)
plot(t,sin(t),’-r’)
title(′x=sin(t)′);
xlabel(′t′)
ylabel(′sin(t)′)
subplot(2,2,2);
plot(t,cos(t),′-b′)
title(′y=cos(t)′);
xlabel(′t′)
ylabel(′cos(t)′)
subplot(2,2,3);
[X,Y]=meshgrid(0:0.1:2*pi);
Z=sin(X).*cos(Y);
surf(X,Y,Z)
zlabel(′sin(t)cos(t)′)
title(′z=sin(t)*cos(t)′);
subplot(2,2,4)
mesh(X,Y,Z)
xlabel(′sin(t)′)
ylabel(′cos(t)′)
zlabel(′sin(t)cos(t)′)
title(′z=sin(t)*cos(t)′);

% 5.简单动画
% (1)质点运动轨迹的显示
% comet函数
% 格式:
% comet(x)
% comet(x,y)
% comet(x,y,p)
% 说明:
% comet(x)显示质点绕向量x运动
% comet(x,y)显示质点绕向量x与y运动
% comet(x,y,p),其中为p轨迹尾巴的长度,缺省值为P=0.1
% 示例:模拟上抛运动
% V(x)=100cos(pi/4)
% V(y)=100sin(pi/4)
% x = V(x)*t
% y = V(y)*t - 1/2*g*t*t
%
vx = 100*cos(1/4*pi);
vy = 100*sin(1/4*pi);
t = 0:0.001:15;
x = vx*t;
y = vy*t-9.8*t.^2/2;
comet(x,y)

% (2)以电影播放的方式显示动画
% 保存想要产生动画的图片,存储为一系列各种类型的二维、三维图,再像放电影的方式按次序播放出来
% 步骤:
% 由getframe函数将当前的图片抓取为电影的画面,存入矩阵中。
% 再由movie函数将矩阵 中的画面以动画显示出来
% 格式: m(j) = getframe
% movie(m)
% 示例:演示驻波
a=50;w=0.2;v=3.64;
for t=1:100
x=0:1:80*pi;
y=a*cos(w*(t-x./v)+pi/2);
z=a*cos(w*(t+x./v));
plot(x,y,x,z,x,y+z);
axis([0 80*pi -100 100]);
grid on;
m(t)=getframe;
end
movie(m)

%
%
% MATLAB的工作模式:
% 1.指令驱动模式
% 通常matlab以指令驱动模式工作,即在matlab 窗口下当用户输入单行指令时,matlab 立即处理这条指令,并显示结果,这就是matlab命令行模式。
% ?命令行模式操作时,matlab窗口只允许一次执行一行上的一个或几个语句
% 在matlab窗口输入数据和命令进行计算时,当处理复杂问题和大量数据时是不方便的。
% 命令行方式程序可读性差,而且不能存储,对于复杂的问题,应编写成能存储的程序文件。
%
% 2.m文件模式
% m文件是matlab所特有的使用该语言编写的磁盘文件。
% 将matlab语句构成的程序存储成以m为扩展名的文件,然后再执行该程序文件,这种工作模式称为程序文件模式。
% 程序文件不能在指令窗口下建立,因为指令窗口只允许一次执行一行上的一个或几个语句。
%
% M文件有两类:
% 1.脚本文件(Script File) : 独立的m文件
% 2.函数文件(Function File): 可调制的m文件
% 脚本文件是一串matlab命令的集合,完成制定的功能,变量空间是workspace;
% 脚本文件既不接受输入参数也不返回输出参数,脚本文件实际上是一串指令的集合,与在命令窗口逐行执行文件中的所有指令,其结果是一样的。
% 函数文件是一个黑箱,根据输入作出输出,变量空间是独立的函数变量空间,在函数运行完成后关闭。
%
% M文件建立:
% 1.命令行输入edit
% 2.新建文件
%
% 1.脚本文件包括两部分:注释部分和程序部分
% a.m
% %绘制sin(x)*cos(x)图形
% x=0:0.1:4*pi;
% y=sin(x).*cos(x);
% plot(x,y)
% %注释部分,用以说明函数的作用及有关内容,使用help命令时显示。只显示程序中的第一句注释语。
%
% 2.函数文件:指能够接受并输出参数的m文件
% 函数文件的格式:
% function 输出形参表=函数名(输入形参表)
% 注释说明语句段
% 程序语句段
% 说明:
% 1.函数m文件第一行必须以单词function作为引导词。
% 2.函数文件的文件命名规则与变量相同,必须是函数名 .m。
% 3.当输出形参多于一个时,须用方括号括起来,以矩阵形式表示。
% 4.程序中的变量均为局部变量,不保存在工作空间中。其变量只在函数运行期间有效。
% 示例:
% peri_area.m
function [l s]=peri_area(r)
%计算圆的周长和面积
l=2*pi*r; %周长
s=pi*r.^2; %面积
[a,b] = per_area(3);
%
% MATLAB的程序结构:(与C语言类似)
% 1.顺序结构
% 2.选择结构
% 3.循环结构
%
% 一.顺序结构
% (1)数据的输入
% input 函数
% 调用格式
% A = input(‘提示信息‘.‘选项‘)
% 说明:选项用于控制输入数据格式,‘s’允许输入字符串。
% 除了字符串输入需要加‘s’,其余的一律不加(如果不加想要输入字符串,则需要‘’包裹输入的字符串)
%
% (2)数据的输出
% disp 函数
% 调用格式
% disp(输出项)
% 说明:
% 输出项可以是变量或字符串。
% 输出变量时,不显示变量名。
%
% sprintf函数
% 调用格式
% S = sprintf(格式,输出项)
% 说明:
% 输出格式控制同c语言
% 输出项可以是变量或表达式
% 示例:输出格式控制
x=pi;m=2;
sprintf(‘x1=%f x2=%4.3f\nx3=%10.9f‘, (1+sqrt(5))/2,x,pi)
sprintf(‘%d’,round(x)) %不显示小数点后面的数字
sprintf(‘%s‘,‘hello‘)
S=sprintf(‘The array is %dx%d.‘,m,3) ;
disp(S)

% 重点再次回顾:
%解一元二次方程
%a b c方程系数
%x1 x2为方程的根
disp(‘pls input a b c‘)
a=input(‘a=‘);
b=input(‘b=‘);
c=input(‘c=‘);
d=b*b-4*a*c;
x1=(-b+sqrt(d))/(2*a);
x2=(-b-sqrt(d))/(2*a);
disp(‘x1 x2 is‘)
disp([x1,x2])
disp(sprintf(‘x1=%6.5f\nx2=%6.5f‘,x1,x2))

% (3)程序的暂停
% pause函数
% 调用格式
% pause(延迟秒数)
% 说明:
% 如果省略延迟时间,直接使用pause,则将暂停程序,直到用户按任一键后程序继续执行。
% 若要强行中止程序的运行可使用Ctrl+C命令。
%
% 二。选择结构
% (1)if语句
if 条件1
语句组1
elseif 条件2
语句组2
……
elseif 条件m
语句组m
else
语句组n
end
% (2)switch语句
switch 表达式
case 表达式1
语句组1
case 表达式2
语句组2
……
case 表达式m
语句组m
otherwise
语句组n
end
% try语句
try
语句组1
catch
语句组2
end

% 示例
A=[1,2,3;4,5,6]; B=[7,8,9;10,11,12];
try
C=A*B;
catch
C=A.*B;
end
C
erro=lasterr%显示出错原因

% 三。循环结构
% (1)for语句
% 格式
% for 循环变量=表达式1:表达式2:表达式3
% 循环体语句
% end
% 说明:
% 表达式1的值为循环变量的初值,表达式2的值为步长,表达式3的值为循环变量的终值。
% 步长为1时,表达式2可以省略。
% (2)while语句
% while(条件)
% 循环体语句
% end
% 说明:若条件成立,则执行循环体语句,执行后再判断条件是否成立,如果不成立则跳出循环
%
% break和continues语句
% break:用于终止循环的执行。
% continue: 直接进行下一次循环
% 示例:
%计算100~200之间第一个能被21整除的整数
for n=100:200
if rem(n,21)~=0 %rem(x,y) :计算x./y的余数
continue
end
break
end
n

% 一、MATLAB文件的打开和关闭
% (1)文件的打开
% 格式: fid = fopen(文件名,打开方式)
% 说明:
% 其中文件名用字符串形式,表示待打开的文件。
% 常见的打开方式有:
% ′r′表示对打开的文件读数据;
% ′r+′表示读写;
% ′w′删除已经存在的文件内容或建立一个新文件,并打开文件写;
% ′w+′删除已经存在的文件内容或建立一个新文件,读写;
% ′a′表示在打开的文件末尾添加数据。
% fid用于存储文件句柄值,句柄值用来标识该数据文件,其它函数可以利用它对该数据文件进行操作。
%
% (2)文件的关闭
% 格式 sta = fclose(fid)
% 说明:该函数关闭fid所表示的文件。
% sta表示关闭文件操作的返回代码,若关闭成功,返回0,否则返回–1。
%
% 二、文件的读写操作
% 文件数据格式有两种形式,一是二进制文件,二是文本文件。对不同类型的文件读写是不同的。
%
% (1)二进制文件的读写操作
% 读二进制文件
% 格式 [A,COUNT]=fread(fid,size, precision)
% 说明:
% A用于存放读取的数据。
% COUNT返回所读取的数据元素个数
% fid为文件句柄
% precision代表读写数据的类型:′int′ ′float′ ′char′等
% size为可选项,若不选用则读取整个文件内容,若选用则它的值可以是下列值: (
% 1) N表示读取 N个元素到一个列向量。
% (2) [M,N]表示读数据到M×N的矩阵中,数据按列存放。
%
% 写二进制文件
% 格式: COUNT=fwrite (fid, A, precision)
% 说明:
% COUNT返回所写的数据元素个数。
% fid为文件句柄。
% A用来存放写入文件的数据,
% precision用于控制所写数据的类型,其形式与fread函数相同
%
% (2)文本文件的读写操作
% 读文本文件
% 格式 A = fscanf(fid,format)
% [A,COUNT] = fscanf(fid,format,size)
% 说明:
% A用以存放读取的数据。
% COUNT返回所读取的数据元素个数。fid为文件句柄。
% format用以控制读取的数据格式,由%加上格式符组成,常见的格式符有%d,%f,%e,%c,%s等。
% size为可选项,决定矩阵A中数据的排列形式。
%
% 写文本文件
% 格式: COUNT= fprintf(fid, format, A)
% 说明:
% A存放要写入文件的数据。
% 先按format指定的格式将数据矩阵A格式化,然后写入到fid所指定的文件。格式符与fscanf函数相同。
%
% 示例:文件的操作
clear
x=0:0.1:2*pi;
fid=fopen(‘example.mat‘,‘w+‘)
count=fwrite(fid,x,‘float‘)
fclose(fid)
fid1=fopen(‘example.mat‘,‘r+‘)
t=fread(fid1,‘float‘);
y=[t,sin(t)];
plot(t,y(:,2))
fclose(fid1)
%
%
% 二、函数的调用
% 格式
% [输出实参表]=函数名(输入实参表)
% 说明:
% 函数调用时各实参出现的顺序、个数,应与函数定义时形参的顺序、个数一致,否则会出错。
% 函数调用时,先将实参传递给相应的形参,从而实现参数传递,然后再执行函数的功能。
% 示例:利用函数文件,实现直角坐标和极坐标的转换

%定义函数文件tran.m:
function [rho,theta]=tran(x,y)
rho=sqrt(x*x+y*y);
theta=atan(y/x);


%调用tran.m的命令文件exam406.m
x=input(‘Please input x=:‘);
y=input(‘Please input y=:‘);
[rho,theta]=tran(x,y);
rho
theta

%
% 函数参数的可调性
% nargin和nargout
% 说明:
% 在调用函数时,用nargin和nargout分别记录调用该函数时的输入实参和输出实参的个数。
% 只要在自定义函数文件中包含这两个函数,就可以准确地知道该函数文件被调用时的输入输出参数个数,从而决定函数如何进行处理。
% nargin和nargout用法示例:
%定义函数文件sumproduct.m:
function [out1,out2]=sumproduct(x,y,z)
if nargin==0
disp(‘no input arguments‘),return
elseif nargin==1
sum=x;product=x;
elseif nargin==2
sum=x+y;product=x.*y;
elseif nargin==3
sum=x+y+z;product=x.*y.*z;
else
disp(‘too many input arguments‘)
end

if nargout==0
return
elseif nargout==1
out1=sum;
else
out1=sum;out2=product;
end

%调用sumproduct.m
sumproduct
[x,y]=sumproduct(1,2,3)
[x,y]=sumproduct(1)
x=sumproduct(1,2,3)
[x,y]=sumproduct(1,2,3,4)


% 定义全局变量
% globle X
% 说明:
% 说明:
% 全局变量名一般用大写。
% 自定义函数文件中的变量,作用区域仅在函数内部。对于自定义函数中的一些公用变量,可以定义成全局变量,在主程序中统一定义.
% global用法示例:
function y=Ep(x)
%计算弹性势能
global K
y=K*x.^2/2;

function f=F(x)
%计算弹力
global K
f=K*x;

%计算exam408.m
global K
K=0.01;
x=input(‘Please input x=‘);
ep=Ep(x)
f=F(x)

% 来自臧小飞老师在最后留下的话:
% 良好的编程习惯:
% 在动手编程之前,明确程序的目的,设想解决方案,作出初步的流程图。如果程序较大,就要把程序分成几个相对独立的模块,各司其职。一个一个模块解决。
% %后面的内容是程序的注解,要养成注释程序的习惯,在关键的命令行,变量定义处必须要有注释,在整个程序开头有一个总结性的注释。这样便于自己或别人查看和修改程序。
% 变量命名规则要统一, 含义清晰。
% 编辑m文件要注意排版,这可以使程序层次变得更清晰,有序,增加程序的可读性。
% 养成在主程序开头用clear指令清除变量的习惯,以消除工作空间中其它变量对程序运行的影响。但注意在函数或子程序中不要用clear。
% 参数值要集中放在程序的开始部分,以便维护。在语句行之后输入分号使其中间结果不在屏幕上显示,以提高执行速度。
% input指令可以用来输入一些临时的数据;而对于大量参数,则通过建立一个存储参数的子程序,在主程序中用子程序的名称来调用。
% 程序尽量模块化,也就是采用主程序调用子程序的方法,将所有子程序合并在一起来执行全部的操作。
% 设置好MATLAB的工作路径,以便程序运行。
%
% 提升MATLAB的运行效率
% 1. 尽量避免使用循环
% a.尽量用向量化的运算来代替循环操作。
% b.在必须使用多重循环时下,则在循环的外环执行循环次数少的,内环执行循环次数多的。这样可以显著提高速度。
% 2.预分配矩阵空间,即事先确定变量的大小、维数。
% 这一类的函数有zeros、ones等。
% 3. 将耗时的循环调用C或fortrun等低级语言运算
% 4.改用更有效的算法
%
% 最后一道题:输出100-200之间第三个被15整除的数
% x=100:200; n=find(rem(x,15)==0); X=x(n(3))
%
%
% 往年期末考试习题:
% 以两种方式求线性方程组的解:
% x1 + 2*x2 - 7*x3 + 6*x4 = 0
% 2*x2 + x3 + x4 = 6
% 4*x2 - 7*x3 + 2*x4 = 7
% 2*x2 - 2*x3 + x4 = 4
% X=A\B
% 在一个图形窗口中绘制直线、圆、极坐标和曲面图四个独立的图形,数据自设,要求图形标注完整。
% 实现上抛和平抛运动的动画演示(速度,时间自设)。
% 输出[100,300]之间第3个能被13整除的整数,另将该 数值范围内所有能被13整除的整数输出。
% Created by Andrew Duan 2016.10.17 00:27:52

MATLAB入门笔记