--matlab的使用

四. 函数文件
• 函数文件是另一种形式的M文件，每一个函数文件都 定义一个函数。事实上，MATLAB提供的标准函数大部 分都是由函数文件定义的。
1．函数文件格式 • 函数文件由function语句引导，其格式为： • function 输出形参表=函数名(输入形参 表) • 注释说明部分 • 函数体
二、矩阵与数组

1 创建 1）A=[3,4,5;sin(pi/6) 4^2 9] 2) 单位阵eye(m,n) rand(m,n) zeros(m,n)
3) 向量A=[1:0.1:3]
2 矩阵的操作
1）提取：a) A=rand(3,4)； A2=A(2,:)
b) A3=A(1:2,[1,3,4])
二．数据的输入输出
• input函数：用于向计算机输入一个参数。
• 调用格式： A=input(提示信息，选项)；
• 【例2】 求一元二次方程a2 +bx+c=0的根。
• a=input('a=?'); • b=input('b=?'); • c=input('c=?'); • d=b*b-4*a*c; • x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)] • 将该程序以aa.m文件存盘，然后运行aa.m文件。
提取12行，134列
2 矩阵的操作
2）运算：a) 矩阵加法 A+B，A+3 b) 矩阵乘法 A*B 向量乘法 A.*B c) 除法 inv(A) 左除右除’/’ ’\’
d) 乘方 A^2
三 标量函数 向量函数 矩阵函数
1 标量函数 1）三角函数 sin cos asin acot 2) 标量函数 sqrt exp log log10 abs 3) Round floor ceil sign 取整 +整 –整 符号
四. 函数文件
【例11】 利用函数的递归调用，求n！。
• • • • • • • • • • • • function f=factor(n) if n<=1 f=1; else f=factor(n-1)*n; end return; %返回 在 命 令 文 件 main2.m 中 调 用 函 数 文 件 factor.m： for i=1:10 fac(i)=factor(i); end fac
University of Science and Technology of China
L2
L1
a c
w2
w1
w1 w2 L sin( a c) sin c
华北科技学院
DEPARTMENT OF MATHEMATICS
University of Science and Technology of China
• 注：其中函数名的命名规则与变量名相同。输入形参 为函数的输入参数，输出形参为函数的输出参数。当 输出形参多于1个时，则应该用方括号括起来。
四. 函数文件
【例10】 利用函数文件，实现直角坐标(x,y)与 极坐标(γ,θ)之间的转换。
已知转换公式为：
极坐标的矢径： γ= 极坐标的幅角： θ= • 函数文件tran.m： • function [gama,theta]=tran(x,y) • gama=sqrt(x*x+y*y); • theta=atan(y/x); • 调 用 tran.m 的 命 令 文 件 main1.m： • x=input('Please input x=:'); • y=input('Please input y=:'); • [gam,the]=tran(x,y); • gam • the
三．程序设计
(二)循环结构
• 2． while语句 • 格式为： • while (条件) • 循环体语句 • end n=1; while(prod(1:n)<99^99) %prod 是所有元素乘 n=n+1; end n
三．程序设计
(二)循环结构
3． 循环的嵌套 • 如果一个循环结构的循环体又包括一个循 环结构，就称为循环的嵌套，或称为多重循 环结构。 • 多重循环的嵌套层数可以是任意的。可以 按照嵌套层数，分别叫做二重循环、三重循 环等。处于内部的循环叫作内循环，处于外 部的循环叫作外循环。
• 运行： • 请输入三角形的三条边：[4 5 6] 9.9216
三．程序设计
(一)选择结构
【例5】 输入一个字符，若为大写字母，则输出其后继 字符，若为小写字母，则输出其前导字符，若为其他字 符则原样输出。
• c=input('','s'); • if c>='A' & c<='Z' • disp(setstr(abs(c)+1)); • elseif c>='a'& c<='z' • disp(setstr(abs(c)-1)); • else • disp(c); • end
• • • • • a=1:9; b=[11,12,13;14,15,16;17,18,19]; c=a;a=b;b=c; a b 在MATLAB的命令窗口中输入e31，将会执行该命令文件。
【例2】 求一元二次方程ax2 +bx+c=0的根。
• • •
funcion y=qiugen(a,b,c) d=b*b-4*a*c; y=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];
三．程序设计
（二）循环结构
• 实现循环结构的语句：for语句和while语句。
• 1． for语句：
• 格式： for 循环变量=表达式1:表达式2:表达式3 • 循环体语句 • end
• 注：其中表达式1的值为循环变量的初值，表达式2的值为 步长，表达式3的值为循环变量的终值。步长为1时，表达 式2可以省略。
方法一： A=123*pi/180; L=inline(‘9/sin(pi-2.1468-c)+7/sin(c)’); fplot(L,[0.4,0.5]);grid on 方法二：x=fminbnd(L,0.4,0.5);L(x) 方法三：fminbnd(L,0.4,0.5,optimset(’Display’,’iter’))
三．程序设计
(二)循环结构 • n=10;%求10!的表达式 • a=1; • for i=1:n • a=a*i; • end • a ------------------------------------------------• n=4; %求等比 • for i=1:n • b(i)=2^i; • end • b
三．程序设计
(一)选择结构
【例4】 输入三角形的三条边，求面积。
• function y=mianji(A)
• if A(1)+A(2)>A(3) & A(1)+A(3)>A(2) & A(2)+A(3)>A(1) • p=(A(1)+A(2)+A(3))/2; • s=sqrt(p*(p-A(1))*(p-A(2))*(p-A(3))); • disp(s); • else • disp('不能构成一个三角形。') • end • s
程序运行结果是： fac = Columns 1 through 6 1 2 6 24 120 720 Columns 7 through 10 504 40320 362880 3628800
MatLab 的使用
一、基本操作与简单语句输入

求 y=a b+a b(a=3,b=4) 的值。

语句表达一般式： 变量=表达式
一、基本操作与简单语句输入

>>a=3;b=4; >>y=a*b+a/b

注意：“;”有无的区别！
一、基本操作与简单语句输入
练习
4
a 5sin c (a 64, c 5)
方法
华北科技学院
DEPARTMENT OF MATHEMATICS
University of Science and Technology of China
MatLab 的介绍
例、用简短命令计算并绘制在0x6范围内的 sin(2x)、sinx2、sin2x。 x=0:0.01:6; y1=sin(2*x);y2=sin(x.^2);y3=(sin(x)).^2; plot(x,y1,x, y2,x, y3)
原始图像lena.bmp
第一次猫脸变换
第 2次
第 3次
第 4次
第 5次
程序：需要一个取模运算





I=imread('f:/lena.bmp'); for i=1:512 for j=1:512 a=[1,1;1,2]*[i;j]; b=mod(a(1),512)+1;c=mod(a(2),512)+1; II(i,j)=I(b,c); end end imshow(II)
1 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.8 -1 0 1 2 3 4 5 6
Arnold(猫脸变换)置乱 : 设I为一幅灰度图像,也可 以看成一个二维数 组,I(i,j)是它的每个元素,令 i' 1 1 i II(i,j)=I(i’,j’), = , j' 1 2 j II是变换后的图像(imread,imshow):
三 标量函数 向量函数 矩阵函数