matlab例题

例7 输入x,y的值，并将它们的值互换后输出。

程序如下：

x=input('Input x please.');

y=input('Input y please.');

z=x;

x=y;

y=z;

disp(x);

disp(y);

例8 求一元二次方程ax2 +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)];

disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);

例9 计算分段函数的值。

程序如下：

x=input('请输入x的值:');

if x<=0

y= (x+sqrt(pi))/exp(2);

else

y=log(x+sqrt(1+x*x))/2;

end

例10 输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出。

c=input('请输入一个字符','s');

if c>='A' & c<='Z'

disp(setstr(abs(c)+abs('a')-abs('A')));

elseif c>='a'& c<='z'

disp(setstr(abs(c)- abs('a')+abs('A')));

elseif c>='0'& c<='9'

disp(abs(c)-abs('0'));

else

disp(c);

end

例11 某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)：price<200 没有折扣

200≤price<500 3%折扣

500≤price<1000 5%折扣

1000≤price<2500 8%折扣

2500≤price<5000 10%折扣

5000≤price 14%折扣

输入所售商品的价格，求其实际销售价格。

程序如下：

price=input('请输入商品价格');

switch fix(price/100)

case {0,1} %价格小于200

rate=0;

case {2,3,4} %价格大于等于200但小于500

rate=3/100;

case num2cell(5:9) %价格大于等于500但小于1000

rate=5/100;

case num2cell(10:24) %价格大于等于1000但小于2500

rate=8/100;

case num2cell(25:49) %价格大于等于2500但小于5000

rate=10/100;

otherwise %价格大于等于5000

rate=14/100;

end

price=price*(1-rate) %输出商品实际销售价格

例12 矩阵乘法运算要求两矩阵的维数相容，否则会出错。先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。

程序如下：

A=[1,2,3;4,5,6]; B=[7,8,9;10,11,12];

try

C=A*B;

catch

C=A.*B;

end

C

lasterr

例13 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。

程序如下：

for m=100:999

m1=fix(m/100); %求m的百位数字

m2=rem(fix(m/10),10); %求m的十位数字

m3=rem(m,10); %求m的个位数字rem求余数

if m==m1*m1*m1+m2*m2*m2+m3*m3*m3

disp(m)

end

end

例14 已知，当n=100时，求的值。

程序如下：

程序如下:

y=0;

n=100;

for k=1:n

y=y+1/(2*k-1);

end

y

例16 从键盘输入若干个数，当输入0时结束输入，求

这些数的平均值和它们之和。

程序如下：

sum=0;

cnt=0;

val=input('Enter a number (end in 0):');

while (val~=0)

sum=sum+val;

cnt=cnt+1;

val=input('Enter a number (end in 0):');

end

if (cnt > 0)

sum

mean=sum/cnt

end

例17 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求[1,500]之间的全部完数。

for m=1:500

s=0;

for k=1:m/2

if rem(m,k)==0 %rem求m/k的余数

s=s+k;

end

end

if m==s

disp(m);

end

end

例18 求[100，200]之间第一个能被21整除的整数。

程序如下：

for n=100:200

if rem(n,21)~=0

continue

end

break

end

n

例19：函数的递归调用

编写函数nj计算自然数n的阶乘

function y=nj(n)

%y=nj(n)计算n的阶乘

%输入参数n是自然数

%输出参数y是n的阶乘

if(n<0)|(floor(n)~=n)%如果输入参数不是自然数

error('输入参数必须是自然数')%则显示出错信息

end

if n==0

y=1;

elseif n==1

y=1;

else

y=n*nj(n-1);

end

例20 猜数游戏。首先由计算机产生[1,100]之间的随机整数，然后由用户猜测所产生的随机数。根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“Y ou won”，同时退出游戏。用户最多可以猜7次。

程序如下：

%Play the game of guess the number

x=fix(100*rand); %a random number calculated by the computer

n=7;

test=1;

for k=1:7

numb=int2str(n);

disp(['Y ou have a right to ',numb,' guesses'])

disp(['A guess is a number between 0 and 100'])

guess=input('Enter your guess:');

if guess

disp('Low')

elseif guess>x

disp('High')

else

disp('Y ou won')

test=0;

break;

end

n=n-1;

end

if test==1

disp('Y ou lost')

end

例21 ：递推与递归

核反应堆中有A和B两种粒子，每秒中内1个A粒子分裂成

3个B粒子，而1个B粒子分裂成1个A粒子和2个B粒子。

若在时刻t=0反应堆中只有1个A粒子，问t=20时反应堆中

有多少个A粒子？多少个B粒子？共有多少个粒子？

function [outf,outg,outfg]=fg(x)

%函数[outf,outg,outfg]=fg(x)计算A粒子与B粒子的个数

%输入参数x是自然数

%输出参数outf是A粒子个数

%输出参数outg是B粒子个数

%输出参数outfg是A粒子和B粒子的总数

%主函数

[m,n]=size(x);%输入参数x的行数和列数

if (m~=1)|(n~=1)|(x<0)|(round(x)~=x) %如果输入参数不是自然数

y=[]; %输出为空

disp('输入参数应为自然数'); %显示出错

elseif x>=0

outf=f(x);

outg=g(x);

outfg=outf+outg;

end

function y=f(x);

%子函数，计算A粒子

if x==0

y=1;

elseif x>0

y=g(x-1);

end

function y=g(x);

%子函数，计算B粒子

if x==0

y=0;

elseif x>0

y=3*f(x-1)+2*g(x-1);

end

例4-3 在0≤x≤2 区间内，绘制曲线

y=2e-0.5xcos(4πx)

程序如下：

x=0:pi/100:2*pi;

y=2*exp(-0.5*x).*cos(4*pi*x);

plot(x,y)

例4-6 用不同标度在同一坐标内绘制曲线y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)。程序如下：

x=0:pi/100:2*pi;

y1=0.2*exp(-0.5*x).*cos(4*pi*x);

y2=2*exp(-0.5*x).*cos(pi*x);

plotyy(x,y1,x,y2

例4-7 采用图形保持，在同一坐标内绘制曲线y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)。

程序如下：

x=0:pi/100:2*pi;

y1=0.2*exp(-0.5*x).*cos(4*pi*x);

plot(x,y1)

hold on

y2=2*exp(-0.5*x).*cos(pi*x);

plot(x,y2);

hold off

例4-8 在同一坐标内，分别用不同线型和颜色绘制曲线y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)，标记两曲线交叉点。

程序如下：

x=linspace(0,2*pi,1000);

y1=0.2*exp(-0.5*x).*cos(4*pi*x);

y2=2*exp(-0.5*x).*cos(pi*x);

k=find(abs(y1-y2)<1e-2); %查找y1与y2相等点(近似相等)的下标

x1=x(k); %取y1与y2相等点的x坐标

y3=0.2*exp(-0.5*x1).*cos(4*pi*x1); %求y1与y2值相等点的y坐标

plot(x,y1,x,y2,…k:?,x1,y3,…bp?); %y1缺省蓝色实线，y2黑色

点线，y3蓝色五角星

例4-9 在0≤x≤2 区间内，绘制曲线y1=2e-0.5x和y2=cos(4πx)，并给图形添加图形标注。

程序如下：

x=0:pi/100:2*pi;

y1=2*exp(-0.5*x);

y2=cos(4*pi*x);

plot(x,y1,x,y2)

title('x from 0 to 2{\pi}'); %加图形标题

xlabel('Variable X'); %加X轴说明

ylabel('Variable Y'); %加Y轴说明

text(0.8,1.5,'曲线y1=2e^{-0.5x}'); %在指定位置添加图形说明

text(2.5,1.1,'曲线y2=cos(4{\pi}x)');

legend(…y1?,… y2?) %加图例

例4-10 在同一坐标中，可以绘制3个同心圆，并加坐标控制。

程序如下：

t=0:0.01:2*pi;

x=exp(i*t);

y=[x;2*x;3*x]';

plot(y)

grid on; %加网格线

box on; %加坐标边框

axis equal %坐标轴采用等刻度

例4-11 用fplot函数绘制f(x)=cos(tan(πx))的曲线。

命令如下：

fplot('cos(tan(pi*x))',[ 0,1],1e-4)

例4-13 绘制y=10x2的对数坐标图并与直角

坐标图进行比较。

x=1:10;

y=10*(x.^2);

subplot(2,2,1),plot(x,y),title('直角坐标')

subplot(2,2,2),loglog(x,y),title('双对数')

subplot(2,2,3),semilogx(x,y),title('x轴对数')

subplot(2,2,4),semilogy(x,y),title('y轴对数')

例4-12 绘制r=sin(t)cos(t)的极坐标图，并标记数据点。

程序如下：

t=0:pi/50:2*pi;

r=sin(t).*cos(t);

polar(t,r,'-*');

例4-14 分别以条形图、阶梯图、杆图和填充图形式绘制曲线y=2sin(x)。

程序如下：

x=0:pi/10:2*pi;

y=2*sin(x);

subplot(2,2,1);bar(x,y,'g');

title('bar(x,y,''g'')');axis([0,7,-2,2]);

subplot(2,2,2);stairs(x,y,'b');

title('stairs(x,y,''b'')');axis([0,7,-2,2]);

subplot(2,2,3);stem(x,y,'k');

title('stem(x,y,''k'')');axis([0,7,-2,2]);

subplot(2,2,4);fill(x,y,'y');

title('fill(x,y,''y'')');axis([0,7,-2,2]);

例4-15 绘制图形：

(1) 某企业全年各季度的产值(单位：万元)分别为：2347,1827,2043,3025，试用饼图作统计分析。

(2) 绘制复数的相量图：7+2.9i、2-3i和-1.5-6i。

程序如下：

subplot(1,2,1);

pie([2347,1827,2043,3025]);

title('饼图');

legend('一季度','二季度','三季度','四季度');

subplot(1,2,2);

compass([7+2.9i,2-3i,-1.5-6i]);

title('相量图');

例4-17 绘制三维曲线。

程序如下：

t=0:pi/100:20*pi;

x=sin(t);

y=cos(t);

z=t.*sin(t).*cos(t);

plot3(x,y,z);

title('Line in 3-D Space');

xlabel('X');ylabel('Y');

zlabel('Z');

grid on;

例4-18 绘制三维曲面图z=sin(x+sin(y))-x/10。

程序如下：

[x,y]=meshgrid(0:0.25:4*pi);

z=sin(x+sin(y))-x/10;

mesh(x,y,z);

axis([0 4*pi 0 4*pi -2.5 1]);

例4-19 在xy平面内选择区域[-8,8]×[-8,8]，绘制4种三维曲面图。程序如下：

[x,y]=meshgrid(-8:0.5:8);

z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);

subplot(2,2,1);

mesh(x,y,z);title('mesh(x,y,z)')

subplot(2,2,2);

meshc(x,y,z);title('meshc(x,y,z)')

subplot(2,2,3);

meshz(x,y,z)title('meshz(x,y,z)')

subplot(2,2,4);

surf(x,y,z);title('surf(x,y,z)')

例4-20 绘制标准三维曲面图形。

程序如下：

t=0:pi/20:2*pi;

[x,y,z]= cylinder(2+sin(t),30);

subplot(2,2,1);surf(x,y,z);

[x,y,z]=sphere;

subplot(2,2,2);surf(x,y,z);

[x,y,z]=peaks(30);

subplot(2,1,2);surf(x,y,z);

例4-21 绘制三维图形：

(1) 绘制魔方阵的三维条形图。

(2) 以三维杆图形式绘制曲线y=2sin(x)。

(3) 已知x=[2347,1827,2043,3025]，绘制饼图。

(4) 用随机的顶点坐标值画出五个黄色三角形。

程序如下：

subplot(2,2,1);

bar3(magic(4))

subplot(2,2,2);

y=2*sin(0:pi/10:2*pi);

stem3(y);

subplot(2,2,3);

pie3([2347,1827,2043,3025]);

subplot(2,2,4);

fill3(rand(3,5),rand(3,5),rand(3,5), 'y' )

例4-22 绘制多峰函数的瀑布图和等高线图。

程序如下：

subplot(1,2,1);

[X,Y,Z]=peaks(30);

waterfall(X,Y ,Z)

xlabel('X-axis'),ylabel('Y -axis'),zlabel('Z-axis');

subplot(1,2,2);

contour3(X,Y ,Z,12,'k'); %其中12代表高度的等级数

xlabel('X-axis'),ylabel('Y -axis'),zlabel('Z-axis');

例4-23 绘制椭球面

u=-2*pi:0.1:2*pi;

v=-2*pi:0.1:2*pi;

ezmesh('3*cos(u)*sin(v)','2*cos(u)*cos(v)','sin(u)',

[-2*pi,2*pi,-2*pi,2*pi])

例6-31 绘制了peaks 函数曲面并且将它绕z 轴旋转。

程序如下

[X,Y,Z]=peaks(30);

surf(X,Y,Z)

axis([-3,3,-3,3,-10,10])

axis off;

shading interp;

colormap(hot);

m=moviein(20); %建立一个20列大矩阵

for i=1:20

view(-37.5+24*(i-1),30) %改变视点

m(:,i)=getframe; %将图形保存到m 矩阵

end

movie(m,2); %播放画面2次

读取图像信息INFO=imfinfo('FILENAME.FMT')

【例】按256灰度级显示

3cos sin 2cos cos sin x u v

y u v

z u =??=??=?

I=i m r e a d('m o o n.t i f');

i m s h o w(I,256)

c o l o r b a r

【例】按8灰度级显示

I=i m r e a d('m o o n.t i f');

i m s h o w(I,8)

c o l o r b a r

【例】按最大灰度范围显示I=i m r e a d('m o o n.t i f');

i m s h o w(I,[])

c o l o r b a r

【例】按指定灰度范围显示I=i m r e a d('m o o n.t i f');

i m s h o w(I,[64,128])

c o l o r b a r

【例】双精度灰度图像

b w1=z e r o s(100,100);

b w1(2:2:98,2:2:98)=1;

i m s h o w(b w1);

w h o s b w1

【例】8位的灰度图像

b w1=z e r o s(100,100);

b w1(2:2:98,2:2:98)=1;

b w2=u i n t8(b w1);

i m s h o w(b w2,[]);

c o l o r b a r

w h o s b w2

【例】二值图像

b w1=z e r o s(100,100);

b w1(2:2:98,2:2:98)=1;

b w2=u i n t8(b w1);

b w3=b w2~=0;

i m s h o w(b w3);

w h o s b w3

【例】显示二值图像

B W=i m r e a d('c i r c l e s.t i f');

i m s h o w(B W);

【例】取反显示二值图像

B W=i m r e a d('c i r c l e s.t i f');

i m s h o w(~B W);

【例】显示二值图像

c l f

B W=i m r e a d('c i r c l e s.t i f');

i m s h o w(B W,[100;001]);

【例】显示索引图像

[X,M A P]=i m r e a d('c a n o e.t i f');

i m s h o w(X,M A P);

c o l o r b a r

【例】显示伪彩色图像

I=i m r e a d('m o o n.t i f');

i m s h o w(I,j e t(256))

c o l o r b a r

显示真彩色图像。

【例】

r g b=i m r e a d('f l o w e r s.t i f');

i m s h o w(r g b);

例如：显示多帧图像数组I的第七帧

命令为：i m s h o w(I(:,:,:,7))

c l e a r a l l

c l c

c l f

m r i=u i n t8(z e r o s(128,128,1,27));

%初始化一个27帧的图像数组

f o r f r a m e=1:27

[m r i(:,:,:,f r a m e),m a p]=i m r e a d('m r i.t i f');

%把对应的每一帧读入内存中

e n d

i m s h o w(m r i(:,:,:,7),m a p)%显示图像中的第七帧

显示多帧图像示例：

c l e a r a l l

c l c

c l f

m r i=u i n t8(z e r o s(128,128,1,27));

%初始化一个27帧的图像数组

f o r f r a m e=1:27

[m r i(:,:,:,f r a m e),m a p]=i m r e a d('m r i.t i f');

%把对应的每一帧读入内存中

e n d

m o n t a g e(m r i,m a p)%显示图像中的所有帧

显示动画

i m m o v i e函数可以从一个多帧图像数组中创建m a t l a b动画

例如：

c l e a r a l l

c l c

c l f

m r i=u i n t8(z e r o s(128,128,1,27));

%初始化一个27帧的图像数组

f o r f r a m e=1:27

[m r i(:,:,:,f r a m e),m a p]=i m r e a d('m r i.t i f');

%把对应的每一帧读入内存中

e n d

mov=immovie(mri,map);%创建动画帧结构

movie(mov)%播放动画

多图显示多个调色板的图像。

【例】

l o a d t r e e s;

[x2,m a p2]=i m r e a d('f o r e s t.t i f');

s u b p l o t(2,1,1),s u b i m a g e(X,m a p);c o l o r b a r

subplot(2,1,2),subimage(x2,map2);colorbar

【例】将t i f图像保存为j p g图像

[x,m a p]=i m r e a d('c a n o e.t i f');

i m w r i t e(x,m a p,'c a n o e.j p g','J P G','Q u a l i t y',75)

1、索引图像

假设：X8－表示8位的索引图像，X16－表示16位的索引图像，X64－表示双精度的索引图像，则有：

X64=d o u b l e(X8)+1－u i n t8格式转d o u b l e格式

X64=d o u b l e(X16)+1－u i n t16格式转d o u b l e格式

X8=u i n t8(X64-1)－d o u b l e格式转u i n t8格式

X16=u i n t16(X64-1)－d o u b l e格式转u i n t16格式

X8=u i n t8(X16)－u i n t16格式转u i n t8格式

X16=u i n t16(X8)－u i n t8格式转u i n t16格式

2、灰度图像

G=m a t2g r a y(A,[A m i n,A m a x])

A m i n A m a x默认为A的最小和最大值返回是个D o u b l e类型的值

假设：I8－表示8位的灰度图像，I16－表示16位的灰度图像，I64－表示双精度的灰度图像，则有：

I64=d o u b l e(I8/255)－u i n t8格式转d o u b l e格式

I64=d o u b l e(I16/65535)－u i n t16格式转d o u b l e格式

I8=u i n t8(r o u n d(I64*255))－d o u b l e格式转u i n t8格式

I16=u i n t16(r o u n d(I64*65535))－d o u b l e格式转u i n t16格式

I8=u i n t8(r o u n d(I16/257))－u i n t16格式转u i n t8格式

I16=u i n t16(r o u n d(I8*257))－u i n t8格式转u i n t16格式

3、真彩色图像

假设：R G B8－表示8位的真彩图像，R G B16－表示16位的真彩图像，R G B64－表示双精度的真彩图像，则有：

R G B64=d o u b l e(R G B8/255)－u i n t8格式转d o u b l e格式

R G B64=d o u b l e(R G B16/65535)－u i n t16格式转d o u b l e格式

R G B8=u i n t8(r o u n d(R G B64*255))－d o u b l e格式转u i n t8格式

R G B16=u i n t16(r o u n d(R G B64*65535))－d o u b l e格式转u i n t16格式

R G B8=u i n t8(r o u n d(R G B16/257))－u i n t16格式转u i n t8格式

R G B16=u i n t16(r o u n d(R G B8*257))－u i n t8格式转u i n t16格式

4、二值图像

G=i m2b w(f,T)T默认为0.5

假设B W8表示8位二值图像，B W64表示双精度二值图像，则有：

B W64=d o u b l e(B W8)u i n t8格式转d o u b l e格式

B W8=u i n t8(B W64)~=0d o u b l e格式转u i n t8格式

1、真彩图像→索引图像

【格式】X=d i t h e r(R G B,m a p)

【例】

C L F

R G B=i m r e a d('f l o w e r s.t i f');

m a p=j e t(256);

X=d i t h e r(R G B,m a p);

s u b p l o t(1,2,1);s u b i m a g e(R G B);t i t l e('真彩图')

s u b p l o t(1,2,2);s u b i m a g e(X,m a p);t i t l e('索引图')

2、索引图像→真彩图像

【格式】R G B=i n d2r g b(X,m a p)

【例】

C L F,l o a d t r e e s;

R G B=i n d2r g b(X,m a p);

s u b p l o t(1,2,1);s u b i m a g e(X,m a p);t i t l e('索引图')

s u b p l o t(1,2,2);s u b i m a g e(R G B);t i t l e('真彩图')

3、真彩图像→灰度图像

【格式】I=r g b2g r a y(R G B)

【例】

C L F,R G B=i m r e a d('f l o w e r s.t i f');

I=r g b2g r a y(R G B);

s u b p l o t(1,2,1);s u b i m a g e(R G B);t i t l e('真彩图')

s u b p l o t(1,2,2);s u b i m a g e(I);t i t l e('灰度图')

4、真彩图像→二值图像

【格式】B W=i m2b w(R G B,l e v e l)

【例】

C L F,R G B=i m r e a d('f l o w e r s.t i f');

B W=i m2b w(R G B,0.5);

s u b p l o t(1,2,1);s u b i m a g e(R G B);t i t l e('真彩图')

s u b p l o t(1,2,2);s u b i m a g e(B W);t i t l e('二值图')

5、索引图像→灰度图像

【格式1】I=i n d2g r a y(X,m a p)

【格式2】N e w m a p=r g b2g r a y(m a p)

【例】

C L F,l o a d t r e e s;

I=i n d2g r a y(X,m a p);

N e w m a p=r g b2g r a y(m a p);

s u b p l o t(2,2,1);s u b i m a g e(X,m a p);t i t l e('索引图')

s u b p l o t(2,2,3);s u b i m a g e(I);t i t l e('格式1灰度图')

s u b p l o t(2,2,4);s u b i m a g e(X,N e w m a p);t i t l e('格式2灰度图') 6、灰度图像→索引图像

【格式1】1、[X,m a p]=g r a y2i n d(I,n)

2、X=g r a y s l i c e(I,n)

3、X=g r a y s l i c e(I,v)

C L F,I=i m r e a d('r i c e.t i f');

[X1,m a p1]=g r a y2i n d(I,16);

X2=g r a y s l i c e(I,8);

X3=g r a y s l i c e(I,255*[00.210.230.260.300.350.61.0]');

s u b p l o t(2,2,1);s u b i m a g e(I);t i t l e('灰度图')

s u b p l o t(2,2,2);s u b i m a g e(X1,m a p1);t i t l e('16灰度级图')

s u b p l o t(2,2,3);s u b i m a g e(X2,h o t(8));t i t l e('均匀量化图')

s u b p l o t(2,2,4);s u b i m a g e(X3,j e t(8));t i t l e('非均匀量化图') 7、灰度图像→二值图像

【格式一】B W=d i t h e r(I)；

【格式二】B W=i m2b w(I,l e v e l)

【例】

C L F,I=i m r e a d('s a t u r n.t i f');

B W1=d i t h e r(I);

B W2=i m2b w(I,0.5);

s u b p l o t(2,2,1);s u b i m a g e(I);t i t l e('灰度图')

s u b p l o t(2,2,3);s u b i m a g e(B W1);t i t l e('抖动二值化')

s u b p l o t(2,2,4);s u b i m a g e(B W2);t i t l e('阈值二值化')

8、索引图像→二值图像

【格式】B W=i m2b w(X,m a p,l e v e l)

【输入】X可以是d o u b l e、u i n t16或u i n t8类型

【输出】B W为u i n t8类型

【例】

C L F,l o a d t r e e s;

B W=i m2b w(X,m a p,0.5);

s u b p l o t(1,2,1);s u b i m a g e(X,m a p);t i t l e('索引图')

s u b p l o t(1,2,2);s u b i m a g e(B W);t i t l e('二值图')

9、数据矩阵→灰度图像

【格式】1、I=m a t2g r a y(A,[m a x,m i n])

2、I=m a t2g r a y(A)

【例】

C L F;I=i m r e a d('r i c e.t i f');

A=f i l t e r2(f s p e c i a l('s o b e l'),I);

J=m a t2g r a y(A);

s u b p l o t(1,2,1);s u b i m a g e(I);t i t l e('原图')

s u b p l o t(1,2,2);s u b i m a g e(J);t i t l e('转换图')

I=i m r e a d('r i c e.p n g');

J=i m r e a d('c a m e r a m a n.t i f');

K=i m a d d(I,J,'u i n t16');

i m s h o w(K,[])

注意类型处理

>> X=uint8([255 0 75; 44 225 100]);

>> Y=uint8([50 50 50; 50 50 50]);

>> Z=imadd(X,Y)

结果：

Z =

255 50 125

94 255 150

两图相加

I=imread('rice.png');

J=imread('cameraman.tif');

K=imadd(I,J,'uint16');

imshow(K,[])

图加数

RGB=imread(‘peppers.png’);

RGB2=imadd(RGB,150);

Subplot(121);imshow(RGB);

Subplot(122);imshow(RGB2);

加噪声

?I=imread('eight.tif');

?>> J=imnoise(I,'salt & pepper',0.02);

?>> imshow(I)

?>> figure,imshow(J)

减法

?J=imread('cameraman.tif');

?I=uint8(filter2(fspecial('gaussian'),J));

?K=imabsdiff(J,I);

?imshow(K

乘法

?I=imread('moon.tif');

?I16=uint16(I);

?J=immultiply(I16,I16);

?imshow(I);

?figure,imshow(J)

除法

?I=imread('rice.png');

?>> J=imdivide(I,2);

?>> imshow(I)

?>> figure(2),imshow(J)

例：计算两幅图像的均值

>> I=imread('rice.png');

>> J=imread('cameraman.tif');

>> k=imdivide(imadd(I,J), 2);

>> imshow(k

>> I=imread('rice.png');

>> J=imread('cameraman.tif'); >> L=imlincomb(.5,I,.5,J);

>> imshow(L)

补运算

>> bw=imread('text.png');

>> bw2=imcomplement(bw);

>> imshow(bw);

>> figure,imshow(bw2)

补或非

clear all

clc

I=imread('rice.png');

J=imdivide(I,2);

K1=bitand(I,J);

imshow(K1)

title('图像的位与运算')

K2=bitcmp(I,8);

figure(2)

imshow(K2);

title('图像的位补运算')

K3=bitor(I,J);

figure(3)

imshow(K3);

title('图像的位或运算')

K4=bitxor(I,J);

figure(4)

imshow(K4);

title('图像的位异或运算')

K5=bitshift(I,2);

figure(5)

imshow(K5);

title('图像的位移位运算')

A=magic(5);

x=[19.5,23.5];

y=[8.0 12.0];

Image(A,?XDAT A?,x, ?YDAT A?,y,)；

axis image, colormap(jet(25))

图像缩放

B=imresize(A,M,method)

?A:原图像矩阵

?M:放大倍数。

?M>1放大

?0

?Method:插值方法

clear all

clc

close all

I=imread('peppers.png');

imshow(I);

I1=imresize(I,2,'bilinear');

figure,imshow(I1);

I2=imresize(I,2,'bicubic');

figure,imshow(I2);

I3=imresize(I,0.3,'nearest');

figure,imshow(I3);

图像旋转

clear all

clc

close all

I=imread('peppers.png');

imshow(I);

I1=imrotate(I,30,'nearest');

figure,imshow(I1);

I2=imrotate(I,30,'bilinear');

figure,imshow(I2);

I3=imrotate(I,80,'bicubic');

figure,imshow(I3);

I4=imrotate(I,30,'bicubic','crop');

figure,imshow(I4);

图像剪切

I = imread(‘circuit.tif');

I2 = imcrop(I,[75 68 130 112]);

figure,imshow(I)

figure,imshow(I2)

例10.1 利用函数impixel显示像素信息imshow canoe.tif

vals=impixel

例10.2 利用函数improfile绘制线段像素分布图

?imshow peppers.png

?improfle

分段线性变换

例9.3 对于原图pout.tif，将其小于30的灰度值不变，将30到150的灰度值拉伸到30到200，同时压缩150到255的灰度值到200与255间。

I=imread('pout.tif');

imshow(I);

I=double(I);

[M,N]=size(I);

for k = 1:M

for l = 1:N

if I(k,l)<=30

I(k,l)=I(k,l);

elseif I(k,l)<=150

I(k,l)=(200-30)/(150-30)*( I(k,l)-30)+30;

else

I(k,l)=(255-200)/(155-150)*( I(k,l)-150)+200;

end

end

end

figure(2);

imshow(uint8(I));

imadjust函数可以实现图像的灰度变换

clc;

I=imread('pout.tif');

imshow(I);

J=imadjust(I,[0.3 0.7],[0 1],1 );

figure;

imshow(J);

J=imadjust(I,[0.3 0.7],[0 1],0.5);

figure;

imshow(J);

J=imadjust(I,[0.3 0.7],[0 1],1.5);

figure;

imshow(J);

以下程序代码示例就是用imadjust函数将输入图像的灰度范围从[0 51]映射到[128 255]，输入图像如图所示：

I=imread('cameraman.tif');

J=imadjust(I,[0 0.2],[0.5 1]);

imshow(I)

figure,imshow(J

直方图均衡化

I=imread('tire.tif');

imshow(I)

figure;imhist(I);

J=histeq(I,64);

figure;

imshow(J);

figure;

imhist(J);

J=histeq(I,32);

figure;

imshow(J);

利用adapthisteq函数来调整图像的对比度问题I=imread('pout.tif');

J=adapthisteq(I);

imshow(I)

figure,imshow(J)

imhist(J,64)

锐化（拉普拉斯）滤波器

例3

?f=imread('moon.tif');

?imshow(f)

?w4=fspecial('laplacian',0);

?w8=[1,1,1;1,-8,1;1,1,1];

?f=im2double(f);

?g4=f-imfilter(f,w4,'replicate');

?g8=f-imfilter(f,w8,'replicate');

?imshow(f)

?figure,imshow(g4)

?figure,imshow(g8)

锐化（拉普拉斯）滤波器例4

?f=imread('moon.tif');

?imshow(f)

?w4=fspecial('laplacian',0);

?w8=[1,1,1;1,-8,1;1,1,1];

?f=im2double(f);

?g4=f-imfilter(f,w4,'replicate');

?g8=f-imfilter(f,w8,'replicate');

?imshow(f)

?figure,imshow(g4)

?figure,imshow(g8)