matlab简单程序

imread：读图像存入矩阵中

rgb2gray：RGB彩图转灰度图

fun是一个函数，作用是将返回一个与输入矩阵x一样尺寸的矩阵，其中每个元素都是输入矩阵的x标准差

blkproc：对图像分块处理，本程序中表示分成32×32的块，每个块应用fun函数

通俗点讲，最后效果是把输入彩图转灰度图并马赛克化。

这个net是个结构体类型,其中结构体中的变量LW又是个元胞数组。

net.LW{2,1}=net.LW{2,1}*0.01这句话是把结构体net中的lw元胞数组的第2行第一个元素变为原来的0.01倍

下面一句意思类似。

matlab元胞数组：http://www.3800hk.com/Article/jxdz/matlab/matlabjc/2005-08-06/Article_35366.html

matlab结构体：http://hi.baidu.com/%C6%BF%D6%D0%B5%C4%C5%AE%CE%D7/blog/item/91939438af235ef43a87ce08.html

%%(1)题

I=imread('abc.jpg')

II=rgb2gray(I)

imshow(II)

for i=1:size(II,1)

for j=1:size(II,2)

if (0<=II(i,j) &&II(i,j)<=50)

III(i,j)=2

elseif (50<II(i,j) &&II(i,j)<=100)

III(i,j)=8

elseif (100<II(i,j) &&II(i,j)<=150)

III(i,j)=16

elseif (150<II(i,j) &&II(i,j)<=200)

III(i,j)=64

elseif(200<II(i,j) &&II(i,j)<=250)

III(i,j)=128

end

end

end

figure

imshow(III)

%%(2)题

x=linspace(0,2*pi,500)

for i=1:length(x)

if (x(i)<=pi)

y1(i)=1

else

y1(i)=0

end

if (x(i)<=pi)

y2(i)=x(i)/pi

else

y2(i)=(2*pi-x(i))/pi

end

y3(i)=sin(x(i))

end

figureplot(x,y1)

figureplot(x,y2)

figureplot(x,y3)

figure

subplot(1,3,1)plot(x,y1)

subplot(1,3,2)plot(x,y2)

subplot(1,3,3)plot(x,y3)

figureplot(x,y1,x,y2,x,y3)legend('方波','三角波','正弦')

%%(3)题

eps=0.001

xl=3fl=xl^3-xl^2-6*xl-7 %不支持函数就直接写。

xr=4fr=xr^3-xr^2-6*xr-7

while (1)

xm=(xl+xr)/2fm=xm^3-xm^2-6*xm-7

if (abs(fm)<eps)

break

elseif (fm*fl<0)

xr=xmfr=fm

else

xl=xmfl=fm

end

end

x=xm

f(x)

---